Overview
Comment: | Fixed ae for non-sorted code. |
---|---|
Downloads: | Tarball | ZIP archive | SQL archive |
Timelines: | family | ancestors | descendants | both | origin/master | trunk |
Files: | files | file ages | folders |
SHA3-256: |
d420d498e123342890194091d602bba3 |
User & Date: | geraint@users.sourceforge.net on 2004-08-29 16:04:44 |
Other Links: | branch diff | manifest | tags |
Context
2004-08-29
| ||
16:09:49 |
Added option -no-reduce
Code sorting can now be used with Reduce (via sese.r) so that existing CRs | |
16:04:44 | Fixed ae for non-sorted code. check-in: d420d498e1 user: geraint@users.sourceforge.net tags: origin/master, trunk | |
14:39:23 | Added rule to create dae from sorted equations (sese) if sorting is requested. check-in: f98d517043 user: geraint@users.sourceforge.net tags: origin/master, trunk | |
Changes
Modified mttroot/mtt/bin/trans/make_ode2odes from [674a5c4a87] to [b57eec8f4f].
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | #! /bin/sh ###################################### ##### Model Transformation Tools ##### ###################################### ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ ## 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 | > > > | 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 | # 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") 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 ode=sese odeo=sese ;; "no" | *) 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 ode=sese odeo=sese ;; "no" | *) ode=ode odeo=odeo ;; esac algorithm="mtt_euler(x,dx,ddt,$Nx,open_switches)" ;; esac | > > > > > | 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 |
︙ | ︙ |