Overview
| Comment: | Writes algebraic equation Jacobian _aej.r (not used yet). |
|---|---|
| Downloads: | Tarball | ZIP archive | SQL archive |
| Timelines: | family | ancestors | descendants | both | origin/numerical-algebraic-solution | trunk |
| Files: | files | file ages | folders |
| SHA3-256: |
f057909af3127b388e68ecd5a2f09a83 |
| User & Date: | geraint@users.sourceforge.net on 2001-06-26 00:55:48 |
| Other Links: | branch diff | manifest | tags |
Context
|
2001-06-26
| ||
| 00:57:21 | Prints more useful name. check-in: 528fe3d379 user: geraint@users.sourceforge.net tags: origin/numerical-algebraic-solution, trunk | |
| 00:55:48 | Writes algebraic equation Jacobian _aej.r (not used yet). check-in: f057909af3 user: geraint@users.sourceforge.net tags: origin/numerical-algebraic-solution, trunk | |
|
2001-06-25
| ||
| 23:28:29 | Generic mtt_rate and mtt_output - allows method independent calls. check-in: 43a0e0bde2 user: geraint@users.sourceforge.net tags: origin/numerical-algebraic-solution, trunk | |
Changes
Modified mttroot/mtt/bin/trans/dae2cse_r from [7e2f3c3197] to [75e9fca884].
| ︙ | |||
11 12 13 14 15 16 17 18 19 20 21 22 23 24 | 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 | + + + | # Copyright (c) P.J.Gawthrop 1991, 1992, 1994. ############################################################### ## 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. ## sys_ae.cc written for unsolved inputs. ## Solution of equations using hybrd from MINPACK (as used by Octave fsolve). |
| ︙ | |||
303 304 305 306 307 308 309 310 311 312 313 314 315 316 | 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 | + + + + + + + + + + + + + + + + |
WRITE "%File: $1_ae.r";
FOR i := 1:MTTNyz DO
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
write "MATRIX MTTEdx(", MTTNx, ",", 1, ")$";
write "MATRIX MTTE(", MTTNx, ",", MTTNx, ")$";
|
| ︙ |