1
2
3
4
5
6
7
8
9
10
11
12
13
14
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
+
+
+
|
function [bonds,components] = rbg2abg(rbonds,rstrokes,rcomponents,rports,infofile)
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% 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
|
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
|
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
|
+
+
+
-
+
+
+
+
+
+
-
+
+
+
-
+
+
-
-
-
+
+
+
-
-
-
+
+
+
-
-
-
-
-
-
-
-
-
+
-
-
-
+
+
-
-
+
+
+
+
-
+
-
+
+
+
-
-
-
+
+
+
+
-
+
-
+
-
-
-
-
|
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
near_bond = adjbond(rports(i,1:2),arrow_end,other_end);
signed_bond = near_bond(1)*sign(1.5-near_bond(2))
port_near_bond(i,:) = [near_bond, rports(i,3)];
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!
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
% Are the component ports numbered? (either they all are or none are)
k=0;
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 = bond_list(j); e = bond_end(j);
% Find all ports on this bond
[port_indices,m] = getindex(port_near_bond(:,1),b);
b = signed_bond_list(j);
% Find the port on component end of bond (if any)
[port_index,m] = getindex(port_near_bond(:,1),b);
% Now find the one at this end - if any
port_index = 0;
for l=1:m
if port_near_bond(port_indices(l),2)==e
port_index = port_indices(l);
break;
end;
end;
if port_index>0
% and put the corresponding number in the list
port_number(k,1) = port_near_bond(port_index,3);
end;
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
if sum(sorted==[1:n]')~=n % The there is something wrong
if sum(sorted==[1:n]')~=n % There is something wrong
mtt_info(port_error,fnum);
mtt_info(sprintf('\t it must have ports from 1:%1.0f\n', n), fnum);
mtt_info(sprintf('\t it must have ports from 1:%1.0f\n', n), ...
fnum);
end;
else
elseif k~=0
mtt_info(port_error,fnum);
mtt_info(sprintf('\t it must have 0 or %1.0f ports\n', n), fnum);
end;
end;
% direction is 1 if arrow at component else -1
direction = -sign(index(:,2)-1.5*one);
signed_bond_list = bond_list.*direction;
|