Index: mttroot/mtt/bin/trans/p/svbksb.p ================================================================== --- mttroot/mtt/bin/trans/p/svbksb.p +++ mttroot/mtt/bin/trans/p/svbksb.p @@ -1,14 +1,18 @@ PROCEDURE svbksb(u: glmpbynp; w: glnparray; v: glnpbynp; - m,n,mp,np: integer; b: glmparray; VAR x: glnparray); + m,n: integer; b: glmparray; VAR x: glnparray); { % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % $Id$ % % $Log$ +% % Revision 1.2 1998/08/12 11:09:54 peterg +% % Taken from NR share library nrpas13 as SVBKSB.PAS +% % and renamed svbksb.p +% % % % Revision 1.1 1998/08/12 11:09:02 peterg % % Initial revision % % % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% } Index: mttroot/mtt/bin/trans/p/svdcmp.p ================================================================== --- mttroot/mtt/bin/trans/p/svdcmp.p +++ mttroot/mtt/bin/trans/p/svdcmp.p @@ -1,15 +1,20 @@ -PROCEDURE svdcmp(VAR a: glmpbynp; m,n,mp,np: integer; +PROCEDURE svdcmp(VAR a: glmpbynp; m,n: integer; VAR w: glnparray; VAR v: glnpbynp); { % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % $Id$ % % $Log$ +% % Revision 1.3 1998/08/12 11:08:03 peterg +% % Put in pythag routine to compute z = sqrt(y^2 + z^2) (as in book) +% % Save loop index l as ll when jumping from loop - using l itself is +% % undefined +% % % % Revision 1.2 1998/08/12 11:05:33 peterg % % Taken from NR share library nrpas13 as SVDCMP.PAS % % and renamed svdcmp.p % % % % Revision 1.1 1998/08/12 11:03:57 peterg