Steady-state problem operator with an equation system. More...

#include <EquationSystemProblemOperator.h>

Inheritance diagram for Moose::MFEM::EquationSystemProblemOperator:

Public Member Functions
	EquationSystemProblemOperator (MFEMProblemData &problem)

void	SetGridFunctions () override

void	Init (mfem::BlockVector &X) override

virtual void	Solve (mfem::Vector &X) override

Moose::MFEM::EquationSystem *	GetEquationSystem () const override
	Returns a pointer to the operator's equation system. More...

void	Mult (const mfem::Vector &, mfem::Vector &) const override

virtual void	SetTestVariablesFromTrueVectors ()

virtual void	SetTrialVariablesFromTrueVectors ()

Public Attributes
mfem::Array< int >	_block_true_offsets

mfem::BlockVector	_true_x

mfem::BlockVector	_true_rhs

mfem::OperatorHandle	_equation_system_operator

Protected Attributes
MFEMProblemData &	_problem

std::vector< std::string >	_test_var_names

std::vector< mfem::ParGridFunction * >	_test_variables

std::vector< std::string >	_trial_var_names

std::vector< mfem::ParGridFunction * >	_trial_variables

Private Attributes
std::shared_ptr< Moose::MFEM::EquationSystem >	_equation_system {nullptr}

Detailed Description

Steady-state problem operator with an equation system.

Definition at line 19 of file EquationSystemProblemOperator.h.

Constructor & Destructor Documentation

◆ EquationSystemProblemOperator()

Moose::MFEM::EquationSystemProblemOperator::EquationSystemProblemOperator ( MFEMProblemData & problem )

inline

Definition at line 22 of file EquationSystemProblemOperator.h.

     : ProblemOperator(problem), _equation_system(problem.eqn_system)
   {
   }

Member Function Documentation

◆ GetEquationSystem()

Moose::MFEM::EquationSystem* Moose::MFEM::EquationSystemProblemOperator::GetEquationSystem ( ) const

inlineoverridevirtual

Returns a pointer to the operator's equation system.

Implements Moose::MFEM::EquationSystemInterface.

Definition at line 31 of file EquationSystemProblemOperator.h.

Referenced by Init(), SetGridFunctions(), and Solve().

   {
     if (!_equation_system)
     {
       MFEM_ABORT("No equation system has been added to ProblemOperator.");
     }
 
     return _equation_system.get();
   }

◆ Init()

void Moose::MFEM::EquationSystemProblemOperator::Init ( mfem::BlockVector & X )

overridevirtual

Reimplemented from Moose::MFEM::ProblemOperatorInterface.

Definition at line 24 of file EquationSystemProblemOperator.C.

 {
   ProblemOperator::Init(X);
 
   GetEquationSystem()->BuildEquationSystem();
 }

◆ Mult()

void Moose::MFEM::ProblemOperator::Mult	(	const mfem::Vector &	,
		mfem::Vector &
	)		const

inlineoverrideinherited

Definition at line 27 of file ProblemOperator.h.

27 {}

◆ SetGridFunctions()

void Moose::MFEM::EquationSystemProblemOperator::SetGridFunctions ( )

overridevirtual

Reimplemented from Moose::MFEM::ProblemOperatorInterface.

Definition at line 17 of file EquationSystemProblemOperator.C.

 {
   _trial_var_names = GetEquationSystem()->TrialVarNames();
   ProblemOperator::SetGridFunctions();
 }

◆ SetTestVariablesFromTrueVectors()

void Moose::MFEM::ProblemOperatorInterface::SetTestVariablesFromTrueVectors ( )

virtualinherited

Definition at line 47 of file ProblemOperatorInterface.C.

Referenced by Moose::MFEM::TimeDomainEquationSystemProblemOperator::ImplicitSolve().

 {
   for (unsigned int ind = 0; ind < _test_variables.size(); ++ind)
   {
     _test_variables.at(ind)->SetFromTrueVector();
   }
 }

◆ SetTrialVariablesFromTrueVectors()

void Moose::MFEM::ProblemOperatorInterface::SetTrialVariablesFromTrueVectors ( )

virtualinherited

Definition at line 56 of file ProblemOperatorInterface.C.

Referenced by Moose::MFEM::TimeDomainEquationSystemProblemOperator::ImplicitSolve().

 {
   for (unsigned int ind = 0; ind < _trial_variables.size(); ++ind)
   {
     _trial_variables.at(ind)->SetFromTrueVector();
   }
 }

◆ Solve()

void Moose::MFEM::EquationSystemProblemOperator::Solve ( mfem::Vector & X )

overridevirtual

Reimplemented from Moose::MFEM::ProblemOperator.

Definition at line 32 of file EquationSystemProblemOperator.C.

 {
   GetEquationSystem()->BuildJacobian(_true_x, _true_rhs);
 
   if (_problem.jacobian_solver->isLOR() && _equation_system->_test_var_names.size() > 1)
     mooseError("LOR solve is only supported for single-variable systems");
 
   _problem.jacobian_solver->updateSolver(
       *_equation_system->_blfs.Get(_equation_system->_test_var_names.at(0)),
       _equation_system->_ess_tdof_lists.at(0));
 
   _problem.nonlinear_solver->SetSolver(_problem.jacobian_solver->getSolver());
   _problem.nonlinear_solver->SetOperator(*GetEquationSystem());
   _problem.nonlinear_solver->Mult(_true_rhs, _true_x);
 
   GetEquationSystem()->RecoverFEMSolution(_true_x, _problem.gridfunctions);
 }