Test iterating over boundary and skeleton elements only when L2 field is presents on the domain.
static char help[] =
"...\n\n";
};
int main(
int argc,
char *argv[]) {
try {
moab::Core moab_core;
moab::Interface &moab = moab_core;
DMType dm_name = "DMMOFEM";
{
int count_fe;
int count_side_fe;
op_side_fe->doWorkRhsHook = [&](
DataOperator *op_ptr,
int side,
<< "Side element name [ " << count_side_fe << " ] "
<< domain_op->getFEName();
++count_side_fe;
};
auto side_fe = boost::make_shared<DomainEle>(m_field);
side_fe->getOpPtrVector().push_back(op_side_fe);
MOFEM_LOG(
"SELF", Sev::verbose) <<
"Element name [ " << count_fe
<< " ] " << bdy_op->getFEName();
++count_fe;
};
op_bdy_fe->doWorkRhsHook = do_work_rhs;
op_skeleton_fe->doWorkRhsHook = do_work_rhs;
count_fe = 0;
count_side_fe = 0;
MOFEM_LOG(
"SELF", Sev::inform) <<
"Number of elements " << count_fe;
<< "Number of side elements " << count_side_fe;
if (count_fe != 16)
"Wrong numbers of FEs");
if (count_side_fe != 24)
"Wrong numbers of side FEs");
}
}
return 0;
}
DomainEle::UserDataOperator DomainEleOp
Finire element operator type.
ElementsAndOps< SPACE_DIM >::DomainEle DomainEle
ElementsAndOps< SPACE_DIM >::BoundaryEle BoundaryEle
BoundaryEle::UserDataOperator BoundaryEleOp
#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_deque< UserDataOperator > & getOpBoundaryRhsPipeline()
Get the Op Boundary Rhs Pipeline object.
boost::ptr_deque< UserDataOperator > & getOpSkeletonRhsPipeline()
Get the Op Skeleton Rhs Pipeline object.
#define MOFEM_LOG(channel, severity)
Log.
implementation of Data Operators for Forces and Sources
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.
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 loadFile(const std::string options, const std::string mesh_file_name, LoadFileFunc loadFunc=defaultLoadFileFunc)
Load mesh file.
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 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.