10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
# Copyright (c) P.J.Gawthrop, 1991, 1994.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.6 2000/04/06 11:27:56 peterg
## Removed the in _cr.r
##
## Revision 1.5 1999/12/08 02:04:46 peterg
## Removed bug - uj := MTTU(j,1); not commented out
##
## Revision 1.4 1999/11/22 23:49:50 peterg
|
>
>
>
|
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
# Copyright (c) P.J.Gawthrop, 1991, 1994.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.7 2000/04/07 19:11:59 peterg
## reduce --> $SYMBOLIC
##
## Revision 1.6 2000/04/06 11:27:56 peterg
## Removed the in _cr.r
##
## Revision 1.5 1999/12/08 02:04:46 peterg
## Removed bug - uj := MTTU(j,1); not commented out
##
## Revision 1.4 1999/11/22 23:49:50 peterg
|
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
OFF Echo;
OFF Nat;
% Get rid of the old mttx and u - now use mkid instead
clear MTTx, MTTu;
% Find MTTA : the A matrix
matrix MTTA(MTTNx,MTTNx);
FOR j := 1:MTTNx DO
BEGIN
%xj := MTTX(j,1);
xj := mkid(MTTx,j);
FOR i := 1:MTTNx DO
MTTA(i,j) := df(MTTEdx(i,1), xj, 1);
END;
% Find MTTB : the B matrix
matrix MTTB(MTTNx,MTTNu);
FOR j := 1:MTTNu DO
BEGIN
%uj := MTTU(j,1);
uj := mkid(MTTu,j);
FOR i := 1:MTTNx DO
MTTB(i,j) := df(MTTEdx(i,1), uj, 1);
END;
% Find MTTC : the C matrix
matrix MTTC(MTTNy,MTTNx);
FOR i := 1:MTTNy DO
FOR j := 1:MTTNx DO
BEGIN
%xj := MTTX(j,1);
xj := mkid(MTTx,j);
MTTC(i,j) := df(MTTY(i,1), xj, 1);
END;
% Find MTTD : the D matrix
matrix MTTD(MTTNy,MTTNu);
FOR i := 1:MTTNy DO
FOR j := 1:MTTNu DO
BEGIN
%uj := MTTU(j,1);
uj := mkid(MTTu,j);
MTTD(i,j) := df(MTTY(i,1), uj, 1);
END;
%Substitute the ss values
in "$1_sspar.r";
%Create the output file
OUT "$1_$2.r";
|
>
>
>
>
>
>
>
>
>
>
>
>
>
|
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
|
OFF Echo;
OFF Nat;
% Get rid of the old mttx and u - now use mkid instead
clear MTTx, MTTu;
IF (MTTNx>0) THEN
BEGIN
% Find MTTA : the A matrix
matrix MTTA(MTTNx,MTTNx);
FOR j := 1:MTTNx DO
BEGIN
%xj := MTTX(j,1);
xj := mkid(MTTx,j);
FOR i := 1:MTTNx DO
MTTA(i,j) := df(MTTEdx(i,1), xj, 1);
END;
END;
IF (MTTNx>0) AND (MTTNu>0) THEN
BEGIN
% Find MTTB : the B matrix
matrix MTTB(MTTNx,MTTNu);
FOR j := 1:MTTNu DO
BEGIN
%uj := MTTU(j,1);
uj := mkid(MTTu,j);
FOR i := 1:MTTNx DO
MTTB(i,j) := df(MTTEdx(i,1), uj, 1);
END;
END;
IF (MTTNx>0) AND (MTTNy>0) THEN
BEGIN
% Find MTTC : the C matrix
matrix MTTC(MTTNy,MTTNx);
FOR i := 1:MTTNy DO
FOR j := 1:MTTNx DO
BEGIN
%xj := MTTX(j,1);
xj := mkid(MTTx,j);
MTTC(i,j) := df(MTTY(i,1), xj, 1);
END;
END;
IF (MTTNy>0) AND (MTTNu>0) THEN
BEGIN
% Find MTTD : the D matrix
matrix MTTD(MTTNy,MTTNu);
FOR i := 1:MTTNy DO
FOR j := 1:MTTNu DO
BEGIN
%uj := MTTU(j,1);
uj := mkid(MTTu,j);
MTTD(i,j) := df(MTTY(i,1), uj, 1);
END;
END;
%Substitute the ss values
in "$1_sspar.r";
%Create the output file
OUT "$1_$2.r";
|