## -*-octave-*-
function [bonds] = \
rbg2ibg(name,rbonds,rstrokes,rcomponents,port_coord,port_name, \
infofile, errorfile)
rbonds
rstrokes
rcomponents
port_coord
port_name
## Default to no components
components = [0];
## Xfig scaling factor
scale = 1200.0/2.54546;
## Rotation matrix
rot = [0, -1; 1, 0];
## Find number of strokes
[n_strokes,columns] = size(rstrokes);
if ((columns ~= 4) & (n_strokes > 0))
error('Incorrect rstrokes matrix: must have 4 columns');
endif
## Find number of bonds
[n_bonds,columns] = size(rbonds);
if ((columns ~= 6) & (n_bonds > 0))
error('Incorrect rbonds matrix: must have 6 columns');
endif
## Find number of components
[n_components,columns] = size(rcomponents);
## Find number of ports referred to
[n_ports,columns] = size(port_coord);
## If port_name empty, make it empty string
if (length(port_name) == 0)
port_name="";
endif
## Determine coordinates of arrow end of bond and other end
other_end_1 = rbonds(:,1:2);
arrow_end = rbonds(:,3:4);
other_end_2 = rbonds(:,5:6);
distance_1 = length2d(other_end_1 - arrow_end);
distance_2 = length2d(other_end_2 - arrow_end);
which_end = (distance_1 > distance_2) * [1, 1];
one = ones(size(which_end));
other_end = which_end .* other_end_1 + (one - which_end) .* \
other_end_2;
arrow_vector = (which_end .* other_end_2 + (one - which_end) .* \
other_end_1) - arrow_end;
## Locate bond end nearest each port
## col 1 of port_near_bond contains a signed bond number (+ for arrow
## end)
## col 2 of port_near_bond contains the corresponding port index
port_bond = [];
for i = 1:n_ports
near_bond = adjbond(port_coord(i,1:2), arrow_end, other_end);
[rows,cols] = size(near_bond);
if (rows > 1)
error(sprintf ...
("A port is near to more than one bond at coordinates \
%g,%g %s\n", ...
port_coord(i,1)/scale, port_coord(i,2)/scale, \
deblank(port_name(i,:))));
endif
## The (signed) bond corresponding to the ith port label
port_bond(i) = near_bond(1) * sign(1.5 - near_bond(2));
endfor
port_bond
## Now have (signed) bond (port_bond(i)) corresponding to the ith port
## Create inverse mapping
for i = 1:n_bonds
eval(sprintf('bond_port_head%i = "[]"', i))
eval(sprintf('bond_port_tail%i = "[]"', i))
endfor
for i = 1:n_ports
if (port_bond(i) > 0)
eval(sprintf('bond_port_head%i = "%s"', port_bond(i), \
deblank(port_name(i,:))))
else
eval(sprintf('bond_port_tail%i = "%s"', -port_bond(i), \
deblank(port_name(i,:))))
endif
endfor
## Locate the components at the ends of each bond
## col 1 of comp_near_bond contains component nearest to the arrow end
## col 2 of comp_near_bond contains component nearest other end
for i = 1:n_bonds
comp_near_bond(i,:) = adjcomp(arrow_end(i,:), other_end(i,:), \
rcomponents);
endfor
comp_near_bond
## Deduce causality from the strokes (if any)
causality = zeros(n_bonds,2);
if (n_strokes > 0)
## Find location of centre and ends of stroke
stroke_end_1 = [rstrokes(:,1), rstrokes(:,2)];
stroke_end_2 = [rstrokes(:,3), rstrokes(:,4)];
stroke_centre = (stroke_end_1 + stroke_end_2)/2;
stroke_vector = (stroke_end_1 - stroke_end_2);
stroke_length = length2d(stroke_vector);
for i = 1:n_strokes
stroke = [stroke_centre(i,:)
stroke_end_1(i,:)
stroke_end_2(i,:)];
## Find the nearest bond end
[index,distance] = adjbond(stroke(1,:),arrow_end,other_end);
if (distance > (2 * stroke_length(i)))
info = sprintf('Stroke at (%4.3f,%4.3f) is %4.3f away from the nearest bond\n', ...
stroke(1,1)/scale, stroke(1,2)/scale, \
distance/scale);
endif
## Bond end coordinates
j = index(1,1);
which_end = (index(1,2) == 1);
bond_end = arrow_end(j,:) * which_end + other_end(j,:) * (1 - \
which_end);
## Now decide which bit of the stroke is nearest
stroke_index = adjbond(bond_end, stroke, zeros(size(stroke)));
if (stroke_index(1) == 1) # uni-causal stroke
causality(j,1:2) = (2 * which_end - 1) * [1, 1];
else # bi-causal stroke
stroke_direction = stroke(1,:) - stroke(stroke_index(1),:);
flow_side = stroke_direction * arrow_vector(j,:)' > 0;
causality(j,1+flow_side) = 2 * which_end - 1;
endif
endfor
endif
causality
## Write data
for i = 1:n_bonds
[hc_type, hc_name] = eval([name, '_cmp(comp_near_bond(i,1))']);
[tc_type, tc_name] = eval([name, '_cmp(comp_near_bond(i,2))']);
## components
eval(sprintf("bonds.bond%i.head.component\t= '%s:%s'", i, hc_type, \
hc_name));
eval(sprintf("bonds.bond%i.tail.component\t= '%s:%s'", i, tc_type, \
tc_name));
## ports
eval(sprintf("bonds.bond%i.head.ports\t= bond_port_head%i", i, i));
eval(sprintf("bonds.bond%i.tail.ports\t= bond_port_tail%i", i, i));
## causality
if (causality(i,1) == 1)
effort_causality = "head"
elseif (causality(i,1) == -1)
effort_causality = "tail"
else
effort_causality = "none"
endif
if (causality(i,2) == 1)
flow_causality = "tail"
elseif (causality(i,2) == -1)
flow_causality = "head"
else
flow_causality = "none"
endif
eval(sprintf("bonds.bond%i.causality.effort\t= '%s'", i, effort_causality));
eval(sprintf("bonds.bond%i.causality.flow\t= '%s'", i, flow_causality));
endfor
endfunction