Index: mttroot/mtt/bin/trans/p/mtt_sparse.p ================================================================== --- mttroot/mtt/bin/trans/p/mtt_sparse.p +++ mttroot/mtt/bin/trans/p/mtt_sparse.p @@ -1,15 +1,18 @@ -PROCEDURE mtt_sparse( b : StateVector; -n : integer; - VAR x : StateVector; +PROCEDURE mtt_sparse( b : glnarray; + n : integer; + VAR x : glnarray; VAR rsq : real); {* ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.1 1998/08/13 14:58:35 peterg +## Initial revision +## ## Revision 1.1 1998/08/13 08:40:40 peterg ## Initial revision ## ############################################################### *} @@ -55,11 +58,11 @@ anum := anum+g[j]*h[j]; aden := aden+sqr(xi[j]) END; IF (aden = 0.0) THEN BEGIN writeln('pause in routine SPARSE'); - writeln('very singular matrix'); readln + writeln('very singular matrix'); {***readln ***} END; anum := anum/aden; FOR j := 1 TO n DO BEGIN xi[j] := x[j]; x[j] := x[j]+anum*h[j] @@ -91,8 +94,10 @@ FOR j := 1 TO n DO BEGIN g[j] := -xi[j]; h[j] := g[j]+gam*h[j] END END; - writeln('pause in routine SPARSE'); - writeln('too many iterations'); readln; -99: END; + {***writeln('pause in routine SPARSE'); + writeln('too many iterations'); readln; **} + +99:{**writeln("---",iter,rsq);**} +END;