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
|
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
|
+
+
+
+
-
+
-
+
-
+
|
% Copyright (c) P.J. Gawthrop, 1996.
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% %% $Id$
% %% $Log$
% %% Revision 1.3 1996/12/05 10:02:45 peterg
% %% Explicit generation of strong-causality equations -- oneequation not
% %% used now.
% %%
% %% Revision 1.2 1996/08/08 18:09:10 peter
% %% Changed to finenumber format.
% %%
% %% Revision 1.1 1996/08/08 16:38:50 peter
% %% Initial revision
% %%
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Same causality as junction
column = (3-jun)/2;
[causing_bond,n,other_bonds,m] = getindex(bonds(:,column),jun);
[causing_bond,n,other_bonds,m] = getindex(bonds(:,column),jun)
for i=other_bonds'
fprintf(filenumber, '%s \t:=\t%s;\n', ...
varname(name, bond_number(i),jun), ...
varname(name, bond_number(causing_bond),jun) ...
);
end;
% Opposite causality to junction
column = (3+jun)/2;
[caused_bond,n,other_bonds,m] = getindex(bonds(:,column),jun);
[caused_bond,n,other_bonds,m] = getindex(bonds(:,column),jun)
fprintf(filenumber, '%s\t:= \n', varname(name, bond_number(caused_bond),-jun));
for i=other_bonds'
term_sign = -direction(caused_bond,column)*direction(i,column);
term_sign = -direction(caused_bond,column)*direction(i,column)
fprintf(filenumber, '\t\t%s %s\n', sign2name(term_sign), ...
varname(name, bond_number(i),-jun));
end;
fprintf(filenumber, ';\n');
|