v0.10.0
Files | Classes | Functions
PipelineManager interface
MoFEM

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 interfaceMake basic interface, to speedup problem setup and analysts. See discussion here link to google groups
 

Classes

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

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...
 
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::createTS (SmartPetscObj< DM > dm=nullptr)
 Create TS (time) 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>
Examples
helmholtz.cpp, and lesson3_poisson.cpp.

Definition at line 68 of file PipelineManager.cpp.

68  {
69  Interface &m_field = cOre;
70  Simple *simple_interface = m_field.getInterface<Simple>();
71 
72  auto copy_dm_struture = [&](auto simple_dm) {
73  MPI_Comm comm;
74  CHKERR PetscObjectGetComm(getPetscObject(simple_dm.get()), &comm);
75  DMType type;
76  CHKERR DMGetType(simple_dm, &type);
77  dm = createSmartDM(comm, type);
78  CHKERR DMMoFEMDuplicateDMCtx(simple_interface->getDM(), dm);
79  return dm;
80  };
81 
82  if (!dm)
83  dm = copy_dm_struture(simple_interface->getDM());
84  else
85  dm = copy_dm_struture(dm);
86 
87  auto set_dm_section = [&](auto dm) {
88  MoFEMFunctionBegin;
89  PetscSection section;
90  CHKERR m_field.getInterface<ISManager>()->sectionCreate(
91  simple_interface->getProblemName(), &section);
92  CHKERR DMSetDefaultSection(dm, section);
93  CHKERR DMSetDefaultGlobalSection(dm, section);
94  CHKERR PetscSectionDestroy(&section);
95  MoFEMFunctionReturn(0);
96  };
97  CHKERR set_dm_section(dm);
98 
99  boost::shared_ptr<FEMethod> null;
100 
101  // Add element to calculate lhs of stiff part
102  if (feDomainLhs)
103  CHKERR DMMoFEMKSPSetComputeOperators(
104  dm, simple_interface->getDomainFEName(), feDomainLhs, null, null);
105  if (feBoundaryLhs)
106  CHKERR DMMoFEMKSPSetComputeOperators(
107  dm, simple_interface->getBoundaryFEName(), feBoundaryLhs, null, null);
108  if (feSkeletonLhs)
109  CHKERR DMMoFEMKSPSetComputeOperators(
110  dm, simple_interface->getSkeletonFEName(), feSkeletonLhs, null, null);
111 
112  // Add element to calculate rhs of stiff part
113  if (feDomainRhs)
114  CHKERR DMMoFEMKSPSetComputeRHS(dm, simple_interface->getDomainFEName(),
115  feDomainRhs, null, null);
116  if (feBoundaryRhs)
117  CHKERR DMMoFEMKSPSetComputeRHS(dm, simple_interface->getBoundaryFEName(),
118  feBoundaryRhs, null, null);
119  if (feSkeletonRhs)
120  CHKERR DMMoFEMKSPSetComputeRHS(dm, simple_interface->getSkeletonFEName(),
121  feSkeletonRhs, null, null);
122 
123  auto ksp = MoFEM::createKSP(m_field.get_comm());
124  CHKERR KSPSetDM(ksp, dm);
125  return ksp;
126 }
MoFEM::DMMoFEMKSPSetComputeRHS
PetscErrorCode DMMoFEMKSPSetComputeRHS(DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
Set compute operator for KSP solver via sub-matrix and IS.
Definition: DMMMoFEM.cpp:556
MoFEM::PipelineManager::feBoundaryLhs
boost::shared_ptr< FEMethod > feBoundaryLhs
Element to assemble LHS side by integrating boundary.
Definition: PipelineManager.hpp:193
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:485
MoFEM::PipelineManager::feDomainLhs
boost::shared_ptr< FEMethod > feDomainLhs
Element to assemble LHS side by integrating domain.
Definition: PipelineManager.hpp:189
MoFEM::PipelineManager::feSkeletonRhs
boost::shared_ptr< FEMethod > feSkeletonRhs
Element to assemble RHS side by integrating skeleton.
Definition: PipelineManager.hpp:195
MoFEM::DMMoFEMDuplicateDMCtx
PetscErrorCode DMMoFEMDuplicateDMCtx(DM dm, DM dm_duplicate)
Duplicate internal data struture.
Definition: DMMMoFEM.cpp:160
MoFEM::getPetscObject
PetscObject getPetscObject(T obj)
Definition: AuxPETSc.hpp:23
MoFEM::DMMoFEMKSPSetComputeOperators
PetscErrorCode DMMoFEMKSPSetComputeOperators(DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
Set KSP operators and push mofem finite element methods.
Definition: DMMMoFEM.cpp:597
MoFEM::PipelineManager::feDomainRhs
boost::shared_ptr< FEMethod > feDomainRhs
Element to assemble RHS side by integrating domain.
Definition: PipelineManager.hpp:187
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(const MOFEMuuid &uuid, IFACE *&iface) const
Get interface by uuid and return reference to pointer of interface.
Definition: UnknownInterface.hpp:130
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:604
MoFEM::createKSP
auto createKSP
Definition: AuxPETSc.hpp:289
MoFEM::Interface
DeprecatedCoreInterface Interface
Definition: Interface.hpp:1943
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:415
MoFEM::createSmartDM
auto createSmartDM
Creates smart DM object.
Definition: AuxPETSc.hpp:174
MoFEM::PipelineManager::feBoundaryRhs
boost::shared_ptr< FEMethod > feBoundaryRhs
Element to assemble RHS side by integrating boundary.
Definition: PipelineManager.hpp:191
MoFEM::PipelineManager::feSkeletonLhs
boost::shared_ptr< FEMethod > feSkeletonLhs
Element to assemble LHS side by integrating skeleton.
Definition: PipelineManager.hpp:197

◆ createSNES()

SmartPetscObj< SNES > MoFEM::PipelineManager::createSNES ( SmartPetscObj< DM >  dm = nullptr)

Create SNES (nonlinear) solver.

Parameters
dm
Returns
SmartPetscObj<SNES>

Definition at line 128 of file PipelineManager.cpp.

128  {
129  Interface &m_field = cOre;
130  Simple *simple_interface = m_field.getInterface<Simple>();
131 
132  auto copy_dm_struture = [&](auto simple_dm) {
133  MPI_Comm comm;
134  CHKERR PetscObjectGetComm(getPetscObject(simple_dm.get()), &comm);
135  DMType type;
136  CHKERR DMGetType(simple_dm, &type);
137  dm = createSmartDM(comm, type);
138  CHKERR DMMoFEMDuplicateDMCtx(simple_interface->getDM(), dm);
139  return dm;
140  };
141 
142  if (!dm)
143  dm = copy_dm_struture(simple_interface->getDM());
144  else
145  dm = copy_dm_struture(dm);
146 
147  auto set_dm_section = [&](auto dm) {
148  MoFEMFunctionBegin;
149  PetscSection section;
150  CHKERR m_field.getInterface<ISManager>()->sectionCreate(
151  simple_interface->getProblemName(), &section);
152  CHKERR DMSetDefaultSection(dm, section);
153  CHKERR DMSetDefaultGlobalSection(dm, section);
154  CHKERR PetscSectionDestroy(&section);
155  MoFEMFunctionReturn(0);
156  };
157  CHKERR set_dm_section(dm);
158 
159  boost::shared_ptr<FEMethod> null;
160 
161  // Add element to calculate lhs of stiff part
162  if (feDomainLhs)
163  CHKERR DMMoFEMSNESSetJacobian(dm, simple_interface->getDomainFEName(),
164  feDomainLhs, null, null);
165  if (feBoundaryLhs)
166  CHKERR DMMoFEMSNESSetJacobian(dm, simple_interface->getBoundaryFEName(),
167  feBoundaryLhs, null, null);
168  if (feSkeletonLhs)
169  CHKERR DMMoFEMSNESSetJacobian(dm, simple_interface->getSkeletonFEName(),
170  feSkeletonLhs, null, null);
171 
172  // Add element to calculate rhs of stiff part
173  if (feDomainRhs)
174  CHKERR DMMoFEMSNESSetFunction(dm, simple_interface->getDomainFEName(),
175  feDomainRhs, null, null);
176  if (feBoundaryRhs)
177  CHKERR DMMoFEMSNESSetFunction(dm, simple_interface->getBoundaryFEName(),
178  feBoundaryRhs, null, null);
179  if (feSkeletonRhs)
180  CHKERR DMMoFEMSNESSetFunction(dm, simple_interface->getSkeletonFEName(),
181  feSkeletonRhs, null, null);
182 
183  auto snes = MoFEM::createSNES(m_field.get_comm());
184  CHKERR SNESSetDM(snes, dm);
185  return snes;
186 }
MoFEM::DMMoFEMSNESSetFunction
PetscErrorCode DMMoFEMSNESSetFunction(DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
set SNES residual evaluation function
Definition: DMMMoFEM.cpp:637
MoFEM::PipelineManager::feBoundaryLhs
boost::shared_ptr< FEMethod > feBoundaryLhs
Element to assemble LHS side by integrating boundary.
Definition: PipelineManager.hpp:193
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:485
MoFEM::PipelineManager::feDomainLhs
boost::shared_ptr< FEMethod > feDomainLhs
Element to assemble LHS side by integrating domain.
Definition: PipelineManager.hpp:189
MoFEM::PipelineManager::feSkeletonRhs
boost::shared_ptr< FEMethod > feSkeletonRhs
Element to assemble RHS side by integrating skeleton.
Definition: PipelineManager.hpp:195
MoFEM::DMMoFEMDuplicateDMCtx
PetscErrorCode DMMoFEMDuplicateDMCtx(DM dm, DM dm_duplicate)
Duplicate internal data struture.
Definition: DMMMoFEM.cpp:160
MoFEM::createSNES
auto createSNES
Definition: AuxPETSc.hpp:282
MoFEM::getPetscObject
PetscObject getPetscObject(T obj)
Definition: AuxPETSc.hpp:23
MoFEM::PipelineManager::feDomainRhs
boost::shared_ptr< FEMethod > feDomainRhs
Element to assemble RHS side by integrating domain.
Definition: PipelineManager.hpp:187
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(const MOFEMuuid &uuid, IFACE *&iface) const
Get interface by uuid and return reference to pointer of interface.
Definition: UnknownInterface.hpp:130
MoFEM::DMMoFEMSNESSetJacobian
PetscErrorCode DMMoFEMSNESSetJacobian(DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
set SNES Jacobian evaluation function
Definition: DMMMoFEM.cpp:678
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:604
MoFEM::Interface
DeprecatedCoreInterface Interface
Definition: Interface.hpp:1943
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:415
MoFEM::createSmartDM
auto createSmartDM
Creates smart DM object.
Definition: AuxPETSc.hpp:174
MoFEM::PipelineManager::feBoundaryRhs
boost::shared_ptr< FEMethod > feBoundaryRhs
Element to assemble RHS side by integrating boundary.
Definition: PipelineManager.hpp:191
MoFEM::PipelineManager::feSkeletonLhs
boost::shared_ptr< FEMethod > feSkeletonLhs
Element to assemble LHS side by integrating skeleton.
Definition: PipelineManager.hpp:197

◆ createTS()

SmartPetscObj< TS > MoFEM::PipelineManager::createTS ( SmartPetscObj< DM >  dm = nullptr)

Create TS (time) solver.

Parameters
dm
Returns
SmartPetscObj<TS>
Examples
lesson6_radiation.cpp.

Definition at line 188 of file PipelineManager.cpp.

188  {
189  Interface &m_field = cOre;
190  Simple *simple_interface = m_field.getInterface<Simple>();
191 
192  auto copy_dm_struture = [&](auto simple_dm) {
193  MPI_Comm comm;
194  CHKERR PetscObjectGetComm(getPetscObject(simple_dm.get()), &comm);
195  DMType type;
196  CHKERR DMGetType(simple_dm, &type);
197  dm = createSmartDM(comm, type);
198  CHKERR DMMoFEMDuplicateDMCtx(simple_interface->getDM(), dm);
199  return dm;
200  };
201 
202  if (!dm)
203  dm = copy_dm_struture(simple_interface->getDM());
204  else
205  dm = copy_dm_struture(dm);
206 
207  auto set_dm_section = [&](auto dm) {
208  MoFEMFunctionBegin;
209  PetscSection section;
210  CHKERR m_field.getInterface<ISManager>()->sectionCreate(
211  simple_interface->getProblemName(), &section);
212  CHKERR DMSetDefaultSection(dm, section);
213  CHKERR DMSetDefaultGlobalSection(dm, section);
214  CHKERR PetscSectionDestroy(&section);
215  MoFEMFunctionReturn(0);
216  };
217  CHKERR set_dm_section(dm);
218 
219  boost::shared_ptr<FEMethod> null;
220 
221  // Add element to calculate lhs of stiff part
222  if (feDomainLhs)
223  CHKERR DMMoFEMTSSetIJacobian(dm, simple_interface->getDomainFEName(),
224  feDomainLhs, null, null);
225  if (feBoundaryLhs)
226  CHKERR DMMoFEMTSSetIJacobian(dm, simple_interface->getBoundaryFEName(),
227  feBoundaryLhs, null, null);
228  if (feSkeletonLhs)
229  CHKERR DMMoFEMTSSetIJacobian(dm, simple_interface->getSkeletonFEName(),
230  feSkeletonLhs, null, null);
231 
232  // Add element to calculate rhs of stiff part
233  if (feDomainRhs)
234  CHKERR DMMoFEMTSSetIFunction(dm, simple_interface->getDomainFEName(),
235  feDomainRhs, null, null);
236  if (feBoundaryRhs)
237  CHKERR DMMoFEMTSSetIFunction(dm, simple_interface->getBoundaryFEName(),
238  feBoundaryRhs, null, null);
239  if (feSkeletonRhs)
240  CHKERR DMMoFEMTSSetIFunction(dm, simple_interface->getSkeletonFEName(),
241  feSkeletonRhs, null, null);
242 
243  // Note: More cases for explit, and implicit time ingeration cases can be
244  // implemented here.
245 
246  auto ts = MoFEM::createTS(m_field.get_comm());
247  CHKERR TSSetDM(ts, dm);
248  return ts;
249 }
MoFEM::PipelineManager::feBoundaryLhs
boost::shared_ptr< FEMethod > feBoundaryLhs
Element to assemble LHS side by integrating boundary.
Definition: PipelineManager.hpp:193
MoFEM::createTS
auto createTS
Definition: AuxPETSc.hpp:275
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:485
MoFEM::PipelineManager::feDomainLhs
boost::shared_ptr< FEMethod > feDomainLhs
Element to assemble LHS side by integrating domain.
Definition: PipelineManager.hpp:189
MoFEM::PipelineManager::feSkeletonRhs
boost::shared_ptr< FEMethod > feSkeletonRhs
Element to assemble RHS side by integrating skeleton.
Definition: PipelineManager.hpp:195
MoFEM::DMMoFEMDuplicateDMCtx
PetscErrorCode DMMoFEMDuplicateDMCtx(DM dm, DM dm_duplicate)
Duplicate internal data struture.
Definition: DMMMoFEM.cpp:160
MoFEM::getPetscObject
PetscObject getPetscObject(T obj)
Definition: AuxPETSc.hpp:23
MoFEM::DMMoFEMTSSetIJacobian
PetscErrorCode DMMoFEMTSSetIJacobian(DM dm, const std::string fe_name, boost::shared_ptr< MoFEM::FEMethod > method, boost::shared_ptr< MoFEM::BasicMethod > pre_only, boost::shared_ptr< MoFEM::BasicMethod > post_only)
set TS Jacobian evaluation function
Definition: DMMMoFEM.cpp:772
MoFEM::PipelineManager::feDomainRhs
boost::shared_ptr< FEMethod > feDomainRhs
Element to assemble RHS side by integrating domain.
Definition: PipelineManager.hpp:187
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(const MOFEMuuid &uuid, IFACE *&iface) const
Get interface by uuid and return reference to pointer of interface.
Definition: UnknownInterface.hpp:130
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:604
MoFEM::Interface
DeprecatedCoreInterface Interface
Definition: Interface.hpp:1943
MoFEM::DMMoFEMTSSetIFunction
PetscErrorCode DMMoFEMTSSetIFunction(DM dm, const char fe_name[], MoFEM::FEMethod *method, MoFEM::BasicMethod *pre_only, MoFEM::BasicMethod *post_only)
set TS implicit function evaluation function
Definition: DMMMoFEM.cpp:719
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:415
MoFEM::createSmartDM
auto createSmartDM
Creates smart DM object.
Definition: AuxPETSc.hpp:174
MoFEM::PipelineManager::feBoundaryRhs
boost::shared_ptr< FEMethod > feBoundaryRhs
Element to assemble RHS side by integrating boundary.
Definition: PipelineManager.hpp:191
MoFEM::PipelineManager::feSkeletonLhs
boost::shared_ptr< FEMethod > feSkeletonLhs
Element to assemble LHS side by integrating skeleton.
Definition: PipelineManager.hpp:197

◆ 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>&
Examples
helmholtz.cpp, lesson3_poisson.cpp, lesson6_radiation.cpp, and lesson8_contact.cpp.

Definition at line 490 of file PipelineManager.hpp.

490  {
491  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
492  createBoundaryFEPipeline<DIM>(feBoundaryLhs))
493  ->getOpPtrVector();
494 }
MoFEM::PipelineManager::feBoundaryLhs
boost::shared_ptr< FEMethod > feBoundaryLhs
Element to assemble LHS side by integrating boundary.
Definition: PipelineManager.hpp:193

◆ 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>&
Examples
helmholtz.cpp, lesson3_poisson.cpp, lesson6_radiation.cpp, lesson7_plastic.cpp, and lesson8_contact.cpp.

Definition at line 514 of file PipelineManager.hpp.

514  {
515  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
516  createBoundaryFEPipeline<DIM>(feBoundaryRhs))
517  ->getOpPtrVector();
518 }
MoFEM::PipelineManager::feBoundaryRhs
boost::shared_ptr< FEMethod > feBoundaryRhs
Element to assemble RHS side by integrating boundary.
Definition: PipelineManager.hpp:191

◆ 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>&
Examples
helmholtz.cpp, lesson3_poisson.cpp, lesson6_radiation.cpp, lesson7_plastic.cpp, and lesson8_contact.cpp.

Definition at line 442 of file PipelineManager.hpp.

442  {
443  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
444  createDomainFEPipeline<DIM>(feDomainLhs))
445  ->getOpPtrVector();
446 }
MoFEM::PipelineManager::feDomainLhs
boost::shared_ptr< FEMethod > feDomainLhs
Element to assemble LHS side by integrating domain.
Definition: PipelineManager.hpp:189

◆ 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>&
Examples
hello_world.cpp, lesson3_poisson.cpp, lesson6_radiation.cpp, lesson7_plastic.cpp, and lesson8_contact.cpp.

Definition at line 466 of file PipelineManager.hpp.

466  {
467  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
468  createDomainFEPipeline<DIM>(feDomainRhs))
469  ->getOpPtrVector();
470 }
MoFEM::PipelineManager::feDomainRhs
boost::shared_ptr< FEMethod > feDomainRhs
Element to assemble RHS side by integrating domain.
Definition: PipelineManager.hpp:187

◆ 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 538 of file PipelineManager.hpp.

538  {
539  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
540  createBoundaryFEPipeline<DIM>(feSkeletonRhs))
541  ->getOpPtrVector();
542 }
MoFEM::PipelineManager::feSkeletonRhs
boost::shared_ptr< FEMethod > feSkeletonRhs
Element to assemble RHS side by integrating skeleton.
Definition: PipelineManager.hpp:195

◆ 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 562 of file PipelineManager.hpp.

562  {
563  return boost::dynamic_pointer_cast<ForcesAndSourcesCore>(
564  createBoundaryFEPipeline<DIM>(feSkeletonLhs))
565  ->getOpPtrVector();
566 }
MoFEM::PipelineManager::feSkeletonLhs
boost::shared_ptr< FEMethod > feSkeletonLhs
Element to assemble LHS side by integrating skeleton.
Definition: PipelineManager.hpp:197

◆ loopFiniteElements()

MoFEMErrorCode MoFEM::PipelineManager::loopFiniteElements ( SmartPetscObj< DM >  dm = nullptr)

Iterate finite elements.

Returns
MoFEMErrorCode
Examples
helmholtz.cpp, lesson3_poisson.cpp, and lesson6_radiation.cpp.

Definition at line 36 of file PipelineManager.cpp.

36  {
37  MoFEMFunctionBegin;
38  Interface &m_field = cOre;
39  Simple *simple_interface = m_field.getInterface<Simple>();
40  if (!dm)
41  dm = simple_interface->getDM();
42 
43  // Add element to calculate lhs of stiff part
44  if (feDomainLhs)
45  CHKERR DMoFEMLoopFiniteElements(dm, simple_interface->getDomainFEName(),
46  feDomainLhs);
47  if (feBoundaryLhs)
48  CHKERR DMoFEMLoopFiniteElements(dm, simple_interface->getBoundaryFEName(),
49  feBoundaryLhs);
50  if (feSkeletonLhs)
51  CHKERR DMoFEMLoopFiniteElements(dm, simple_interface->getSkeletonFEName(),
52  feSkeletonLhs);
53 
54  // Add element to calculate rhs of stiff part
55  if (feDomainRhs)
56  CHKERR DMoFEMLoopFiniteElements(dm, simple_interface->getDomainFEName(),
57  feDomainRhs);
58  if (feBoundaryRhs)
59  CHKERR DMoFEMLoopFiniteElements(dm, simple_interface->getBoundaryFEName(),
60  feBoundaryRhs);
61  if (feSkeletonRhs)
62  CHKERR DMoFEMLoopFiniteElements(dm, simple_interface->getSkeletonFEName(),
63  feSkeletonRhs);
64 
65  MoFEMFunctionReturn(0);
66 }
MoFEM::PipelineManager::feBoundaryLhs
boost::shared_ptr< FEMethod > feBoundaryLhs
Element to assemble LHS side by integrating boundary.
Definition: PipelineManager.hpp:193
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:485
MoFEM::PipelineManager::feDomainLhs
boost::shared_ptr< FEMethod > feDomainLhs
Element to assemble LHS side by integrating domain.
Definition: PipelineManager.hpp:189
MoFEM::PipelineManager::feSkeletonRhs
boost::shared_ptr< FEMethod > feSkeletonRhs
Element to assemble RHS side by integrating skeleton.
Definition: PipelineManager.hpp:195
MoFEM::DMoFEMLoopFiniteElements
PetscErrorCode DMoFEMLoopFiniteElements(DM dm, const char fe_name[], MoFEM::FEMethod *method)
Executes FEMethod for finite elements in DM.
Definition: DMMMoFEM.cpp:507
MoFEM::PipelineManager::feDomainRhs
boost::shared_ptr< FEMethod > feDomainRhs
Element to assemble RHS side by integrating domain.
Definition: PipelineManager.hpp:187
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(const MOFEMuuid &uuid, IFACE *&iface) const
Get interface by uuid and return reference to pointer of interface.
Definition: UnknownInterface.hpp:130
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:604
MoFEM::Interface
DeprecatedCoreInterface Interface
Definition: Interface.hpp:1943
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:415
MoFEM::PipelineManager::feBoundaryRhs
boost::shared_ptr< FEMethod > feBoundaryRhs
Element to assemble RHS side by integrating boundary.
Definition: PipelineManager.hpp:191
MoFEM::PipelineManager::feSkeletonLhs
boost::shared_ptr< FEMethod > feSkeletonLhs
Element to assemble LHS side by integrating skeleton.
Definition: PipelineManager.hpp:197
Generated on Sun Sep 27 2020 23:34:55 for MoFEM by Doxygen 1.8.15 and hosted at University of Glasgow