MeshsetsManager

Interface for meshsets with entities with data and boundary conditions. More...

Collaboration diagram for MeshsetsManager:

Files
file	MeshsetsManager.hpp
	MeshsetsManager interface.

Classes
struct	MoFEM::MeshsetsManager
	Interface for managing meshsets containing materials and boundary conditions. More...

Macros
#define	_IT_CUBITMESHSETS_FOR_LOOP_(MESHSET_MANAGER, IT)
	Iterator that loops over all the Cubit MeshSets in a moFEM field. More...
#define	_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(MESHSET_MANAGER, CUBITBCTYPE, IT)
	Iterator that loops over a specific Cubit MeshSet in a moFEM field. More...
#define	_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(MESHSET_MANAGER, CUBITBCTYPE, IT)
	Iterator that loops over a specific Cubit MeshSet having a particular BC meshset in a moFEM field. More...
#define	_IT_CUBITMESHSETS_BY_NAME_FOR_LOOP_(MESHSET_MANAGER, NAME, IT)
	Iterator that loops over Cubit BlockSet having a particular name. More...

Functions
CubitMeshSet_multiIndex::iterator	MoFEM::MeshsetsManager::getBegin () const
	get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT) More...
CubitMeshSet_multiIndex::iterator	MoFEM::MeshsetsManager::getEnd () const
	get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT) More...
CubitMeshsetByType::iterator	MoFEM::MeshsetsManager::getBegin (const unsigned int cubit_bc_type) const
	get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT) More...
CubitMeshsetByType::iterator	MoFEM::MeshsetsManager::getEnd (const unsigned int cubit_bc_type) const
	get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT) More...
CubitMeshsetByMask::iterator	MoFEM::MeshsetsManager::getBySetTypeBegin (const unsigned int cubit_bc_type) const
	get end iterator of cubit meshset of given type (instead you can use IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT) More...
CubitMeshsetByMask::iterator	MoFEM::MeshsetsManager::getBySetTypeEnd (const unsigned int cubit_bc_type) const
	get end iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT) More...
CubitMeshsetByName::iterator	MoFEM::MeshsetsManager::getBegin (const std::string &name) const
	get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT) More...
CubitMeshsetByName::iterator	MoFEM::MeshsetsManager::getEnd (const std::string &name) const
	get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT) More...
bool	MoFEM::MeshsetsManager::checkMeshset (const int ms_id, const CubitBCType cubit_bc_type) const
	check for CUBIT Id and CUBIT type More...
MoFEMErrorCode	MoFEM::MeshsetsManager::addMeshset (const CubitBCType cubit_bc_type, const int ms_id, const std::string name="")
	add cubit meshset More...
MoFEMErrorCode	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const int ms_id, const CubitBCType cubit_bc_type, const CubitMeshSets **cubit_meshset_ptr) const
	get cubit meshset More...
const CubitMeshSets *	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const int ms_id, const CubitBCType cubit_bc_type) const
	get cubit meshset More...
MoFEMErrorCode	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const CubitBCType cubit_bc_type, std::vector< const CubitMeshSets * > &vec_ptr) const
	Get vector of pointer by bc type. More...
std::vector< const CubitMeshSets * >	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const CubitBCType cubit_bc_type) const
	Get vector of pointer by bc type. More...
MoFEMErrorCode	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const string name, const CubitMeshSets **cubit_meshset_ptr) const
	get cubit meshset More...
MoFEMErrorCode	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const std::regex reg_exp_name, std::vector< const CubitMeshSets * > &vec_ptr) const
	Get vector of pointer to blocksets with name satisfying regular expression. More...
std::vector< const CubitMeshSets * >	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const std::regex reg_exp_name) const
	Get vector of pointer to blocksets with name satisfying regular expression. More...
MoFEMErrorCode	MoFEM::MeshsetsManager::getEntitiesByDimension (const int ms_id, const unsigned int cubit_bc_type, const int dimension, Range &entities, const bool recursive=true) const
	get entities from CUBIT/meshset of a particular entity dimension More...
MoFEMErrorCode	MoFEM::MeshsetsManager::getEntitiesByDimension (const int ms_id, const unsigned int cubit_bc_type, Range &entities, const bool recursive=true) const
	get entities related to CUBIT/meshset, More...
MoFEMErrorCode	MoFEM::MeshsetsManager::getMeshset (const int ms_id, const unsigned int cubit_bc_type, EntityHandle &meshset) const
	get meshset from CUBIT Id and CUBIT type More...
MoFEMErrorCode	MoFEM::MeshsetsManager::getMeshsetsByType (const unsigned int cubit_bc_type, Range &meshsets) const
	get all CUBIT meshsets by CUBIT type More...

Detailed Description

Interface for meshsets with entities with data and boundary conditions.

Macro Definition Documentation

_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_

#define _IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_	(		MESHSET_MANAGER,
			CUBITBCTYPE,
			IT
	)

Value:

  CubitMeshsetByType::iterator IT =                                            \
      MESHSET_MANAGER.get_meshsets_manager_ptr()->getBegin(CUBITBCTYPE);       \
  IT != MESHSET_MANAGER.get_meshsets_manager_ptr()->getEnd(CUBITBCTYPE);       \
  IT++

Iterator that loops over a specific Cubit MeshSet in a moFEM field.

Parameters

mField	moFEM Field
see	CubitBC (NODESET, SIDESET or BLOCKSET and more)
iterator

Examples

ElasticityMixedFormulation.hpp, HookeElement.cpp, add_cubit_meshsets.cpp, bone_adaptation.cpp, elasticity.cpp, field_blas.cpp, nonlinear_dynamics.cpp, and simple_elasticity.cpp.

Definition at line 49 of file MeshsetsManager.hpp.

_IT_CUBITMESHSETS_BY_NAME_FOR_LOOP_

#define _IT_CUBITMESHSETS_BY_NAME_FOR_LOOP_	(		MESHSET_MANAGER,
			NAME,
			IT
	)

Value:

  CubitMeshsetByName::iterator IT =                                            \
      MESHSET_MANAGER.get_meshsets_manager_ptr()->getBegin(NAME);              \
  IT != MESHSET_MANAGER.get_meshsets_manager_ptr()->getEnd(NAME);              \
  IT++

Iterator that loops over Cubit BlockSet having a particular name.

Parameters

MESHSET_MANAGER	meshset manager (works as well with Interface)
NAME	name
IT	iterator

Example:

for(_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(mField,"SOME_BLOCK_NAME",it) {
...
  } 

Definition at line 94 of file MeshsetsManager.hpp.

_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_

#define _IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_	(		MESHSET_MANAGER,
			CUBITBCTYPE,
			IT
	)

Value:

  CubitMeshsetByMask::iterator IT =                                            \
      MESHSET_MANAGER.get_meshsets_manager_ptr()->getBySetTypeBegin(           \
          CUBITBCTYPE);                                                        \
  IT != MESHSET_MANAGER.get_meshsets_manager_ptr()->getBySetTypeEnd(           \
            CUBITBCTYPE);                                                      \
  IT++

Iterator that loops over a specific Cubit MeshSet having a particular BC meshset in a moFEM field.

Parameters

MESHSET_MANAGER	meshset manager (works as well with Interface)
see	CubitBC (NODESET, SIDESET or BLOCKSET and more)
iterator

Example:

for(_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(mField,NODESET|DISPLACEMENTSET,it) {
...
  } 

Examples

MagneticElement.hpp, Remodeling.cpp, UnsaturatedFlow.hpp, add_cubit_meshsets.cpp, cell_forces.cpp, continuity_check_on_contact_prism_side_ele.cpp, cubit_bc_test.cpp, elasticity.cpp, forces_and_sources_testing_flat_prism_element.cpp, helmholtz.cpp, hertz_surface.cpp, mesh_cut.cpp, mesh_insert_interface_atom.cpp, mortar_contact_thermal.cpp, navier_stokes.cpp, photon_diffusion.cpp, simple_contact.cpp, simple_contact_thermal.cpp, test_jacobian_of_simple_contact_element.cpp, unsaturated_transport.cpp, and wavy_surface.cpp.

Definition at line 71 of file MeshsetsManager.hpp.

_IT_CUBITMESHSETS_FOR_LOOP_

#define _IT_CUBITMESHSETS_FOR_LOOP_	(		MESHSET_MANAGER,
			IT
	)

Value:

  CubitMeshSet_multiIndex::iterator IT =                                       \
      MESHSET_MANAGER.get_meshsets_manager_ptr()->getBegin();                  \
  IT != MESHSET_MANAGER.get_meshsets_manager_ptr()->getEnd();                  \
  IT++

Iterator that loops over all the Cubit MeshSets in a moFEM field.

Parameters

MESHSET_MANAGER	meshset manager (works as well with Interface)
iterator

Examples

add_cubit_meshsets.cpp, bone_adaptation.cpp, continuity_check_on_contact_prism_side_ele.cpp, forces_and_sources_testing_flat_prism_element.cpp, mesh_insert_interface_atom.cpp, mortar_contact_thermal.cpp, simple_contact.cpp, simple_contact_thermal.cpp, split_sideset.cpp, and unsaturated_transport.cpp.

Definition at line 34 of file MeshsetsManager.hpp.

Function Documentation

addMeshset()

MoFEMErrorCode MoFEM::MeshsetsManager::addMeshset	(	const CubitBCType	cubit_bc_type,
		const int	ms_id,
		const std::string	name = ""
	)

add cubit meshset

Parameters

see	CubitBC (NODESET, SIDESET or BLOCKSET and more)
ms_id	id of the BLOCKSET/SIDESET/BLOCKSET
name	of set

Examples

add_cubit_meshsets.cpp, and mesh_cut.cpp.

Definition at line 312 of file MeshsetsManager.cpp.

                                                                 {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
  if (checkMeshset(ms_id, cubit_bc_type)) {
    SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "such cubit meshset is already there", ms_id);
  }
 
  CubitMeshSets cmeshset(moab, cubit_bc_type, ms_id);
  if ((cmeshset.cubitBcType & CubitBCType(NODESET | SIDESET | BLOCKSET))
          .any()) {
    auto p = cubitMeshsets.insert(cmeshset);
    if (!p.second) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
              "meshset not inserted");
    }
    if (name.size() > 0) {
      bool success =
          cubitMeshsets.modify(p.first, CubitMeshSets_change_name(moab, name));
      if (!success) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
                "name to cubit meshset can not be set");
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

checkMeshset()

bool MoFEM::MeshsetsManager::checkMeshset	(	const int	ms_id,
		const CubitBCType	cubit_bc_type
	)		const

check for CUBIT Id and CUBIT type

Todo:

All cubit interface functions should be outsourced to dedicated interface

Parameters

ms_id	id of the BLOCKSET/SIDESET/BLOCKSET: from CUBIT
see	CubitBC (NODESET, SIDESET or BLOCKSET and more)

Examples

heat_equation.cpp, mesh_cut.cpp, mesh_smoothing.cpp, and reaction_diffusion.cpp.

Definition at line 287 of file MeshsetsManager.cpp.

                                                                          {
  auto miit =
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().find(
          boost::make_tuple(ms_id, cubit_bc_type.to_ulong()));
  if (miit !=
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().end()) {
    return true;
  }
  return false;
}

getBegin() [1/3]

CubitMeshSet_multiIndex::iterator MoFEM::MeshsetsManager::getBegin ( ) const

inline

get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT)

for(IT_CUBITMESHSETS_FOR_LOOP(mField,it) { ... }

Definition at line 243 of file MeshsetsManager.hpp.

                                                          {
    return cubitMeshsets.begin();
  }

getBegin() [2/3]

CubitMeshsetByName::iterator MoFEM::MeshsetsManager::getBegin ( const std::string &  name ) const

inline

get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT)

for(IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP(mField,NODESET|DISPLACEMENTSET,it) { ... }

Parameters

type	of meshset (NODESET, SIDESET or BLOCKSET and more)

Definition at line 351 of file MeshsetsManager.hpp.

                                                                          {
    return cubitMeshsets.get<CubitMeshsets_name>().lower_bound(name);
  }

getBegin() [3/3]

CubitMeshsetByType::iterator MoFEM::MeshsetsManager::getBegin ( const unsigned int  cubit_bc_type ) const

inline

get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT)

for(IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(mField,NODESET|DISPLACEMENTSET,it) { ... }

Parameters

type	of meshset (NODESET, SIDESET or BLOCKSET and more)

Definition at line 276 of file MeshsetsManager.hpp.

                                                   {
    return cubitMeshsets.get<CubitMeshsetType_mi_tag>().lower_bound(
        cubit_bc_type);
  }

getBySetTypeBegin()

CubitMeshsetByMask::iterator MoFEM::MeshsetsManager::getBySetTypeBegin ( const unsigned int  cubit_bc_type ) const

inline

get end iterator of cubit meshset of given type (instead you can use IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT)

for(IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(mField,NODESET,it) { ... }

Parameters

type	of meshset (NODESET, SIDESET or BLOCKSET and more)

Definition at line 314 of file MeshsetsManager.hpp.

                                                            {
    return cubitMeshsets.get<CubitMeshsetMaskedType_mi_tag>().lower_bound(
        cubit_bc_type);
  }

getBySetTypeEnd()

CubitMeshsetByMask::iterator MoFEM::MeshsetsManager::getBySetTypeEnd ( const unsigned int  cubit_bc_type ) const

inline

get end iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT)

for(IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(mField,NODESET,it) { ... }

Parameters

type	of meshset (NODESET, SIDESET or BLOCKSET and more)

Definition at line 332 of file MeshsetsManager.hpp.

                                                          {
    return cubitMeshsets.get<CubitMeshsetMaskedType_mi_tag>().upper_bound(
        cubit_bc_type);
  }

getCubitMeshsetPtr() [1/7]

std::vector< const CubitMeshSets * > MoFEM::MeshsetsManager::getCubitMeshsetPtr

(

const CubitBCType

cubit_bc_type

)

const

Get vector of pointer by bc type.

Returns

MoFEMErrorCode

Definition at line 569 of file MeshsetsManager.cpp.

                                                                         {
  std::vector<const CubitMeshSets *> vec_ptr;
  CHK_MOAB_THROW(getCubitMeshsetPtr(cubit_bc_type, vec_ptr),
                 "Error in getting meshsets by name");
  return vec_ptr;
}

getCubitMeshsetPtr() [2/7]

MoFEMErrorCode MoFEM::MeshsetsManager::getCubitMeshsetPtr	(	const CubitBCType	cubit_bc_type,
		std::vector< const CubitMeshSets * > &	vec_ptr
	)		const

Get vector of pointer by bc type.

Parameters

std::vector<const

CubitMeshSets *>

Returns

MoFEMErrorCode

Definition at line 551 of file MeshsetsManager.cpp.

                                                     {
  MoFEMFunctionBegin;
  for (auto &c : cubitMeshsets) {
    if ((c.getBcType() & cubit_bc_type) == cubit_bc_type) {
      vec_ptr.push_back(&c);
    }
  }
  std::sort(vec_ptr.begin(), vec_ptr.end(),
            [](const CubitMeshSets *a, const CubitMeshSets *b) {
              return a->getMeshsetId() < b->getMeshsetId();
            });
 
  MoFEMFunctionReturn(0);
}

getCubitMeshsetPtr() [3/7]

const CubitMeshSets * MoFEM::MeshsetsManager::getCubitMeshsetPtr	(	const int	ms_id,
		const CubitBCType	cubit_bc_type
	)		const

get cubit meshset

Definition at line 524 of file MeshsetsManager.cpp.

                                                                           {
  const CubitMeshSets *cubit_meshset_ptr;
  CHK_THROW_MESSAGE(
      getCubitMeshsetPtr(ms_id, cubit_bc_type, &cubit_meshset_ptr),
      "Get not get meshset");
  return cubit_meshset_ptr;
}

getCubitMeshsetPtr() [4/7]

MoFEMErrorCode MoFEM::MeshsetsManager::getCubitMeshsetPtr	(	const int	ms_id,
		const CubitBCType	cubit_bc_type,
		const CubitMeshSets **	cubit_meshset_ptr
	)		const

get cubit meshset

Examples

EshelbianPlasticity.cpp, dynamic_first_order_con_law.cpp, plastic.cpp, and thermo_elastic.cpp.

Definition at line 505 of file MeshsetsManager.cpp.

                                                   {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  auto miit =
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().find(
          boost::make_tuple(ms_id, cubit_bc_type.to_ulong()));
  if (miit !=
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().end()) {
    *cubit_meshset_ptr = &*miit;
  } else {
    SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "msId = %d is not there", ms_id);
  }
  MoFEMFunctionReturn(0);
}

getCubitMeshsetPtr() [5/7]

std::vector< const CubitMeshSets * > MoFEM::MeshsetsManager::getCubitMeshsetPtr

(

const std::regex

reg_exp_name

)

const

Get vector of pointer to blocksets with name satisfying regular expression.

for (auto m :
     m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(std::regex(
 
         (boost::format("%s(.*)") % block_name).str()
 
             ))
 
) {
  MOFEM_LOG("WORLD", Sev::inform) << m->getName();
 
  Rants ents;
  m_field.get_moab().get_entities_by_dimension(
      m->getMeshset(), 0, ents, false);
 
  auto print_vets = [](boost::shared_ptr<FieldEntity> ent_ptr) {
    MoFEMFunctionBegin;
    if(!(ent_ptr->getPStatus() & PSTATUS_NOT_OWNED)) {
      MOFEM_LOG("SYNC", Sev::inform) << ents_ptr->getName() << ": " <<
ent_ptr->getEntFieldData();
    }
    MoFEMFunctionReturn(0);
  };
 
  CHKERR m_field.getInterface<FieldBlas>()->fieldLambdaOnEntities(
      print_vets, "VELOCITIES", &ents);
  CHKERR m_field.getInterface<FieldBlas>()->fieldLambdaOnEntities(
      print_vets, "DISPLACEMENTS", &ents);
  MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::ifrom);
}

Parameters

reg_exp_name

Returns

std::vector<const CubitMeshSets *>

Definition at line 596 of file MeshsetsManager.cpp.

                                                                     {
  std::vector<const CubitMeshSets *> vec_ptr;
  CHK_MOAB_THROW(getCubitMeshsetPtr(reg_exp_name, vec_ptr),
                 "Error in getting meshsets by name");
  return vec_ptr;
}

getCubitMeshsetPtr() [6/7]

MoFEMErrorCode MoFEM::MeshsetsManager::getCubitMeshsetPtr	(	const std::regex	reg_exp_name,
		std::vector< const CubitMeshSets * > &	vec_ptr
	)		const

Get vector of pointer to blocksets with name satisfying regular expression.

Parameters

reg_exp_name
std::vector<const	CubitMeshSets *>

Returns

MoFEMErrorCode

Definition at line 576 of file MeshsetsManager.cpp.

                                                     {
  MoFEMFunctionBegin;
  auto r =
      cubitMeshsets.get<CubitMeshsetMaskedType_mi_tag>().equal_range(BLOCKSET);
  for (; r.first != r.second; ++r.first) {
    const auto name = r.first->getName();
    if (std::regex_match(name, reg_exp_name)) {
      vec_ptr.push_back(&*r.first);
    }
  }
  std::sort(vec_ptr.begin(), vec_ptr.end(),
            [](const CubitMeshSets *a, const CubitMeshSets *b) {
              return a->getMeshsetId() < b->getMeshsetId();
            });
  MoFEMFunctionReturn(0);
}

getCubitMeshsetPtr() [7/7]

MoFEMErrorCode MoFEM::MeshsetsManager::getCubitMeshsetPtr	(	const string	name,
		const CubitMeshSets **	cubit_meshset_ptr
	)		const

get cubit meshset

Definition at line 533 of file MeshsetsManager.cpp.

                                                                      {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  auto miit = cubitMeshsets.get<CubitMeshsets_name>().lower_bound(name);
  auto hi_miit = cubitMeshsets.get<CubitMeshsets_name>().upper_bound(name);
  if (std::distance(miit, hi_miit) == 0) {
    SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "meshset name <%s> is not there", name.c_str());
  }
  if (std::distance(miit, hi_miit) > 1) {
    SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "more that one meshser of that name <%s>", name.c_str());
  }
  *cubit_meshset_ptr = &*miit;
  MoFEMFunctionReturn(0);
}

getEnd() [1/3]

CubitMeshSet_multiIndex::iterator MoFEM::MeshsetsManager::getEnd ( ) const

inline

get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT)

for(IT_CUBITMESHSETS_FOR_LOOP(mField,it) { ... }

Definition at line 257 of file MeshsetsManager.hpp.

                                                 {
    return cubitMeshsets.end();
  }

getEnd() [2/3]

CubitMeshsetByName::iterator MoFEM::MeshsetsManager::getEnd ( const std::string &  name ) const

inline

get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT)

for(IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP(mField,NODESET|DISPLACEMENTSET,it) { ... }

Parameters

type	of meshset (NODESET, SIDESET or BLOCKSET and more)

Definition at line 369 of file MeshsetsManager.hpp.

                                                                        {
    return cubitMeshsets.get<CubitMeshsets_name>().upper_bound(name);
  }

getEnd() [3/3]

CubitMeshsetByType::iterator MoFEM::MeshsetsManager::getEnd ( const unsigned int  cubit_bc_type ) const

inline

get begin iterator of cubit mehset of given type (instead you can use IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(MFIELD,CUBITBCTYPE,IT)

for(IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP(mField,NODESET|DISPLACEMENTSET,it) { ... }

Parameters

type	of meshset (NODESET, SIDESET or BLOCKSET and more)

Definition at line 296 of file MeshsetsManager.hpp.

                                                 {
    return cubitMeshsets.get<CubitMeshsetType_mi_tag>().upper_bound(
        cubit_bc_type);
  }

getEntitiesByDimension() [1/2]

MoFEMErrorCode MoFEM::MeshsetsManager::getEntitiesByDimension	(	const int	ms_id,
		const unsigned int	cubit_bc_type,
		const int	dimension,
		Range &	entities,
		const bool	recursive = true
	)		const

get entities from CUBIT/meshset of a particular entity dimension

Nodeset can contain nodes, edges, triangles and tets. This applies to other meshsets too. The nodeset's meshset contain the nodes in the MIDDLE of the surface or volume which is done by default in Cubit, Hence if all nodes on a particular nodeset are required, one should get all triangles or tetrahedrons for which the nodeset was create in Cubit, and get all the connectivities of tris/tets.

Parameters

ms_id	id of the BLOCKSET/SIDESET/BLOCKSET: from CUBIT
see	CubitBC (NODESET, SIDESET or BLOCKSET and more)
dimensions	(0 - Nodes, 1 - Edges, 2 - Faces, 3 - Volume(tetrahedral))
Range	containing the retrieved entities
recursive	If true, meshsets containing meshsets are queried recursively. Returns the contents of meshsets, but not the meshsets themselves if true.

Examples

mesh_cut.cpp, and prism_elements_from_surface.cpp.

Definition at line 603 of file MeshsetsManager.cpp.

                                                 {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
  auto miit =
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().find(
          boost::make_tuple(msId, cubit_bc_type));
  if (miit !=
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().end()) {
    CHKERR miit->getMeshsetIdEntitiesByDimension(moab, dimension, entities,
                                                 recursive);
  } else {
    SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "msId = %d is not there", msId);
  }
  MoFEMFunctionReturn(0);
}

getEntitiesByDimension() [2/2]

MoFEMErrorCode MoFEM::MeshsetsManager::getEntitiesByDimension	(	const int	ms_id,
		const unsigned int	cubit_bc_type,
		Range &	entities,
		const bool	recursive = true
	)		const

get entities related to CUBIT/meshset,

NODESET will get Vertices only, even if the NODESET contains edges, tris and tets SIDESET will get Tris, BLOCKSET will get Tets, DISPLACEMENTSET and FORCESET are stored in NODESET, PRESSURESET is stored in Sideset.

Parameters

ms_id	id of the BLOCKSET/SIDESET/BLOCKSET: from CUBIT
see	CubitBC (NODESET, SIDESET or BLOCKSET and more)
Range	containing the retrieved entities related to the
recursive	If true, meshsets containing meshsets are queried recursively. Returns the contents of meshsets, but not the meshsets themselves if true.

Definition at line 623 of file MeshsetsManager.cpp.

                                {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
  auto miit =
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().find(
          boost::make_tuple(ms_id, cubit_bc_type));
  if (miit !=
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().end()) {
    CHKERR miit->getMeshsetIdEntitiesByDimension(moab, entities, recursive);
  } else {
    SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "ms_id = %d is not there", ms_id);
  }
  MoFEMFunctionReturn(0);
}

getMeshset()

MoFEMErrorCode MoFEM::MeshsetsManager::getMeshset	(	const int	ms_id,
		const unsigned int	cubit_bc_type,
		EntityHandle &	meshset
	)		const

get meshset from CUBIT Id and CUBIT type

Parameters

ms_id	id of the BLOCKSET/SIDESET/BLOCKSET: from CUBIT
see	CubitBC (NODESET, SIDESET or BLOCKSET and more)
meshset	where to store the retrieved entities

Examples

free_surface.cpp, and mesh_cut.cpp.

Definition at line 642 of file MeshsetsManager.cpp.

                                                                        {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  auto miit =
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().find(
          boost::make_tuple(ms_id, cubit_bc_type));
  if (miit !=
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().end()) {
    meshset = miit->meshset;
  } else {
    SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "ms_id = %d is not there", ms_id);
  }
  MoFEMFunctionReturn(0);
}

getMeshsetsByType()

MoFEMErrorCode MoFEM::MeshsetsManager::getMeshsetsByType	(	const unsigned int	cubit_bc_type,
		Range &	meshsets
	)		const

get all CUBIT meshsets by CUBIT type

Parameters

see	CubitBC (NODESET, SIDESET or BLOCKSET and more).
meshsets	is range of meshsets

Definition at line 676 of file MeshsetsManager.cpp.

                                                          {
  MoFEMFunctionBegin;
  auto miit =
      cubitMeshsets.get<CubitMeshsetType_mi_tag>().lower_bound(cubit_bc_type);
  auto hi_miit =
      cubitMeshsets.get<CubitMeshsetType_mi_tag>().upper_bound(cubit_bc_type);
  for (; miit != hi_miit; miit++) {
    meshsets.insert(miit->meshset);
  }
  MoFEMFunctionReturn(0);
}