Overview
Comment:Fixed errors when MTTNx=0.
Downloads: Tarball | ZIP archive | SQL archive
Timelines: family | ancestors | descendants | both | origin/numerical-algebraic-solution | trunk
Files: files | file ages | folders
SHA3-256: fc9f65477ab709300892e452591252621eec6c5c631080f27ff8c901a2594cbe
User & Date: geraint@users.sourceforge.net on 2001-05-05 20:50:16
Other Links: branch diff | manifest | tags
Context
2001-05-09
00:19:22
Fixed EOF error when MTTNYZ=0. check-in: 0c40464f88 user: geraint@users.sourceforge.net tags: origin/numerical-algebraic-solution, trunk
2001-05-05
20:50:16
Fixed errors when MTTNx=0. check-in: fc9f65477a user: geraint@users.sourceforge.net tags: origin/numerical-algebraic-solution, trunk
2001-05-04
04:07:24
Numerical solution of algebraic equations.
sys_ae.cc written for unsolved inputs.
Solution of equations using hybrd from MINPACK (as used by Octave fsolve).
check-in: 2640c1bcbb user: geraint@users.sourceforge.net tags: origin/numerical-algebraic-solution, trunk
Changes

Modified mttroot/mtt/bin/trans/dae2cse_r from [c759129236] to [c5316418a2].

11
12
13
14
15
16
17





18
19
20
21
22
23
24
# Copyright (c) P.J.Gawthrop 1991, 1992, 1994.

###############################################################
## Version control history
###############################################################
## $Id$
## $Log$





## Revision 1.15  2001/03/19 02:28:52  geraint
## Branch merge: merging-ode2odes-exe back to MAIN.
##
## Revision 1.14.2.1  2001/03/19 00:29:08  geraint
## Parse switches (-A) before calling def2write_r.
## Update $1_def.* instead of removing.
##







>
>
>
>
>







11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
# Copyright (c) P.J.Gawthrop 1991, 1992, 1994.

###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## 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).
##
## Revision 1.15  2001/03/19 02:28:52  geraint
## Branch merge: merging-ode2odes-exe back to MAIN.
##
## Revision 1.14.2.1  2001/03/19 00:29:08  geraint
## Parse switches (-A) before calling def2write_r.
## Update $1_def.* instead of removing.
##
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303



304
305

306
307
308


309
310
311
312
313
314
315
316
317
318
319
320
321
322
  MTTNYz := 0;
END; % IF MTTNyz>0 and $solve

IF (MTTNyz>0) THEN % not $solve (or perhaps solution failed?)
BEGIN
    OUT "$1_ae.r";
    WRITE "MATRIX MTTyz(",MTTNyz,",1)";
    WRITE "%file $1_ae.r";
    FOR i := 1:MTTNyz DO
	WRITE "MTTyz(",i,",1) := ",MTTyz(i,1);
    WRITE ";END;";
    SHUT "$1_ae.r";
END; % if MTTNyz>0 (and !$solve)

% Create the matrix declarations
OUT "$1_cse.r1";



write "MATRIX MTTEdx(", MTTNx, ",", 1, ")$"; 
write "MATRIX MTTE(", MTTNx, ",", MTTNx, ")$"; 

SHUT "$1_cse.r1";

OUT "$1_csex.r1";


write "MATRIX MTTEdx(", MTTNx, ",", 1, ")$"; 
SHUT "$1_csex.r1";

IF MTTNy>0 THEN
BEGIN
  OUT "$1_cseo.r1";
  write "MATRIX MTTY(", MTTNy, ",", MTTNx, ")$"; 
  SHUT "$1_cseo.r1";
END;

%%Create the _cse.r file
OUT "$1_cse.r2";
write "%File: $1_cse.r";
in ("$1_cse_write.r");







|








>
>
>
|
|
>



>
>
|





|







293
294
295
296
297
298
299
300
301
302
303
304
305
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
  MTTNYz := 0;
END; % IF MTTNyz>0 and $solve

IF (MTTNyz>0) THEN % not $solve (or perhaps solution failed?)
BEGIN
    OUT "$1_ae.r";
    WRITE "MATRIX MTTyz(",MTTNyz,",1)";
    WRITE "%File: $1_ae.r";
    FOR i := 1:MTTNyz DO
	WRITE "MTTyz(",i,",1) := ",MTTyz(i,1);
    WRITE ";END;";
    SHUT "$1_ae.r";
END; % if MTTNyz>0 (and !$solve)

% Create the matrix declarations
OUT "$1_cse.r1";
write "%";
IF (MTTNx > 0) THEN
BEGIN
    write "MATRIX MTTEdx(", MTTNx, ",", 1, ")$"; 
    write "MATRIX MTTE(", MTTNx, ",", MTTNx, ")$"; 
END;
SHUT "$1_cse.r1";

OUT "$1_csex.r1";
write "%File:$1_csex.r1";
IF (MTTNx > 0) THEN
    write "MATRIX MTTEdx(", MTTNx, ",", 1, ")$"; 
SHUT "$1_csex.r1";

IF MTTNy>0 THEN
BEGIN
  OUT "$1_cseo.r1";
  write "MATRIX MTTY(", MTTNy, ",", 1, ")$"; 
  SHUT "$1_cseo.r1";
END;

%%Create the _cse.r file
OUT "$1_cse.r2";
write "%File: $1_cse.r";
in ("$1_cse_write.r");


MTT: Model Transformation Tools
GitHub | SourceHut | Sourceforge | Fossil RSS ]