1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
+
+
+
|
#! /bin/sh
######################################
##### Model Transformation Tools #####
######################################
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.88 2004/08/29 13:15:28 geraint
## Uses sys_sae instead of sys_ae if sorted equations are being used.
##
## Revision 1.87 2004/08/29 01:46:56 geraint
## Added rules to create ode2odes for sorted system: sesx and sesy.
##
## Revision 1.86 2004/08/29 00:19:49 geraint
## Defaults to noAlgebraicSolver.
##
## Revision 1.85 2003/06/25 12:46:06 gawthrop
|
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
|
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
|
+
+
+
+
+
|
# Find system constants
Nx=`mtt_getsize $sys x` # States
Nu=`mtt_getsize $sys u` # Inputs
Ny=`mtt_getsize $sys y` # Outputs
case "$method" in
"implicit")
ae=ae
ode=csex
odeo=cseo
algorithm="mtt_implicit(x,dx,AA,AAx,ddt,$Nx,open_switches)"
;;
"dassl")
case "$sorted_equations" in
"make") # used by sese generated by make
ae=sae
ode=sesx
odeo=sesy
;;
"seqn") # shouldn't be here unless mtt has changed
ae=ae
ode=sese
odeo=sese
;;
"no" | *)
ae=ae
ode=ode
odeo=odeo
;;
esac
algorithm="mtt_dassl(x,u,t,par,dx,ddt,MTTNX,MTTNYZ,open_switches)"
;;
"sorted_euler") # used by sese generated from seqn
algorithm="mtt_euler(x,dx,ddt,$Nx,open_switches)"
;;
"euler" | "rk4" | *)
case "$sorted_equations" in
"make") # used by sese generated by make
ae=sae
ode=sesx
odeo=sesy
;;
"seqn") # shouldn't be here unless mtt has changed
ae=ae
ode=sese
odeo=sese
;;
"no" | *)
ae=ae
ode=ode
odeo=odeo
;;
esac
algorithm="mtt_euler(x,dx,ddt,$Nx,open_switches)"
;;
esac
|