22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
+
+
+
|
% Copyright (c) P.J.Gawthrop, 1996.
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.6 1997/08/19 09:41:47 peterg
% %% Some debugging lines added.
% %%
% %% Revision 1.5 1997/05/16 07:33:45 peterg
% %% Now checks to see if sub system is a simple component before
% %% recursion.
% %% 0 --> zero
% %% 1 --> one
% %%
% %% Revision 1.4 1996/12/07 21:34:52 peterg
|
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
|
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
|
-
+
|
filenum = fopen(fig_name, 'a');
% Get the raw and the processed bonds
eval(['[rbonds,rstrokes,rcomponents] = ', system_type, '_rbg;']);
eval(['[bonds] = ', system_type, '_abg;']);
% Original number of bonds
[n_bonds,junk] = size(bonds);
[n_bonds,junk] = size(rbonds);
% Get the causal bonds
eval(['[cbonds,status]=', full_name, '_cbg;']);
% Check sizes
[N_components,Columns] = size(rcomponents);
if (Columns ~= 13)
|
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
|
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
|
-
-
+
+
|
unit_arrow_vector = arrow_vector./(length2d(arrow_vector)*[1 1]);
bond_vector = (arrow_end - other_end);
unit_bond_vector = bond_vector./(length2d(bond_vector)*[1 1]);
unit_stroke_vector = (rot*unit_bond_vector')';
% Get indices of bonds with changed causality -- but ignore the extra bonds
% due to vector bond expansion
changed_e = bonds(:,1)~=cbonds(1:n_bonds,1)
changed_f = bonds(:,2)~=cbonds(1:n_bonds,2)
changed_e = bonds(1:n_bonds,1)~=cbonds(1:n_bonds,1)
changed_f = bonds(1:n_bonds,2)~=cbonds(1:n_bonds,2)
changed = changed_e|changed_f
index_e = getindex(changed_e,1)'
index_f = getindex(changed_f,1)'
index = getindex(changed,1)';
% Print the new strokes in fig format
if index(1,1)>0
|