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.
#define	_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(MESHSET_MANAGER, CUBITBCTYPE, IT)
	Iterator that loops over a specific Cubit MeshSet in a moFEM field.
#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.
#define	_IT_CUBITMESHSETS_BY_NAME_FOR_LOOP_(MESHSET_MANAGER, NAME, IT)
	Iterator that loops over Cubit BlockSet having a particular name.

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)
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)
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)
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)
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)
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)
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)
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)
bool	MoFEM::MeshsetsManager::checkMeshset (const int ms_id, const CubitBCType cubit_bc_type) const
	check for CUBIT Id and CUBIT type
MoFEMErrorCode	MoFEM::MeshsetsManager::addMeshset (const CubitBCType cubit_bc_type, const int ms_id, const std::string name="")
	add cubit meshset
MoFEMErrorCode	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const int ms_id, const CubitBCType cubit_bc_type, const CubitMeshSets **cubit_meshset_ptr) const
	get cubit meshset
const CubitMeshSets *	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const int ms_id, const CubitBCType cubit_bc_type) const
	get cubit meshset
MoFEMErrorCode	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const CubitBCType cubit_bc_type, std::vector< const CubitMeshSets * > &vec_ptr) const
	Get vector of pointer by bc type.
std::vector< const CubitMeshSets * >	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const CubitBCType cubit_bc_type) const
	Get vector of pointer by bc type.
MoFEMErrorCode	MoFEM::MeshsetsManager::getCubitMeshsetPtr (const string name, const CubitMeshSets **cubit_meshset_ptr) const
	get cubit meshset
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.
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.
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
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,
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
MoFEMErrorCode	MoFEM::MeshsetsManager::getMeshsetsByType (const unsigned int cubit_bc_type, Range &meshsets) const
	get all CUBIT meshsets by CUBIT type

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
	)

#include <src/interfaces/MeshsetsManager.hpp>

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.

                       : public UnknownInterface {
 
  MoFEMErrorCode query_interface(boost::typeindex::type_index type_index,
                                 UnknownInterface **iface) const;
 
  MoFEM::Core &cOre;
  MeshsetsManager(const MoFEM::Core &core);
  virtual ~MeshsetsManager() = default;
 
  /**
   * \brief get tags handlers used on meshsets
 
   * On meshsets range of tags in set. Depending on tag type and data on that
   * tag type of meshset could be determined. This function get hanldes to
   * tags.
   *
   * Most of the tags are followinf convention used by MoAB or Cubit and other
   * meshing softwares, f.e. gmesh.
 
   */
  MoFEMErrorCode getTags(int verb = -1);
 
  static MoFEMErrorCode getTags(moab::Interface &moab, Tag &nsTag, Tag &ssTag,
                                Tag &nsTag_data, Tag &ssTag_data, Tag &bhTag,
                                Tag &bhTag_header, int verb = -1);
 
  /**
   * \brief get tag handle used to store "id" of NODESET
   */
  inline Tag get_nsTag() const { return nsTag; }
 
  /**
   * \brief get tag handle used to store "id" of SIDESET
   */
  inline Tag get_ssTag() const { return ssTag; }
 
  /**
   * \brief get tag handle used to store boundary data on NODESET
   */
  inline Tag get_nsTag_data() const { return nsTag_data; }
 
  /**
   * \brief get tag handle used to store boundary data on SIDESET
   */
  inline Tag get_ssTag_data() const { return ssTag_data; }
 
  /**
   * \brief get tag handle used to store "id" of BLOCKSET
   */
  inline Tag get_bhTag() const { return bhTag; }
 
  /**
   * \brief get tag handle used to store of block set header (Used by Cubit)
   */
  inline Tag get_bhTag_header() const { return bhTag_header; }
 
  /**
   * \brief return pointer to meshset manager
   */
  MeshsetsManager *get_meshsets_manager_ptr() { return this; }
 
  /**
   * \brief return pointer to meshset manager
   */
  const MeshsetsManager *get_meshsets_manager_ptr() const { return this; }
 
  /**
   * \brief clear multi-index container
   * @return error code
   */
  MoFEMErrorCode clearMap();
 
  /**
   * \brier initialize container form data on mesh
   * @return [description]
   */
  MoFEMErrorCode initialiseDatabaseFromMesh(int verb = DEFAULT_VERBOSITY);
 
  /**
   * @brief Read meshsets
   *
   * @param verb
   * @return MoFEMErrorCode
   */
  MoFEMErrorCode readMeshsets(int verb = DEFAULT_VERBOSITY);
 
  /**
   * @brief Broadcast meshsets
   *
   * @param verb
   * @return MoFEMErrorCode
   */
  MoFEMErrorCode broadcastMeshsets(int verb = DEFAULT_VERBOSITY);
 
  template <class CUBIT_BC_DATA_TYPE>
  MoFEMErrorCode printBcSet(CUBIT_BC_DATA_TYPE &data,
                            unsigned long int type) const;
 
  /**
   * \brief print meshsets with displacement boundary conditions data
   * structure
   */
  MoFEMErrorCode printDisplacementSet() const;
 
  /**
   * \brief print meshsets with pressure boundary conditions data structure
   */
  MoFEMErrorCode printPressureSet() const;
 
  /**
   * \brief print meshsets with force boundary conditions data structure
   */
  MoFEMErrorCode printForceSet() const;
 
  /**
   * \brief print meshsets with temperature boundary conditions data structure
   */
  MoFEMErrorCode printTemperatureSet() const;
 
  /**
   * \brief print meshsets with heat flux boundary conditions data structure
   */
  MoFEMErrorCode printHeatFluxSet() const;
 
  /**
   * \brief print meshsets with material data structure set on it
   */
  MoFEMErrorCode printMaterialsSet() const;
 
  inline CubitMeshSet_multiIndex &getMeshsetsMultindex() {
    return cubitMeshsets;
  }
 
  /**
   * \ingroup mofem_meshset_mng
   * \brief 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) {
   *    ...
   * }
   *
   */
  inline CubitMeshSet_multiIndex::iterator getBegin() const {
    return cubitMeshsets.begin();
  }
 
  /**
   * \ingroup mofem_meshset_mng
   * \brief 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) {
   *    ...
   * }
   *
   */
  CubitMeshSet_multiIndex::iterator getEnd() const {
    return cubitMeshsets.end();
  }
 
  /**
    * \brief get begin iterator of cubit mehset of given type (instead you can
    use _IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(MFIELD,CUBITBCTYPE,IT)
    * \ingroup mofem_meshset_mng
 
    *
    *
    for(_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(mField,NODESET|DISPLACEMENTSET,it)
    {
    *   ...
    * }
    *
    * \param  type of meshset (NODESET, SIDESET or BLOCKSET and more)
    */
  inline CubitMeshsetByType::iterator
  getBegin(const unsigned int cubit_bc_type) const {
    return cubitMeshsets.get<CubitMeshsetType_mi_tag>().lower_bound(
        cubit_bc_type);
  }
 
  /**
    * \brief get begin iterator of cubit mehset of given type (instead you can
    use _IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(MFIELD,CUBITBCTYPE,IT)
    * \ingroup mofem_meshset_mng
 
    *
    *
    for(_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(mField,NODESET|DISPLACEMENTSET,it)
    {
    *   ...
    * }
    *
    * \param  type of meshset (NODESET, SIDESET or BLOCKSET and more)
    */
  inline CubitMeshsetByType::iterator
  getEnd(const unsigned int cubit_bc_type) const {
    return cubitMeshsets.get<CubitMeshsetType_mi_tag>().upper_bound(
        cubit_bc_type);
  }
 
  /**
    * \brief get end iterator of cubit meshset of given type (instead you can
    use _IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(MFIELD,CUBITBCTYPE,IT)
    * \ingroup mofem_meshset_mng
 
    *
    * for(_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(mField,NODESET,it) {
    *   ...
    * }
    *
    * \param  type of meshset (NODESET, SIDESET or BLOCKSET and more)
    */
  inline CubitMeshsetByMask::iterator
  getBySetTypeBegin(const unsigned int cubit_bc_type) const {
    return cubitMeshsets.get<CubitMeshsetMaskedType_mi_tag>().lower_bound(
        cubit_bc_type);
  }
 
  /**
    * \brief get end iterator of cubit mehset of given type (instead you can
    use _IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(MFIELD,CUBITBCTYPE,IT)
    * \ingroup mofem_meshset_mng
 
    *
    * for(_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(mField,NODESET,it) {
    *   ...
    * }
    *
    * \param  type of meshset (NODESET, SIDESET or BLOCKSET and more)
    */
  inline CubitMeshsetByMask::iterator
  getBySetTypeEnd(const unsigned int cubit_bc_type) const {
    return cubitMeshsets.get<CubitMeshsetMaskedType_mi_tag>().upper_bound(
        cubit_bc_type);
  }
 
  /**
    * \brief get begin iterator of cubit mehset of given type (instead you can
    use _IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(MFIELD,CUBITBCTYPE,IT)
    * \ingroup mofem_meshset_mng
 
    *
    *
    for(_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(mField,NODESET|DISPLACEMENTSET,it)
    {
    *   ...
    * }
    *
    * \param  type of meshset (NODESET, SIDESET or BLOCKSET and more)
    */
  inline CubitMeshsetByName::iterator getBegin(const std::string &name) const {
    return cubitMeshsets.get<CubitMeshsets_name>().lower_bound(name);
  }
 
  /**
    * \brief get begin iterator of cubit mehset of given type (instead you can
    use _IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(MFIELD,CUBITBCTYPE,IT)
    * \ingroup mofem_meshset_mng
 
    *
    *
    for(_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(mField,NODESET|DISPLACEMENTSET,it)
    {
    *   ...
    * }
    *
    * \param  type of meshset (NODESET, SIDESET or BLOCKSET and more)
    */
  inline CubitMeshsetByName::iterator getEnd(const std::string &name) const {
    return cubitMeshsets.get<CubitMeshsets_name>().upper_bound(name);
  }
 
  /**
    * \brief check for CUBIT Id and CUBIT type
    * \ingroup mofem_meshset_mng
 
    \todo All cubit interface functions should be outsourced to dedicated
    interface
 
    * \param ms_id id of the BLOCKSET/SIDESET/BLOCKSET: from CUBIT
    * \param  see CubitBC (NODESET, SIDESET or BLOCKSET and more)
    */
  bool checkMeshset(const int ms_id, const CubitBCType cubit_bc_type) const;
 
  /**
   * \brief check if meshset of given name exist
   * @param  name name of meshset
   * @return      error code
   */
  bool checkMeshset(const string name,
                    int *const number_of_meshsets_ptr = NULL) const;
 
  /**
    * \brief add cubit meshset
    * \ingroup mofem_meshset_mng
 
    *
    * \param see CubitBC (NODESET, SIDESET or BLOCKSET and more)
    * \param ms_id id of the BLOCKSET/SIDESET/BLOCKSET
    * \param name of set
 
    */
  MoFEMErrorCode addMeshset(const CubitBCType cubit_bc_type, const int ms_id,
                            const std::string name = "");
 
  /**
   * \brief add entities to cubit meshset
   * @param  cubit_bc_type type of meshset, f.e. NODESET, SIDESET or BLOCKSET
   * @param  ms_id         id of meshset
   * @param  ents          entities to add
   * @return               error code
   */
  MoFEMErrorCode addEntitiesToMeshset(const CubitBCType cubit_bc_type,
                                      const int ms_id, const Range &ents);
 
  /**
   * \brief add entities to cubit meshset
   * @param  cubit_bc_type type of meshset, f.e. NODESET, SIDESET or BLOCKSET
   * @param  ms_id         id of meshset
   * @param  ents          pointer to entities array
   * @param  nb_ents       number of entities in array
   * @return               error code
   */
  MoFEMErrorCode addEntitiesToMeshset(const CubitBCType cubit_bc_type,
                                      const int ms_id, const EntityHandle *ents,
                                      const int nb_ents);
 
  /**
   * \brief set attributes to cubit meshset
   * @param  cubit_bc_type type of meshset, see CubitBC, i.e. BLOCKSET,
   * NODESET, SIDESET
   * @param  ms_id         id of meshset
   * @param  attributes    attributes
   * @param  name          set name to blockset
   * @return               error code
   */
  MoFEMErrorCode setAtributes(const CubitBCType cubit_bc_type, const int ms_id,
                              const std::vector<double> &attributes,
                              const std::string name = "");
 
  /**
   * \brief set (material) data structure to cubit meshset
   * @param  cubit_bc_type type of meshset, see CubitBC, i.e. BLOCKSET,
   * NODESET, SIDESET
   * @param  ms_id         id of meshset
   * @param  attributes    attributes
   * @return               error code
   */
  MoFEMErrorCode setAtributesByDataStructure(const CubitBCType cubit_bc_type,
                                             const int ms_id,
                                             const GenericAttributeData &data,
                                             const std::string name = "");
 
  /**
   * \brief set boundary data structure to meshset
   * @param  cubit_bc_type type of meshset, see CubitBC, i.e. BLOCKSET,
   * NODESET, SIDESET
   * @param  ms_id         id of meshset
   * @param  data          data structure
   * @return               error code
   */
  MoFEMErrorCode setBcData(const CubitBCType cubit_bc_type, const int ms_id,
                           const GenericCubitBcData &data);
 
  /**
    * \brief delete cubit meshset
    * \ingroup mopfem_bc
 
    *
    * \param  see CubitBC (NODESET, SIDESET or BLOCKSET and more)
    * \param ms_id id of the BLOCKSET/SIDESET/BLOCKSET: from CUBIT
    *
    */
  MoFEMErrorCode deleteMeshset(const CubitBCType cubit_bc_type, const int ms_id,
                               const MoFEMTypes bh = MF_EXIST);
 
  /**
   * \brief get cubit meshset
   * \ingroup mofem_meshset_mng
   *
   *
   */
  MoFEMErrorCode
  getCubitMeshsetPtr(const int ms_id, const CubitBCType cubit_bc_type,
                     const CubitMeshSets **cubit_meshset_ptr) const;
 
  /**
   * \brief get cubit meshset
   * \ingroup mofem_meshset_mng
   *
   *
   */
  const CubitMeshSets *
  getCubitMeshsetPtr(const int ms_id, const CubitBCType cubit_bc_type) const;
 
  /**
   * @brief Get vector of pointer by bc type
   *
   * \ingroup mofem_meshset_mng
   *
   * @param std::vector<const CubitMeshSets *>
   * @return MoFEMErrorCode
   */
  MoFEMErrorCode
  getCubitMeshsetPtr(const CubitBCType cubit_bc_type,
                     std::vector<const CubitMeshSets *> &vec_ptr) const;
 
/**
   * @brief Get vector of pointer by bc type
   *
   * \ingroup mofem_meshset_mng
   *
   * @return MoFEMErrorCode
   */
  std::vector<const CubitMeshSets *>
  getCubitMeshsetPtr(const CubitBCType cubit_bc_type) const;
 
  /**
   * \brief get cubit meshset
   *
   * \ingroup mofem_meshset_mng
   */
  MoFEMErrorCode
  getCubitMeshsetPtr(const string name,
                     const CubitMeshSets **cubit_meshset_ptr) const;
 
  /**
   * @brief Get vector of pointer to blocksets with name satisfying regular
   * expression.
   *
   * \ingroup mofem_meshset_mng
   *
   * @param reg_exp_name
   * @param std::vector<const CubitMeshSets *>
   * @return MoFEMErrorCode
   */
  MoFEMErrorCode
  getCubitMeshsetPtr(const std::regex reg_exp_name,
                     std::vector<const CubitMeshSets *> &vec_ptr) const;
 
  /**
   * @brief Get vector of pointer to blocksets with name satisfying regular
   * expression.
   *
   * \ingroup mofem_meshset_mng
   *
   * \code
  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);
  }
   * \endcode
   *
   * @param reg_exp_name
   * @return std::vector<const CubitMeshSets *>
   */
  std::vector<const CubitMeshSets *>
  getCubitMeshsetPtr(const std::regex reg_exp_name) const;
 
  /**
    * \brief get entities from CUBIT/meshset of a particular entity dimension
    * \ingroup mofem_meshset_mng
 
    * 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.
 
    * \param ms_id id of the BLOCKSET/SIDESET/BLOCKSET: from CUBIT
    * \param  see CubitBC (NODESET, SIDESET or BLOCKSET and more)
    * \param dimensions (0 - Nodes, 1 - Edges, 2 - Faces, 3 -
    Volume(tetrahedral))
    * \param Range containing the retrieved entities
    * \param recursive If true, meshsets containing meshsets are queried
    recursively. Returns the contents of meshsets, but not the meshsets
    themselves if true.
    */
  MoFEMErrorCode getEntitiesByDimension(const int ms_id,
                                        const unsigned int cubit_bc_type,
                                        const int dimension, Range &entities,
                                        const bool recursive = true) const;
 
  /**
    * \brief get entities related to CUBIT/meshset,
    * \ingroup mofem_meshset_mng
 
    * 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.
 
    * \param ms_id id of the BLOCKSET/SIDESET/BLOCKSET: from CUBIT
    * \param  see CubitBC (NODESET, SIDESET or BLOCKSET and more)
    * \param Range containing the retrieved entities related to the
    * \param recursive If true, meshsets containing meshsets are queried
    recursively.  Returns the contents of meshsets, but not the meshsets
    themselves if true.
    */
  MoFEMErrorCode getEntitiesByDimension(const int ms_id,
                                        const unsigned int cubit_bc_type,
                                        Range &entities,
                                        const bool recursive = true) const;
 
  /**
   * \ingroup mofem_meshset_mng
   * \brief get meshset from CUBIT Id and CUBIT type
   *
   * \param ms_id id of the BLOCKSET/SIDESET/BLOCKSET: from CUBIT
   * \param  see CubitBC (NODESET, SIDESET or BLOCKSET and more)
   * \param meshset where to store the retrieved entities
   */
  MoFEMErrorCode getMeshset(const int ms_id, const unsigned int cubit_bc_type,
                            EntityHandle &meshset) const;
 
  /**
   * @brief Check if meshset constains entities
   *
   * @param ms_id
   * @param cubit_bc_type
   * @param entities
   * @param num_entities
   * @param operation_type
   * @return true
   * @return false
   */
  bool checkIfMeshsetContainsEntities(
      const int ms_id, const unsigned int cubit_bc_type,
      const EntityHandle *entities, int num_entities,
      const int operation_type = moab::Interface::INTERSECT);
 
  /**
   * \ingroup mofem_meshset_mng
   * \brief get all CUBIT meshsets by CUBIT type
   *
   * \param  see CubitBC (NODESET, SIDESET or BLOCKSET and more).
   * \param meshsets is range of meshsets
   */
  MoFEMErrorCode getMeshsetsByType(const unsigned int cubit_bc_type,
                                   Range &meshsets) const;
 
  /**
   * \brief add blocksets reading config file
 
   Example of config file
   \code
 
   [block_1001]
 
   # Example applying attributes to blockset
 
   id=2001
   add=BLOCKSET
   user1=1.0  # attribute value 1
   user2=2.0  # you can set up to 10 attributes (if needed could be easily
   extended to more, let us know)
   user3=3.0
 
   [block_1002]
 
   # Example applying material block (isotropic elastic material)
 
   id=2002
   add=BLOCKSET
   name=MAT_ELASTIC
   young=10
   poisson=0.25
   thermalexpansion=0
 
   [block_1003]
 
   # Example applying displacement constrains
 
   id=2003
   add=NODESET
 
   # Each flag means that boundary consition on displacements is set.
   disp_flag1=1 # Setting constrains in x- direction
   disp_flag2=1 # Setting constrains in y- direction
   disp_flag3=1 # Setting constrains in z- direction
   disp_flag4=1 # Setting constrains on rotation over x- axis
   disp_flag5=1 # Setting constrains on rotation over y- axis
   disp_flag6=1 # Setting constrains on rotation over z-axis
   disp_ux=1 # value of disp in x- direction
   disp_uy=2
   disp_uz=3
   disp_rx=4 # value of rotation in y-direction
   disp_ry=5
   disp_rz=6
 
   # Note above values could be interpreted differently if needed.
 
   [block_1004]
 
   # Example applying force boundary conditions
 
   id=2004
   add=NODESET
   force_magnitude=1
   moment_magnitude=1
   force_fx=1
   force_fy=1
   force_fz=1
   moment_mx=1
   moment_my=1
   moment_mz=1
 
   [block_1005]
 
   # Example applying pressure boundary conditions
 
   id=2005
   add=SIDESET
   pressure_flag2=1       # 0: Pressure is interpreted as pure pressure 1:
   pressure is interpreted as total force
   pressure_magnitude=1
 
   # Example applying temperature boundary conditions
 
   [block_1006]
 
   id=2006
   add=NODESET
   temperature_flag1=1        # 0: N/A, 1: temperature value applied
   temperature_t=1
 
   [block_1007]
 
   id=2007
   add=SIDESET
   heatflux_flag1=1        # 0: N/A, 1: heat flux applied
   heatflux_magnitude=1
 
   [block_1008]
 
   # Example applying material block (isotropic thermal material)
 
   id=2008
   add=BLOCKSET
   name=MAT_THERMAL # Hast to be set for Thermal Mat
   conductivity=1
   capacity=1
 
   [block_1009]
 
   # Example applying interface
   id=2009
   add=SIDESET
   interface_type=1
 
   [block_1010]
 
   # Example applying material block for interface element
 
   id=2010
   add=BLOCKSET
   name=MAT_INTERF
   interface_alpha = 1
   interface_beta = 0
   interface_ft = 1
   interface_Gf = 1
 
 
   [block_1009]
 
   # Example applying material block (isotropic trans iso material)
 
   id=2011
   add=BLOCKSET
   name=MAT_ELASTIC_TRANS_ISO
   Youngp=1
   Youngz=2
   Poissonp=3
   Poissonpz=4
   Shearzp=5
 
   [SET_ATTR_foo]
 
   # Example set atttributes to block name "foo"
   number_of_attributes=3
   user1=1
   user2=2
   user3=3
 
   \endcode
 
   * @param  file_name config file
   * @return           error code
 
   */
  MoFEMErrorCode setMeshsetFromFile(const string file_name,
                                    const bool clean_file_options = true);
 
  /**
   * \brief get name of config file from line command '-meshsets_config'
   * @return error code
 
   Option is "-meshsets_config file_name.cfg"
 
   */
  MoFEMErrorCode setMeshsetFromFile();
 
  /**
   * @brief save cubit meshset entities on the moab mesh
   *
   * @param ms_id id of the cubit meshset (NODESET  SIDESET  BLOCKSET)
   * @param cubit_bc_type type of a cubit mesheset
   * @param file_name optional name for the file
   * @param file_type optional file type for moab (VTK MOAB)
   * @param options optional parameters for moab writer (PARALLEL=WRITE_PART)
   * @return MoFEMErrorCode
   */
  MoFEMErrorCode
  saveMeshsetToFile(const int ms_id, const unsigned int cubit_bc_type,
                    const std::string file_name = "out_meshset.vtk",
                    const std::string file_type = "VTK",
                    const std::string options = "") const;
 
  /**
   * @brief save cubit meshset entities on the moab mesh
   *
   * @param ms_id id of the cubit meshset
   * @param cubit_bc_type type of a cubit mesheset (NODESET  SIDESET  BLOCKSET)
   * @param dim dimension of the entities
   * @param file_name optional name for the file
   * @param file_type optional file type for moab (VTK MOAB)
   * @param options optional parameters for moab writer (PARALLEL=WRITE_PART)
   * @return MoFEMErrorCode
   */
  MoFEMErrorCode saveMeshsetToFile(
      const int ms_id, const unsigned int cubit_bc_type, const int dim,
      const std::string file_name = "out_meshset.vtk",
      const bool recursive = false, const std::string file_type = "VTK",
      const std::string options = "") const;
 
  /**
   * \brief Get config file options, use with care
   * @return error code
   */
  inline boost::shared_ptr<boost::program_options::options_description> &
  getConfigFileOptionsPtr() {
    return configFileOptionsPtr;
  }
 
  /**
   * @brief Update all blolsets, sidesets and node sets
   * 
   * @param bit 
   * @return MoFEMErrorCode 
   */
  MoFEMErrorCode updateAllMeshsetsByEntitiesChildren(const BitRefLevel &bit);
 
  static bool broadcastMeshsetsOn; ///< if true meshsets are synchronised
                                   ///< between processors
 
  static void sortMeshsets(std::vector<const CubitMeshSets *> &vec_ptr);
 
protected:
  Tag nsTag;
  Tag ssTag;
  Tag nsTag_data;
  Tag ssTag_data;
  Tag bhTag;
  Tag bhTag_header;
 
  // cubit
  CubitMeshSet_multiIndex cubitMeshsets; ///< cubit meshsets
  boost::shared_ptr<boost::program_options::options_description>
      configFileOptionsPtr; ///< config file options
 
};
 
template <class CUBIT_BC_DATA_TYPE>
MoFEMErrorCode MeshsetsManager::printBcSet(CUBIT_BC_DATA_TYPE &data,
                                           unsigned long int type) const {
  MoFEMFunctionBegin;
  const MoFEM::Interface &m_field = cOre;
  const moab::Interface &moab = m_field.get_moab();
  for (auto it : getCubitMeshsetPtr(type)) {
    CHKERR it->getBcDataStructure(data);
    MOFEM_LOG("MeshsetMngWorld", Sev::inform) << *it;
    MOFEM_LOG("MeshsetMngWorld", Sev::inform) << data;
    MOFEM_LOG("MeshsetMngWorld", Sev::inform) << "name " << it->getName();
    for (EntityType t = MBVERTEX; t != MBENTITYSET; ++t) {
      int nb;
      CHKERR moab.get_number_entities_by_type(it->meshset, t, nb, true);
      if (nb > 0) {
        MOFEM_LOG("MeshsetMngSync", Sev::inform)
            << "msId " << it->getMeshsetId() << " number of "
            << moab::CN::EntityTypeName(t) << " " << nb;
      }
    }
    MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::inform);
  }
  MoFEMFunctionReturn(0);
}
 
} // namespace MoFEM
 
#endif //__MESHSETSMANAGER_HPP__
 
/**
 * \defgroup mofem_meshset_mng MeshsetsManager
 * \brief Interface for meshsets with entities with data and boundary conditions
 *
 * \ingroup mofem
 **/

_IT_CUBITMESHSETS_BY_NAME_FOR_LOOP_

#define _IT_CUBITMESHSETS_BY_NAME_FOR_LOOP_	(		MESHSET_MANAGER,
			NAME,
			IT
	)

#include <src/interfaces/MeshsetsManager.hpp>

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
	)

#include <src/interfaces/MeshsetsManager.hpp>

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, continuity_check_on_contact_prism_side_ele.cpp, cubit_bc_test.cpp, elasticity.cpp, forces_and_sources_testing_flat_prism_element.cpp, helmholtz.cpp, mesh_cut.cpp, mesh_insert_interface_atom.cpp, navier_stokes.cpp, photon_diffusion.cpp, and unsaturated_transport.cpp.

Definition at line 71 of file MeshsetsManager.hpp.

_IT_CUBITMESHSETS_FOR_LOOP_

#define _IT_CUBITMESHSETS_FOR_LOOP_	(		MESHSET_MANAGER,
			IT
	)

#include <src/interfaces/MeshsetsManager.hpp>

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, 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 = ""
	)

#include <src/interfaces/MeshsetsManager.hpp>

add cubit meshset

Parameters

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

Examples

EshelbianPlasticity.cpp, add_cubit_meshsets.cpp, and mesh_cut.cpp.

Definition at line 396 of file MeshsetsManager.cpp.

                                                                 {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
  if (checkMeshset(ms_id, cubit_bc_type)) {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "such cubit meshset is already there id  = %d", 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

#include <src/interfaces/MeshsetsManager.hpp>

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 371 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

#include <src/interfaces/MeshsetsManager.hpp>

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 247 of file MeshsetsManager.hpp.

                                                          {
    return cubitMeshsets.begin();
  }

getBegin() [2/3]

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

inline

#include <src/interfaces/MeshsetsManager.hpp>

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 355 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

#include <src/interfaces/MeshsetsManager.hpp>

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 280 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

#include <src/interfaces/MeshsetsManager.hpp>

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 318 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

#include <src/interfaces/MeshsetsManager.hpp>

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 336 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

#include <src/interfaces/MeshsetsManager.hpp>

Get vector of pointer by bc type.

Returns

MoFEMErrorCode

Definition at line 649 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

#include <src/interfaces/MeshsetsManager.hpp>

Get vector of pointer by bc type.

Parameters

std::vector<const

CubitMeshSets *>

Returns

MoFEMErrorCode

Definition at line 635 of file MeshsetsManager.cpp.

                                                     {
  MoFEMFunctionBegin;
  auto r = cubitMeshsets.get<CubitMeshsetType_mi_tag>().equal_range(
      cubit_bc_type.to_ulong());
  for (auto it = r.first; it != r.second; ++it) {
    vec_ptr.push_back(&*it);
  }
  sortMeshsets(vec_ptr);
  MoFEMFunctionReturn(0);
}

getCubitMeshsetPtr() [3/7]

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

#include <src/interfaces/MeshsetsManager.hpp>

get cubit meshset

Definition at line 608 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

#include <src/interfaces/MeshsetsManager.hpp>

get cubit meshset

Examples

EshelbianPlasticity.cpp, mofem/tutorials/vec-7/adjoint.cpp, seepage.cpp, and thermo_elastic.cpp.

Definition at line 589 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 {
    SETERRQ(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

#include <src/interfaces/MeshsetsManager.hpp>

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 673 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

#include <src/interfaces/MeshsetsManager.hpp>

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

Parameters

reg_exp_name
std::vector<const	CubitMeshSets *>

Returns

MoFEMErrorCode

Definition at line 656 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);
    }
  }
  sortMeshsets(vec_ptr);
  MoFEMFunctionReturn(0);
}

getCubitMeshsetPtr() [7/7]

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

#include <src/interfaces/MeshsetsManager.hpp>

get cubit meshset

Definition at line 617 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) {
    SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "meshset name <%s> is not there", name.c_str());
  }
  if (std::distance(miit, hi_miit) > 1) {
    SETERRQ(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

#include <src/interfaces/MeshsetsManager.hpp>

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 261 of file MeshsetsManager.hpp.

                                                 {
    return cubitMeshsets.end();
  }

getEnd() [2/3]

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

inline

#include <src/interfaces/MeshsetsManager.hpp>

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 373 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

#include <src/interfaces/MeshsetsManager.hpp>

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 300 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

#include <src/interfaces/MeshsetsManager.hpp>

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 680 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 {
    SETERRQ(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

#include <src/interfaces/MeshsetsManager.hpp>

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 700 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 {
    SETERRQ(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

#include <src/interfaces/MeshsetsManager.hpp>

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 719 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 {
    SETERRQ(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

#include <src/interfaces/MeshsetsManager.hpp>

get all CUBIT meshsets by CUBIT type

Parameters

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

Definition at line 753 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);
}