11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
+
+
+
|
# Copyright (c) P.J.Gawthrop 1998
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.8 1998/07/26 19:38:17 peterg
## Replaced t0..t9 by mtt_temp0..
##
## Revision 1.7 1998/07/25 20:06:23 peterg
## Does the mtt_update function
##
## Revision 1.6 1998/07/25 16:59:40 peterg
## Give other procedures local i,j vars.
##
## Revision 1.5 1998/07/25 15:06:17 peterg
|
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
|
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
|
-
+
|
if (doing_header==1){
if (rep=="ode2odes"){
printf("VAR ");
for (k=1;k<i;k++) printf("%s,",global[k])
printf("%s : REAL;\n", global[i])
printf("VAR "); for (k=1;k<=Nx;k++) printf("mttx%i,", k)
for (k=1;k<=Nu;k++) printf("mttu%i,", k)
for (k=0;k<Nt;k++) printf("mtt_temp%i,", k)
for (k=1;k<=Nt;k++) printf("mtt_temp%i,", k)
printf("t%i : REAL;\n",Nt)
printf("{%s sign.p}\n",inc,sys)
printf("{%s mtt_update.p}\n",inc,sys)
printf("{%s mtt_write.p}\n",inc,sys)
printf("{%s %s_simpar.p}\n",inc,sys)
printf("{%s %s_numpar.p}\n",inc,sys)
printf("{%s %s_state.p}\n",inc,sys)
|
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
|
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
|
-
-
+
+
|
}
else
if (rep=="ode2odes")
print "END."
else
print "END;"
}' sys=$sys rep=$rep Nx=$Nx Nu=$Nu |\
sed 's/(\([ijk0-9]*\))/\[\1\]/g' |\
sed 's/(\(MTT[ijk0-9]*\))/\[\1\]/g' \
sed 's/(\([ijk0-9,]*\))/\[\1\]/g' |\
sed 's/(\(MTT[ijk0-9,]*\))/\[\1\]/g' \
>> $Filename
# p2c doesn't like mixed case filenames!
if [ "$Filename" != "$filename" ]; then
echo Creating $filename
cp -f $Filename $filename
fi
|