diff --git a/src/problem/Formulation.h b/src/problem/Formulation.h index 2828273f94872aedaeaee89cc71df66f58ac0225..149ce6281cf9d9e152d01ebd6cbc41181913a3ee 100644 --- a/src/problem/Formulation.h +++ b/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) - { - } - };*/ + }; diff --git a/src/problem/Solvers.cpp b/src/problem/Solvers.cpp index 1ebfd843616eed42f7218599ad3370248d5c3ce3..0736ce95b27e8737fca8808567a99eed3dd69b70 100644 --- a/src/problem/Solvers.cpp +++ b/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);