MFEMLinearSolverBase.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#ifdef MOOSE_MFEM_ENABLED

#pragma once

#include "MFEMSolverBase.h"

class MFEMProblemSolve;

namespace Moose::MFEM
{
class LinearSolverBase : public SolverBase
{
public:
  static InputParameters validParams();

  LinearSolverBase(const InputParameters & parameters);

  template <typename T>
  void SetPreconditioner(T & solver);

  LinearSolverBase * GetPreconditioner() { return _preconditioner.get(); };

  virtual void SetupLOR(mfem::ParBilinearForm & /*a*/, mfem::Array<int> & /*tdofs*/) {}

  virtual void SetOperator(mfem::OperatorHandle & op);

  bool IsLOR() const { return _lor || (_preconditioner && _preconditioner->IsLOR()); }

  virtual void Solve() { mooseError("'solve' method not used in this solver type."); }

protected:
  virtual void CheckSpectralEquivalence(mfem::ParBilinearForm & blf) const;

  bool _lor;

  std::shared_ptr<LinearSolverBase> _preconditioner;

private:
  friend class ::MFEMProblemSolve;
};
} // namespace Moose::MFEM

#endif