28static char help[] =
"...\n\n";
47int main(
int argc,
char *argv[]) {
63 DMType dm_name =
"DMMOFEM";
94 op_side_fe->doWorkRhsHook = [&](
DataOperator *op_ptr,
int side,
97 auto domain_op =
static_cast<DomainEleOp *
>(op_ptr);
101 <<
"Side element name [ " << count_side_fe <<
" ] "
102 << domain_op->getFEName();
110 auto side_fe = boost::make_shared<DomainEle>(m_field);
111 side_fe->getOpPtrVector().push_back(op_side_fe);
120 MOFEM_LOG(
"SELF", Sev::verbose) <<
"Element name [ " << count_fe
121 <<
" ] " << bdy_op->getFEName();
132 op_bdy_fe->doWorkRhsHook = do_work_rhs;
135 op_skeleton_fe->doWorkRhsHook = do_work_rhs;
145 MOFEM_LOG(
"SELF", Sev::inform) <<
"Number of elements " << count_fe;
147 <<
"Number of side elements " << count_side_fe;
151 "Wrong numbers of FEs");
152 if (count_side_fe != 24)
154 "Wrong numbers of side FEs");
ForcesAndSourcesCore::UserDataOperator UserDataOperator
MoFEM::FaceElementForcesAndSourcesCore FaceEle
ElementsAndOps< SPACE_DIM >::DomainEle DomainEle
ElementsAndOps< SPACE_DIM >::BoundaryEle BoundaryEle
MoFEM::EdgeElementForcesAndSourcesCore EdgeEle
#define CATCH_ERRORS
Catch errors.
@ AINSWORTH_LEGENDRE_BASE
Ainsworth Cole (Legendre) approx. base .
@ L2
field with C-1 continuity
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
@ MOFEM_ATOM_TEST_INVALID
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define CHKERR
Inline error check.
PetscErrorCode DMRegister_MoFEM(const char sname[])
Register MoFEM problem.
MoFEMErrorCode loopFiniteElements(SmartPetscObj< DM > dm=nullptr)
Iterate finite elements.
boost::ptr_vector< UserDataOperator > & getOpSkeletonRhsPipeline()
Get the Op Skeleton Rhs Pipeline object.
boost::ptr_vector< UserDataOperator > & getOpBoundaryRhsPipeline()
Get the Op Boundary Rhs Pipeline object.
#define MOFEM_LOG(channel, severity)
Log.
implementation of Data Operators for Forces and Sources
DeprecatedCoreInterface Interface
int main(int argc, char *argv[])
EntitiesFieldData::EntData EntData
DomainEle::UserDataOperator DomainEleOp
BoundaryEle::UserDataOperator BoundaryEleOp
static MoFEMErrorCode Initialize(int *argc, char ***args, const char file[], const char help[])
Initializes the MoFEM database PETSc, MOAB and MPI.
static MoFEMErrorCode Finalize()
Checks for options to be called at the conclusion of the program.
base operator to do operations at Gauss Pt. level
Deprecated interface functions.
default operator for EDGE element
Data on single entity (This is passed as argument to DataOperator::doWork)
default operator for TRI element
@ OPSPACE
operator do Work is execute on space data
PipelineManager interface.
Simple interface for fast problem set-up.
bool & getAddBoundaryFE()
Get the addSkeletonFE.
MoFEMErrorCode addDomainField(const std::string &name, const FieldSpace space, const FieldApproximationBase base, const FieldCoefficientsNumber nb_of_coefficients, const TagType tag_type=MB_TAG_SPARSE, const enum MoFEMTypes bh=MF_ZERO, int verb=-1)
Add field on domain.
MoFEMErrorCode getOptions()
get options
MoFEMErrorCode loadFile(const std::string options, const std::string mesh_file_name)
Load mesh file.
MoFEMErrorCode setFieldOrder(const std::string field_name, const int order, const Range *ents=NULL)
Set field order.
bool & getAddSkeletonFE()
Get the addSkeletonFE.
MoFEMErrorCode setUp(const PetscBool is_partitioned=PETSC_TRUE)
Setup problem.
const std::string getDomainFEName() const
Get the Domain FE Name.
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.