v0.14.0
Loading...
Searching...
No Matches
Classes | Typedefs | Functions | Variables
boundary_marker.cpp File Reference
#include <MoFEM.hpp>

Go to the source code of this file.

Classes

struct  OpFace
 

Typedefs

using FaceEle = MoFEM::FaceElementForcesAndSourcesCore
 
using FaceEleOp = FaceEle::UserDataOperator
 

Functions

int main (int argc, char *argv[])
 

Variables

static char help [] = "...\n\n"
 

Typedef Documentation

◆ FaceEle

Examples
boundary_marker.cpp.

Definition at line 16 of file boundary_marker.cpp.

◆ FaceEleOp

Examples
boundary_marker.cpp.

Definition at line 17 of file boundary_marker.cpp.

Function Documentation

◆ main()

int main ( int argc,
char * argv[] )

Definition at line 74 of file boundary_marker.cpp.

74 {
75
76 MoFEM::Core::Initialize(&argc, &argv, (char *)0, help);
77
78 try {
79
80 moab::Core mb_instance;
81 moab::Interface &moab = mb_instance;
82
83 // Create MoFEM instance
84 MoFEM::Core core(moab);
85 MoFEM::Interface &m_field = core;
86
87 int order = 4;
88 CHKERR PetscOptionsGetInt(PETSC_NULL, "", "-order", &order, PETSC_NULL);
89
90 DMType dm_name = "DMMOFEM";
91 CHKERR DMRegister_MoFEM(dm_name);
92
93 auto *simple_interface = m_field.getInterface<Simple>();
94 CHKERR simple_interface->getOptions();
95 CHKERR simple_interface->loadFile("");
96 CHKERR simple_interface->addDomainField("FIELD1", H1,
98 CHKERR simple_interface->setFieldOrder("FIELD1", order);
99 CHKERR simple_interface->setUp();
100
101 auto get_ents_on_mesh_skin = [&]() {
102 Range faces;
103 CHKERR m_field.get_moab().get_entities_by_type(0, MBTRI, faces);
104 Skinner skin(&m_field.get_moab());
105 Range skin_edges;
106 CHKERR skin.find_skin(0, faces, false, skin_edges);
107 Range skin_verts;
108 CHKERR moab.get_connectivity(skin_edges, skin_verts, true);
109 skin_edges.merge(skin_verts);
110 return skin_edges;
111 };
112
113 auto mark_boundary_dofs = [&](Range &skin_edges) {
114 auto problem_manager = m_field.getInterface<ProblemsManager>();
115 std::vector<unsigned char> marker;
116 problem_manager->markDofs(simple_interface->getProblemName(), ROW,
117 ProblemsManager::OR, skin_edges, marker);
118 // Unset coefficient 1, e.g. u_y
119 problem_manager->modifyMarkDofs(simple_interface->getProblemName(), ROW,
120 "FIELD1", 1, 1, ProblemsManager::AND, 0,
121 marker);
122 return marker;
123 };
124
125 auto skin_ents = get_ents_on_mesh_skin();
126
127 // Get global local vector of marked DOFs. Is global, since is set for all
128 // DOFs on processor. Is local since only DOFs on processor are in the
129 // vector. To access DOFs use local indices.
130 auto marker = mark_boundary_dofs(skin_ents);
131
132 boost::shared_ptr<FaceEle> fe(new FaceEle(m_field));
133 fe->getOpPtrVector().push_back(new OpFace(skin_ents, marker));
134
135 auto dm = simple_interface->getDM();
136 CHKERR DMoFEMLoopFiniteElements(dm, simple_interface->getDomainFEName(),
137 fe);
138 }
140
142}
static char help[]
MoFEM::FaceElementForcesAndSourcesCore FaceEle
@ ROW
#define CATCH_ERRORS
Catch errors.
@ AINSWORTH_LEGENDRE_BASE
Ainsworth Cole (Legendre) approx. base .
Definition definitions.h:60
@ H1
continuous field
Definition definitions.h:85
#define CHKERR
Inline error check.
constexpr int order
PetscErrorCode DMRegister_MoFEM(const char sname[])
Register MoFEM problem.
Definition DMMoFEM.cpp:47
PetscErrorCode DMoFEMLoopFiniteElements(DM dm, const char fe_name[], MoFEM::FEMethod *method, CacheTupleWeakPtr cache_ptr=CacheTupleSharedPtr())
Executes FEMethod for finite elements in DM.
Definition DMMoFEM.cpp:572
auto marker
set bit to marker
PetscErrorCode PetscOptionsGetInt(PetscOptions *, const char pre[], const char name[], PetscInt *ivalue, PetscBool *set)
virtual moab::Interface & get_moab()=0
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
Deprecated interface functions.
Problem manager is used to build and partition problems.
MoFEMErrorCode markDofs(const std::string problem_name, RowColData rc, const enum MarkOP op, const Range ents, std::vector< unsigned char > &marker) const
Create vector with marked indices.
Simple interface for fast problem set-up.
Definition Simple.hpp:27
MoFEMErrorCode getOptions()
get options
Definition Simple.cpp:180
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.

Variable Documentation

◆ help

char help[] = "...\n\n"
static

Definition at line 14 of file boundary_marker.cpp.