allow other Krylov solvers (bcgs, bcgsl, ...)

src/problem/Formulation.cpp

@@ -234,34 +234,22 @@ namespace gmshddm
         } while(_relativeResidual.back() > tolerance && i < iterMax);
         VecDestroy(&W);
       }
-      else if(solver == "gmres") { // gmres
+      else { // gmres, bcgs, bcgsl, ...
         _IA = true;
         KSP ksp;
         PC pc;
         KSPCreate(PETSC_COMM_SELF, &ksp);
         KSPSetFromOptions(ksp);
         KSPSetOperators(ksp, A, A);
-        KSPSetType(ksp, "gmres");
+        KSPSetType(ksp, solver.c_str());
         KSPMonitorSet(ksp, kspPrint< T_Scalar >, this, PETSC_NULL);
         KSPSetTolerances(ksp, tolerance, PETSC_DEFAULT, PETSC_DEFAULT, iterMax);
         KSPGetPC(ksp, &pc);
         PCSetType(pc, PCNONE);
+        if(solver == "gmres")
           KSPGMRESSetRestart(ksp, 10000);
-        KSPSolve(ksp, _rhs.getPetsc(), Y);
-      }
-      else if(solver == "bicgl") { // bicgl
-        _IA = true;
-        KSP ksp;
-        PC pc;
-        KSPCreate(PETSC_COMM_SELF, &ksp);
-        KSPSetFromOptions(ksp);
-        KSPSetOperators(ksp, A, A);
-        KSPSetType(ksp, "bcgsl");
-        KSPMonitorSet(ksp, kspPrint< T_Scalar >, this, PETSC_NULL);
-        KSPSetTolerances(ksp, tolerance, PETSC_DEFAULT, PETSC_DEFAULT, iterMax);
-        KSPGetPC(ksp, &pc);
-        PCSetType(pc, PCNONE);
-        // KSPBCGSLSetEll(ksp, 2);
+        else if(solver == "bcgsl")
+          KSPBCGSLSetEll(ksp, 6);
         KSPSolve(ksp, _rhs.getPetsc(), Y);
       }