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

[GMRES2stage.git] / code / ex15.c

diff --git a/code/ex15.c b/code/ex15.c
index 0e0a8c467b785a3ebbaae3f4135057fa2abf1390..94c93659b6acb2e0ca63c874f44ed14e79f649c4 100644 (file)
--- a/code/ex15.c
+++ b/code/ex15.c
@@ -1,5 +1,5 @@
 
-// /home/couturie/work/petsc-3.5.1/arch-linux2-c-debug/bin/mpirun -np 4    ./ex15 -m 400 -n 400
+// /home/couturie/work/petsc-3.5.1/arch-linux2-c-debug/bin/mpirun -np 4    ./ex15 -m 400 -n 400 -ksp_type fgmres
 
 
 
@@ -232,7 +232,7 @@ int KrylovMinimizeLSQR(Mat A, Vec b, Vec x) {
 
   PetscScalar  alpha, beta;
   PetscReal norm=20, Eprecision=1e-8, tol=1e-40;     
-  PetscInt giter=0, ColS=12, col=0, Emaxiter=50000000, iter=0, iterations=15, Iiter=0;
+  PetscInt giter=0, ColS=12, col=0, Emaxiter=50000000, iter=0, iterations=20, Iiter=0;
   PetscErrorCode ierr;
   PetscScalar T1, T2;
   KSP ksp;
@@ -419,14 +419,14 @@ int KrylovMinimizeLSQR(Mat A, Vec b, Vec x) {
         VecAYPX(d,-thet,v);       //d = (v - thet * d);
         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;
@@ -630,7 +630,7 @@ int main(int argc,char **args)
     VecDuplicate(b,&x2);
     VecDuplicate(b,&sol);
     
-    KrylovMinimize(A, b, x2);
+    KrylovMinimizeLSQR(A, b, x2);