Managing BitRefLevels. More...

#include <src/interfaces/BitRefManager.hpp>

Inheritance diagram for MoFEM::BitRefManager:

Collaboration diagram for MoFEM::BitRefManager:

Public Member Functions
MoFEMErrorCode	query_interface (boost::typeindex::type_index type_index, UnknownInterface **iface) const

	BitRefManager (const MoFEM::Core &core)

virtual	~BitRefManager ()=default
	Destructor. More...

Public Member Functions inherited from MoFEM::UnknownInterface
template<class IFACE >
MoFEMErrorCode	registerInterface (bool error_if_registration_failed=true)
	Register 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 ()=default

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, Range *adj_ents_ptr=nullptr) const
	add entities to database and set bit ref level More...

MoFEMErrorCode	setElementsBitRefLevel (const Range &ents, const BitRefLevel bit=BitRefLevel(), int verb=QUIET) const
	add entities to database and set bit ref level More...

MoFEMErrorCode	setEntitiesBitRefLevel (const Range &ents, const BitRefLevel bit=BitRefLevel(), int verb=QUIET) const
	add entities to database and set bit ref level More...

MoFEMErrorCode	setFieldEntitiesBitRefLevel (const std::string field_name, const BitRefLevel bit=BitRefLevel(), int verb=QUIET) const
	Set the bit ref level to entities in the field meshset. More...

MoFEMErrorCode	setBitRefLevelByDim (const EntityHandle meshset, const int dim, const BitRefLevel bit, int verb=QUIET) const
	Set the Bit Ref Level By Dim object. More...

MoFEMErrorCode	setBitRefLevelByType (const EntityHandle meshset, const EntityType type, const BitRefLevel bit, int verb=QUIET) const
	Set the Bit Ref Level By Type object. More...

MoFEMErrorCode	setBitLevelToMeshset (const EntityHandle meshset, const BitRefLevel bit, int verb=0) const

MoFEMErrorCode	addToDatabaseBitRefLevelByType (const EntityType type, const BitRefLevel bit, const BitRefLevel mask=BitRefLevel().set(), int verb=QUIET) const
	Add entities which exist in MoAB database, and have set appropiate BitRef level tag, to multi-indices in MoFEM. More...

MoFEMErrorCode	addToDatabaseBitRefLevelByDim (const int dim, const BitRefLevel bit, const BitRefLevel mask=BitRefLevel().set(), int verb=QUIET) const
	Add entities which exist in MoAB database, and have set appropiate BitRef level tag, to multi-indices in MoFEM. More...

MoFEMErrorCode	lambdaBitRefLevel (boost::function< void(EntityHandle ent, BitRefLevel &bit)> fun) const
	Process bit ref level by lambda function. More...

MoFEMErrorCode	lambdaBitRefLevel (const Range &ents, boost::function< void(EntityHandle ent, BitRefLevel &bit)> fun) const
	Process bit ref level by lambda function. More...

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=QUIET) const
	Set nth bit ref level. More...

MoFEMErrorCode	setNthBitRefLevel (const int n, const bool b, int verb=QUIET) const
	Set nth bit ref level to all entities in database. More...

MoFEMErrorCode	shiftLeftBitRef (const int shift, const BitRefLevel mask=BitRefLevel().set(), int verb=DEFAULT_VERBOSITY) const
	left shift bit ref level More...

MoFEMErrorCode	shiftRightBitRef (const int shift, const BitRefLevel mask=BitRefLevel().set(), int verb=DEFAULT_VERBOSITY, MoFEMTypes mf=MF_ZERO) 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 int verb=QUIET) 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. 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. 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. 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. 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 entities. 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	updateRangeByChildren (const Range &parent, Range &child, MoFEMTypes bh=MF_ZERO)
	Update range by childrens. More...

MoFEMErrorCode	updateRangeByParent (const Range &child_ents, Range &parent_ents, MoFEMTypes bh=MF_ZERO)
	Update range by parents. More...

DEPRECATED MoFEMErrorCode	updateRange (const Range &parent, Range &child, MoFEMTypes bh=MF_ZERO)

Writing files
MoFEMErrorCode	writeBitLevel (const BitRefLevel bit, const BitRefLevel mask, const char file_name, const char file_type, const char *options, const bool check_for_empty=true) const
	Write bit ref level to file. More...

MoFEMErrorCode	writeBitLevelByDim (const BitRefLevel bit, const BitRefLevel mask, const int dim, const char file_name, const char file_type, const char *options, const bool check_for_empty=true) const
	Write bit ref level to file. More...

MoFEMErrorCode	writeBitLevelByType (const BitRefLevel bit, const BitRefLevel mask, const EntityType type, const char file_name, const char file_type, const char *options, const bool check_for_empty=true) const
	Write bit ref level to file. More...

MoFEMErrorCode	writeEntitiesNotInDatabase (const char file_name, const char file_type, const char *options, const bool check_for_empty=true) const
	Write ents not in database. More...

MoFEMErrorCode	writeEntitiesAllBitLevelsByType (const BitRefLevel mask, const EntityType type, const char file_name, const char file_type, const char *options)
	Write all entities by bit levels and type. More...

static MoFEMErrorCode	fixTagSize (moab::Interface &moab, bool *changed=nullptr)
	Fix tag size when BITREFLEVEL_SIZE of core library is different than file BITREFLEVEL_SIZE. More...

Get tag handles to data on the mesh
Tag	get_th_RefParentHandle () const

Tag	get_th_RefBitLevel () const

Additional Inherited Members
Static Public Member Functions inherited from MoFEM::UnknownInterface
static MoFEMErrorCode	getLibVersion (Version &version)
	Get library version. More...

static MoFEMErrorCode	getFileVersion (moab::Interface &moab, Version &version)
	Get database major version. More...

static MoFEMErrorCode	setFileVersion (moab::Interface &moab, Version version=Version(MoFEM_VERSION_MAJOR, MoFEM_VERSION_MINOR, MoFEM_VERSION_BUILD))
	Get database major version. More...

static MoFEMErrorCode	getInterfaceVersion (Version &version)
	Get database major version. More...

Detailed Description

Managing BitRefLevels.

Examples: bernstein_bezier_generate_base.cpp, build_large_problem.cpp, build_problems.cpp, child_and_parent.cpp, continuity_check_on_contact_prism_side_ele.cpp, continuity_check_on_skeleton_3d.cpp, dm_build_partitioned_mesh.cpp, dm_create_subdm.cpp, dm_partitioned_no_field.cpp, elasticity.cpp, elasticity_mixed_formulation.cpp, ep.cpp, field_blas.cpp, field_to_vertices.cpp, forces_and_sources_testing_edge_element.cpp, forces_and_sources_testing_flat_prism_element.cpp, forces_and_sources_testing_users_base.cpp, free_surface.cpp, hanging_node_approx.cpp, hcurl_divergence_operator_2d.cpp, level_set.cpp, mesh_cut.cpp, mesh_insert_interface_atom.cpp, mesh_smoothing.cpp, minimal_surface_area.cpp, navier_stokes.cpp, node_merge.cpp, nonlinear_dynamics.cpp, plot_base.cpp, prism_elements_from_surface.cpp, prism_polynomial_approximation.cpp, quad_polynomial_approximation.cpp, Remodeling.cpp, remove_entities_from_problem.cpp, remove_entities_from_problem_not_partitioned.cpp, simple_contact.cpp, simple_contact_thermal.cpp, split_sideset.cpp, test_jacobian_of_simple_contact_element.cpp, and unsaturated_transport.cpp.

Definition at line 21 of file BitRefManager.hpp.

Constructor & Destructor Documentation

◆ BitRefManager()

MoFEM::BitRefManager::BitRefManager ( const MoFEM::Core & core )

Definition at line 16 of file BitRefManager.cpp.

     : cOre(const_cast<MoFEM::Core &>(core)), dEbug(false) {
  
   if (!LogManager::checkIfChannelExist("BitRefSelf")) {
     auto core_log = logging::core::get();
     core_log->add_sink(
         LogManager::createSink(LogManager::getStrmSelf(), "BitRefSelf"));
     LogManager::setLog("BitRefSelf");
     core_log->add_sink(
         LogManager::createSink(LogManager::getStrmWorld(), "BitRefWorld"));
     LogManager::setLog("BitRefWorld");
     core_log->add_sink(
         LogManager::createSink(LogManager::getStrmSync(), "BitRefSync"));
     LogManager::setLog("BitRefSync");
     MOFEM_LOG_TAG("BitRefSelf", "BitRefManager");
     MOFEM_LOG_TAG("BitRefWorld", "BitRefManager");
     MOFEM_LOG_TAG("BitRefSync", "BitRefManager");
   }
  
   MOFEM_LOG_FUNCTION();
   MOFEM_LOG("BitRefWorld", Sev::noisy) << "BitRefManager interface created";
 }

◆ ~BitRefManager()

virtual MoFEM::BitRefManager::~BitRefManager ( )

virtualdefault

Destructor.

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

Examples: hanging_node_approx.cpp.

Definition at line 746 of file BitRefManager.cpp.

                                                                        {
   MoFEM::Interface &m_field = cOre;
   moab::Interface &moab(m_field.get_moab());
   MoFEMFunctionBeginHot;
  
   std::vector<EntityHandle> ents_vec;
   ents_vec.reserve(ents.size());
  
   std::vector<BitRefLevel *> tags_bits_ptr_vec(ents.size());
  
   Range swap_ents;
   auto hint = swap_ents.begin();
  
   for (Range::pair_iterator p_eit = ents.pair_begin(); p_eit != ents.pair_end();
        ++p_eit) {
  
     EntityHandle f = p_eit->first;
     const EntityHandle s = p_eit->second;
  
     // get bits on entities
     rval = moab.tag_get_by_ptr(cOre.get_th_RefBitLevel(), Range(f, s),
                                (const void **)(&*tags_bits_ptr_vec.begin()));
  
     if (rval == MB_SUCCESS) {
  
       auto bit_it = tags_bits_ptr_vec.begin();
  
       auto check = [&bit, &mask](const auto &entity_bit) -> bool {
         return
  
             (entity_bit & bit).any() &&
  
             ((entity_bit & mask) == entity_bit);
       };
  
       while (f != s + 1) {
  
         while (f != s + 1 && !check(**bit_it)) {
           ++bit_it;
           ++f;
         }
  
         if (f != s + 1) {
  
           const EntityHandle start = f;
  
           while (f != (s + 1) && check(**bit_it)) {
             ++bit_it;
             ++f;
           };
  
           hint = swap_ents.insert(hint, start, f - 1);
         }
       }
     }
   }
  
   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 908 of file BitRefManager.cpp.

                                                                            {
   MoFEM::Interface &m_field = cOre;
   auto ref_ents_ptr = m_field.get_ref_ents();
   MoFEMFunctionBeginHot;
   auto eit = ents.begin();
   for (; eit != ents.end();) {
     auto 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);
 }

◆ fixTagSize()

MoFEMErrorCode MoFEM::BitRefManager::fixTagSize	(	moab::Interface &	moab,
		bool *	changed = `nullptr`
	)

static

Fix tag size when BITREFLEVEL_SIZE of core library is different than file BITREFLEVEL_SIZE.

Returns: MoFEMErrorCode

Definition at line 1232 of file BitRefManager.cpp.

                                                                            {
   MoFEMFunctionBegin;
   MOFEM_LOG_CHANNEL("WORLD");
  
   if (changes)
     *changes = false;
  
   if (Tag th = 0; moab.tag_get_handle("_RefBitLevel", th) == MB_SUCCESS) {
  
     MOFEM_TAG_AND_LOG("WORLD", Sev::verbose, "BitRefManager") << "Tag found";
  
     auto get_old_tag = [&](auto &&name) {
       Tag th;
       CHK_MOAB_THROW(moab.tag_get_handle(name, th),
                      "bit ref level handle does not exist");
       return th;
     };
  
     auto get_new_tag = [&](auto &&name, auto &&def_val) {
       Tag th;
       CHK_MOAB_THROW(moab.tag_get_handle(
                          name, sizeof(BitRefLevel), MB_TYPE_OPAQUE, th,
                          MB_TAG_CREAT | MB_TAG_BYTES | MB_TAG_SPARSE, &def_val),
                      "can not create tag");
       return th;
     };
  
     int length;
     CHKERR moab.tag_get_length(get_old_tag("_RefBitLevel"), length);
  
     if (sizeof(BitRefLevel) != length) {
  
       if(changes)
         *changes = true;
  
       MOFEM_TAG_AND_LOG("WORLD", Sev::verbose, "BitRefManager")
           << "Fixing tag length";
  
       Range all_ents;
       CHKERR moab.get_entities_by_type(0, MBENTITYSET, all_ents, true);
       CHKERR moab.get_entities_by_handle(0, all_ents, true);
  
       auto process_tag = [&](auto &&name, auto &&def_val) {
         MoFEMFunctionBegin;
         auto tag_old = get_old_tag(name);
         auto get_bits = [&]() {
           std::vector<BitRefLevel> bits;
           bits.reserve(all_ents.size());
           auto it_bit = bits.begin();
           for (auto e : all_ents) {
             const void *data;
             int data_size;
             CHKERR moab.tag_get_by_ptr(tag_old, &e, 1, (const void **)&data,
                                        &data_size);
             bcopy(
                 data, &*it_bit,
                 std::min(sizeof(BitRefLevel), static_cast<size_t>(data_size)));
             ++it_bit;
           }
           return bits;
         };
         auto bits = get_bits();
         CHKERR moab.tag_delete(tag_old);
         auto tag_new = get_new_tag(name, def_val);
         auto it_bit = bits.begin();
         for (auto e : all_ents) {
           if (it_bit->any()) {
             CHKERR moab.tag_set_data(tag_new, &e, 1, &*it_bit);
           }
           ++it_bit;
         }
         MoFEMFunctionReturn(0);
       };
  
       CHKERR process_tag("_RefBitLevel", BitRefLevel() = 0);
       CHKERR process_tag("_RefBitLevelMask", BitRefLevel().set());
     }
   }
  
   MOFEM_LOG_CHANNEL("WORLD");
   MoFEMFunctionReturn(0);
 }

◆ get_th_RefBitLevel()

Tag MoFEM::BitRefManager::get_th_RefBitLevel ( ) const

inline

Definition at line 709 of file BitRefManager.hpp.

709 { return cOre.get_th_RefBitLevel(); }

◆ get_th_RefParentHandle()

Tag MoFEM::BitRefManager::get_th_RefParentHandle ( ) const

inline

Definition at line 706 of file BitRefManager.hpp.

                                             {
     return cOre.get_th_RefParentHandle();
   }

◆ getAllEntitiesNotInDatabase()

MoFEMErrorCode MoFEM::BitRefManager::getAllEntitiesNotInDatabase ( Range & ents ) const

Get all entities not in database.

Parameters

ents

Returns: MoFEMErrorCode

Definition at line 898 of file BitRefManager.cpp.

                                                                            {
   MoFEM::Interface &m_field = cOre;
   moab::Interface &moab = m_field.get_moab();
   MoFEMFunctionBegin;
   CHKERR moab.get_entities_by_handle(0, ents, false);
   ents = subtract(ents, ents.subset_by_type(MBENTITYSET));
   CHKERR filterEntitiesNotInDatabase(ents);
   MoFEMFunctionReturn(0);
 }

◆ getEntitiesByParentType()

MoFEMErrorCode MoFEM::BitRefManager::getEntitiesByParentType	(	const BitRefLevel	bit,
		const BitRefLevel	mask,
		const EntityType	type,
		Range &	ents,
		const int	verb = `QUIET`
	)		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 873 of file BitRefManager.cpp.

                                                                             {
   MoFEM::Interface &m_field = cOre;
   auto ref_ents_ptr = m_field.get_ref_ents();
   MoFEMFunctionBegin;
   auto &ref_ents = ref_ents_ptr->get<Ent_Ent_mi_tag>();
   auto it = ref_ents.lower_bound(get_id_for_min_type(type));
   auto hi_it = ref_ents.upper_bound(get_id_for_max_type(type));
   std::vector<EntityHandle> ents_vec;
   ents_vec.reserve(std::distance(it, hi_it));
   for (; it != hi_it; it++) {
     const BitRefLevel &ent_bit = it->get()->getBitRefLevel();
     if ((ent_bit & mask) == ent_bit && (ent_bit & bit).any())
       ents_vec.emplace_back(it->get()->getEnt());
   }
   ents.insert_list(ents_vec.begin(), ents_vec.end());
   if (verb > NOISY)
     MOFEM_LOG("BitRefSelf", Sev::noisy)
         << "getEntitiesByParentType: " << ents << endl;
   MoFEMFunctionReturn(0);
 }

◆ query_interface()

MoFEMErrorCode MoFEM::BitRefManager::query_interface	(	boost::typeindex::type_index	type_index,
		UnknownInterface **	iface
	)		const

virtual

Implements MoFEM::UnknownInterface.

Definition at line 9 of file BitRefManager.cpp.

                                                                {
   MoFEMFunctionBeginHot;
   *iface = const_cast<BitRefManager *>(this);
   MoFEMFunctionReturnHot(0);
 }

◆ writeBitLevel()

MoFEMErrorCode MoFEM::BitRefManager::writeBitLevel	(	const BitRefLevel	bit,
		const BitRefLevel	mask,
		const char *	file_name,
		const char *	file_type,
		const char *	options,
		const bool	check_for_empty = `true`
	)		const

Write bit ref level to file.

Parameters

bit	bit ref level
mask	mask of bit ref level
dim	dimension
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 634 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 getEntitiesByRefLevel(bit, mask, meshset);
   int nb_ents;
   CHKERR moab.get_number_entities_by_handle(meshset, nb_ents, true);
   if (check_for_empty && !nb_ents) {
     MOFEM_LOG("SELF", Sev::warning)
         << "No entities to save < " << file_name << " > in writeBitLevel";
     MoFEMFunctionReturnHot(0);
   }
  
   CHKERR moab.write_file(file_name, file_type, options, &meshset, 1);
   CHKERR moab.delete_entities(&meshset, 1);
   MoFEMFunctionReturn(0);
 }

◆ writeBitLevelByDim()

MoFEMErrorCode MoFEM::BitRefManager::writeBitLevelByDim	(	const BitRefLevel	bit,
		const BitRefLevel	mask,
		const int	dim,
		const char *	file_name,
		const char *	file_type,
		const char *	options,
		const bool	check_for_empty = `true`
	)		const

Write bit ref level to file.

Parameters

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

Returns: error code

Examples: hanging_node_approx.cpp.

Definition at line 660 of file BitRefManager.cpp.

                                                                     {
   MoFEM::Interface &m_field = cOre;
   moab::Interface &moab(m_field.get_moab());
   MoFEMFunctionBegin;
   Range ents;
   CHKERR getEntitiesByDimAndRefLevel(bit, mask, dim, ents);
   if (check_for_empty && ents.empty()) {
     MOFEM_LOG("SELF", Sev::warning)
         << "No entities to save < " << file_name << " > in writeBitLevelByDim";
     MoFEMFunctionReturnHot(0);
   }
   EntityHandle meshset;
   CHKERR moab.create_meshset(MESHSET_SET, meshset);
   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);
 }

◆ 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 bool	check_for_empty = `true`
	)		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 682 of file BitRefManager.cpp.

                                       {
   MoFEM::Interface &m_field = cOre;
   moab::Interface &moab(m_field.get_moab());
   MoFEMFunctionBegin;
   Range ents;
   CHKERR getEntitiesByTypeAndRefLevel(bit, mask, type, ents);
   if (check_for_empty && ents.empty()) {
     MOFEM_LOG("SELF", Sev::warning)
         << "No entities to save < " << file_name << " > in writeBitLevelByType";
     MoFEMFunctionReturnHot(0);
   }
   EntityHandle meshset;
   CHKERR moab.create_meshset(MESHSET_SET, meshset);
   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);
 }

◆ writeEntitiesAllBitLevelsByType()

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

Write all entities by bit levels and type.

Parameters

mask
type
file_name
file_type
options

Returns: MoFEMErrorCode

Definition at line 722 of file BitRefManager.cpp.

                                                 {
   MoFEMFunctionBegin;
   for (int ll = 0; ll != BITREFLEVEL_SIZE; ++ll) {
     std::string name = boost::lexical_cast<std::string>(ll) + "_" + file_name;
     CHKERR writeBitLevelByType(BitRefLevel().set(ll), mask, type, name.c_str(),
                                file_type, options, true);
   }
   MoFEMFunctionReturn(0);
 }

◆ writeEntitiesNotInDatabase()

MoFEMErrorCode MoFEM::BitRefManager::writeEntitiesNotInDatabase	(	const char *	file_name,
		const char *	file_type,
		const char *	options,
		const bool	check_for_empty = `true`
	)		const

Write ents not in database.

Parameters

file_name
file_type	for example "VTK"
options
check_for_empty

Returns: MoFEMErrorCode

Definition at line 704 of file BitRefManager.cpp.

                                       {
   MoFEM::Interface &m_field = cOre;
   moab::Interface &moab(m_field.get_moab());
   MoFEMFunctionBegin;
   EntityHandle meshset;
   Range ents;
   CHKERR getAllEntitiesNotInDatabase(ents);
   if (check_for_empty && ents.empty())
     MoFEMFunctionReturnHot(0);
   CHKERR moab.create_meshset(MESHSET_SET, meshset);
   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 26 of file BitRefManager.hpp.

◆ dEbug

bool MoFEM::BitRefManager::dEbug

Definition at line 27 of file BitRefManager.hpp.

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

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

Public Member Functions

Public Attributes

Setting and shifting bits

Entity handlers by bit ref level

Get adjacencies bit ref level

Update meshsets and ranges by children

Writing files

Get tag handles to data on the mesh

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ BitRefManager()

◆ ~BitRefManager()

Member Function Documentation

◆ filterEntitiesByRefLevel()

◆ filterEntitiesNotInDatabase()

◆ fixTagSize()

◆ get_th_RefBitLevel()

◆ get_th_RefParentHandle()

◆ getAllEntitiesNotInDatabase()

◆ getEntitiesByParentType()

◆ query_interface()

◆ writeBitLevel()

◆ writeBitLevelByDim()

◆ writeBitLevelByType()

◆ writeEntitiesAllBitLevelsByType()

◆ writeEntitiesNotInDatabase()

Member Data Documentation

◆ cOre

◆ dEbug