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.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).
##
## Revision 1.15.2.3 2001/05/09 00:19:22 geraint
|
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
336
|
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
336
337
338
339
|
-
+
-
+
|
OUT "$1_ae.r";
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";
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;
WRITE "MTTyzj(",i,",",j,") := ",didj," +0";
END;
END;
WRITE ";END;";
SHUT "$1_aej.r";
% Create the matrix declarations
OUT "$1_cse.r1";
|