PipelineManager interface

Implementation of basic interface for rapid problem implementation. More...

Collaboration diagram for PipelineManager interface:

Files
file	PipelineManager.cpp
	Implementation of basic interface.
file	PipelineManager.hpp
	Header file for basic interface.

Classes
struct	MoFEM::PipelineManager
	PipelineManager interface. More...

Functions
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpDomainLhsPipeline ()
	Get the Op Domain Lhs Pipeline object.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpDomainRhsPipeline ()
	Get the Op Domain Rhs Pipeline object.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpBoundaryLhsPipeline ()
	Get the Op Boundary Lhs Pipeline object.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpBoundaryRhsPipeline ()
	Get the Op Boundary Rhs Pipeline object.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpDomainPostProcPipeline ()
	Get the Op Domain PostProc Pipeline object.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpBoundaryPostProcPipeline ()
	Get the Op Boundary PostProc Pipeline object.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpEvaluationPipeline ()
	Get the Op Domain Evaluation Pipeline object.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpSkeletonLhsPipeline ()
	Get the Op Skeleton Lhs Pipeline object.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpSkeletonRhsPipeline ()
	Get the Op Skeleton Rhs Pipeline object.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpDomainExplicitRhsPipeline ()
	Get the Op Domain Rhs Pipeline object for implicit-explicit G term.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpBoundaryExplicitRhsPipeline ()
	Get the Op Boundary Rhs Pipeline object for implicit-explicit G term.
template<int DIM = -1>
boost::ptr_deque< UserDataOperator > &	MoFEM::PipelineManager::getOpSkeletonExplicitRhsPipeline ()
	Get the Op Skeleton Rhs Pipeline object for implicit-explicit G term.
MoFEMErrorCode	MoFEM::PipelineManager::loopFiniteElements (SmartPetscObj< DM > dm=nullptr)
	Iterate finite elements.
MoFEMErrorCode	MoFEM::PipelineManager::loopFiniteElementsPostProc (SmartPetscObj< DM > dm=nullptr)
	Iterate postprocessing finite elements.
MoFEMErrorCode	MoFEM::PipelineManager::loopDomainPostProc (SmartPetscObj< DM > dm=nullptr)
	Iterate domain postprocessing finite elements only.
MoFEMErrorCode	MoFEM::PipelineManager::loopBoundaryPostProc (SmartPetscObj< DM > dm=nullptr)
	Iterate boundary postprocessing finite elements only.
MoFEMErrorCode	MoFEM::PipelineManager::loopFiniteElementsEvaluation (SmartPetscObj< DM > dm=nullptr)
	Iterate evaluation finite elements.
SmartPetscObj< KSP >	MoFEM::PipelineManager::createKSP (SmartPetscObj< DM > dm=nullptr)
	Create KSP (linear) solver.
SmartPetscObj< SNES >	MoFEM::PipelineManager::createSNES (SmartPetscObj< DM > dm=nullptr)
	Create SNES (nonlinear) solver.
SmartPetscObj< TS >	MoFEM::PipelineManager::createTSEX (SmartPetscObj< DM > dm=nullptr)
	Create TS (time) explicit solver.
SmartPetscObj< TS >	MoFEM::PipelineManager::createTSIM (SmartPetscObj< DM > dm=nullptr)
	Create TS (time) implicit solver.
SmartPetscObj< TS >	MoFEM::PipelineManager::createTSIM2 (SmartPetscObj< DM > dm=nullptr)
	Create TS (time) solver for second order equation in time.
SmartPetscObj< TS >	MoFEM::PipelineManager::createTSIMEX (SmartPetscObj< DM > dm=nullptr)
	Create TS (time) implicit-explicit solver.

Detailed Description

Implementation of basic interface for rapid problem implementation.

Function Documentation

createKSP()

SmartPetscObj< KSP > MoFEM::PipelineManager::createKSP

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Create KSP (linear) solver.

Parameters

dm	Optional DM object (default: nullptr)

Returns

SmartPetscObj<KSP> Smart pointer to PETSc KSP linear solver

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/tutorials/adv-6_dg_projection/between_meshes_dg_projection.cpp, mofem/tutorials/clx-0_acoustics/helmholtz.cpp, mofem/tutorials/mix-0_mixed_poisson/mixed_poisson.cpp, mofem/tutorials/mix-1_light_intensity_equation/phase.cpp, mofem/tutorials/scl-0_least_squares/approximation.cpp, and mofem/tutorials/scl-1_poisson_2d_homogeneous/poisson_2d_homogeneous.cpp.

Definition at line 187 of file PipelineManager.cpp.

                                                                  {
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
 
  auto copy_dm_struture = [&](auto simple_dm) {
    MPI_Comm comm;
    CHKERR PetscObjectGetComm(getPetscObject(simple_dm.get()), &comm);
    DMType type;
    CHKERR DMGetType(simple_dm, &type);
    dm = createDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
 
  const MoFEM::Problem *prb_ptr;
  CHKERR DMMoFEMGetProblemPtr(dm, &prb_ptr);
 
  auto set_dm_section = [&](auto dm) {
    MoFEMFunctionBegin;
    auto section =
        m_field.getInterface<ISManager>()->sectionCreate(prb_ptr->getName());
    CHKERR DMSetSection(dm, section);
    MoFEMFunctionReturn(0);
  };
  CHKERR set_dm_section(dm);
 
  boost::shared_ptr<FEMethod> null;
 
  getDMKspCtx(dm)->clearLoops();
 
  // Add element to calculate lhs of stiff part
  if (feDomainLhs)
    CHKERR DMMoFEMKSPSetComputeOperators(dm, simple->getDomainFEName(),
                                         feDomainLhs, null, null);
  if (feBoundaryLhs)
    CHKERR DMMoFEMKSPSetComputeOperators(dm, simple->getBoundaryFEName(),
                                         feBoundaryLhs, null, null);
  if (feSkeletonLhs)
    CHKERR DMMoFEMKSPSetComputeOperators(dm, simple->getSkeletonFEName(),
                                         feSkeletonLhs, null, null);
  if (feMeshsetLhs)
    CHKERR DMMoFEMKSPSetComputeOperators(dm, simple->getMeshsetFEName(),
                                         feMeshsetLhs, null, null);
 
  // Add element to calculate rhs of stiff part
  if (feDomainRhs)
    CHKERR DMMoFEMKSPSetComputeRHS(dm, simple->getDomainFEName(), feDomainRhs,
                                   null, null);
  if (feBoundaryRhs)
    CHKERR DMMoFEMKSPSetComputeRHS(dm, simple->getBoundaryFEName(),
                                   feBoundaryRhs, null, null);
  if (feSkeletonRhs)
    CHKERR DMMoFEMKSPSetComputeRHS(dm, simple->getSkeletonFEName(),
                                   feSkeletonRhs, null, null);
  if (feMeshsetRhs)
    CHKERR DMMoFEMKSPSetComputeRHS(dm, simple->getMeshsetFEName(), feMeshsetRhs,
                                   null, null);
 
  auto ksp = MoFEM::createKSP(m_field.get_comm());
  CHKERR KSPSetDM(ksp, dm);
  return ksp;
}

createSNES()

SmartPetscObj< SNES > MoFEM::PipelineManager::createSNES

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Create SNES (nonlinear) solver.

Parameters

dm	Optional DM object (default: nullptr)

Returns

SmartPetscObj<SNES> Smart pointer to PETSc SNES nonlinear solver

Definition at line 253 of file PipelineManager.cpp.

                                                                    {
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
 
  auto copy_dm_struture = [&](auto simple_dm) {
    MPI_Comm comm;
    CHKERR PetscObjectGetComm(getPetscObject(simple_dm.get()), &comm);
    DMType type;
    CHKERR DMGetType(simple_dm, &type);
    dm = createDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
 
  const MoFEM::Problem *prb_ptr;
  CHKERR DMMoFEMGetProblemPtr(dm, &prb_ptr);
 
  auto set_dm_section = [&](auto dm) {
    MoFEMFunctionBegin;
    auto section =
        m_field.getInterface<ISManager>()->sectionCreate(prb_ptr->getName());
    CHKERR DMSetSection(dm, section);
    MoFEMFunctionReturn(0);
  };
  CHKERR set_dm_section(dm);
 
  getDMSnesCtx(dm)->clearLoops();
 
  boost::shared_ptr<FEMethod> null;
 
  // Add element to calculate lhs of stiff part
  if (feDomainLhs)
    CHKERR DMMoFEMSNESSetJacobian(dm, simple->getDomainFEName(), feDomainLhs,
                                  null, null);
  if (feBoundaryLhs)
    CHKERR DMMoFEMSNESSetJacobian(dm, simple->getBoundaryFEName(),
                                  feBoundaryLhs, null, null);
  if (feSkeletonLhs)
    CHKERR DMMoFEMSNESSetJacobian(dm, simple->getSkeletonFEName(),
                                  feSkeletonLhs, null, null);
  if (feMeshsetLhs)
    CHKERR DMMoFEMSNESSetJacobian(dm, simple->getMeshsetFEName(), feMeshsetLhs,
                                  null, null);
 
  // Add element to calculate rhs of stiff part
  if (feDomainRhs)
    CHKERR DMMoFEMSNESSetFunction(dm, simple->getDomainFEName(), feDomainRhs,
                                  null, null);
  if (feBoundaryRhs)
    CHKERR DMMoFEMSNESSetFunction(dm, simple->getBoundaryFEName(),
                                  feBoundaryRhs, null, null);
  if (feSkeletonRhs)
    CHKERR DMMoFEMSNESSetFunction(dm, simple->getSkeletonFEName(),
                                  feSkeletonRhs, null, null);
  if (feMeshsetRhs)
    CHKERR DMMoFEMSNESSetFunction(dm, simple->getMeshsetFEName(), feMeshsetRhs,
                                  null, null);
 
  auto snes = MoFEM::createSNES(m_field.get_comm());
  CHKERR SNESSetDM(snes, dm);
  return snes;
}

createTSEX()

SmartPetscObj< TS > MoFEM::PipelineManager::createTSEX

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Create TS (time) explicit solver.

Parameters

dm	Optional DM object (default: nullptr)

Returns

SmartPetscObj<TS> Smart pointer to PETSc TS explicit time solver

Definition at line 342 of file PipelineManager.cpp.

                                                                  {
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
 
  auto copy_dm_struture = [&](auto simple_dm) {
    MPI_Comm comm;
    CHKERR PetscObjectGetComm(getPetscObject(simple_dm.get()), &comm);
    DMType type;
    CHKERR DMGetType(simple_dm, &type);
    dm = createDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
 
  const MoFEM::Problem *prb_ptr;
  CHKERR DMMoFEMGetProblemPtr(dm, &prb_ptr);
 
  auto set_dm_section = [&](auto dm) {
    MoFEMFunctionBegin;
    auto section =
        m_field.getInterface<ISManager>()->sectionCreate(prb_ptr->getName());
    CHKERR DMSetSection(dm, section);
    MoFEMFunctionReturn(0);
  };
  CHKERR set_dm_section(dm);
 
  boost::shared_ptr<FEMethod> null;
 
  getDMTsCtx(dm)->clearLoops();
 
  // Add element to calculate rhs of slow part
  if (feDomainExplicitRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getDomainFEName(),
                                   feDomainExplicitRhs, null, null);
  if (feBoundaryExplicitRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getBoundaryFEName(),
                                   feBoundaryExplicitRhs, null, null);
  if (feSkeletonExplicitRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getSkeletonFEName(),
                                   feSkeletonExplicitRhs, null, null);
  if (feMeshsetRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getMeshsetFEName(),
                                   feMeshsetExplicitRhs, null, null);
 
  // Note: More cases for explict, and implicit time interation cases can be
  // implemented here.
 
  auto ts = MoFEM::createTS(m_field.get_comm());
  CHKERR TSSetDM(ts, dm);
  return ts;
}

createTSIM()

SmartPetscObj< TS > MoFEM::PipelineManager::createTSIM

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Create TS (time) implicit solver.

Parameters

dm	Optional DM object (default: nullptr)

Returns

SmartPetscObj<TS> Smart pointer to PETSc TS implicit time solver

Examples

mofem/tutorials/adv-5_poroelasticity/seepage.cpp, and mofem/tutorials/scl-8_heat_radiation/radiation.cpp.

Definition at line 397 of file PipelineManager.cpp.

                                                                  {
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
 
  auto copy_dm_struture = [&](auto simple_dm) {
    MPI_Comm comm;
    CHKERR PetscObjectGetComm(getPetscObject(simple_dm.get()), &comm);
    DMType type;
    CHKERR DMGetType(simple_dm, &type);
    dm = createDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
 
  const MoFEM::Problem *prb_ptr;
  CHKERR DMMoFEMGetProblemPtr(dm, &prb_ptr);
 
  auto set_dm_section = [&](auto dm) {
    MoFEMFunctionBegin;
    auto section =
        m_field.getInterface<ISManager>()->sectionCreate(prb_ptr->getName());
    CHKERR DMSetSection(dm, section);
    MoFEMFunctionReturn(0);
  };
  CHKERR set_dm_section(dm);
 
  boost::shared_ptr<FEMethod> null;
 
  // Add element to calculate lhs of stiff part
  if (feDomainLhs)
    CHKERR DMMoFEMTSSetIJacobian(dm, simple->getDomainFEName(), feDomainLhs,
                                 null, null);
  if (feBoundaryLhs)
    CHKERR DMMoFEMTSSetIJacobian(dm, simple->getBoundaryFEName(), feBoundaryLhs,
                                 null, null);
  if (feSkeletonLhs)
    CHKERR DMMoFEMTSSetIJacobian(dm, simple->getSkeletonFEName(), feSkeletonLhs,
                                 null, null);
  if (feMeshsetLhs)
    CHKERR DMMoFEMTSSetIJacobian(dm, simple->getMeshsetFEName(), feMeshsetLhs,
                                 null, null);
 
  // Add element to calculate rhs of stiff part
  if (feDomainRhs)
    CHKERR DMMoFEMTSSetIFunction(dm, simple->getDomainFEName(), feDomainRhs,
                                 null, null);
  if (feBoundaryRhs)
    CHKERR DMMoFEMTSSetIFunction(dm, simple->getBoundaryFEName(), feBoundaryRhs,
                                 null, null);
  if (feSkeletonRhs)
    CHKERR DMMoFEMTSSetIFunction(dm, simple->getSkeletonFEName(), feSkeletonRhs,
                                 null, null);
  if (feMeshsetRhs)
    CHKERR DMMoFEMTSSetIFunction(dm, simple->getMeshsetFEName(), feMeshsetRhs,
                                 null, null);
 
  // Note: More cases for explict, and implicit time interation cases can be
  // implemented here.
 
  auto ts = MoFEM::createTS(m_field.get_comm());
  CHKERR TSSetDM(ts, dm);
  return ts;
}

createTSIM2()

SmartPetscObj< TS > MoFEM::PipelineManager::createTSIM2

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Create TS (time) solver for second order equation in time.

Parameters

dm	Optional DM object (default: nullptr)

Returns

SmartPetscObj<TS> Smart pointer to PETSc TS second-order time solver

Definition at line 464 of file PipelineManager.cpp.

                                                                   {
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
 
  auto copy_dm_struture = [&](auto simple_dm) {
    MPI_Comm comm;
    CHKERR PetscObjectGetComm(getPetscObject(simple_dm.get()), &comm);
    DMType type;
    CHKERR DMGetType(simple_dm, &type);
    dm = createDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
 
  const MoFEM::Problem *prb_ptr;
  CHKERR DMMoFEMGetProblemPtr(dm, &prb_ptr);
 
  auto set_dm_section = [&](auto dm) {
    MoFEMFunctionBegin;
    auto section =
        m_field.getInterface<ISManager>()->sectionCreate(prb_ptr->getName());
    CHKERR DMSetSection(dm, section);
    MoFEMFunctionReturn(0);
  };
  CHKERR set_dm_section(dm);
 
  boost::shared_ptr<FEMethod> null;
 
  // Add element to calculate lhs of stiff part
  if (feDomainLhs)
    CHKERR DMMoFEMTSSetI2Jacobian(dm, simple->getDomainFEName(), feDomainLhs,
                                  null, null);
  if (feBoundaryLhs)
    CHKERR DMMoFEMTSSetI2Jacobian(dm, simple->getBoundaryFEName(),
                                  feBoundaryLhs, null, null);
  if (feSkeletonLhs)
    CHKERR DMMoFEMTSSetI2Jacobian(dm, simple->getSkeletonFEName(),
                                  feSkeletonLhs, null, null);
  if (feMeshsetLhs)
    CHKERR DMMoFEMTSSetI2Jacobian(dm, simple->getMeshsetFEName(), feMeshsetLhs,
                                  null, null);
 
  // Add element to calculate rhs of stiff part
  if (feDomainRhs)
    CHKERR DMMoFEMTSSetI2Function(dm, simple->getDomainFEName(), feDomainRhs,
                                  null, null);
  if (feBoundaryRhs)
    CHKERR DMMoFEMTSSetI2Function(dm, simple->getBoundaryFEName(),
                                  feBoundaryRhs, null, null);
  if (feSkeletonRhs)
    CHKERR DMMoFEMTSSetI2Function(dm, simple->getSkeletonFEName(),
                                  feSkeletonRhs, null, null);
  if (feMeshsetRhs)
    CHKERR DMMoFEMTSSetI2Function(dm, simple->getMeshsetFEName(), feMeshsetRhs,
                                  null, null);
 
  // Note: More cases for explict, and implicit time interation cases can be
  // implemented here.
 
  auto ts = MoFEM::createTS(m_field.get_comm());
  CHKERR TSSetDM(ts, dm);
  return ts;
}

createTSIMEX()

SmartPetscObj< TS > MoFEM::PipelineManager::createTSIMEX

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Create TS (time) implicit-explicit solver.

Parameters

dm	Optional DM object (default: nullptr)

Returns

SmartPetscObj<TS> Smart pointer to PETSc TS implicit-explicit time solver

Definition at line 531 of file PipelineManager.cpp.

                                                                    {
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
 
  auto copy_dm_struture = [&](auto simple_dm) {
    MPI_Comm comm;
    CHKERR PetscObjectGetComm(getPetscObject(simple_dm.get()), &comm);
    DMType type;
    CHKERR DMGetType(simple_dm, &type);
    dm = createDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
 
  const MoFEM::Problem *prb_ptr;
  CHKERR DMMoFEMGetProblemPtr(dm, &prb_ptr);
 
  auto set_dm_section = [&](auto dm) {
    MoFEMFunctionBegin;
    auto section =
        m_field.getInterface<ISManager>()->sectionCreate(prb_ptr->getName());
    CHKERR DMSetSection(dm, section);
    MoFEMFunctionReturn(0);
  };
  CHKERR set_dm_section(dm);
 
  boost::shared_ptr<FEMethod> null;
 
  // Add element to calculate lhs of stiff part
  if (feDomainLhs)
    CHKERR DMMoFEMTSSetIJacobian(dm, simple->getDomainFEName(), feDomainLhs,
                                 null, null);
  if (feBoundaryLhs)
    CHKERR DMMoFEMTSSetIJacobian(dm, simple->getBoundaryFEName(), feBoundaryLhs,
                                 null, null);
  if (feSkeletonLhs)
    CHKERR DMMoFEMTSSetIJacobian(dm, simple->getSkeletonFEName(), feSkeletonLhs,
                                 null, null);
  if (feMeshsetLhs)
    CHKERR DMMoFEMTSSetIJacobian(dm, simple->getMeshsetFEName(), feMeshsetLhs,
                                 null, null);
 
  // Add element to calculate rhs of stiff part
  if (feDomainRhs)
    CHKERR DMMoFEMTSSetIFunction(dm, simple->getDomainFEName(), feDomainRhs,
                                 null, null);
  if (feBoundaryRhs)
    CHKERR DMMoFEMTSSetIFunction(dm, simple->getBoundaryFEName(), feBoundaryRhs,
                                 null, null);
  if (feSkeletonRhs)
    CHKERR DMMoFEMTSSetIFunction(dm, simple->getSkeletonFEName(), feSkeletonRhs,
                                 null, null);
  if (feMeshsetRhs)
    CHKERR DMMoFEMTSSetIFunction(dm, simple->getMeshsetFEName(), feMeshsetRhs,
                                 null, null);
 
  // Add element to calculate rhs of stiff part
  if (feDomainExplicitRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getDomainFEName(),
                                   feDomainExplicitRhs, null, null);
  if (feBoundaryExplicitRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getBoundaryFEName(),
                                   feBoundaryExplicitRhs, null, null);
  if (feSkeletonExplicitRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getSkeletonFEName(),
                                   feSkeletonExplicitRhs, null, null);
  if (feMeshsetExplicitRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getMeshsetFEName(),
                                   feMeshsetExplicitRhs, null, null);
 
  // Note: More cases for explict, and implicit time interation cases can be
  // implemented here.
 
  auto ts = MoFEM::createTS(m_field.get_comm());
  CHKERR TSSetDM(ts, dm);
  return ts;
}

getOpBoundaryExplicitRhsPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpBoundaryExplicitRhsPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Boundary Rhs Pipeline object for implicit-explicit G term.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Definition at line 1505 of file PipelineManager.hpp.

                                                  {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createBoundaryFEPipeline<DIM>(feBoundaryExplicitRhs))
      ->getOpPtrVector();
}

getOpBoundaryLhsPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpBoundaryLhsPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Boundary Lhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Examples

mofem/tutorials/clx-0_acoustics/helmholtz.cpp, and mofem/tutorials/scl-8_heat_radiation/radiation.cpp.

Definition at line 1313 of file PipelineManager.hpp.

                                          {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createBoundaryFEPipeline<DIM>(feBoundaryLhs))
      ->getOpPtrVector();
}

getOpBoundaryPostProcPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpBoundaryPostProcPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Boundary PostProc Pipeline object.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Examples

mofem/tutorials/scl-8_heat_radiation/radiation.cpp.

Definition at line 1385 of file PipelineManager.hpp.

                                               {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createBoundaryFEPipeline<DIM>(feBoundaryPostProc))
      ->getOpPtrVector();
}

getOpBoundaryRhsPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpBoundaryRhsPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Boundary Rhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/tutorials/clx-0_acoustics/helmholtz.cpp, and mofem/tutorials/scl-8_heat_radiation/radiation.cpp.

Definition at line 1337 of file PipelineManager.hpp.

                                          {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createBoundaryFEPipeline<DIM>(feBoundaryRhs))
      ->getOpPtrVector();
}

getOpDomainExplicitRhsPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpDomainExplicitRhsPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Domain Rhs Pipeline object for implicit-explicit G term.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Definition at line 1481 of file PipelineManager.hpp.

                                                {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createDomainFEPipeline<DIM>(feDomainExplicitRhs))
      ->getOpPtrVector();
}

getOpDomainLhsPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpDomainLhsPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Domain Lhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/tutorials/adv-3_level_set/level_set.cpp, mofem/tutorials/adv-6_dg_projection/between_meshes_dg_projection.cpp, mofem/tutorials/clx-0_acoustics/helmholtz.cpp, mofem/tutorials/mix-0_mixed_poisson/mixed_poisson.cpp, mofem/tutorials/scl-0_least_squares/approximation.cpp, mofem/tutorials/scl-8_heat_radiation/radiation.cpp, mofem/tutorials/scl-9_heat_method/heat_method.cpp, and mofem/tutorials/vec-4_shallow_wave/approx_sphere.cpp.

Definition at line 1265 of file PipelineManager.hpp.

                                        {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createDomainFEPipeline<DIM>(feDomainLhs))
      ->getOpPtrVector();
}

getOpDomainPostProcPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpDomainPostProcPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Domain PostProc Pipeline object.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Examples

mofem/tutorials/mix-0_mixed_poisson/mixed_poisson.cpp, and mofem/tutorials/scl-0_least_squares/approximation.cpp.

Definition at line 1361 of file PipelineManager.hpp.

                                             {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createDomainFEPipeline<DIM>(feDomainPostProc))
      ->getOpPtrVector();
}

getOpDomainRhsPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpDomainRhsPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Domain Rhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/tutorials/adv-6_dg_projection/between_meshes_dg_projection.cpp, mofem/tutorials/fun-0_hello_world/hello_world.cpp, mofem/tutorials/fun-1_integration/integration.cpp, mofem/tutorials/mix-0_mixed_poisson/mixed_poisson.cpp, mofem/tutorials/scl-0_least_squares/approximation.cpp, mofem/tutorials/scl-8_heat_radiation/radiation.cpp, mofem/tutorials/scl-9_heat_method/heat_method.cpp, and mofem/tutorials/vec-4_shallow_wave/approx_sphere.cpp.

Definition at line 1289 of file PipelineManager.hpp.

                                        {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createDomainFEPipeline<DIM>(feDomainRhs))
      ->getOpPtrVector();
}

getOpEvaluationPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpEvaluationPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Domain Evaluation Pipeline object.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Examples

mofem/atom_tests/higher_derivatives.cpp.

Definition at line 1409 of file PipelineManager.hpp.

                                         {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createDomainFEPipeline<DIM>(feEvaluation))
      ->getOpPtrVector();
}

getOpSkeletonExplicitRhsPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpSkeletonExplicitRhsPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Skeleton Rhs Pipeline object for implicit-explicit G term.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Definition at line 1529 of file PipelineManager.hpp.

                                                  {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createBoundaryFEPipeline<DIM>(feSkeletonExplicitRhs))
      ->getOpPtrVector();
}

getOpSkeletonLhsPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpSkeletonLhsPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Skeleton Lhs Pipeline object.

Returns

boost::ptr_deque<UserDataOperator>&

Examples

mofem/tutorials/adv-3_level_set/level_set.cpp.

Definition at line 1433 of file PipelineManager.hpp.

                                          {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createBoundaryFEPipeline<DIM>(feSkeletonLhs))
      ->getOpPtrVector();
}

getOpSkeletonRhsPipeline()

template<int DIM>

boost::ptr_deque< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpSkeletonRhsPipeline ( )

inline

#include <src/interfaces/PipelineManager.hpp>

Get the Op Skeleton Rhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_deque<UserDataOperator>&

Examples

mofem/atom_tests/simple_l2_only.cpp.

Definition at line 1457 of file PipelineManager.hpp.

                                          {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createBoundaryFEPipeline<DIM>(feSkeletonRhs))
      ->getOpPtrVector();
}

loopBoundaryPostProc()

MoFEMErrorCode MoFEM::PipelineManager::loopBoundaryPostProc

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Iterate boundary postprocessing finite elements only.

Parameters

dm	Optional DM object (default: nullptr)

Returns

MoFEMErrorCode Error code (0 on success)

Definition at line 147 of file PipelineManager.cpp.

                                                                         {
  MoFEMFunctionBegin;
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
  if (!dm)
    dm = simple->getDM();
  if (!feBoundaryPostProc) {
    MOFEM_LOG("WORLD", Sev::warning)
        << "PipelineManager::loopBoundaryPostProc called with no active finite "
           "element pipeline (boundary postprocessing is null)";
  }
 
  if (feBoundaryPostProc)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getBoundaryFEName(),
                                    feBoundaryPostProc);
 
  MoFEMFunctionReturn(0);
}

loopDomainPostProc()

MoFEMErrorCode MoFEM::PipelineManager::loopDomainPostProc

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Iterate domain postprocessing finite elements only.

Parameters

dm	Optional DM object (default: nullptr)

Returns

MoFEMErrorCode Error code (0 on success)

Definition at line 128 of file PipelineManager.cpp.

                                                                       {
  MoFEMFunctionBegin;
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
  if (!dm)
    dm = simple->getDM();
  if (!feDomainPostProc) {
    MOFEM_LOG("WORLD", Sev::warning)
        << "PipelineManager::loopDomainPostProc called with no active finite "
           "element pipeline (domain postprocessing is null)";
  }
 
  if (feDomainPostProc)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getDomainFEName(),
                                    feDomainPostProc);
 
  MoFEMFunctionReturn(0);
}

loopFiniteElements()

MoFEMErrorCode MoFEM::PipelineManager::loopFiniteElements

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Iterate finite elements.

Parameters

dm	Optional DM object (default: nullptr)

Returns

MoFEMErrorCode Error code (0 on success)

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/simple_l2_only.cpp, and mofem/tutorials/fun-1_integration/integration.cpp.

Definition at line 63 of file PipelineManager.cpp.

                                                                       {
  MoFEMFunctionBegin;
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
  if (!dm)
    dm = simple->getDM();
  if (!feDomainLhs && !feBoundaryLhs && !feSkeletonLhs && !feMeshsetLhs &&
      !feDomainRhs && !feBoundaryRhs && !feSkeletonRhs && !feMeshsetRhs) {
    MOFEM_LOG("WORLD", Sev::warning)
        << "PipelineManager::loopFiniteElements called with no active finite "
           "element pipelines (all are null)";
  }
 
  // Add element to calculate lhs of stiff part
  if (feDomainLhs)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getDomainFEName(), feDomainLhs);
  if (feBoundaryLhs)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getBoundaryFEName(),
                                    feBoundaryLhs);
  if (feSkeletonLhs)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getSkeletonFEName(),
                                    feSkeletonLhs);
  if (feMeshsetLhs)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getMeshsetFEName(),
                                    feMeshsetLhs);
 
  // Add element to calculate rhs of stiff part
  if (feDomainRhs)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getDomainFEName(), feDomainRhs);
  if (feBoundaryRhs)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getBoundaryFEName(),
                                    feBoundaryRhs);
  if (feSkeletonRhs)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getSkeletonFEName(),
                                    feSkeletonRhs);
  if (feMeshsetRhs)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getMeshsetFEName(),
                                    feMeshsetRhs);
 
  MoFEMFunctionReturn(0);
}

loopFiniteElementsEvaluation()

MoFEMErrorCode MoFEM::PipelineManager::loopFiniteElementsEvaluation

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Iterate evaluation finite elements.

Parameters

dm	Optional DM object (default: nullptr)

Returns

MoFEMErrorCode Error code (0 on success)

Examples

mofem/atom_tests/higher_derivatives.cpp.

Definition at line 167 of file PipelineManager.cpp.

                                                                  {
  MoFEMFunctionBegin;
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
  if (!dm)
    dm = simple->getDM();
  if (!feEvaluation) {
    MOFEM_LOG("WORLD", Sev::warning)
        << "PipelineManager::loopFiniteElementsEvaluation called with no "
           "active "
           "finite element pipelines (all are null)";
  }
 
  if (feEvaluation)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getDomainFEName(),
                                    feEvaluation);
 
  MoFEMFunctionReturn(0);
}

loopFiniteElementsPostProc()

MoFEMErrorCode MoFEM::PipelineManager::loopFiniteElementsPostProc

(

SmartPetscObj< DM >

dm = nullptr

)

#include <src/interfaces/PipelineManager.hpp>

Iterate postprocessing finite elements.

Parameters

dm	Optional DM object (default: nullptr)

Returns

MoFEMErrorCode Error code (0 on success)

Examples

mofem/tutorials/clx-0_acoustics/helmholtz.cpp, mofem/tutorials/fun-2_plot_base/plot_base.cpp, mofem/tutorials/mix-0_mixed_poisson/mixed_poisson.cpp, mofem/tutorials/mix-1_light_intensity_equation/phase.cpp, mofem/tutorials/scl-0_least_squares/approximation.cpp, mofem/tutorials/scl-1_poisson_2d_homogeneous/poisson_2d_homogeneous.cpp, mofem/tutorials/scl-8_heat_radiation/radiation.cpp, mofem/tutorials/scl-9_heat_method/heat_method.cpp, and mofem/tutorials/vec-6_dg_kirchhoff_plate/plate.cpp.

Definition at line 106 of file PipelineManager.cpp.

                                                                {
  MoFEMFunctionBegin;
  Interface &m_field = cOre;
  Simple *simple = m_field.getInterface<Simple>();
  if (!dm)
    dm = simple->getDM();
  if (!feDomainPostProc && !feBoundaryPostProc) {
    MOFEM_LOG("WORLD", Sev::warning)
        << "PipelineManager::loopFiniteElementsPostProc called with no active "
           "finite element pipelines (all are null)";
  }
 
  if (feDomainPostProc)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getDomainFEName(),
                                    feDomainPostProc);
  if (feBoundaryPostProc)
    CHKERR DMoFEMLoopFiniteElements(dm, simple->getBoundaryFEName(),
                                    feBoundaryPostProc);
 
  MoFEMFunctionReturn(0);
}