Index: mttroot/mtt/bin/trans/dae2cse_r ================================================================== --- mttroot/mtt/bin/trans/dae2cse_r +++ mttroot/mtt/bin/trans/dae2cse_r @@ -13,10 +13,13 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.15.2.3 2001/05/09 00:19:22 geraint +## Fixed EOF error when MTTNYZ=0. +## ## Revision 1.15.2.2 2001/05/05 20:50:16 geraint ## Fixed errors when MTTNx=0. ## ## Revision 1.15.2.1 2001/05/04 04:07:24 geraint ## Numerical solution of algebraic equations. @@ -305,10 +308,26 @@ WRITE "MTTyz(",i,",1) := ",MTTyz(i,1); END; % if MTTNyz>0 (and !$solve) WRITE ";END;"; SHUT "$1_ae.r"; +OUT "$1_aej.r"; +IF (MTTNyz>0) THEN % as above +BEGIN + WRITE "MATRIX MTTyzj(",MTTNyz,",",MTTNyz,")"; + WRITE "%File: $1_aej.r"; + FOR i := 1:MTTNyz DO + FOR j := 1:MTTNyz DO + BEGIN + didj := df(MTTyz(i,1),mkid('mttui,j)); + IF (didj NEQ 0) THEN + WRITE "MTTyzj(",i,",",j,") := ",didj; + END; +END; +WRITE ";END;"; +SHUT "$1_aej.r"; + % Create the matrix declarations OUT "$1_cse.r1"; write "%"; IF (MTTNx > 0) THEN BEGIN