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.37 2000/08/01 12:25:06 peterg +## Now includes euler +## ## Revision 1.36 2000/05/19 17:48:16 peterg ## Argument to state ## ## Revision 1.35 2000/05/18 18:59:40 peterg ## Removed the First time stuff @@ -136,11 +139,11 @@ Nx=`mtt_getsize $Sys x` # States Nu=`mtt_getsize $Sys u` # Inputs Ny=`mtt_getsize $Sys y` # Inputs if [ "$method" = "implicit" ]; then - ode=cse + ode=csex odeo=cseo algorithm="mtt_implicit(x,dx,AA,AAx,ddt,$Nx,open_switches)" else ode=ode odeo=odeo @@ -174,12 +177,16 @@ ## Initialise t = 0.0; ddt = simpar.dt/simpar.stepfactor; ilast = round(simpar.last/ddt)+1; # Total number of steps -## Following remove due to p2c bug -## [u] = zero_input(1); # Zero the input +## Following removed due to p2c bug +## [u] = zero_input($Nu); # Zero the input + for MTTi=1:$Nu + u(MTTi) = 0; + endfor; + mttj = 0; for it = 1:ilast #Integration loop [y] = $1_$odeo(x,u,t,par);# Output [u] = $1_input(t,x,y); # Input if mttj==0