Index: mttroot/mtt/bin/trans/p/mtt_solve_svd.p ================================================================== --- mttroot/mtt/bin/trans/p/mtt_solve_svd.p +++ mttroot/mtt/bin/trans/p/mtt_solve_svd.p @@ -1,10 +1,9 @@ -PROCEDURE mtt_solve_svd(VAR x : StateVector; - A : StateMatrix; - VAR B : StateVector; - n : integer; - Small : real); +PROCEDURE mtt_solve_svd(VAR x : StateVector; + A : StateMatrix; + VAR B : StateVector; + n : integer); {Solve linear equation via svd} { @@ -11,10 +10,14 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % $Id$ % % $Log$ +% % Revision 1.2 1998/08/17 12:40:25 peterg +% % Renamed mtt_solve to mtt_solve_svd +% % -- now also have mtt_solve_lu +% % % % Revision 1.1 1998/08/17 12:38:52 peterg % % Initial revision % % % % Revision 1.2 1998/08/14 12:09:13 peterg % % A passed by value - its destroyed by SVDcm @@ -23,10 +26,12 @@ % % Initial revision % % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% } +CONST + Small = 1e-6; (* Criteroion for small svd *) VAR i : integer; wmax,wmin : real; w : StateVector ; v : StateMatrix;