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::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.
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/clx-0/helmholtz.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/mix-1/phase.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, and mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp.

Definition at line 99 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 165 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 254 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/seepage.cpp, and mofem/tutorials/scl-8/radiation.cpp.

Definition at line 309 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 376 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 443 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 1245 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/helmholtz.cpp, and mofem/tutorials/scl-8/radiation.cpp.

Definition at line 1125 of file PipelineManager.hpp.

                                          {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createBoundaryFEPipeline<DIM>(feBoundaryLhs))
      ->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/helmholtz.cpp, and mofem/tutorials/scl-8/radiation.cpp.

Definition at line 1149 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 1221 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.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-8/radiation.cpp, mofem/tutorials/scl-9/heat_method.cpp, and mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 1077 of file PipelineManager.hpp.

                                        {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createDomainFEPipeline<DIM>(feDomainLhs))
      ->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/fun-0/hello_world.cpp, mofem/tutorials/fun-1/integration.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-8/radiation.cpp, mofem/tutorials/scl-9/heat_method.cpp, and mofem/tutorials/vec-4/approx_sphere.cpp.

Definition at line 1101 of file PipelineManager.hpp.

                                        {
  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
             createDomainFEPipeline<DIM>(feDomainRhs))
      ->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 1269 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.cpp.

Definition at line 1173 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 1197 of file PipelineManager.hpp.

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

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/higher_derivatives.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/fun-1/integration.cpp, mofem/tutorials/fun-2/plot_base.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/mix-1/phase.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-8/radiation.cpp, and mofem/tutorials/scl-9/heat_method.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();
 
  // 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);
}