MFEMHypreLOBPCG.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 "MFEMEigensolverBase.h"

class MFEMHypreLOBPCG : public Moose::MFEM::EigensolverBase
{
public:
  static InputParameters validParams();

  MFEMHypreLOBPCG(const InputParameters & parameters);

  virtual void SetOperator(mfem::OperatorHandle & op) override
  {
    if (_preconditioner)
      _preconditioner->SetOperator(op);
    _eigensolver->SetOperator(*op);
  }

  virtual void SetMassMatrix(mfem::OperatorHandle & mass) override
  {
    _eigensolver->SetMassMatrix(*mass);
  }

  virtual void Solve() override { _eigensolver->Solve(); }

  virtual void GetEigenvalues(mfem::Array<mfem::real_t> & eigenvalues) const override
  {
    _eigensolver->GetEigenvalues(eigenvalues);
  }

  virtual const mfem::HypreParVector & GetEigenvector(int index) const override
  {
    return _eigensolver->GetEigenvector(index);
  }

protected:
  virtual void ConstructSolver() override;

  std::unique_ptr<mfem::HypreLOBPCG> _eigensolver;
};

#endif