Index: mttroot/mtt/bin/trans/make_ode2odes ================================================================== --- mttroot/mtt/bin/trans/make_ode2odes +++ mttroot/mtt/bin/trans/make_ode2odes @@ -7,10 +7,14 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.28 1999/12/08 05:56:52 peterg +## Reordered the writing of the input and output. +## Note that last value now discarded. +## ## Revision 1.27 1999/11/15 22:47:53 peterg ## Generates method-specific code. ## ## Revision 1.26 1999/10/20 01:31:43 peterg ## *** empty log message *** @@ -148,11 +152,11 @@ fi cat << EOF >> $1_ode2odes.m for MTTit = 1:MTTilast #Integration loop - [MTTy] = $1_$odeo(MTTx,MTTu,MTTt); # Output + [MTTy] = $1_$odeo(MTTx,MTTu,MTTt,MTTpar); # Output [MTTu] = $1_input(MTTt,MTTx,MTTy); # Input mtt_write(MTTt,MTTx,MTTy,$Nx,$Ny); # Write it out if $Nx>0 # Dont if no states EOF @@ -159,11 +163,11 @@ if [ "$method" = "euler" ]; then cat << EOF >> $1_ode2odes.m # if mttmethod==1 # Euler for MTTjt = 1:mttstepfactor - [MTTdx] = $1_$ode(MTTx,MTTu,MTTt); # State derivative + [MTTdx] = $1_$ode(MTTx,MTTu,MTTt,MTTpar); # State derivative [MTTopen] = $1_switchopen(MTTx); # Open switches [MTTx] = mtt_euler(MTTx,MTTdx,MTTddt,$Nx,MTTopen); # Euler update MTTt = MTTt + MTTddt; endfor; # endif;