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.
294 auto get_nb_of_ref_levels_from_options = [
this] {
297 PetscBool flg = PETSC_TRUE;
299 "-my_max_post_proc_ref_level", &max_level, &flg);
306 const int max_level = get_nb_of_ref_levels_from_options();
309 moab::Interface &moab_ref = core_ref;
311 auto create_reference_element = [&moab_ref]() {
313 const double base_coords[] = {0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1};
315 for (
int nn = 0; nn < 4; nn++) {
316 CHKERR moab_ref.create_vertex(&base_coords[3 * nn], nodes[nn]);
319 CHKERR moab_ref.create_element(MBTET, nodes, 4, tet);
323 MoFEM::Core m_core_ref(moab_ref, PETSC_COMM_SELF, -2);
326 auto refine_ref_tetrahedron = [
this, &m_field_ref, max_level]() {
332 for (
int ll = 0; ll != max_level; ++ll) {
333 PetscPrintf(T::mField.get_comm(),
"Refine Level %d\n", ll);
336 ->getEntitiesByTypeAndRefLevel(
BitRefLevel().set(ll),
340 ->getEntitiesByTypeAndRefLevel(
BitRefLevel().set(ll),
343 MeshRefinement *m_ref;
345 CHKERR m_ref->add_vertices_in_the_middle_of_edges(
352 auto get_ref_gauss_pts_and_shape_functions = [
this, max_level, &moab_ref,
355 for (
int ll = 0; ll != max_level + 1; ++ll) {
358 ->getEntitiesByTypeAndRefLevel(
BitRefLevel().set(ll),
362 CHKERR moab_ref.create_meshset(MESHSET_SET, meshset);
363 CHKERR moab_ref.add_entities(meshset, tets);
364 CHKERR moab_ref.convert_entities(meshset,
true,
false,
false);
365 CHKERR moab_ref.delete_entities(&meshset, 1);
368 CHKERR moab_ref.get_connectivity(tets, elem_nodes,
false);
371 gauss_pts.resize(elem_nodes.size(), 4,
false);
372 std::map<EntityHandle, int> little_map;
373 Range::iterator nit = elem_nodes.begin();
374 for (
int gg = 0; nit != elem_nodes.end(); nit++, gg++) {
375 CHKERR moab_ref.get_coords(&*nit, 1, &gauss_pts(gg, 0));
376 little_map[*nit] = gg;
378 gauss_pts = trans(gauss_pts);
381 Range::iterator tit = tets.begin();
382 for (
int tt = 0; tit != tets.end(); ++tit, ++tt) {
385 CHKERR moab_ref.get_connectivity(*tit, conn, num_nodes,
false);
389 ref_tets.resize(tets.size(), num_nodes);
391 for (
int nn = 0; nn != num_nodes; ++nn) {
392 ref_tets(tt, nn) = little_map[conn[nn]];
397 shape_functions.resize(elem_nodes.size(), 4);
399 &gauss_pts(1, 0), &gauss_pts(2, 0),
409 CHKERR create_reference_element();
410 CHKERR refine_ref_tetrahedron();
411 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 >