10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
+
+
+
|
# Copyright (c) P.J.Gawthrop 1991, 1994, 1995, 1996
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.2 1996/08/16 08:57:34 peter
## Reads $1_input.
##
## Revision 1.1 1996/08/15 16:46:06 peter
## Initial revision
##
###############################################################
# WARNING this is experimental!!
|
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
|
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
|
-
+
|
MTTdx := MTTdx;
MTTGx := MTTGx;
OFF fort;
write "for i=1:", MTTNy, ";;";
write "for i=1:", MTTNx, ";;";
write " residual(i) = MTTdX(i)-MTTdX0(i);";
write "end;";
IF MTTNz>0 THEN
BEGIN
write "dxz = MTTGx*MTTdx";
write "for i=1:", MTTNz, ";;";
|