Create adjacency to parent skeleton elements. More...

#include <src/finite_elements/MeshProjectionDataOperators.hpp>

Inheritance diagram for MoFEM::ParentFiniteElementAdjacencyFunctionSkeleton< DIM >:

Collaboration diagram for MoFEM::ParentFiniteElementAdjacencyFunctionSkeleton< DIM >:

Public Member Functions
	ParentFiniteElementAdjacencyFunctionSkeleton (BitRefLevel bit_parent, BitRefLevel bit_parent_mask, BitRefLevel bit_ent, BitRefLevel bit_ent_mask)

MoFEMErrorCode	operator() (moab::Interface &moab, const Field &field, const EntFiniteElement &fe, std::vector< EntityHandle > &adjacency)
	Function setting adjacencies to DOFs of parent element. More...

Public Member Functions inherited from MoFEM::ParentFiniteElementAdjacencyFunction< DIM >
	ParentFiniteElementAdjacencyFunction (BitRefLevel bit_parent, BitRefLevel bit_parent_mask, BitRefLevel bit_ent, BitRefLevel bit_ent_mask)

MoFEMErrorCode	operator() (moab::Interface &moab, const Field &field, const EntFiniteElement &fe, std::vector< EntityHandle > &adjacency)
	Function setting adjacencies to DOFs of parent element. More...

Protected Attributes
std::vector< EntityHandle >	brideAdjacencyEdge

Protected Attributes inherited from MoFEM::ParentFiniteElementAdjacencyFunction< DIM >
BitRefLevel	bitParent

BitRefLevel	bitParentMask

BitRefLevel	bitEnt

BitRefLevel	bitEntMask

std::vector< EntityHandle >	adjTmp

Additional Inherited Members
Protected Member Functions inherited from MoFEM::ParentFiniteElementAdjacencyFunction< DIM >
MoFEMErrorCode	getParent (EntityHandle fe, std::vector< EntityHandle > &parents, moab::Interface &moab, Tag th_parent_handle, Tag th_bit_level)

MoFEMErrorCode	getParentsAdjacencies (const Field &field, moab::Interface &moab, std::vector< EntityHandle > &parents, std::vector< EntityHandle > &adjacency)

MoFEMErrorCode	getDefaultAdjacencies (moab::Interface &moab, const Field &field, const EntFiniteElement &fe, std::vector< EntityHandle > &adjacency)

Detailed Description

template<int DIM>
struct MoFEM::ParentFiniteElementAdjacencyFunctionSkeleton< DIM >

Create adjacency to parent skeleton elements.

That class is used during entity finite element construction.

Template Parameters

DIM	dimension of parent element

Definition at line 313 of file MeshProjectionDataOperators.hpp.

Constructor & Destructor Documentation

◆ ParentFiniteElementAdjacencyFunctionSkeleton()

template<int DIM>

MoFEM::ParentFiniteElementAdjacencyFunctionSkeleton< DIM >::ParentFiniteElementAdjacencyFunctionSkeleton	(	BitRefLevel	bit_parent,
		BitRefLevel	bit_parent_mask,
		BitRefLevel	bit_ent,
		BitRefLevel	bit_ent_mask
	)

inline

Definition at line 316 of file MeshProjectionDataOperators.hpp.

320 : ParentFiniteElementAdjacencyFunction<DIM>(bit_parent, bit_parent_mask,

321 bit_ent, bit_ent_mask) {}

Member Function Documentation

◆ operator()()

template<int DIM>

MoFEMErrorCode MoFEM::ParentFiniteElementAdjacencyFunctionSkeleton< DIM >::operator()	(	moab::Interface &	moab,
		const Field &	field,
		const EntFiniteElement &	fe,
		std::vector< EntityHandle > &	adjacency
	)

inline

Function setting adjacencies to DOFs of parent element.

Note: elements form child, see dofs from parent, so DOFs located on adjacencies of parent entity has adjacent to dofs of child.

Template Parameters

DIM	dimension of the element entity

Parameters

moab
field
fe
adjacency

Returns: MoFEMErrorCode

Definition at line 336 of file MeshProjectionDataOperators.hpp.

                                                                 {
     MoFEMFunctionBegin;
  
     adjacency.clear();
     CHKERR this->getDefaultAdjacencies(moab, field, fe, adjacency);
  
     if (field.getSpace() != NOFIELD) {
  
       const auto fe_ent = fe.getEnt();
       brideAdjacencyEdge.clear();
       CHKERR moab.get_adjacencies(&fe_ent, 1, DIM + 1, false,
                                   brideAdjacencyEdge);
  
       std::vector<EntityHandle> parents;
  
       if (this->bitParent.any()) {
         auto basic_entity_data_ptr = fe.getBasicDataPtr();
         auto th_parent_handle = basic_entity_data_ptr->th_RefParentHandle;
         auto th_bit_level = basic_entity_data_ptr->th_RefBitLevel;
         parents.reserve(BITREFLEVEL_SIZE);
         for (auto bridge_fe : brideAdjacencyEdge) {
           CHKERR this->getParent(bridge_fe, parents, moab, th_parent_handle,
                                  th_bit_level);
         };
         parents.insert(parents.end(), brideAdjacencyEdge.begin(),
                        brideAdjacencyEdge.end());
       } else {
         parents.swap(brideAdjacencyEdge);
       }
  
       CHKERR this->getParentsAdjacencies(field, moab, parents, adjacency);
  
       std::sort(adjacency.begin(), adjacency.end());
       auto it = std::unique(adjacency.begin(), adjacency.end());
       adjacency.resize(std::distance(adjacency.begin(), it));
  
       for (auto e : adjacency) {
         auto &side_table = fe.getSideNumberTable();
         if (side_table.find(e) == side_table.end())
           const_cast<SideNumber_multiIndex &>(side_table)
               .insert(
                   boost::shared_ptr<SideNumber>(new SideNumber(e, -1, 0, 0)));
       }
     }
  
     MoFEMFunctionReturn(0);
   }