v0.14.0
Loading...
Searching...
No Matches
simple_l2_only.cpp
Go to the documentation of this file.
1/**
2 * \file simple_l2_only.cpp
3 * \ingroup mofem_simple_interface
4 * \example simple_l2_only.cpp
5 *
6 * Test iterating over boundary and skeleton elements only when L2 field is
7 * presents on the domain.
8 *
9 */
10
11
12
13#include <MoFEM.hpp>
14using namespace MoFEM;
15
16static char help[] = "...\n\n";
17
18constexpr int SPACE_DIM = 2;
19
20template <int DIM> struct ElementsAndOps {};
21
22template <> struct ElementsAndOps<2> {
27};
28
34
35int main(int argc, char *argv[]) {
36
37 // initialize petsc
38 MoFEM::Core::Initialize(&argc, &argv, (char *)0, help);
39
40 try {
41
42 // Create MoAB database
43 moab::Core moab_core;
44 moab::Interface &moab = moab_core;
45
46 // Create MoFEM database and link it to MoAB
47 MoFEM::Core mofem_core(moab);
48 MoFEM::Interface &m_field = mofem_core;
49
50 // Register DM Manager
51 DMType dm_name = "DMMOFEM";
52 CHKERR DMRegister_MoFEM(dm_name);
53
54 // Simple interface
55 Simple *simple_interface;
56 CHKERR m_field.getInterface(simple_interface);
57 {
58
59 // get options from command line
60 CHKERR simple_interface->getOptions();
61 // load mesh file
62 CHKERR simple_interface->loadFile();
63
64 CHKERR simple_interface->addDomainField("FIELD", L2,
66
67 simple_interface->getAddBoundaryFE() = true;
68 simple_interface->getAddSkeletonFE() = true;
69
70 // set fields order
71 CHKERR simple_interface->setFieldOrder("FIELD", 1);
72 // setup problem
73 CHKERR simple_interface->setUp();
74
75 int count_fe;
76 int count_side_fe;
77
78 PipelineManager *pipeline_mng = m_field.getInterface<PipelineManager>();
79
80 // Create OP for side FE
81 auto op_side_fe = new DomainEleOp(NOSPACE, DomainEleOp::OPSPACE);
82 op_side_fe->doWorkRhsHook = [&](DataOperator *op_ptr, int side,
83 EntityType type,
85 auto domain_op = static_cast<DomainEleOp *>(op_ptr);
87
88 MOFEM_LOG("SELF", Sev::verbose)
89 << "Side element name [ " << count_side_fe << " ] "
90 << domain_op->getFEName();
91
92 ++count_side_fe;
93
95 };
96
97 // Create side FE
98 auto side_fe = boost::make_shared<DomainEle>(m_field);
99 side_fe->getOpPtrVector().push_back(op_side_fe);
100
101 // Create boundary FE OP
102
103 auto do_work_rhs = [&](DataOperator *op_ptr, int side, EntityType type,
105 auto bdy_op = static_cast<BoundaryEleOp *>(op_ptr);
107
108 MOFEM_LOG("SELF", Sev::verbose) << "Element name [ " << count_fe
109 << " ] " << bdy_op->getFEName();
110
111 CHKERR bdy_op->loopSide(simple_interface->getDomainFEName(),
112 side_fe.get(), SPACE_DIM);
113
114 ++count_fe;
115
117 };
118
119 auto op_bdy_fe = new BoundaryEleOp(NOSPACE, DomainEleOp::OPSPACE);
120 op_bdy_fe->doWorkRhsHook = do_work_rhs;
121
122 auto op_skeleton_fe = new BoundaryEleOp(NOSPACE, DomainEleOp::OPSPACE);
123 op_skeleton_fe->doWorkRhsHook = do_work_rhs;
124
125 // Count boundary
126 count_fe = 0;
127 count_side_fe = 0;
128
129 pipeline_mng->getOpBoundaryRhsPipeline().push_back(op_bdy_fe);
130 pipeline_mng->getOpSkeletonRhsPipeline().push_back(op_skeleton_fe);
131 pipeline_mng->loopFiniteElements();
132
133 MOFEM_LOG("SELF", Sev::inform) << "Number of elements " << count_fe;
134 MOFEM_LOG("SELF", Sev::inform)
135 << "Number of side elements " << count_side_fe;
136
137 if (count_fe != 16)
138 SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
139 "Wrong numbers of FEs");
140 if (count_side_fe != 24)
141 SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
142 "Wrong numbers of side FEs");
143 }
144 }
146
147 // finish work cleaning memory, getting statistics, etc.
149
150 return 0;
151}
int main()
ElementsAndOps< SPACE_DIM >::DomainEle DomainEle
ElementsAndOps< SPACE_DIM >::BoundaryEle BoundaryEle
#define CATCH_ERRORS
Catch errors.
@ AINSWORTH_LEGENDRE_BASE
Ainsworth Cole (Legendre) approx. base .
Definition definitions.h:60
@ L2
field with C-1 continuity
Definition definitions.h:88
@ NOSPACE
Definition definitions.h:83
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
@ MOFEM_ATOM_TEST_INVALID
Definition definitions.h:40
#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.
Definition DMMoFEM.cpp:47
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
Definition Common.hpp:10
static char help[]
constexpr int SPACE_DIM
DomainEle::UserDataOperator DomainEleOp
BoundaryEle::UserDataOperator BoundaryEleOp
Core (interface) class.
Definition Core.hpp:82
static MoFEMErrorCode Initialize(int *argc, char ***args, const char file[], const char help[])
Initializes the MoFEM database PETSc, MOAB and MPI.
Definition Core.cpp:72
static MoFEMErrorCode Finalize()
Checks for options to be called at the conclusion of the program.
Definition Core.cpp:112
base operator to do operations at Gauss Pt. level
Deprecated interface functions.
Data on single entity (This is passed as argument to DataOperator::doWork)
@ OPSPACE
operator do Work is execute on space data
PipelineManager interface.
MoFEM::FaceElementForcesAndSourcesCore FaceEle
MoFEM::EdgeElementForcesAndSourcesCore EdgeEle
Simple interface for fast problem set-up.
Definition Simple.hpp:27
bool & getAddBoundaryFE()
Get the addSkeletonFE.
Definition Simple.hpp:436
MoFEMErrorCode loadFile(const std::string options, const std::string mesh_file_name, LoadFileFunc loadFunc=defaultLoadFileFunc)
Load mesh file.
Definition Simple.cpp:194
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.
Definition Simple.cpp:264
MoFEMErrorCode getOptions()
get options
Definition Simple.cpp:180
MoFEMErrorCode setFieldOrder(const std::string field_name, const int order, const Range *ents=NULL)
Set field order.
Definition Simple.cpp:473
bool & getAddSkeletonFE()
Get the addSkeletonFE.
Definition Simple.hpp:425
MoFEMErrorCode setUp(const PetscBool is_partitioned=PETSC_TRUE)
Setup problem.
Definition Simple.cpp:611
const std::string getDomainFEName() const
Get the Domain FE Name.
Definition Simple.hpp:341
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.