cleanup

src/problem/Formulation.h

@@ -95,7 +95,7 @@ namespace gmshddm
       friend int MatVectProduct< T_Scalar >(Mat A, Vec X, Vec Y);
       friend int MatVectProductOverride< T_Scalar >(Mat A, Vec X, Vec Y, bool IA);
 
-      //*
+      
       // T_Scalar, PetscScalar, T_Interger
       template< class T_Scalar1, class T_Scalar2, class T_Integer >
       struct PetscInterface {
@@ -115,26 +115,7 @@ namespace gmshddm
           delete[] array;
         }
       };
-      //*/
-      /*
-      template< class T_Scalar1, class T_Integer >
-      struct PetscInterface< T_Scalar1, T_Scalar1, T_Integer > {
-        static int kspPrint(KSP ksp, T_Integer it, gmshfem::scalar::Precision< T_Scalar1 > rnorm, void *mctx)
-        {
-          Formulation< T_Scalar1 > *formulation = static_cast< Formulation< T_Scalar1 > * >(mctx);
-          //formulation->_runIterationOperations(it, rnorm);
-          return 0;
-        }
       
-        static const T_Scalar1 *arrayInterface(const T_Scalar1 *array, const T_Integer size)
-        {
-          return array;
-        }
-
-        static void freeArray(const T_Scalar1 *array)
-        {
-        }
-      };*/
     };
 
 

src/problem/Solvers.cpp

@@ -53,7 +53,6 @@ namespace gmshddm
     {
       _absoluteResidual.push_back(rnorm);
       _relativeResidual.push_back(rnorm / _absoluteResidual[0]);
-      //_numberOfIterations = iteration;
 
       const unsigned int MPI_Rank = mpi::getMPIRank();
       if(MPI_Rank == 0) {
@@ -125,7 +124,6 @@ namespace gmshddm
             solverMod.insert(*it);
         }
       }
-      gmshfem::msg::print << "solverBase: " << solverBase << "solverName : " << _solverName << gmshfem::msg::endl;
       KSP ksp;
       PC pc;
       KSPCreate(PETSC_COMM_WORLD, &ksp);