GMRES.hh

Go to the documentation of this file.

#pragma once
 
// external libraries
#include <deal.II/lac/precondition.h>
#include <deal.II/lac/solver_gmres.h>
 
// DiFfRG
#include <DiFfRG/timestepping/linear_solver/abstract_linear_solver.hh>
 
namespace DiFfRG
{
  template <typename SparseMatrixType, typename VectorType>
  class GMRES : public AbstractLinearSolver<SparseMatrixType, VectorType>
  {
  public:
    GMRES() : matrix(nullptr) {}
 
    void init(const SparseMatrixType &matrix) { this->matrix = &matrix; }
 
    bool invert() { return false; }
 
    int solve(const VectorType &src, VectorType &dst, const double tol)
    {
      if (!matrix) throw std::runtime_error("GMRES::solve: matrix not initialized");
      dealii::SolverControl solver_control(std::max<std::size_t>(1000, src.size() / 10), tol);
      dealii::SolverGMRES<VectorType> solver(solver_control);
 
      preconditioner.initialize(*matrix, 1.0);
      solver.solve(*matrix, dst, src, preconditioner);
 
      return solver_control.last_step();
    }
 
  private:
    const SparseMatrixType *matrix;
    dealii::PreconditionJacobi<SparseMatrixType> preconditioner;
  };
} // namespace DiFfRG