function [bonds,components] = rbg2abg(rbonds,rstrokes,rcomponents,rports,infofile)
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.10 1997/03/17 13:45:42 peterg
% %% Added more error info.
% %%
% %% Revision 1.9 1996/12/31 11:25:57 peterg
% %% Clearer error messages for incorrect ports.
% %%
% %% Revision 1.8 1996/12/04 21:52:39 peterg
% %% Now uses fopen.
% %%
% %% Revision 1.7 1996/11/01 18:01:57 peterg
% %% Rationalised port ordering.
% %% Fixed port bug.
% %%
% %% Revision 1.6 1996/08/25 08:27:14 peter
% %% Now checks ports correctely - I hope.
% %%
% %% Revision 1.5 1996/08/24 19:21:26 peter
% %% More specific error messages.
% %%
% %% Revision 1.4 1996/08/24 18:00:33 peter
% %% Fixed bug with finding ports.
% %%
% %% Revision 1.3 1996/08/09 08:26:35 peter
% %% Cosmetic tidy up.
% %%
% %% Revision 1.2 1996/08/04 18:37:57 peter
% %% Fixed no causal strokes bug.
% %%
% %% Revision 1.1 1996/08/04 18:30:14 peter
% %% Initial revision
% %%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if nargin<5
infofile='stdout';
else
fnum = fopen(infofile, 'w');
end;
% 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');
end;
% Find number of bonds
[n_bonds,columns] = size(rbonds);
if (columns ~= 6)&(n_bonds>0)
error('Incorrect rbonds matrix: must have 6 columns');
end;
% Find number of components
[n_components,columns] = size(rcomponents);
% Find number of ports
[n_ports,columns] = size(rports);
% Determine coordinates of the arrow end of the bond and the 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 the bond end nearest to 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 number
for i = 1:n_ports
[i rports(i,1:2)/scale rports(i,3)]
near_bond = adjbond(rports(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\n", ...
rports(i,1)/scale, rports(i,2)/scale));
end;
signed_bond = near_bond(1)*sign(1.5-near_bond(2));
port_near_bond(i,:) = [signed_bond, rports(i,3)];
end;
% Locate the components at the ends of each bond
% col 1 of comp_near_bond contain the component nearest to the arrow end
% col 2 of comp_near_bond contain the component nearest to the other end
for i = 1:n_bonds
comp_near_bond(i,:) = adjcomp(arrow_end(i,:),other_end(i,:),rcomponents);
end;
% Produce a list of bonds on each component - sorted if explicit port numbers
for i = 1:n_components
% There are n bonds on this component with corresponding index
[index,n] = getindex(comp_near_bond,i);
% Error message in case we need it!
port_error = sprintf(...
'Component at (%1.3f,%1.3f) has inconsistent port numbers', ...
rcomponents(i,1)/scale, rcomponents(i,2)/scale);
if index(1,1) ~= 0 % Then its a genuine component
% Create the signed list of bonds on this component
one = ones(n,1);
bond_list = index(:,1); % bond at component
bond_end = index(:,2); % which end of bond at component?
direction = -sign(bond_end-1.5*one);
signed_bond_list = bond_list.*direction;
% Default sort of bonds (ie no change)
sort_index = [1:n]';
if n_ports>0 % then there are some numbered ports
% so find those associated with the bonds on this component.
k=0; port_number=[];
for j = 1:n
b = signed_bond_list(j);
% Find the port label on component end of bond (if any)
[port_index,m] = getindex(port_near_bond(:,1),b);
if m==1
k=k+1;
port_number(k,1) = port_near_bond(port_index,2);
end;
end;
% Must have a lable for each port or non at all
if k==n % Then a full set of port numbers
[sorted,sort_index]=sort(port_number);
%Check that all the numbers (1..n) are there
if sum(sorted==[1:n]')~=n % There is something wrong
mtt_info(port_error,fnum);
mtt_info(sprintf('\t it must have consecutive ports labels from 1:%1.0f\n', n), ...
fnum);
port_str = sprintf('\t the port labels are: ');
for kk=1:k
port_str=sprintf('%s %1.0f', port_str, port_number(kk));
end;
mtt_info(port_str, fnum);
end;
else
if k~=0
mtt_info(port_error,fnum);
mtt_info(sprintf('\t it must have 0 or %1.0f ports\n', n), fnum);
port_str = sprintf('\t the port labels are: ');
for kk=1:k
port_str=sprintf('%s %1.0f', port_str, port_number(kk));
end;
mtt_info(port_str, fnum);
end;
end;
end;
end;
% Write out bond list sorted by port number (if any)
for j = 1:length(sort_index)
jj = sort_index(j);
components(i,j) = signed_bond_list(jj);
end;
end;
% Deduce causality from the strokes (if any).
causality = zeros(n_bonds,2);
if n_strokes>0
% Find out 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);
% Deduce bond causality from the strokes
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);
end;
% 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 % bicausal stroke
% Find out whether stroke is on flow side of bond
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;
end;
end;
end;
bonds = causality;
fclose(fnum);