10 #ifdef MOOSE_MFEM_ENABLED 24 mooseError(
"Eigenproblems are only supported in single-variable systems");
44 mfem::Array<mfem::real_t> eigenvalues;
53 for (
int i = 0; i < eigenvalues.Size(); ++i)
54 gridfunctions.
Get(trial_var_name + sep + std::to_string(i))
void BuildEquationSystemOperator()
Add kernels/bcs and assemble the linear part of the equation system.
void BuildEigenproblemJacobian(mfem::BlockVector &trueX, mfem::OperatorHandle &massRHS)
Build eigenproblem system, with essential boundary conditions accounted for.
MFEMProblemData & _problem_data
std::vector< std::string > _trial_var_names
Vector of names of state gridfunctions used in formulation, ordered by appearance in block vector dur...
void mooseError(Args &&... args)
Emit an error message with the given stringified, concatenated args and terminate the application...
Base class for eigensolvers.
EigenproblemEquationSystem * GetEquationSystem() const override
Returns a pointer to the operator's equation system.
virtual void setMassMatrix(mfem::OperatorHandle &mass)=0
Sets the mass matrix for the eigensolver in derived classes.
std::unique_ptr< T_DEST, T_DELETER > dynamic_pointer_cast(std::unique_ptr< T_SRC, T_DELETER > &src)
These are reworked from https://stackoverflow.com/a/11003103.
virtual void getEigenvalues(mfem::Array< mfem::real_t > &eigenvalues) const =0
Retrieves the computed eigenvalues.
T * Get(const std::string &field_name) const
Returns a non-owning pointer to the field. This is guaranteed to return a non-null pointer...
virtual void Solve() override
mfem::BlockVector _true_x
virtual void BuildEquationSystem()
Build all forms comprising this EquationSystem.
void RecoverEigenproblemSolution(Moose::MFEM::GridFunctions &gridfunctions, MFEMEigensolverBase *eigensolver)
virtual const mfem::HypreParVector & getEigenvector(int index) const =0
Retrieves the computed eigenvector corresponding to the given index.
mfem::OperatorHandle _mass_rhs
std::string mode_separator
Utilities for converting between vector(s) of libMesh Points and MFEM Vector(s).
Moose::MFEM::GridFunctions gridfunctions
std::shared_ptr< MFEMSolverBase > jacobian_solver