1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
-
-
-
-
+
+
+
+
+
+
+
+
-
-
+
+
-
-
+
+
|
PROCEDURE mtt_solve(VAR x : StateVector;
A : StateMatrix;
VAR B : StateVector;
n : integer);
PROCEDURE mtt_solve(VAR x : StateVector;
A : StateMatrix;
B : StateVector;
n : integer);
{
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % Version control history
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% % $Id$
% % $Log$
% % Revision 1.3 1999/10/26 23:44:36 peterg
% % Now packages up mtt_solve_lud OR mtt_solve_svd
% %
% % Revision 1.2 1998/08/14 12:09:13 peterg
% % A passed by value - its destroyed by SVDcm
% %
% % Revision 1.1 1998/08/13 08:51:57 peterg
% % Initial revision
% %
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
}
VAR
i : integer;
wmax,wmin : real;
w : StateVector ;
v : StateMatrix;
(*$I mtt_solve_lud.p *)
(* I mtt_solve_svd.p *)
(*$I $MTTPATH/trans/p/mtt_solve_lud.p *)
(* I $MTTPATH/trans/p/mtt_solve_svd.p *)
BEGIN{mtt_solve}
mtt_solve_lud(xsub,AAsub,BB,Nx);
(*** mtt_solve_svd(xsub,AAsub,BB,Nx); ***)
mtt_solve_lud(x,A,B,N);
(*** mtt_solve_svd(x,A,B,N); ***)
END{mtt_solve};
|