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.30 2000/04/18 11:24:19 peterg +## Removed _numpar. +## ## Revision 1.29 2000/04/07 19:10:57 peterg ## *** empty log message *** ## ## Revision 1.28 1999/12/08 05:56:52 peterg ## Reordered the writing of the input and output. @@ -134,12 +137,12 @@ #Initialise MTTt = 0.0; [MTTu] = zero_input($Nu); # Zero the input - [MTTx] = $1_state; # Read in initial state +[MTTpar] = $1_numpar; # Read in initial state #[mttAA] = zero_matrix($Nx); # Zero the A matrix #[mttAAx] = zero_vector($Nx); # Zero the AAx vector #if $Nx>0 @@ -179,11 +182,12 @@ if [ "$method" = "implicit" ]; then cat << EOF >> $1_ode2odes.m # if mttmethod==2 # Implicit [MTTdx] = $1_cse(MTTx,MTTu,MTTt,MTTpar); # State derivative - [mttAA,mttAAx] = $1_smx(MTTx,MTTu,mttdt,MTTpar); # (I-Adt) and (I-Adt)x + [mttAA] = $1_smxa(MTTx,MTTu,mttdt,MTTpar); # (I-Adt) and (I-Adt)x + [mttAAx] = $1_smxax(MTTx,MTTu,mttdt,MTTpar); # (I-Adt) and (I-Adt)x [MTTopen] = $1_switchopen(MTTx); # Open switches [MTTx] = mtt_implicit(MTTx,MTTdx,mttAA,mttAAx,mttdt,$Nx,MTTopen); # Implicit update MTTt = MTTt + mttdt; # endif; EOF