modif de LSQR qui fonctionne quand meme mieux (un poil moins rapide que CGLS quand...

[GMRES2stage.git] / code / ex45.c

diff --git a/code/ex45.c b/code/ex45.c
index a44a2d09174416a6ae6d2a08056a78d5770b4126..11fd5c28340202ee669d6fea574e220a1d8a771d 100644 (file)
--- a/code/ex45.c
+++ b/code/ex45.c
@@ -408,14 +408,6 @@ int KrylovMinimizeLSQR(Mat A, Vec b, Vec x) {
         VecAYPX(d,1/rho,zero_short);     //d=d/ rho;
 
 
-        if (normar/(norma*normr) <= tol) { // check for convergence in min{|b-A*x|}
-          break;
-        }
-        if (normr <= tolb) {           // check for convergence in A*x=b
-          break;
-        }
-
-
         VecAXPY(x_lsqr,phi,d);     // x_lsqr=x_lsqr+phi*d
         normr = abs(s) * normr;
         MatMultTranspose(AS, u, vt); //vt=A'*u;