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.

Functions
template<int DIM = -1>
boost::ptr_vector< UserDataOperator > &	MoFEM::PipelineManager::getOpDomainLhsPipeline ()
	Get the Op Domain Lhs Pipeline object. More...
template<int DIM = -1>
boost::ptr_vector< UserDataOperator > &	MoFEM::PipelineManager::getOpDomainRhsPipeline ()
	Get the Op Domain Rhs Pipeline object. More...
template<int DIM = -1>
boost::ptr_vector< UserDataOperator > &	MoFEM::PipelineManager::getOpBoundaryLhsPipeline ()
	Get the Op Boundary Lhs Pipeline object. More...
template<int DIM = -1>
boost::ptr_vector< UserDataOperator > &	MoFEM::PipelineManager::getOpBoundaryRhsPipeline ()
	Get the Op Boundary Rhs Pipeline object. More...
template<int DIM = -1>
boost::ptr_vector< UserDataOperator > &	MoFEM::PipelineManager::getOpSkeletonLhsPipeline ()
	Get the Op Skeleton Lhs Pipeline object. More...
template<int DIM = -1>
boost::ptr_vector< UserDataOperator > &	MoFEM::PipelineManager::getOpSkeletonRhsPipeline ()
	Get the Op Skeleton Rhs Pipeline object. More...
template<int DIM = -1>
boost::ptr_vector< UserDataOperator > &	MoFEM::PipelineManager::getOpDomainExplicitRhsPipeline ()
	Get the Op Domain Rhs Pipeline object for implicit-explicit G term. More...
template<int DIM = -1>
boost::ptr_vector< UserDataOperator > &	MoFEM::PipelineManager::getOpBoundaryExplicitRhsPipeline ()
	Get the Op Bondary Rhs Pipeline object for implicit-explicit G term. More...
template<int DIM = -1>
boost::ptr_vector< UserDataOperator > &	MoFEM::PipelineManager::getOpSkeletonExplicitRhsPipeline ()
	Get the Op Skeleton Rhs Pipeline object for implicit-explicit G term. More...
MoFEMErrorCode	MoFEM::PipelineManager::loopFiniteElements (SmartPetscObj< DM > dm=nullptr)
	Iterate finite elements. More...
SmartPetscObj< KSP >	MoFEM::PipelineManager::createKSP (SmartPetscObj< DM > dm=nullptr)
	Create KSP (linear) solver. More...
SmartPetscObj< SNES >	MoFEM::PipelineManager::createSNES (SmartPetscObj< DM > dm=nullptr)
	Create SNES (nonlinear) solver. More...
SmartPetscObj< TS >	MoFEM::PipelineManager::createTSEX (SmartPetscObj< DM > dm=nullptr)
	Create TS (time) explit solver. More...
SmartPetscObj< TS >	MoFEM::PipelineManager::createTSIM (SmartPetscObj< DM > dm=nullptr)
	Create TS (time) implicit solver. More...
SmartPetscObj< TS >	MoFEM::PipelineManager::createTSIM2 (SmartPetscObj< DM > dm=nullptr)
	Create TS (time) solver for second order equation in time. More...
SmartPetscObj< TS >	MoFEM::PipelineManager::createTSIMEX (SmartPetscObj< DM > dm=nullptr)
	Create TS (time) implicit-explicit solver. More...

Detailed Description

Implementation of basic interface for rapid problem implementation.

Function Documentation

createKSP()

SmartPetscObj< KSP > MoFEM::PipelineManager::createKSP

(

SmartPetscObj< DM >

dm = nullptr

)

Create KSP (linear) solver.

Parameters

dm

Returns

SmartPetscObj<KSP>

Definition at line 62 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 = createSmartDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
  else
    dm = copy_dm_struture(dm);
 
  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;
 
  smartGetDMKspCtx(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);
 
  // 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);
 
  auto ksp = MoFEM::createKSP(m_field.get_comm());
  CHKERR KSPSetDM(ksp, dm);
  return ksp;
}

createSNES()

SmartPetscObj< SNES > MoFEM::PipelineManager::createSNES

(

SmartPetscObj< DM >

dm = nullptr

)

Create SNES (nonlinear) solver.

Parameters

dm

Returns

SmartPetscObj<SNES>

Definition at line 124 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 = createSmartDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
  else
    dm = copy_dm_struture(dm);
 
  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);
 
  smartGetDMSnesCtx(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);
 
  // 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);
 
  auto snes = MoFEM::createSNES(m_field.get_comm());
  CHKERR SNESSetDM(snes, dm);
  return snes;
}

createTSEX()

SmartPetscObj< TS > MoFEM::PipelineManager::createTSEX

(

SmartPetscObj< DM >

dm = nullptr

)

Create TS (time) explit solver.

Parameters

dm

Returns

SmartPetscObj<TS>

Definition at line 209 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 = createSmartDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
  else
    dm = copy_dm_struture(dm);
 
  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;
 
  smartGetDMTsCtx(dm)->clearLoops();
 
  // Add element to calculate rhs of stiff part
  if (feDomainRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getDomainFEName(),
                                   feDomainExplicitRhs, null, null);
  if (feBoundaryRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getBoundaryFEName(),
                                   feBoundaryExplicitRhs, null, null);
  if (feSkeletonExplicitRhs)
    CHKERR DMMoFEMTSSetRHSFunction(dm, simple->getSkeletonFEName(),
                                   feSkeletonExplicitRhs, null, null);
 
  // Note: More cases for explit, and implicit time ingeration 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

)

Create TS (time) implicit solver.

Parameters

dm

Returns

SmartPetscObj<TS>

Definition at line 263 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 = createSmartDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
  else
    dm = copy_dm_struture(dm);
 
  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);
 
  // 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);
 
  // Note: More cases for explit, and implicit time ingeration 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

)

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

Parameters

dm

Returns

SmartPetscObj<TS>

Definition at line 326 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 = createSmartDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
  else
    dm = copy_dm_struture(dm);
 
  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);
 
  // 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);
 
  // Note: More cases for explit, and implicit time ingeration 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

)

Create TS (time) implicit-explicit solver.

Parameters

dm

Returns

SmartPetscObj<TS>

Definition at line 389 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 = createSmartDM(comm, type);
    CHKERR DMMoFEMDuplicateDMCtx(simple_dm, dm);
    return dm;
  };
 
  if (!dm)
    dm = copy_dm_struture(simple->getDM());
  else
    dm = copy_dm_struture(dm);
 
  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);
 
  // 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);
 
  // 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);
 
  // Note: More cases for explit, and implicit time ingeration 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_vector< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpBoundaryExplicitRhsPipeline

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

Template Parameters

-1

Returns

boost::ptr_vector<UserDataOperator>&

Definition at line 819 of file PipelineManager.hpp.

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

getOpBoundaryLhsPipeline()

template<int DIM>

boost::ptr_vector< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpBoundaryLhsPipeline

Get the Op Boundary Lhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_vector<UserDataOperator>&

Definition at line 699 of file PipelineManager.hpp.

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

getOpBoundaryRhsPipeline()

template<int DIM>

boost::ptr_vector< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpBoundaryRhsPipeline

Get the Op Boundary Rhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_vector<UserDataOperator>&

Definition at line 723 of file PipelineManager.hpp.

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

getOpDomainExplicitRhsPipeline()

template<int DIM>

boost::ptr_vector< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpDomainExplicitRhsPipeline

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

Template Parameters

-1

Returns

boost::ptr_vector<UserDataOperator>&

Definition at line 795 of file PipelineManager.hpp.

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

getOpDomainLhsPipeline()

template<int DIM>

boost::ptr_vector< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpDomainLhsPipeline

Get the Op Domain Lhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_vector<UserDataOperator>&

Definition at line 651 of file PipelineManager.hpp.

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

getOpDomainRhsPipeline()

template<int DIM>

boost::ptr_vector< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpDomainRhsPipeline

Get the Op Domain Rhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_vector<UserDataOperator>&

Definition at line 675 of file PipelineManager.hpp.

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

getOpSkeletonExplicitRhsPipeline()

template<int DIM>

boost::ptr_vector< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpSkeletonExplicitRhsPipeline

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

Template Parameters

-1

Returns

boost::ptr_vector<UserDataOperator>&

Definition at line 843 of file PipelineManager.hpp.

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

getOpSkeletonLhsPipeline()

template<int DIM>

boost::ptr_vector< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpSkeletonLhsPipeline

Get the Op Skeleton Lhs Pipeline object.

Returns

boost::ptr_vector<UserDataOperator>&

Definition at line 747 of file PipelineManager.hpp.

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

getOpSkeletonRhsPipeline()

template<int DIM>

boost::ptr_vector< PipelineManager::UserDataOperator > & MoFEM::PipelineManager::getOpSkeletonRhsPipeline

Get the Op Skeleton Rhs Pipeline object.

Template Parameters

-1

Returns

boost::ptr_vector<UserDataOperator>&

Definition at line 771 of file PipelineManager.hpp.

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

loopFiniteElements()

MoFEMErrorCode MoFEM::PipelineManager::loopFiniteElements

(

SmartPetscObj< DM >

dm = nullptr

)

Iterate finite elements.

Returns

MoFEMErrorCode

Definition at line 32 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);
 
  // 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);
 
  MoFEMFunctionReturn(0);
}