1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
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
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
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
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
|
-
-
-
-
-
+
+
+
+
+
+
-
-
-
-
-
+
+
+
+
-
-
+
-
-
-
+
+
-
-
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
+
+
+
-
-
-
+
+
+
-
-
-
+
+
+
-
-
-
+
+
+
-
-
-
+
+
+
-
-
-
+
+
+
-
-
+
+
-
-
-
-
-
-
-
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
-
-
+
+
-
-
-
-
-
+
+
+
+
+
-
-
+
+
-
-
-
-
+
+
+
+
-
-
-
-
-
-
+
+
+
+
+
+
-
+
-
-
-
-
-
+
+
+
+
+
-
-
+
+
-
-
-
-
-
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
+
+
+
-
-
-
+
+
+
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
-
-
-
-
+
+
+
+
-
-
-
-
+
+
+
+
-
+
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
-
-
-
-
+
+
+
+
-
-
-
-
-
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
+
+
+
+
+
+
+
-
-
-
-
+
+
+
+
-
-
-
-
+
+
+
+
-
-
-
-
-
+
+
+
+
+
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
+
-
-
-
-
-
-
+
+
+
-
+
-
-
-
-
-
-
-
-
-
-
-
-
|
function cbg2fig(system_name, ...
system_type, full_name, ...
stroke_length, stroke_thickness, stroke_colour, ...
comp_font, comp_colour_u, comp_colour_o)
% cbg2fig - converts causal bg to figure
%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %%%%% Model Transformation Tools %%%%%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Matlab function cbg_m2fig
% Causal bond graph: mfile format to fig file format
% The resultant fig file is the original _abg.fig with
% additional causal strokes superimposed.
% cbg2fig(system_name, ...
% system_type, full_name, ...
% stroke_length, stroke_thickness, stroke_colour, ...
% comp_font, comp_colour_u, comp_colour_o)
% P.J.Gawthrop May 1996
% Copyright (c) P.J.Gawthrop, 1996.
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.10 1998/12/14 15:27:32 peterg
% %% Added separate handling of ports ( poss zero)
% %%
% %% Revision 1.9 1998/08/25 06:43:02 peterg
% %% Revised (partially) for data strucures - needs to include graphic info
% %% in abg file.
% %%
% %% Revision 1.8 1998/04/04 10:45:01 peterg
% %% Don't do strokes on port bonds
% %%
% %% Revision 1.7 1997/08/19 09:49:19 peterg
% %% Modified to take account of the expanded vector bonds. Only displays
% %% causality corresponding to the bond connecting the first element of
% %% the vector ports.
% %%
% %% 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
% %% Tests for null string with strcmp
% %%
% %% Revision 1.3 1996/08/08 15:52:28 peter
% %% Recursive version.
% %% Fails due to octave bug - reported.
% %%
% %% Revision 1.2 1996/08/05 20:15:39 peter
% %% Prepared for recursive version.
% %%
% %% Revision 1.1 1996/08/05 18:12:25 peter
% %% Initial revision
% %%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if nargin<4
stroke_length = 20;
end;
if nargin<4
stroke_length = 20;
end;
if nargin<5
stroke_thickness = 2;
end;
if nargin<5
stroke_thickness = 2;
end;
if nargin<6
stroke_colour = 1; %Blue
end;
if nargin<6
stroke_colour = 1; %Blue
end;
if nargin<7
comp_font = 18; %Helvetica bold
end;
if nargin<7
comp_font = 18; %Helvetica bold
end;
if nargin<8
comp_colour_u = 12; %Green
end;
if nargin<8
comp_colour_u = 12; %Green
end;
if nargin<9
comp_colour_o = 4; %Red
end;
if nargin<9
comp_colour_o = 4; %Red
end;
% Create a back slash '\' character.
bs = setstr(92);
% Create the (full) system name
if strcmp(full_name,'')
full_name = system_name;
system_type = system_name;
else
full_name = [full_name, '_', system_name];
end;
full_name_type = [full_name, '_', system_type];
% $$$ fig_name = [full_name_type, '_cbg.fig'];
fig_name = [full_name, '_cbg.fig']
% Return if initial fig file doesn't exist
if exist(fig_name)~=2
return
end;
full_name_type = [full_name, '_', system_type];
## $$$ fig_name = [full_name_type, '_cbg.fig'];
fig_name = [full_name, '_cbg.fig']
## Return if initial fig file doesn't exist
if exist(fig_name)~=2
return
end;
% Setup file - append to the fig file
filenum = fopen(fig_name, 'a');
% Get the raw and the processed bonds
eval(['[rbonds,rstrokes,rcomponents] = ', system_type, '_rbg;']);
%eval(['[bonds,components,n_ports] = ', system_type, '_abg;']);
eval(["ABG = ", system_type, "_abg;"]);
bonds=ABG.bonds;
% Original number of bonds
[n_bonds,junk] = size(rbonds);
% Get the causal bonds
%eval(['[cbonds,status]=', full_name, '_cbg;']);
eval(["CBG = ", full_name, "_cbg;"]);
cbonds=CBG.bonds;
% Check sizes
[Rows,Columns] = size(rcomponents);
if (Columns ~= 13)
error('Incorrect rcomponents matrix: must have 13 columns');
end;
M_components = Columns;
[N_components, Columns] = size(CBG.subsystemlist);
[N_components, Columns] = size(CBG.subsystemlist);
if struct_contains(CBG,'portlist')
[N_ports, Columns] = size(CBG.portlist);
else
N_ports = 0;
end
if struct_contains(CBG,'portlist')
[N_ports, Columns] = size(CBG.portlist);
else
N_ports = 0;
end
% Rotation matrix
rot = [0 -1; 1 0];
% Determine coordinates of the arrow end of the bond and the other end
% and other geometry
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_barb = which_end.*other_end_2 + (one-which_end).*other_end_1;
arrow_vector = arrow_barb-arrow_end;
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')';
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_barb = which_end.*other_end_2 + (one-which_end).*other_end_1;
arrow_vector = arrow_barb-arrow_end;
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: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;
% Don't do port bonds
#if n_ports>0
# port_bonds = sort(abs(components(1:n_ports,1)));
# changed(port_bonds) = zeros(n_ports,1);
# changed_e(port_bonds) = zeros(n_ports,1);
# changed_f(port_bonds) = zeros(n_ports,1);
#end
index_e = getindex(changed_e,1)'
index_f = getindex(changed_f,1)'
index_ef = getindex(changed,1)';
index_e = getindex(changed_e,1)'
index_f = getindex(changed_f,1)'
index_ef = getindex(changed,1)';
% Print the new strokes in fig format
if index_ef(1,1)>0
for i = index_e % Do the effort stroke - opp. side to arrow
if cbonds(i,1)==1 % Stroke at arrow end
stroke_end_1 = arrow_end(i,:);
else
stroke_end_1 = other_end(i,:);
end;
sig = sign(unit_arrow_vector(i,:)*unit_stroke_vector(i,:)');
stroke_end_2 = stroke_end_1 - stroke_length*sig*unit_stroke_vector(i,:);
if cbonds(i,1)==1 # Stroke at arrow end
stroke_end_1 = arrow_end(i,:);
else
stroke_end_1 = other_end(i,:);
end;
sig = sign(unit_arrow_vector(i,:)*unit_stroke_vector(i,:)');
stroke_end_2 = stroke_end_1 - stroke_length*sig*unit_stroke_vector(i,:);
%print the fig3 format firstline spec.
polyline = 2;
firstline = fig3(polyline,stroke_thickness,stroke_colour);
fprintf(filenum, '%s\n', firstline);
## print the fig3 format firstline spec.
polyline = 2;
firstline = fig3(polyline,stroke_thickness,stroke_colour);
fprintf(filenum, '%s\n', firstline);
fprintf(filenum, ' %4.0f %4.0f %4.0f %4.0f \n', ...
stroke_end_1(1), stroke_end_1(2), ...
stroke_end_2(1), stroke_end_2(2) );
end;
fprintf(filenum, ' %4.0f %4.0f %4.0f %4.0f \n', ...
stroke_end_1(1), stroke_end_1(2), ...
stroke_end_2(1), stroke_end_2(2) );
end;
for i = index_f % Do the flow stroke - same side as arrow
for i = index_f # Do the flow stroke - same side as arrow
if cbonds(i,2)==1 % Stroke at arrow end
stroke_end_1 = arrow_end(i,:);
else
stroke_end_1 = other_end(i,:);
end;
sig = sign(unit_arrow_vector(i,:)*unit_stroke_vector(i,:)');
stroke_end_2 = stroke_end_1 + stroke_length*sig*unit_stroke_vector(i,:);
if cbonds(i,2)==1 # Stroke at arrow end
stroke_end_1 = arrow_end(i,:);
else
stroke_end_1 = other_end(i,:);
end;
sig = sign(unit_arrow_vector(i,:)*unit_stroke_vector(i,:)');
stroke_end_2 = stroke_end_1 + stroke_length*sig*unit_stroke_vector(i,:);
%print the fig3 format firstline spec.
polyline = 2;
firstline = fig3(polyline,stroke_thickness,stroke_colour);
fprintf(filenum, '%s\n', firstline);
## print the fig3 format firstline spec.
polyline = 2;
firstline = fig3(polyline,stroke_thickness,stroke_colour);
fprintf(filenum, '%s\n', firstline);
fprintf(filenum, ' %4.0f %4.0f %4.0f %4.0f \n', ...
stroke_end_1(1), stroke_end_1(2), ...
stroke_end_2(1), stroke_end_2(2) );
end;
end;
fprintf(filenum, ' %4.0f %4.0f %4.0f %4.0f \n', ...
stroke_end_1(1), stroke_end_1(2), ...
stroke_end_2(1), stroke_end_2(2) );
end;
end;
% Print all the components - coloured acording to causality.
for i = 1:N_ports+N_components
if i>N_ports # Subsystem
comp_name = CBG.subsystemlist(i-N_ports,:);
eval(["comp_status = CBG.subsystems.", comp_name, ".status;"]);
else
comp_name = CBG.portlist(i,:);
eval(["comp_status = CBG.ports.", comp_name, ".status;"]);
end
fig_params = rcomponents(i,3:M_components);
coords = rcomponents(i,1:2);
if comp_status==-1 %Then under causal
fig_params(3) = comp_colour_u;
fig_params(6) = comp_font;
end;
fig_params = rcomponents(i,3:M_components);
coords = rcomponents(i,1:2);
if comp_status==-1 # Then under causal
fig_params(3) = comp_colour_u;
fig_params(6) = comp_font;
end;
if comp_status==1 %Then over causal
fig_params(3) = comp_colour_o;
fig_params(6) = comp_font;
end;
if comp_status==1 # Then over causal
fig_params(3) = comp_colour_o;
fig_params(6) = comp_font;
end;
%Now print the component in fig format
eval(['[comp_type,comp_name] = ', system_type, '_cmp(i)']);
if index(comp_name,"mtt")==1 # Its a dummy name
## Now print the component in fig format
eval(['[comp_type,comp_name] = ', system_type, '_cmp(i)']);
if index(comp_name,"mtt")==1 # Its a dummy name
typename = comp_type; # just show type
else
typename = [comp_type,":",comp_name];
endif;
Terminator = [bs, '001'];
for j = 1:length(fig_params)
fprintf(filenum, '%1.0f ', fig_params(j));
end;
Terminator = [bs, '001'];
for j = 1:length(fig_params)
fprintf(filenum, '%1.0f ', fig_params(j));
end;
fprintf(filenum, '%1.0f %1.0f ', coords(1), coords(2));
fprintf(filenum, '%s%s\n', typename, Terminator);
# % If it's a subsystem (ie not a component), do the fig file for that as
# % well
# if comp_type=='0'
# comp_type='zero';
# endif
# if comp_type=='1'
# comp_type='one';
# endif
# if (exist([comp_type,'_cause'])==0)
# cbg2fig(comp_name, ...
# comp_type, full_name, ...
# stroke_length, stroke_thickness, stroke_colour, ...
# comp_font, comp_colour_u, comp_colour_o);
# endif
fprintf(filenum, '%1.0f %1.0f ', coords(1), coords(2));
fprintf(filenum, '%s%s\n', typename, Terminator);
# % If it's a subsystem (ie not a component), do the fig file for that as
# % well
# if comp_type=='0'
# comp_type='zero';
# endif
# if comp_type=='1'
# comp_type='one';
# endif
# if (exist([comp_type,'_cause'])==0)
# cbg2fig(comp_name, ...
# comp_type, full_name, ...
# stroke_length, stroke_thickness, stroke_colour, ...
# comp_font, comp_colour_u, comp_colour_o);
# endif
endfor
endfor
% Close the file
fclose(filenum);
return
endfunction
|