Moose::MFEM::TimeDomainEquationSystemProblemOperator Class Reference

Problem operator for time-dependent problems with an equation system. More...

#include <TimeDomainEquationSystemProblemOperator.h>

Inheritance diagram for Moose::MFEM::TimeDomainEquationSystemProblemOperator:

Public Member Functions
	TimeDomainEquationSystemProblemOperator (MFEMProblemData &problem)
void	SetGridFunctions () override
void	Init (mfem::BlockVector &X) override
void	ImplicitSolve (const double dt, const mfem::Vector &X, mfem::Vector &dX_dt) override
Moose::MFEM::TimeDependentEquationSystem *	GetEquationSystem () const override
	Returns a pointer to the operator's equation system. More...
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 Member Functions
void	BuildEquationSystemOperator (double dt)

Protected Attributes
MFEMProblemData &	_problem
	Reference to the current problem. More...
std::vector< std::string >	_test_var_names
	Vector of names of state gridfunctions used in formulation, ordered by appearance in block vector during solve. More...
std::vector< mfem::ParGridFunction * >	_test_variables
const mfem::Vector *	_test_true_vector = nullptr
std::vector< std::string >	_trial_var_names
	Vector of names of state gridfunctions used in formulation, ordered by appearance in block vector during solve. More...
std::vector< mfem::ParGridFunction * >	_trial_variables

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

Detailed Description

Problem operator for time-dependent problems with an equation system.

Definition at line 20 of file TimeDomainEquationSystemProblemOperator.h.

Constructor & Destructor Documentation

TimeDomainEquationSystemProblemOperator()

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

inline

Definition at line 24 of file TimeDomainEquationSystemProblemOperator.h.

    : TimeDomainProblemOperator(problem),
      _equation_system(
          std::dynamic_pointer_cast<Moose::MFEM::TimeDependentEquationSystem>(problem.eqn_system))
  {
  }

Member Function Documentation

BuildEquationSystemOperator()

void Moose::MFEM::TimeDomainEquationSystemProblemOperator::BuildEquationSystemOperator ( double  dt )

protected

Definition at line 60 of file TimeDomainEquationSystemProblemOperator.C.

Referenced by ImplicitSolve().

{
  GetEquationSystem()->SetTimeStep(dt);
  GetEquationSystem()->UpdateEquationSystem();
  GetEquationSystem()->BuildJacobian(_true_x, _true_rhs);
}

GetEquationSystem()

Moose::MFEM::TimeDependentEquationSystem* Moose::MFEM::TimeDomainEquationSystemProblemOperator::GetEquationSystem ( ) const

inlineoverridevirtual

Returns a pointer to the operator's equation system.

Implements Moose::MFEM::EquationSystemInterface.

Definition at line 36 of file TimeDomainEquationSystemProblemOperator.h.

Referenced by BuildEquationSystemOperator(), ImplicitSolve(), Init(), and SetGridFunctions().

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

    return _equation_system.get();
  }

ImplicitSolve()

void Moose::MFEM::TimeDomainEquationSystemProblemOperator::ImplicitSolve ( const double  dt, const mfem::Vector &  X, mfem::Vector &  dX_dt  )

override

Definition at line 32 of file TimeDomainEquationSystemProblemOperator.C.

{
  dX_dt = 0.0;
  SetTestVariablesFromTrueVectors();
  for (unsigned int ind = 0; ind < _trial_variables.size(); ++ind)
  {
    _trial_variables.at(ind)->MakeTRef(
        _trial_variables.at(ind)->ParFESpace(), dX_dt, _block_true_offsets[ind]);
  }
  _problem.coefficients.setTime(GetTime());
  BuildEquationSystemOperator(dt);

  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, dX_dt);
  SetTrialVariablesFromTrueVectors();
}

Init()

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

overridevirtual

Reimplemented from Moose::MFEM::ProblemOperatorInterface.

Definition at line 25 of file TimeDomainEquationSystemProblemOperator.C.

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

SetGridFunctions()

void Moose::MFEM::TimeDomainEquationSystemProblemOperator::SetGridFunctions ( )

overridevirtual

Reimplemented from Moose::MFEM::ProblemOperatorInterface.

Definition at line 17 of file TimeDomainEquationSystemProblemOperator.C.

{
  _test_var_names = GetEquationSystem()->TestVarNames();
  _trial_var_names = GetEquationSystem()->TrialVarNames();
  TimeDomainProblemOperator::SetGridFunctions();
}

SetTestVariablesFromTrueVectors()

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

virtualinherited

Definition at line 57 of file ProblemOperatorInterface.C.

Referenced by ImplicitSolve(), and Moose::MFEM::ProblemOperatorInterface::Init().

{
  for (unsigned int ind = 0; ind < _test_variables.size(); ++ind)
  {
    auto * const test_var = _test_variables.at(ind);

    // We must sync the memory flags from the true true vector to the grid function copy of the true
    // vector
    mooseAssert(_test_true_vector, "The true vector should already have been set");
    test_var->GetTrueVector().SyncMemory(*_test_true_vector);
    test_var->SetFromTrueVector();
  }
}

SetTrialVariablesFromTrueVectors()

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

virtualinherited

Definition at line 72 of file ProblemOperatorInterface.C.

Referenced by ImplicitSolve().

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

Member Data Documentation

_block_true_offsets

mfem::Array<int> Moose::MFEM::ProblemOperatorInterface::_block_true_offsets

inherited

Definition at line 29 of file ProblemOperatorInterface.h.

Referenced by ImplicitSolve(), Moose::MFEM::ProblemOperatorInterface::Init(), Moose::MFEM::ProblemOperator::SetGridFunctions(), Moose::MFEM::ProblemOperatorInterface::SetGridFunctions(), and Moose::MFEM::TimeDomainProblemOperator::SetGridFunctions().

_equation_system

std::shared_ptr<Moose::MFEM::TimeDependentEquationSystem> Moose::MFEM::TimeDomainEquationSystemProblemOperator::_equation_system {nullptr}

private

Definition at line 50 of file TimeDomainEquationSystemProblemOperator.h.

Referenced by GetEquationSystem(), and ImplicitSolve().

_equation_system_operator

mfem::OperatorHandle Moose::MFEM::ProblemOperatorInterface::_equation_system_operator

inherited

Definition at line 32 of file ProblemOperatorInterface.h.

_problem

MFEMProblemData& Moose::MFEM::ProblemOperatorInterface::_problem

protectedinherited

Reference to the current problem.

Definition at line 36 of file ProblemOperatorInterface.h.

Referenced by ImplicitSolve(), Moose::MFEM::ProblemOperatorInterface::SetGridFunctions(), and Moose::MFEM::EquationSystemProblemOperator::Solve().

_test_true_vector

const mfem::Vector* Moose::MFEM::ProblemOperatorInterface::_test_true_vector = nullptr

protectedinherited

Definition at line 42 of file ProblemOperatorInterface.h.

Referenced by Moose::MFEM::ProblemOperatorInterface::Init(), and Moose::MFEM::ProblemOperatorInterface::SetTestVariablesFromTrueVectors().

_test_var_names

std::vector<std::string> Moose::MFEM::ProblemOperatorInterface::_test_var_names

protectedinherited

Vector of names of state gridfunctions used in formulation, ordered by appearance in block vector during solve.

Definition at line 40 of file ProblemOperatorInterface.h.

Referenced by Moose::MFEM::ProblemOperatorInterface::SetGridFunctions(), Moose::MFEM::EquationSystemProblemOperator::SetGridFunctions(), and SetGridFunctions().

_test_variables

std::vector<mfem::ParGridFunction *> Moose::MFEM::ProblemOperatorInterface::_test_variables

protectedinherited

Definition at line 41 of file ProblemOperatorInterface.h.

Referenced by Moose::MFEM::ProblemOperatorInterface::Init(), Moose::MFEM::ProblemOperatorInterface::SetGridFunctions(), and Moose::MFEM::ProblemOperatorInterface::SetTestVariablesFromTrueVectors().

_trial_var_names

std::vector<std::string> Moose::MFEM::ProblemOperatorInterface::_trial_var_names

protectedinherited

Vector of names of state gridfunctions used in formulation, ordered by appearance in block vector during solve.

Definition at line 46 of file ProblemOperatorInterface.h.

Referenced by Moose::MFEM::ProblemOperatorInterface::SetGridFunctions(), Moose::MFEM::EquationSystemProblemOperator::SetGridFunctions(), and SetGridFunctions().

_trial_variables

std::vector<mfem::ParGridFunction *> Moose::MFEM::ProblemOperatorInterface::_trial_variables

protectedinherited

Definition at line 47 of file ProblemOperatorInterface.h.

Referenced by ImplicitSolve(), Moose::MFEM::ProblemOperator::SetGridFunctions(), Moose::MFEM::ProblemOperatorInterface::SetGridFunctions(), Moose::MFEM::TimeDomainProblemOperator::SetGridFunctions(), and Moose::MFEM::ProblemOperatorInterface::SetTrialVariablesFromTrueVectors().

_true_rhs

mfem::BlockVector Moose::MFEM::ProblemOperatorInterface::_true_rhs

inherited

Definition at line 31 of file ProblemOperatorInterface.h.

Referenced by BuildEquationSystemOperator(), ImplicitSolve(), Moose::MFEM::ProblemOperatorInterface::SetGridFunctions(), and Moose::MFEM::EquationSystemProblemOperator::Solve().

_true_x

mfem::BlockVector Moose::MFEM::ProblemOperatorInterface::_true_x

inherited

Definition at line 31 of file ProblemOperatorInterface.h.

Referenced by BuildEquationSystemOperator(), Moose::MFEM::ProblemOperatorInterface::SetGridFunctions(), and Moose::MFEM::EquationSystemProblemOperator::Solve().

---

The documentation for this class was generated from the following files:

• include/mfem/problem_operators/TimeDomainEquationSystemProblemOperator.h
• src/mfem/problem_operators/TimeDomainEquationSystemProblemOperator.C