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.50 2000/05/19 17:46:41 peterg
## New version of state with par argument
##
## Revision 1.49 2000/05/18 09:45:45 peterg
## Fixed missing ;
##
## Revision 1.48 2000/05/18 09:39:38 peterg
## Removed fifth argument from _input
##
## Revision 1.47 2000/05/16 18:56:03 peterg
|
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
|
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
|
-
+
|
echo " glnarray = StateVector;"
echo " glnpbynp = StateMatrix;"
echo " glmpbynp = StateMatrix;"
## echo " IntegrationMethod = 1..4;"
echo ""
echo "VAR"
echo " simpar : SimulationParameters;"
echo " t : REAL;"
echo " t,ddt : REAL;"
echo " x,dx,AAx : StateVector;"
echo " u : InputVector;"
echo " y : OutputVector;"
echo " par : ParameterVector;"
echo " mttpar : ParameterVector;"
echo " mttnpar : INTEGER;"
echo " AA : StateMatrix;"
|
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
|
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
|
-
+
-
+
+
|
{
if (doing_header==1){
if (rep=="ode2odes"){
printf("\n{%s $MTTPATH/trans/p/mtt_write.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_par_update.p}\n",inc)
printf("{%s $MTTPATH/trans/p/sign.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_euler.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_solve.p}\n",inc)
#printf("{%s $MTTPATH/trans/p/mtt_solve.p}\n",inc)
printf("{%s $MTTPATH/trans/p/mtt_implicit.p}\n",inc)
# printf("{%s $MTTPATH/trans/p/zero_matrix.p}\n",inc)
#printf("{%s $MTTPATH/trans/p/zero_matrix.p}\n",inc)
printf("{%s $MTTPATH/trans/p/zero_input.p}\n",inc)
printf("{%s $MTTPATH/trans/p/zero_state.p}\n",inc)
printf("{%s %s_simpar.p}\n",inc,Sys)
printf("{%s %s_numpar.p}\n",inc,Sys)
# printf("{%s $MTTPATH/trans/p/mtt_getargs.p}\n",inc)
printf("{%s %s_state.p}\n",inc,Sys)
printf("{%s %s_input.p}\n",inc,Sys)
if (Method=="euler") {
printf("{%s %s_ode.p}\n",inc,Sys)
printf("{%s %s_odeo.p}\n",inc,Sys)
}
if (Method=="implicit") {
printf("{%s %s_cse.p}\n",inc,Sys)
printf("{%s %s_cseo.p}\n",inc,Sys)
|