Generate reference mesh on single element.
Each element is subdivided on smaller elements, i.e. a reference mesh on single element is created. Nodes of such reference mesh are used as integration points at which field values are calculated and to each node a "moab" tag is attached to store those values.
295 auto get_nb_of_ref_levels_from_options = [
this] {
298 PetscBool flg = PETSC_TRUE;
300 "-my_max_post_proc_ref_level", &max_level, &flg);
307 const int max_level = get_nb_of_ref_levels_from_options();
310 moab::Interface &moab_ref = core_ref;
312 auto create_reference_element = [&moab_ref]() {
314 const double base_coords[] = {0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1};
316 for (
int nn = 0; nn < 4; nn++) {
317 CHKERR moab_ref.create_vertex(&base_coords[3 * nn], nodes[nn]);
320 CHKERR moab_ref.create_element(MBTET, nodes, 4, tet);
324 MoFEM::Core m_core_ref(moab_ref, PETSC_COMM_SELF, -2);
327 auto refine_ref_tetrahedron = [
this, &m_field_ref, max_level]() {
333 for (
int ll = 0; ll != max_level; ++ll) {
334 PetscPrintf(T::mField.get_comm(),
"Refine Level %d\n", ll);
337 ->getEntitiesByTypeAndRefLevel(
BitRefLevel().set(ll),
341 ->getEntitiesByTypeAndRefLevel(
BitRefLevel().set(ll),
344 MeshRefinement *m_ref;
346 CHKERR m_ref->add_vertices_in_the_middel_of_edges(
353 auto get_ref_gauss_pts_and_shape_functions = [
this, max_level, &moab_ref,
356 for (
int ll = 0; ll != max_level + 1; ++ll) {
359 ->getEntitiesByTypeAndRefLevel(
BitRefLevel().set(ll),
363 CHKERR moab_ref.create_meshset(MESHSET_SET, meshset);
364 CHKERR moab_ref.add_entities(meshset, tets);
365 CHKERR moab_ref.convert_entities(meshset,
true,
false,
false);
366 CHKERR moab_ref.delete_entities(&meshset, 1);
369 CHKERR moab_ref.get_connectivity(tets, elem_nodes,
false);
372 gauss_pts.resize(elem_nodes.size(), 4,
false);
373 std::map<EntityHandle, int> little_map;
374 Range::iterator nit = elem_nodes.begin();
375 for (
int gg = 0; nit != elem_nodes.end(); nit++, gg++) {
376 CHKERR moab_ref.get_coords(&*nit, 1, &gauss_pts(gg, 0));
377 little_map[*nit] = gg;
379 gauss_pts = trans(gauss_pts);
382 Range::iterator tit = tets.begin();
383 for (
int tt = 0; tit != tets.end(); ++tit, ++tt) {
386 CHKERR moab_ref.get_connectivity(*tit, conn, num_nodes,
false);
390 ref_tets.resize(tets.size(), num_nodes);
392 for (
int nn = 0; nn != num_nodes; ++nn) {
393 ref_tets(tt, nn) = little_map[conn[nn]];
398 shape_functions.resize(elem_nodes.size(), 4);
400 &gauss_pts(1, 0), &gauss_pts(2, 0),
410 CHKERR create_reference_element();
411 CHKERR refine_ref_tetrahedron();
412 CHKERR get_ref_gauss_pts_and_shape_functions();
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
std::vector< ublas::matrix< int > > levelRefTets
bool tenNodesPostProcTets
MoFEMErrorCode getInterface(const MOFEMuuid &uuid, IFACE *&iface) const
Get interface by uuid and return reference to pointer of interface.
std::vector< MatrixDouble > levelShapeFunctions
std::vector< MatrixDouble > levelGaussPtsOnRefMesh
#define CHKERR
Inline error check.
PetscErrorCode PetscOptionsGetInt(PetscOptions *, const char pre[], const char name[], PetscInt *ivalue, PetscBool *set)
std::bitset< BITREFLEVEL_SIZE > BitRefLevel
Bit structure attached to each entity identifying to what mesh entity is attached.
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Public Member Functions inherited from PostProcTemplateOnRefineMesh< VOLUME_ELEMENT >