diff --git a/Solver/linearSystemGMM.h b/Solver/linearSystemGMM.h
index a611822877d30748edbfe6cada8e5d9eac67663e..40d7a4ba6ba4c49a60c4f3792ce155280de3b2dd 100644
--- a/Solver/linearSystemGMM.h
+++ b/Solver/linearSystemGMM.h
@@ -19,9 +19,9 @@ template <class scalar>
 class linearSystemGmm : public linearSystem<scalar> {
  protected:
   std::vector<scalar> *_x; // the nonLinearSystemGmm has to access to this vector
+  std::vector<scalar> *_b; // idem
  private:
   gmm::row_matrix<gmm::wsvector<scalar> > *_a;
-  std::vector<scalar> *_b;
   double _prec;
   int _noisy, _gmres;
  public:
diff --git a/Solver/linearSystemPETSc.hpp b/Solver/linearSystemPETSc.hpp
index dc66e630410d5f8b6559ef603f243cbfd394b304..2694650dd6c7597bec9607ac1e0ce5f48ee9ae07 100644
--- a/Solver/linearSystemPETSc.hpp
+++ b/Solver/linearSystemPETSc.hpp
@@ -2,13 +2,13 @@
 #include <petsc.h>
 #include <petscksp.h>
 #include "linearSystemPETSc.h"
-static void  _try(int ierr) 
+static void  _try(int ierr)
 {
   CHKERRABORT(PETSC_COMM_WORLD, ierr);
 }
 
 template <class scalar>
-void linearSystemPETSc<scalar>::_kspCreate() 
+void linearSystemPETSc<scalar>::_kspCreate()
 {
   _try(KSPCreate(PETSC_COMM_WORLD, &_ksp));
   PC pc;
@@ -120,7 +120,7 @@ void linearSystemPETSc<scalar>::allocate(int nbRows)
 }
 
 template <class scalar>
-void linearSystemPETSc<scalar>::print() 
+void linearSystemPETSc<scalar>::print()
 {
   _try(MatAssemblyBegin(_a, MAT_FINAL_ASSEMBLY));
   _try(MatAssemblyEnd(_a, MAT_FINAL_ASSEMBLY));
@@ -254,7 +254,7 @@ int linearSystemPETSc<scalar>::systemSolve()
   //Msg::Info("%d iterations", its);
   return 1;
 }
-  
+
 template <class scalar>
 std::vector<scalar> linearSystemPETSc<scalar>::getData()
 {