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
|
else
lang="r"
fi
# Inform user
echo Creating $1_smx.$lang
# Remove the old log file
rm -f ode2smx_r.log
# Use reduce to accomplish the transformation
reduce << EOF >ode2smx_lang.log
lang := $lang;
in "$1_def.r";
in "$1_ode.r";
OFF Echo;
%Procedure to write out the result according to the language used.
PROCEDURE WriteElement(name, i, j, element);
BEGIN
IF (element NEQ 0) THEN
BEGIN
IF (lang = r) THEN
write name, "(", i, ",", j, ") := ", element;
IF (lang = m) THEN
write name, "[", i, ",", j, "] = ", element;
END
END;
%Set up output according to the language.
OFF NAT;
OUT "$1_smx.$lang";
% Set up output according to the language.
IF (lang = r) THEN
BEGIN
write "matrix mtta(", mttnx, ",", mttnx, ");";
write "matrix mttb(", mttnx, ",", mttnu, ");";
write "matrix mttc(", mttny, ",", mttnx, ");";
write "matrix mttd(", mttny, ",", mttnu, ");";
END;
% find MTTA : the A matrix
FOR j := 1:MTTNx DO
BEGIN
xj := MTTX(j,1);
FOR i := 1:MTTNx DO
WriteElement("mtta",i,j,df(MTTdx(i,1), xj, 1));
END;
%% Find MTTB : the B matrix
% FOR j := 1:MTTNu DO
% BEGIN
% uj := MTTU(j,1);
% FOR i := 1:MTTNx DO
% WriteElement("mttb",i,j,df(MTTdx(i,1), uj, 1));
|
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
|
|
>
>
>
>
>
|
|
<
|
>
>
>
>
>
>
>
>
>
>
|
>
>
>
>
>
>
>
>
|
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
|
else
lang="r"
fi
# Inform user
echo Creating $1_smx.$lang
echo Creating $1_smxx.$lang
echo Creating $1_smxtx.$lang
# Remove the old log file
rm -f ode2smx_r.log
rm -f $1_smx.$lang
rm -f $1_smxx.$lang
rm -f $1_smxtx.$lang
# Use reduce to accomplish the transformation
reduce << EOF >ode2smx_lang.log
lang := $lang;
in "$1_def.r";
in "$1_ode.r";
OFF Echo;
% Load the general translator package
LOAD GENTRAN;
GENTRANLANG!* := 'Pascal;
ON GENTRANSEG;
MAXEXPPRINTLEN!* := 80;
TEMPVARNUM!* := 1;
TEMPVARNAME!* := 'mtt_t;
% Matrix output function
in"$MTTPATH/trans/lang_matrix.r";
%Procedure to write out the result according to the language used.
PROCEDURE WriteElement(name, i, j, element);
BEGIN
IF (element NEQ 0) THEN
BEGIN
IF (lang = r) THEN
write name, "(", i, ",", j, ") := ", element;
IF (lang = m) THEN
GENTRAN foo :=: element;
END
END;
%Set up output according to the language.
OFF NAT;
GENTRANOUT "$1_smx.$lang";
% Set up output according to the language.
IF (lang = r) THEN
BEGIN
write "matrix mtta(", mttnx, ",", mttnx, ");";
write "matrix mttb(", mttnx, ",", mttnu, ");";
write "matrix mttc(", mttny, ",", mttnx, ");";
write "matrix mttd(", mttny, ",", mttnu, ");";
END;
% find MTTA : the A matrix
matrix mttAAx(MTTNx,1);
matrix mttAAtx(MTTNx,1);
FOR i := 1:MTTNx DO
BEGIN
xi := MTTX(i,1);
FOR j := 1:MTTNx DO
BEGIN
xj := MTTX(j,1);
a_ij := df(MTTdx(i,1), xj, 1);
a_ji := df(MTTdx(j,1), xi, 1);
IF i=j THEN
BEGIN
aa_ij := 1 - mttdt*a_ij;
aa_ji := aa_ij;
END
ELSE
BEGIN
aa_ij := -mttdt*a_ij;
aa_ji := -mttdt*a_ji;
END;
IF (aa_ij NEQ 0) THEN
GENTRAN mtta(i,j) ::=: aa_ij;
mttAAx(i,1) := mttAAx(i,1) + aa_ij*mkid(mttxx,j);
mttAAtx(i,1) := mttAAtx(i,1) + aa_ji*mkid(mttxx,j);
END;
END;
%% Find MTTB : the B matrix
% FOR j := 1:MTTNu DO
% BEGIN
% uj := MTTU(j,1);
% FOR i := 1:MTTNx DO
% WriteElement("mttb",i,j,df(MTTdx(i,1), uj, 1));
|
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
|
% FOR j := 1:MTTNu DO
% BEGIN
% xj := MTTU(j,1);
% WriteElement("mttd",i,j,df(MTTY(i,1), xj, 1));
% END;
%Shut the output according to the language.
SHUT "$1_smx.$lang";
EOF
if [ "$lang" = "m" ]; then
mv $1_smx.$lang mtt_junk
lang_header $1 smx m 'mttx,mttu' '[mtta]' > $1_smx.m
cat mtt_junk | mtt_p2m >> $1_smx.m
rm -f mtt_junk
fi
# Now invoke the standard error handling.
mtt_error_r ode2smx_lang.log
|
|
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
>
|
>
>
>
>
>
>
>
>
>
>
>
>
|
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
|
% FOR j := 1:MTTNu DO
% BEGIN
% xj := MTTU(j,1);
% WriteElement("mttd",i,j,df(MTTY(i,1), xj, 1));
% END;
%Shut the output according to the language.
GENTRANSHUT "$1_smx.$lang";
% The AAx matrix;
GENTRANOUT "$1_smxx.$lang";
FOR i := 1:MTTNx DO
BEGIN
ax_i := mttAAx(i,1);
IF (ax_i NEQ 0) THEN
GENTRAN mttax(i) ::=: ax_i;
END
GENTRANSHUT "$1_smxx.$lang";
% The AAtx matrix;
GENTRANOUT "$1_smxtx.$lang";
FOR i := 1:MTTNx DO
BEGIN
atx_i := mttAAtx(i,1);
IF (atx_i NEQ 0) THEN
GENTRAN mttax(i) ::=: atx_i;
END
GENTRANSHUT "$1_smxtx.$lang";
EOF
if [ "$lang" = "m" ]; then
mv $1_smx.$lang mtt_junk
lang_header $1 smx m 'mttx,mttu,mttdt' '[mtta]' > $1_smx.m
cat mtt_junk | mtt_p2m >> $1_smx.m
rm -f mtt_junk
mv $1_smxx.$lang mtt_junk
lang_header $1 smxx m 'mttx,mttu,mttxx,mttdt' '[mttax]' > $1_smxx.m
cat mtt_junk | mtt_p2m >> $1_smxx.m
rm -f mtt_junk
mv $1_smxtx.$lang mtt_junk
lang_header $1 smxtx m 'mttx,mttu,mttxx,mttdt' '[mttax]' > $1_smxtx.m
cat mtt_junk | mtt_p2m >> $1_smxtx.m
rm -f mtt_junk
fi
# Now invoke the standard error handling.
mtt_error_r ode2smx_lang.log
|