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.17 2001/07/27 23:29:10 geraint +## Optimises only when requested (-opt). +## ## Revision 1.16 2001/07/13 04:54:04 geraint ## Branch merge: numerical-algebraic-solution back to main. ## ## Revision 1.15.2.4 2001/06/26 00:55:48 geraint ## Writes algebraic equation Jacobian _aej.r (not used yet). @@ -309,11 +312,11 @@ IF (MTTNyz>0) THEN % not $solve (or perhaps solution failed?) BEGIN WRITE "MATRIX MTTyz(",MTTNyz,",1)"; WRITE "%File: $1_ae.r"; FOR i := 1:MTTNyz DO - WRITE "MTTyz(",i,",1) := ",MTTyz(i,1); + WRITE "MTTyz(",i,",1) := ",MTTyz(i,1)," +0"; END; % if MTTNyz>0 (and !$solve) WRITE ";END;"; SHUT "$1_ae.r"; OUT "$1_aej.r"; @@ -324,11 +327,11 @@ 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; + WRITE "MTTyzj(",i,",",j,") := ",didj," +0"; END; END; WRITE ";END;"; SHUT "$1_aej.r";