Index: mttroot/mtt/bin/trans/make_ode2odes ================================================================== --- mttroot/mtt/bin/trans/make_ode2odes +++ mttroot/mtt/bin/trans/make_ode2odes @@ -7,10 +7,13 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.86 2004/08/29 00:19:49 geraint +## Defaults to noAlgebraicSolver. +## ## Revision 1.85 2003/06/25 12:46:06 gawthrop ## Input only changed one per print interval ## No effect if stepfactor=1 ## Fixes bug when _input.m is compiled using -stdin option ## and stepfactor>1 @@ -347,10 +350,16 @@ if [ -n "$4" ]; then algebraic_solver=$4 else algebraic_solver="noAlgebraicSolver" fi + +if [ -n "$5" ]; then + sorted_equations=$5 +else + sorted_equations="no" +fi insertor=\<\< # help emacs sh-mode handle C++ lines echo Creating $filename with $method integration method @@ -368,16 +377,28 @@ "dassl") ode=ode odeo=odeo algorithm="mtt_dassl(x,u,t,par,dx,ddt,MTTNX,MTTNYZ,open_switches)" ;; - "sorted_euler") + "sorted_euler") # used by sese generated from seqn algorithm="mtt_euler(x,dx,ddt,$Nx,open_switches)" ;; "euler" | "rk4" | *) - ode=ode - odeo=odeo + case "$sorted_equations" in + "make") # used by sese generated by make + 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 make_m()