Index: mttroot/mtt/bin/trans/dae2daes_m ================================================================== --- mttroot/mtt/bin/trans/dae2daes_m +++ mttroot/mtt/bin/trans/dae2daes_m @@ -13,10 +13,14 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +# Revision 1.4 1996/08/24 14:27:59 peter +# Global parameters. +# Error handling. +# ## Revision 1.3 1996/08/18 12:00:19 peter ## Unified format of responses. ## ## Revision 1.2 1996/08/16 13:17:57 peter ## Changed default args to include nx+nz @@ -32,51 +36,42 @@ echo Creating $1_daeso.m rm -f dae2daes_m.log rm -f $1_daes.m rm -f $1_daeso.m -if [ "$2" = "" ]; -then - PARAMS='T=[0:0.1:10];x0=zeros(nx+nz,1);dx0=zeros(nx+nz,1);' - echo Using default parameter $PARAMS -else - PARAMS=$2; -fi - -PARAMS="$PARAMS ;" - - -$MATRIX << EOF > dae2daes_m.log 2>mtt_error + +$MATRIX << EOF > dae2daes_m.log 2>mtt_error.txt %Read the parameters $1_numpar; [nx,ny,nu,nz,nyz] = $1_def; + NX = nx+2*nz+nyz; + t=0; %Just in case its in the parameter list - $PARAMS %Defaults if exist('T')==0 T=[0:0.1:10] end; if exist('x0')==0 - x0 = zeros(nx+nz,1); + x0 = zeros(NX,1); end; if exist('dx0')==0 - dx0 = zeros(nx+nz,1); + dx0 = zeros(NX,1); end; [n,m]=size(T); if m>n T=T'; end; -if nx>0 - x = dassl('$1_dae', x0, dx0, T); +if NX>0 + x = dassl("$1_dae", x0, dx0, T); write_matrix([T,x], '$1_daes'); else x = zeros(size(T)); end; @@ -89,22 +84,15 @@ write_matrix([T,y], '$1_daeso'); end; EOF -# Test for errors and print if any -err_length=$(wc -c