v0.13.2
Loading...
Searching...
No Matches
Public Member Functions | List of all members
MoFEM::PostProcGenerateRefMesh< MBTET > Struct Reference

#include <src/post_proc/PostProcBrokenMeshInMoabBase.hpp>

Inheritance diagram for MoFEM::PostProcGenerateRefMesh< MBTET >:
[legend]
Collaboration diagram for MoFEM::PostProcGenerateRefMesh< MBTET >:
[legend]

Public Member Functions

MoFEMErrorCode generateReferenceElementMesh ()
 
 PostProcGenerateRefMeshBase ()
 
- Public Member Functions inherited from MoFEM::PostProcGenerateRefMeshBase
 PostProcGenerateRefMeshBase ()
 
virtual ~PostProcGenerateRefMeshBase ()=default
 
MoFEMErrorCode getOptions ()
 
virtual MoFEMErrorCode generateReferenceElementMesh ()=0
 

Additional Inherited Members

- Public Attributes inherited from MoFEM::PostProcGenerateRefMeshBase
std::vector< MatrixDoublelevelShapeFunctions
 
std::vector< MatrixDoublelevelGaussPtsOnRefMesh
 
std::vector< ublas::matrix< int > > levelRef
 
EntityHandle startingVertEleHandle
 
std::vector< double * > verticesOnEleArrays
 
EntityHandle startingEleHandle
 
EntityHandleeleConn
 
int countEle
 
int countVertEle
 
int nbVertices
 
int nbEles
 
PetscBool hoNodes
 
int defMaxLevel
 
std::string optPrefix
 

Detailed Description

Definition at line 66 of file PostProcBrokenMeshInMoabBase.hpp.

Member Function Documentation

◆ generateReferenceElementMesh()

MoFEMErrorCode MoFEM::PostProcGenerateRefMesh< MBTET >::generateReferenceElementMesh ( )
virtual

Implements MoFEM::PostProcGenerateRefMeshBase.

Definition at line 33 of file PostProcBrokenMeshInMoabBase.cpp.

33 {
35
37
38 const int max_level = defMaxLevel;
39
40 moab::Core core_ref;
41 moab::Interface &moab_ref = core_ref;
42
43 auto create_reference_element = [&moab_ref]() {
45 constexpr double base_coords[] = {0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1};
46 EntityHandle nodes[4];
47 for (int nn = 0; nn < 4; nn++) {
48 CHKERR
49 moab_ref.create_vertex(&base_coords[3 * nn], nodes[nn]);
50 }
51 EntityHandle tet;
52 CHKERR moab_ref.create_element(MBTET, nodes, 4, tet);
54 };
55
56 MoFEM::CoreTmp<-1> m_core_ref(moab_ref, PETSC_COMM_SELF, -2);
57 MoFEM::Interface &m_field_ref = m_core_ref;
58
59 auto refine_ref_tetrahedron = [this, &m_field_ref, max_level]() {
61 // seed ref mofem database by setting bit ref level to reference
62 // tetrahedron
63 CHKERR
64 m_field_ref.getInterface<BitRefManager>()->setBitRefLevelByDim(
65 0, 3, BitRefLevel().set(0));
66 for (int ll = 0; ll != max_level; ++ll) {
67 MOFEM_TAG_AND_LOG_C("WORLD", Sev::noisy, "PostProc", "Refine Level %d",
68 ll);
69 Range edges;
70 CHKERR m_field_ref.getInterface<BitRefManager>()
71 ->getEntitiesByTypeAndRefLevel(BitRefLevel().set(ll),
72 BitRefLevel().set(), MBEDGE, edges);
73 Range tets;
74 CHKERR m_field_ref.getInterface<BitRefManager>()
75 ->getEntitiesByTypeAndRefLevel(BitRefLevel().set(ll),
76 BitRefLevel(ll).set(), MBTET, tets);
77 // refine mesh
78 MeshRefinement *m_ref;
79 CHKERR m_field_ref.getInterface(m_ref);
80 CHKERR m_ref->addVerticesInTheMiddleOfEdges(edges,
81 BitRefLevel().set(ll + 1));
82 CHKERR m_ref->refineTets(tets, BitRefLevel().set(ll + 1));
83 }
85 };
86
87 auto get_ref_gauss_pts_and_shape_functions = [this, max_level, &moab_ref,
88 &m_field_ref]() {
90 for (int ll = 0; ll != max_level + 1; ++ll) {
91 Range tets;
92 CHKERR
93 m_field_ref.getInterface<BitRefManager>()->getEntitiesByTypeAndRefLevel(
94 BitRefLevel().set(ll), BitRefLevel().set(ll), MBTET, tets);
95 if (hoNodes) {
96 EntityHandle meshset;
97 CHKERR moab_ref.create_meshset(MESHSET_SET, meshset);
98 CHKERR moab_ref.add_entities(meshset, tets);
99 CHKERR moab_ref.convert_entities(meshset, true, false, false);
100 CHKERR moab_ref.delete_entities(&meshset, 1);
101 }
102 Range elem_nodes;
103 CHKERR moab_ref.get_connectivity(tets, elem_nodes, false);
104
105 auto &gauss_pts = levelGaussPtsOnRefMesh[ll];
106 gauss_pts.resize(elem_nodes.size(), 4, false);
107 std::map<EntityHandle, int> little_map;
108 Range::iterator nit = elem_nodes.begin();
109 for (int gg = 0; nit != elem_nodes.end(); nit++, gg++) {
110 CHKERR moab_ref.get_coords(&*nit, 1, &gauss_pts(gg, 0));
111 little_map[*nit] = gg;
112 }
113 gauss_pts = trans(gauss_pts);
114
115 auto &ref_tets = levelRef[ll];
116 Range::iterator tit = tets.begin();
117 for (int tt = 0; tit != tets.end(); ++tit, ++tt) {
118 const EntityHandle *conn;
119 int num_nodes;
120 CHKERR moab_ref.get_connectivity(*tit, conn, num_nodes, false);
121 if (tt == 0) {
122 ref_tets.resize(tets.size(), num_nodes);
123 }
124 for (int nn = 0; nn != num_nodes; ++nn) {
125 ref_tets(tt, nn) = little_map[conn[nn]];
126 }
127 }
128
129 auto &shape_functions = levelShapeFunctions[ll];
130 shape_functions.resize(elem_nodes.size(), 4);
131 CHKERR ShapeMBTET(&*shape_functions.data().begin(), &gauss_pts(0, 0),
132 &gauss_pts(1, 0), &gauss_pts(2, 0), elem_nodes.size());
133 }
135 };
136
137 levelRef.resize(max_level + 1);
138 levelGaussPtsOnRefMesh.resize(max_level + 1);
139 levelShapeFunctions.resize(max_level + 1);
140
141 CHKERR create_reference_element();
142 CHKERR refine_ref_tetrahedron();
143 CHKERR get_ref_gauss_pts_and_shape_functions();
144
146}
#define MOFEM_TAG_AND_LOG_C(channel, severity, tag, format,...)
Tag and log in channel.
Definition: LogManager.hpp:367
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:346
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:416
#define CHKERR
Inline error check.
Definition: definitions.h:535
PetscErrorCode ShapeMBTET(double *N, const double *G_X, const double *G_Y, const double *G_Z, int DIM)
calculate shape functions
Definition: fem_tools.c:306
std::bitset< BITREFLEVEL_SIZE > BitRefLevel
Bit structure attached to each entity identifying to what mesh entity is attached.
Definition: Types.hpp:40
Deprecated interface functions.
std::vector< ublas::matrix< int > > levelRef
PetscBool hoNodes
int defMaxLevel
std::vector< MatrixDouble > levelShapeFunctions
std::vector< MatrixDouble > levelGaussPtsOnRefMesh
MoFEMErrorCode getOptions()
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface refernce to pointer of interface.

◆ PostProcGenerateRefMeshBase()

MoFEM::PostProcGenerateRefMeshBase::PostProcGenerateRefMeshBase ( )

The documentation for this struct was generated from the following files: