10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
# Copyright (c) P.J.Gawthrop, 1991, 1994, 1998
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.7 1998/08/27 06:48:20 peterg
## Changed mttxx to mttx - in implicit integration smx is called once
## with x in the linearisation the same as x in the current state
##
## Revision 1.6 1998/08/26 18:25:45 peterg
## SMX now generates both A and Ax
##
|
>
>
>
|
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
# Copyright (c) P.J.Gawthrop, 1991, 1994, 1998
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.8 1998/08/27 07:04:44 peterg
## Cleared out mttx array ans use mkid instead.
##
## Revision 1.7 1998/08/27 06:48:20 peterg
## Changed mttxx to mttx - in implicit integration smx is called once
## with x in the linearisation the same as x in the current state
##
## Revision 1.6 1998/08/26 18:25:45 peterg
## SMX now generates both A and Ax
##
|
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
ELSE
BEGIN
aa_ij := -mttdt*a_ij;
%% aa_ji := -mttdt*a_ji;
END;
IF (aa_ij NEQ 0) THEN
GENTRAN mtta(i,j) ::=: aa_ij;
mttAAx_i := mttAAx_i + aa_ij*mkid(mttx,j);
%%%mttAAtx(i,1) := mttAAtx(i,1) + aa_ji*mkid(mttx,j);
END;
IF (mttAAx_i NEQ 0) THEN
GENTRAN mttax(i) ::=: mttAAx_i;
END;
|
|
|
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
|
ELSE
BEGIN
aa_ij := -mttdt*a_ij;
%% aa_ji := -mttdt*a_ji;
END;
IF (aa_ij NEQ 0) THEN
GENTRAN mtta(i,j) ::=: aa_ij;
mttAAx_i := mttAAx_i + aa_ij*xj;
%%%mttAAtx(i,1) := mttAAtx(i,1) + aa_ji*mkid(mttx,j);
END;
IF (mttAAx_i NEQ 0) THEN
GENTRAN mttax(i) ::=: mttAAx_i;
END;
|