CD_KrylovGMRES.H

Go to the documentation of this file.

/*
 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 */
 
#ifndef CD_KRYLOVGMRES_H
#define CD_KRYLOVGMRES_H
 
// Chombo includes
#include <LevelData.H>
#include <Vector.H>
 
// Our includes
#include <CD_AMRMultigridKrylovOp.H>
#include <CD_NamespaceHeader.H>
 
template <class T>
class KrylovGMRES
{
public:
  Real m_eps = 1.E-10;
 
  int m_maxIter = 50;
 
  int m_restart = 30;
 
  bool m_reorthogonalize = true;
 
  int m_verbosity = -1;
 
  int m_exitStatus = -1;
 
  Real m_residualNorm = -1.0;
 
  int m_iterations = 0;
 
  Real m_residualScale = 1.0;
 
  KrylovGMRES() noexcept = default;
 
  ~KrylovGMRES() noexcept;
 
  void
  define(AMRMultigridKrylovOp<T>* a_op, const Vector<LevelData<T>*>& a_template) noexcept;
 
  void
  undefine() noexcept;
 
  bool
  isDefined() const noexcept;
 
  bool
  solve(Vector<LevelData<T>*>& a_phi, const Vector<LevelData<T>*>& a_rhs) noexcept;
 
protected:
  AMRMultigridKrylovOp<T>* m_op = nullptr;
 
  bool m_isDefined = false;
 
  int m_allocRestart = 0;
 
  Vector<Vector<LevelData<T>*>> m_basis;
 
  Vector<LevelData<T>*> m_w;
 
  Vector<LevelData<T>*> m_z;
 
  Vector<LevelData<T>*> m_c;
 
  Vector<LevelData<T>*> m_r;
 
  Vector<Real> m_H;
 
  Vector<Real> m_cs;
 
  Vector<Real> m_sn;
 
  Vector<Real> m_g;
 
  Vector<Real> m_y;
 
  Vector<Real> m_hcol;
 
  Vector<Real> m_corr;
};
 
#include <CD_NamespaceFooter.H>
 
#include <CD_KrylovGMRESImplem.H>
 
#endif