MoFEM::BitRefManager Struct Reference

Managing BitRefLevels. More...

#include <src/interfaces/BitRefManager.hpp>

Inheritance diagram for MoFEM::BitRefManager:

Collaboration diagram for MoFEM::BitRefManager:

Public Member Functions
MoFEMErrorCode	query_interface (const MOFEMuuid &uuid, UnknownInterface **iface) const
	BitRefManager (const MoFEM::Core &core)
	~BitRefManager ()
	Destructor. More...
Public Member Functions inherited from MoFEM::UnknownInterface
template<class IFACE >
MoFEMErrorCode	registerInterface (const MOFEMuuid &uuid, bool error_if_registration_failed=true)
	Register interface. More...
template<class IFACE , bool VERIFY = false>
MoFEMErrorCode	getInterface (const MOFEMuuid &uuid, IFACE *&iface) const
	Get interface by uuid and return reference to pointer of interface. More...
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE *&iface) const
	Get interface refernce to pointer of interface. More...
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE **const iface) const
	Get interface pointer to pointer of interface. More...
template<class IFACE , typename boost::enable_if< boost::is_pointer< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get interface pointer to pointer of interface. More...
template<class IFACE , typename boost::enable_if< boost::is_reference< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get reference to interface. More...
template<class IFACE >
IFACE *	getInterface () const
	Function returning pointer to interface. More...
virtual	~UnknownInterface ()
virtual MoFEMErrorCode	getLibVersion (Version &version) const
	Get library version. More...
virtual const MoFEMErrorCode	getFileVersion (moab::Interface &moab, Version &version) const
	Get database major version. More...
virtual MoFEMErrorCode	getInterfaceVersion (Version &version) const
	Get database major version. More...
template<>
MoFEMErrorCode	getInterface (const MOFEMuuid &uuid, UnknownInterface *&iface) const

Public Attributes
MoFEM::Core &	cOre
bool	dEbug

Setting and shifting bits
MoFEMErrorCode	setBitRefLevel (const Range &ents, const BitRefLevel &bit, const bool only_tets=true, int verb=0) const
	add entities to database and set bit ref levelThis function set bit ref level, add entries to core database and create ref finite elements. Finite elements are create of entities in function argument, whereas all lower dimension entities are added as a field entities More...
MoFEMErrorCode	setElementsBitRefLevel (const Range &ents, const BitRefLevel bit=BitRefLevel(), int verb=0) const
	add entities to database and set bit ref level More...
MoFEMErrorCode	setEntitiesBitRefLevel (const Range &ents, const BitRefLevel bit=BitRefLevel(), int verb=0) const
	add entities to database and set bit ref level More...
MoFEMErrorCode	addToDatabaseBitRefLevelByType (const EntityType type, const BitRefLevel &bit, int verb=0) const
MoFEMErrorCode	addToDatabaseBitRefLevelByDim (const int dim, const BitRefLevel &bit, int verb=0) const
MoFEMErrorCode	setBitRefLevelByDim (const EntityHandle meshset, const int dim, const BitRefLevel &bit, int verb=0) const
MoFEMErrorCode	setBitRefLevelByType (const EntityHandle meshset, const EntityType type, const BitRefLevel &bit, int verb=0) const
MoFEMErrorCode	setBitLevelToMeshset (const EntityHandle meshset, const BitRefLevel &bit, int verb=0) const
MoFEMErrorCode	addBitRefLevel (const Range &ents, const BitRefLevel &bit, int verb=QUIET) const
	add bit ref level to ref entity More...
MoFEMErrorCode	addBitRefLevelByDim (const EntityHandle meshset, const int dim, const BitRefLevel &bit, int verb=QUIET) const
	add bit ref level by dimension More...
MoFEMErrorCode	setNthBitRefLevel (const Range &ents, const int n, const bool b, int verb=0) const
	Set nth bit ref level. More...
MoFEMErrorCode	setNthBitRefLevel (const int n, const bool b, int verb=0) const
	Set nth bit ref level. More...
MoFEMErrorCode	shiftLeftBitRef (const int shift, const BitRefLevel mask=BitRefLevel().set(), int verb=-1) const
	left shift bit ref levelthis results of deletion of entities on far left side More...
MoFEMErrorCode	shiftRightBitRef (const int shift, const BitRefLevel mask=BitRefLevel().set(), int verb=-1) const
	right shift bit ref level More...

Entity handlers by bit ref level
MoFEMErrorCode	filterEntitiesByRefLevel (const BitRefLevel &bit, const BitRefLevel &mask, Range &ents, int verb=QUIET) const
	filter entities by bit ref level More...
MoFEMErrorCode	getEntitiesByTypeAndRefLevel (const BitRefLevel &bit, const BitRefLevel &mask, const EntityType type, const EntityHandle meshset, int verb=0) const
	add all ents from ref level given by bit to meshset More...
MoFEMErrorCode	getEntitiesByTypeAndRefLevel (const BitRefLevel &bit, const BitRefLevel &mask, const EntityType type, Range &ents, int verb=0) const
	add all ents from ref level given by bit to meshset More...
MoFEMErrorCode	getEntitiesByDimAndRefLevel (const BitRefLevel &bit, const BitRefLevel &mask, const int dim, const EntityHandle meshset, int verb=0) const
	add all ents from ref level given by bit to meshset More...
MoFEMErrorCode	getEntitiesByDimAndRefLevel (const BitRefLevel &bit, const BitRefLevel &mask, const int dim, Range &ents, int verb=0) const
	add all ents from ref level given by bit to meshset More...
MoFEMErrorCode	getEntitiesByRefLevel (const BitRefLevel &bit, const BitRefLevel &mask, const EntityHandle meshset, const int verb=QUIET) const
	add all ents from ref level given by bit to meshset More...
MoFEMErrorCode	getEntitiesByRefLevel (const BitRefLevel &bit, const BitRefLevel &mask, Range &ents, const int verb=QUIET) const
	add all ents from ref level given by bit to meshset More...
MoFEMErrorCode	getEntitiesByParentType (const BitRefLevel &bit, const BitRefLevel &mask, const EntityType type, Range &ents) const
	get entities by bit ref level and type of parent More...
MoFEMErrorCode	getAllEntitiesNotInDatabase (Range &ents) const
	Get all entities not in database. More...
MoFEMErrorCode	filterEntitiesNotInDatabase (Range &ents) const
	Get entities not in database. More...

Get adjacencies bit ref level
virtual MoFEMErrorCode	getAdjacenciesEquality (const EntityHandle from_entity, const int to_dimension, Range &adj_entities) const
	Get the adjacencies associated with a entity to entities of a specified dimension.bit ref level of adjacent entities is equal to bit ref level of adjacent entities. More...
virtual MoFEMErrorCode	getAdjacenciesAny (const EntityHandle from_entity, const int to_dimension, Range &adj_entities) const
	Get the adjacencies associated with a entity to entities of a specified dimension.bit ref level of adjacent entities is any of bit ref level of adjacent entities. More...
virtual MoFEMErrorCode	getAdjacencies (const Problem *problem_ptr, const EntityHandle *from_entities, const int num_entities, const int to_dimension, Range &adj_entities, const int operation_type=moab::Interface::INTERSECT, const int verb=0) const
	Get the adjacencies associated with a entity to entities of a specified dimension.bit ref level of adjacent entities is equal to bit ref level of adjacent entities. More...
virtual MoFEMErrorCode	getAdjacencies (const BitRefLevel &bit, const EntityHandle *from_entities, const int num_entities, const int to_dimension, Range &adj_entities, const int operation_type=moab::Interface::INTERSECT, const int verb=0) const
	Get the adjacencies associated with a entity to entities of a specified dimension.bit ref level of adjacent entities is equal to bit ref level of adjacent entities. More...

Update meshsets and ranges by children
MoFEMErrorCode	updateMeshsetByEntitiesChildren (const EntityHandle parent, const BitRefLevel &parent_bit, const BitRefLevel &parent_mask, const BitRefLevel &child_bit, const BitRefLevel &child_mask, const EntityHandle child, EntityType child_type, const bool recursive=false, int verb=0)
	Get child entities form meshset containing parent entitiesSearch for refined entities of given type whose parent are entities in the parent meshset. It can be used for example to transfer information about boundary conditions to refined mesh or split mesh by interface elements. It is used by function refine_MESHSET, to update MESHSET finite elements. More...
MoFEMErrorCode	updateMeshsetByEntitiesChildren (const EntityHandle parent, const BitRefLevel &child_bit, const EntityHandle child, EntityType child_type, const bool recursive=false, int verb=0)
	Get child entities form meshset containing parent entities. More...
MoFEMErrorCode	updateFieldMeshsetByEntitiesChildren (const BitRefLevel &child_bit, int verb=0)
	update fields meshesets by child entities More...
MoFEMErrorCode	updateFieldMeshsetByEntitiesChildren (const std::string name, const BitRefLevel &child_bit, int verb=0)
	update field meshset by child entities More...
MoFEMErrorCode	updateFiniteElementMeshsetByEntitiesChildren (const std::string name, const BitRefLevel &child_bit, const EntityType fe_ent_type, int verb=0)
	update finite element meshset by child entities More...
MoFEMErrorCode	updateRange (const Range &parent, Range &child)
	Update range by prents. More...

Writting files
MoFEMErrorCode	writeBitLevelByType (const BitRefLevel &bit, const BitRefLevel &mask, const EntityType type, const char *file_name, const char *file_type, const char *options) const
	Write bit ref level to file. More...
MoFEMErrorCode	writeEntitiesNotInDatabase (const char *file_name, const char *file_type, const char *options) const
	Write ents not in database. More...

Additional Inherited Members
Protected Member Functions inherited from MoFEM::UnknownInterface
boost::typeindex::type_index	getClassIdx (const MOFEMuuid &uid) const
	Get type name for interface Id. More...
MOFEMuuid	getUId (const boost::typeindex::type_index &class_idx) const
	Get interface Id for class name. More...

Detailed Description

Managing BitRefLevels.

Examples:

elasticity.cpp, field_to_vertices.cpp, forces_and_sources_testing_edge_element.cpp, mesh_cut.cpp, mesh_smoothing.cpp, minimal_surface_area.cpp, split_sideset.cpp, and unsaturated_transport.cpp.

Definition at line 34 of file BitRefManager.hpp.

Constructor & Destructor Documentation

BitRefManager()

MoFEM::BitRefManager::BitRefManager

(

const MoFEM::Core &

core

)

Definition at line 34 of file BitRefManager.cpp.

    : cOre(const_cast<MoFEM::Core &>(core)), dEbug(false) {}

~BitRefManager()

MoFEM::BitRefManager::~BitRefManager

(

)

Destructor.

Definition at line 36 of file BitRefManager.cpp.

{}

Member Function Documentation

filterEntitiesByRefLevel()

MoFEMErrorCode MoFEM::BitRefManager::filterEntitiesByRefLevel	(	const BitRefLevel &	bit,
		const BitRefLevel &	mask,
		Range &	ents,
		int	verb = QUIET
	)		const

filter entities by bit ref level

Parameters

bit
mask
ents
QUIET

Returns

MoFEMErrorCode

Definition at line 556 of file BitRefManager.cpp.

                                                                       {
  MoFEM::Interface &m_field = cOre;
  moab::Interface &moab(m_field.get_moab());
  MoFEMFunctionBeginHot;

  const BitRefLevel *tag_bit;

  Range swap_ents;
  for (Range::pair_iterator p_eit = ents.pair_begin(); p_eit != ents.pair_end();
       ++p_eit) {
    EntityHandle f = p_eit->first;
    EntityHandle s = p_eit->second;

    EntityHandle ff = 0;
    EntityHandle ss = 0;

    for (; f != (s + 1); ++f) {
      // Get entity bit ref level
      rval = moab.tag_get_by_ptr(cOre.get_th_RefBitLevel(), &f, 1,
                                 (const void **)(&tag_bit));
      if (PetscLikely(rval == MB_SUCCESS)) {
        if ((mask.any() && tag_bit->none()) ||
            ((*tag_bit) & mask) != (*tag_bit) || (((*tag_bit) & bit).none())) {
          // Entity not on BitRefLevel.
          if (ff != 0) {
            // Add sub-range
            swap_ents.insert(ff, ss);
            ff = ss = 0;
          }
        } else {
          if (ff == 0) {
            // Start new sub-range
            ff = ss = f;
          } else {
            // Add entity to sub-range
            ss = f;
          }
        }
      } else {
        if (ff != 0) {
          // Add sub-range
          swap_ents.insert(ff, ss);
        }
        // Start new sub-range
        ff = ss = 0;
      }
    }

    if (ff != 0) {
      swap_ents.insert(ff, ss);
    }
  }

  ents.swap(swap_ents);

  MoFEMFunctionReturnHot(0);
}

filterEntitiesNotInDatabase()

MoFEMErrorCode MoFEM::BitRefManager::filterEntitiesNotInDatabase

(

Range &

ents

)

const

Get entities not in database.

Parameters

ents

Returns

MoFEMErrorCode

Definition at line 722 of file BitRefManager.cpp.

                                                                           {
  MoFEM::Interface &m_field = cOre;
  const RefEntity_multiIndex *ref_ents_ptr;
  MoFEMFunctionBeginHot;
  ierr = m_field.get_ref_ents(&ref_ents_ptr);
  CHKERRG(ierr);
  Range::iterator eit = ents.begin();
  for (; eit != ents.end();) {
    RefEntity_multiIndex::index<Ent_mi_tag>::type::iterator rit;
    rit = ref_ents_ptr->get<Ent_mi_tag>().find(*eit);
    if (rit != ref_ents_ptr->get<Ent_mi_tag>().end()) {
      eit = ents.erase(eit);
    } else {
      eit++;
    }
  }
  MoFEMFunctionReturnHot(0);
}

getAllEntitiesNotInDatabase()

MoFEMErrorCode MoFEM::BitRefManager::getAllEntitiesNotInDatabase

(

Range &

ents

)

const

Get all entities not in database.

Parameters

ents

Returns

MoFEMErrorCode

Definition at line 710 of file BitRefManager.cpp.

                                                                           {
  MoFEM::Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBeginHot;
  rval = moab.get_entities_by_handle(0, ents, false);
  CHKERRQ_MOAB(rval);
  ents = subtract(ents, ents.subset_by_type(MBENTITYSET));
  ierr = filterEntitiesNotInDatabase(ents);
  CHKERRG(ierr);
  MoFEMFunctionReturnHot(0);
}

getEntitiesByParentType()

MoFEMErrorCode MoFEM::BitRefManager::getEntitiesByParentType	(	const BitRefLevel &	bit,
		const BitRefLevel &	mask,
		const EntityType	type,
		Range &	ents
	)		const

get entities by bit ref level and type of parent

Note

Entities have to be added to MoFEM database

Parameters

BitRefLevel	bitLevel
BitRefLevel	mask
type	of parent
ents	returned ents

Returns

error code

Definition at line 684 of file BitRefManager.cpp.

                                                                         {
  MoFEM::Interface &m_field = cOre;
  // moab::Interface &moab = m_field.get_moab();
  const RefEntity_multiIndex *ref_ents_ptr;
  MoFEMFunctionBeginHot;
  ierr = m_field.get_ref_ents(&ref_ents_ptr);
  CHKERRG(ierr);
  typedef RefEntity_multiIndex::index<ParentEntType_mi_tag>::type
      RefEntsByParentType;
  const RefEntsByParentType &ref_ents =
      ref_ents_ptr->get<ParentEntType_mi_tag>();
  RefEntsByParentType::iterator it, hi_it;
  it = ref_ents.lower_bound(type);
  hi_it = ref_ents.upper_bound(type);
  for (; it != hi_it; it++) {
    const BitRefLevel &ent_bit = it->get()->getBitRefLevel();
    if ((ent_bit & mask) == ent_bit && (ent_bit & bit).any()) {
      ents.insert(it->get()->getRefEnt());
    }
  }
  MoFEMFunctionReturnHot(0);
}

query_interface()

MoFEMErrorCode MoFEM::BitRefManager::query_interface ( const MOFEMuuid &  uuid, UnknownInterface **  iface  ) const

virtual

Implements MoFEM::UnknownInterface.

Definition at line 22 of file BitRefManager.cpp.

                                                                              {
  MoFEMFunctionBeginHot;
  *iface = NULL;
  if (uuid == IDD_MOFEMBitRefManager) {
    *iface = const_cast<BitRefManager *>(this);
    MoFEMFunctionReturnHot(0);
  }
  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "unknown interface");
  MoFEMFunctionReturnHot(0);
}

writeBitLevelByType()

MoFEMErrorCode MoFEM::BitRefManager::writeBitLevelByType	(	const BitRefLevel &	bit,
		const BitRefLevel &	mask,
		const EntityType	type,
		const char *	file_name,
		const char *	file_type,
		const char *	options
	)		const

Write bit ref level to file.

Parameters

bit	bit ref level
mask	mask of bit ref level
type	type of entity
file_name	file name (see moab documentation)
file_type	file type (see moab documentation)
options	file options (see moab documentation)

Returns

error code

Definition at line 513 of file BitRefManager.cpp.

                                                                             {
  MoFEM::Interface &m_field = cOre;
  moab::Interface &moab(m_field.get_moab());
  MoFEMFunctionBegin;
  EntityHandle meshset;
  CHKERR moab.create_meshset(MESHSET_SET, meshset);
  CHKERR getEntitiesByTypeAndRefLevel(bit, mask, type, meshset);
  CHKERR moab.write_file(file_name, file_type, options, &meshset, 1);
  CHKERR moab.delete_entities(&meshset, 1);
  MoFEMFunctionReturn(0);
}

writeEntitiesNotInDatabase()

MoFEMErrorCode MoFEM::BitRefManager::writeEntitiesNotInDatabase	(	const char *	file_name,
		const char *	file_type,
		const char *	options
	)		const

Write ents not in database.

Parameters

file_name
file_type	for example "VTK"
options

Returns

MoFEMErrorCode

Definition at line 527 of file BitRefManager.cpp.

                                                                             {
  MoFEM::Interface &m_field = cOre;
  moab::Interface &moab(m_field.get_moab());
  MoFEMFunctionBegin;
  EntityHandle meshset;
  CHKERR moab.create_meshset(MESHSET_SET, meshset);
  Range ents;
  ierr = getAllEntitiesNotInDatabase(ents);
  CHKERR moab.add_entities(meshset, ents);
  CHKERR moab.write_file(file_name, file_type, options, &meshset, 1);
  CHKERR moab.delete_entities(&meshset, 1);
  MoFEMFunctionReturn(0);
}

Member Data Documentation

cOre

MoFEM::Core& MoFEM::BitRefManager::cOre

Definition at line 39 of file BitRefManager.hpp.

dEbug

bool MoFEM::BitRefManager::dEbug

Definition at line 40 of file BitRefManager.hpp.

---

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

• src/interfaces/BitRefManager.hpp
• src/interfaces/impl/BitRefManager.cpp