LinearSystem.cpp

Go to the documentation of this file.

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   This file is part of CEPS.
 
   CEPS is free software: you can redistribute it and/or modify
   it under the terms of the GNU General Public License as published by
   the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.
 
   CEPS is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License
   along with CEPS (see file LICENSE at root of project).
   If not, see <https://www.gnu.org/licenses/>.
 
 
   Copyright 2019-2024 Inria, Universite de Bordeaux
 
   Authors, in alphabetical order:
 
     Pierre-Elliott BECUE, Florian CARO, Yves COUDIERE(*), Andjela DAVIDOVIC, 
     Charlie DOUANLA-LONTSI, Marc FUENTES, Mehdi JUHOOR, Michael LEGUEBE(*), 
     Pauline MIGERDITICHAN, Valentin PANNETIER(*), Nejib ZEMZEMI.
     * : currently active authors
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
#include "linearAlgebra/LinearSystem.hpp"
 
#ifdef CEPS_USE_PETSC
// wrapper
#include "linearAlgebra/PETScKSP.cpp"
#endif // CEPS_USE_PETSC
 
void
LinearSystem::setupWithParameters(InputParameters* params)
{
  // Setup the Linear Solver properties
  setType(params->getString("linear solver type", "GMRES"));
  setTolerances(params->getReal("linear solver relative tolerance"  ,m_relativeTolerance), 
                params->getReal("linear solver absolute tolerance"  ,m_absoluteTolerance),
                params->getReal("linear solver divergence tolerance",m_divergenceTolerance),
                params->getInt ("linear solver max iterations"      ,m_maxIterations));
}
 
CepsSolver&
LinearSystem::getSolver()
{
  return m_solver;
}
 
DMatPtr
LinearSystem::getLhsMatrix()
{
  return m_lhsMatrix;
}
 
DVecPtr 
LinearSystem::getRhsVector()
{
  return m_rhsVector;
}
 
CepsReal
LinearSystem::getRelativeTolerance() const
{
  return m_relativeTolerance;
}
 
CepsReal
LinearSystem::getAbsoluteTolerance() const
{
  return m_absoluteTolerance;
}
 
CepsReal
LinearSystem::getDivergenceTolerance() const
{
  return m_divergenceTolerance;
}