MoFEM::MeshsetsManager Struct Reference

Interface for managing meshsets containing materials and boundary conditions. More...

#include <src/interfaces/MeshsetsManager.hpp>

Inheritance diagram for MoFEM::MeshsetsManager:

Collaboration diagram for MoFEM::MeshsetsManager:

Public Member Functions
MoFEMErrorCode	query_interface (const MOFEMuuid &uuid, UnknownInterface **iface) const
	MeshsetsManager (const MoFEM::Core &core)
	~MeshsetsManager ()
MoFEMErrorCode	getTags (int verb=-1)
	get tags handlers used on meshsets More...
Tag	get_nsTag () const
	get tag handle used to store "id" of NODESET More...
Tag	get_ssTag () const
	get tag handle used to store "id" of SIDESET More...
Tag	get_nsTag_data () const
	get tag handle used to store boundary data on NODESET More...
Tag	get_ssTag_data () const
	get tag handle used to store boundary data on SIDESET More...
Tag	get_bhTag () const
	get tag handle used to store "id" of BLOCKSET More...
Tag	get_bhTag_header () const
	get tag handle used to store of block set header (Used by Cubit) More...
MeshsetsManager *	get_meshsets_manager_ptr ()
	return pointer to meshset manager More...
const MeshsetsManager *	get_meshsets_manager_ptr () const
	return pointer to meshset manager More...
MoFEMErrorCode	clearMap ()
	clear multi-index container More...
MoFEMErrorCode	initialiseDatabaseFromMesh (int verb=0)
template<class CUBIT_BC_DATA_TYPE >
MoFEMErrorCode	printBcSet (CUBIT_BC_DATA_TYPE &data, unsigned long int type) const
MoFEMErrorCode	printDisplacementSet () const
	print meshsets with displacement boundary conditions data structure More...
MoFEMErrorCode	printPressureSet () const
	print meshsets with pressure boundary conditions data structure More...
MoFEMErrorCode	printForceSet () const
	print meshsets with force boundary conditions data structure More...
MoFEMErrorCode	printTemperatureSet () const
	print meshsets with temperature boundary conditions data structure More...
MoFEMErrorCode	printHeatFluxSet () const
	print meshsets with heat flux boundary conditions data structure More...
MoFEMErrorCode	printMaterialsSet () const
	print meshsets with material data structure set on it More...
CubitMeshSet_multiIndex &	getMeshsetsMultindex ()
CubitMeshSet_multiIndex::iterator	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	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	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	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	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	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	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	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	checkMeshset (const int ms_id, const CubitBCType cubit_bc_type) const
	check for CUBIT Id and CUBIT type More...
bool	checkMeshset (const string name, int *const number_of_meshsets_ptr=NULL) const
	check if meshset of given name exist More...
MoFEMErrorCode	addMeshset (const CubitBCType cubit_bc_type, const int ms_id, const std::string name="")
	add cubit meshset More...
MoFEMErrorCode	addEntitiesToMeshset (const CubitBCType cubit_bc_type, const int ms_id, const Range &ents)
	add entities to cubit meshset More...
MoFEMErrorCode	addEntitiesToMeshset (const CubitBCType cubit_bc_type, const int ms_id, const EntityHandle *ents, const int nb_ents)
	add entities to cubit meshset More...
MoFEMErrorCode	setAtributes (const CubitBCType cubit_bc_type, const int ms_id, const std::vector< double > &attributes, const std::string name="")
	set attributes to cubit meshset More...
MoFEMErrorCode	setAtributesByDataStructure (const CubitBCType cubit_bc_type, const int ms_id, const GenericAttributeData &data, const std::string name="")
	set (material) data structure to cubit meshset More...
MoFEMErrorCode	setBcData (const CubitBCType cubit_bc_type, const int ms_id, const GenericCubitBcData &data)
	set boundary data structure to meshset More...
MoFEMErrorCode	deleteMeshset (const CubitBCType cubit_bc_type, const int ms_id, const MoFEMTypes bh=MF_EXIST)
	delete cubit meshset More...
MoFEMErrorCode	getCubitMeshsetPtr (const int ms_id, const CubitBCType cubit_bc_type, const CubitMeshSets **cubit_meshset_ptr)
	get cubit meshset More...
MoFEMErrorCode	getCubitMeshsetPtr (const string name, const CubitMeshSets **cubit_meshset_ptr)
	get cubit meshset More...
MoFEMErrorCode	getEntitiesByDimension (const int ms_id, const unsigned int cubit_bc_type, const int dimension, Range &entities, const bool recursive=true)
	get entities from CUBIT/meshset of a particular entity dimensionNodeset 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. More...
MoFEMErrorCode	getEntitiesByDimension (const int ms_id, const unsigned int cubit_bc_type, Range &entities, const bool recursive=true)
	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. More...
MoFEMErrorCode	getMeshset (const int ms_id, const unsigned int cubit_bc_type, EntityHandle &meshset)
	get meshset from CUBIT Id and CUBIT type More...
MoFEMErrorCode	getMeshsetsByType (const unsigned int cubit_bc_type, Range &meshsets)
	get all CUBIT meshsets by CUBIT type More...
MoFEMErrorCode	setMeshsetFromFile (const string file_name, const bool clean_file_options=true)
	add blocksets reading config file More...
MoFEMErrorCode	setMeshsetFromFile ()
	get name of config file from line command '-meshsets_config' More...
boost::shared_ptr< boost::program_options::options_description > &	getConfigFileOptionsPtr ()
	Get config file options, use with care. More...
MoFEMErrorCode	updateAllMeshsetsByEntitiesChildren (const BitRefLevel &bit)
Public Member Functions inherited from MoFEM::UnknownInterface
template<class IFACE >
MoFEMErrorCode	registerInterface (const MOFEMuuid &uuid, bool error_if_registration_failed=true)
	Register interface. More...
template<class IFACE , bool VERIFY = false>
MoFEMErrorCode	getInterface (const MOFEMuuid &uuid, IFACE *&iface) const
	Get interface by uuid and return reference to pointer of interface. More...
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE *&iface) const
	Get interface refernce to pointer of interface. More...
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE **const iface) const
	Get interface pointer to pointer of interface. More...
template<class IFACE , typename boost::enable_if< boost::is_pointer< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get interface pointer to pointer of interface. More...
template<class IFACE , typename boost::enable_if< boost::is_reference< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get reference to interface. More...
template<class IFACE >
IFACE *	getInterface () const
	Function returning pointer to interface. More...
virtual	~UnknownInterface ()
virtual MoFEMErrorCode	getLibVersion (Version &version) const
	Get library version. More...
virtual const MoFEMErrorCode	getFileVersion (moab::Interface &moab, Version &version) const
	Get database major version. More...
virtual MoFEMErrorCode	getInterfaceVersion (Version &version) const
	Get database major version. More...
template<>
MoFEMErrorCode	getInterface (const MOFEMuuid &uuid, UnknownInterface *&iface) const

Public Attributes
MoFEM::Core &	cOre

Protected Attributes
Tag	nsTag
Tag	ssTag
Tag	nsTag_data
Tag	ssTag_data
Tag	bhTag
Tag	bhTag_header
CubitMeshSet_multiIndex	cubitMeshsets
	cubit meshsets More...
boost::shared_ptr< boost::program_options::options_description >	configFileOptionsPtr
	config file options More...

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

Detailed Description

Interface for managing meshsets containing materials and boundary conditions.

Examples:

cell_forces.cpp, cubit_bc_test.cpp, elasticity.cpp, mesh_cut.cpp, mesh_smoothing.cpp, split_sideset.cpp, and unsaturated_transport.cpp.

Definition at line 118 of file MeshsetsManager.hpp.

Constructor & Destructor Documentation

MeshsetsManager()

MoFEM::MeshsetsManager::MeshsetsManager

(

const MoFEM::Core &

core

)

Definition at line 43 of file MeshsetsManager.cpp.

      : cOre(const_cast<Core &>(core)) {}

~MeshsetsManager()

MoFEM::MeshsetsManager::~MeshsetsManager

(

)

Definition at line 125 of file MeshsetsManager.hpp.

{}

Member Function Documentation

addEntitiesToMeshset() [1/2]

MoFEMErrorCode MoFEM::MeshsetsManager::addEntitiesToMeshset	(	const CubitBCType	cubit_bc_type,
		const int	ms_id,
		const Range &	ents
	)

add entities to cubit meshset

Parameters

cubit_bc_type	type of meshset, f.e. NODESET, SIDESET or BLOCKSET
ms_id	id of meshset
ents	entities to add

Returns

error code

Examples:

mesh_cut.cpp.

Definition at line 274 of file MeshsetsManager.cpp.

                                                                            {
    Interface &m_field = cOre;
    moab::Interface &moab = m_field.get_moab();
    MoFEMFunctionBeginHot;
    CubitMeshSet_multiIndex::index<
        Composite_Cubit_msId_And_MeshSetType_mi_tag>::type::iterator cit =
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>().find(
            boost::make_tuple(ms_id, cubit_bc_type.to_ulong()));
    if (cit ==
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>()
            .end()) {
      SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
               "such cubit meshset is already there", ms_id);
    }
    EntityHandle meshset = cit->getMeshset();
    rval = moab.add_entities(meshset, ents);
    CHKERRQ_MOAB(rval);
    MoFEMFunctionReturnHot(0);
  }

addEntitiesToMeshset() [2/2]

MoFEMErrorCode MoFEM::MeshsetsManager::addEntitiesToMeshset	(	const CubitBCType	cubit_bc_type,
		const int	ms_id,
		const EntityHandle *	ents,
		const int	nb_ents
	)

add entities to cubit meshset

Parameters

cubit_bc_type	type of meshset, f.e. NODESET, SIDESET or BLOCKSET
ms_id	id of meshset
ents	pointer to entities array
nb_ents	number of entities in array

Returns

error code

Definition at line 295 of file MeshsetsManager.cpp.

                                                   {
    Interface &m_field = cOre;
    moab::Interface &moab = m_field.get_moab();
    MoFEMFunctionBeginHot;
    CubitMeshSet_multiIndex::index<
        Composite_Cubit_msId_And_MeshSetType_mi_tag>::type::iterator cit =
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>().find(
            boost::make_tuple(ms_id, cubit_bc_type.to_ulong()));
    if (cit ==
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>()
            .end()) {
      SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
               "such cubit meshset is already there", ms_id);
    }
    EntityHandle meshset = cit->getMeshset();
    rval = moab.add_entities(meshset, ents, nb_ents);
    CHKERRQ_MOAB(rval);
    MoFEMFunctionReturnHot(0);
  }

checkMeshset()

bool MoFEM::MeshsetsManager::checkMeshset	(	const string	name,
		int *const	number_of_meshsets_ptr = NULL
	)		const

check if meshset of given name exist

Parameters

name	name of meshset

Returns

error code

clearMap()

MoFEMErrorCode MoFEM::MeshsetsManager::clearMap

(

)

clear multi-index container

Returns

error code

Definition at line 46 of file MeshsetsManager.cpp.

                                           {
    MoFEMFunctionBeginHot;
    cubitMeshsets.clear();
    MoFEMFunctionReturnHot(0);
  }

deleteMeshset()

MoFEMErrorCode MoFEM::MeshsetsManager::deleteMeshset	(	const CubitBCType	cubit_bc_type,
		const int	ms_id,
		const MoFEMTypes	bh = MF_EXIST
	)

delete cubit meshset

Parameters

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

Definition at line 408 of file MeshsetsManager.cpp.

                                                                                     {
    Interface &m_field = cOre;
    moab::Interface &moab = m_field.get_moab();
    MoFEMFunctionBeginHot;
    CubitMeshSet_multiIndex::index<
        Composite_Cubit_msId_And_MeshSetType_mi_tag>::type::iterator 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()) {
      if (bh & MF_EXIST) {
        SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
                 "meshset not found", ms_id);
      } else {
        MoFEMFunctionReturnHot(0);
      }
    }
    EntityHandle meshset = miit->getMeshset();
    cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>().erase(
        miit);
    rval = moab.delete_entities(&meshset, 1);
    CHKERRQ_MOAB(rval);
    MoFEMFunctionReturnHot(0);
  }

get_bhTag()

Tag MoFEM::MeshsetsManager::get_bhTag

(

)

const

get tag handle used to store "id" of BLOCKSET

Definition at line 163 of file MeshsetsManager.hpp.

{ return bhTag; }

get_bhTag_header()

Tag MoFEM::MeshsetsManager::get_bhTag_header

(

)

const

get tag handle used to store of block set header (Used by Cubit)

Definition at line 168 of file MeshsetsManager.hpp.

{ return bhTag_header; }

get_meshsets_manager_ptr() [1/2]

MeshsetsManager* MoFEM::MeshsetsManager::get_meshsets_manager_ptr

(

)

return pointer to meshset manager

Definition at line 173 of file MeshsetsManager.hpp.

{ return this; }

get_meshsets_manager_ptr() [2/2]

const MeshsetsManager* MoFEM::MeshsetsManager::get_meshsets_manager_ptr

(

)

const

return pointer to meshset manager

Definition at line 178 of file MeshsetsManager.hpp.

{ return this; }

get_nsTag()

Tag MoFEM::MeshsetsManager::get_nsTag

(

)

const

get tag handle used to store "id" of NODESET

Definition at line 143 of file MeshsetsManager.hpp.

{ return nsTag; }

get_nsTag_data()

Tag MoFEM::MeshsetsManager::get_nsTag_data

(

)

const

get tag handle used to store boundary data on NODESET

Definition at line 153 of file MeshsetsManager.hpp.

{ return nsTag_data; }

get_ssTag()

Tag MoFEM::MeshsetsManager::get_ssTag

(

)

const

get tag handle used to store "id" of SIDESET

Definition at line 148 of file MeshsetsManager.hpp.

{ return ssTag; }

get_ssTag_data()

Tag MoFEM::MeshsetsManager::get_ssTag_data

(

)

const

get tag handle used to store boundary data on SIDESET

Definition at line 158 of file MeshsetsManager.hpp.

{ return ssTag_data; }

getConfigFileOptionsPtr()

boost::shared_ptr<boost::program_options::options_description>& MoFEM::MeshsetsManager::getConfigFileOptionsPtr

(

)

Get config file options, use with care.

Returns

error code

Definition at line 739 of file MeshsetsManager.hpp.

                            {
    return configFileOptionsPtr;
  }

getMeshsetsMultindex()

CubitMeshSet_multiIndex& MoFEM::MeshsetsManager::getMeshsetsMultindex

(

)

Examples:

split_sideset.cpp.

Definition at line 264 of file MeshsetsManager.hpp.

                                                         {
    return cubitMeshsets;
  }

getTags()

MoFEMErrorCode MoFEM::MeshsetsManager::getTags

(

int

verb = -1

)

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.

Definition at line 80 of file MeshsetsManager.cpp.

                                                  {
    MoFEMFunctionBeginHot;
    Interface &m_field = cOre;
    moab::Interface &moab = m_field.get_moab();
    int default_val = -1;
    rval =
        moab.tag_get_handle(DIRICHLET_SET_TAG_NAME, 1, MB_TYPE_INTEGER, nsTag,
                            MB_TAG_SPARSE | MB_TAG_CREAT, &default_val);
    CHKERRQ_MOAB(rval);
    rval = moab.tag_get_handle(NEUMANN_SET_TAG_NAME, 1, MB_TYPE_INTEGER, ssTag,
                               MB_TAG_SPARSE | MB_TAG_CREAT, &default_val);
    CHKERRQ_MOAB(rval);
    const int def_bc_data_len = 0;
    std::string tag_name = std::string(DIRICHLET_SET_TAG_NAME) + "__BC_DATA";
    rval = moab.tag_get_handle(
        tag_name.c_str(), def_bc_data_len, MB_TYPE_OPAQUE, nsTag_data,
        MB_TAG_CREAT | MB_TAG_SPARSE | MB_TAG_BYTES | MB_TAG_VARLEN, NULL);
    CHKERRQ_MOAB(rval);
    tag_name = std::string(NEUMANN_SET_TAG_NAME) + "__BC_DATA";
    rval = moab.tag_get_handle(
        tag_name.c_str(), def_bc_data_len, MB_TYPE_OPAQUE, ssTag_data,
        MB_TAG_CREAT | MB_TAG_SPARSE | MB_TAG_BYTES | MB_TAG_VARLEN, NULL);
    CHKERRQ_MOAB(rval);
    rval = moab.tag_get_handle(MATERIAL_SET_TAG_NAME, 1, MB_TYPE_INTEGER, bhTag,
                               MB_TAG_SPARSE | MB_TAG_CREAT, &default_val);
    CHKERRQ_MOAB(rval);
    std::vector<unsigned int> def_uint_zero(3, 0);
    rval = moab.tag_get_handle(
        BLOCK_HEADER, 3 * sizeof(unsigned int), MB_TYPE_INTEGER, bhTag_header,
        MB_TAG_CREAT | MB_TAG_SPARSE | MB_TAG_BYTES, &def_uint_zero[0]);
    CHKERRQ_MOAB(rval);
    Tag block_attribs;
    int def_Block_Attributes_length = 0;
    rval = moab.tag_get_handle(
        BLOCK_ATTRIBUTES, def_Block_Attributes_length, MB_TYPE_DOUBLE,
        block_attribs, MB_TAG_CREAT | MB_TAG_SPARSE | MB_TAG_VARLEN, NULL);
    CHKERRQ_MOAB(rval);
    Tag entity_name_tag;
    rval = moab.tag_get_handle(NAME_TAG_NAME, NAME_TAG_SIZE, MB_TYPE_OPAQUE,
                               entity_name_tag, MB_TAG_SPARSE | MB_TAG_CREAT);
    CHKERRQ_MOAB(rval);
    MoFEMFunctionReturnHot(0);
  }

initialiseDatabaseFromMesh()

MoFEMErrorCode MoFEM::MeshsetsManager::initialiseDatabaseFromMesh

(

int

verb = 0

)

initialize container form data on mesh

Returns

[description]

Definition at line 52 of file MeshsetsManager.cpp.

                                                                     {
    Interface &m_field = cOre;
    moab::Interface &moab = m_field.get_moab();
    MoFEMFunctionBegin;
    Range meshsets;
    CHKERR moab.get_entities_by_type(0, MBENTITYSET, meshsets, false);
    for (Range::iterator mit = meshsets.begin(); mit != meshsets.end(); mit++) {
      // check if meshset is cubit meshset
      CubitMeshSets base_meshset(moab, *mit);
      if ((base_meshset.cubitBcType & CubitBCType(NODESET | SIDESET | BLOCKSET))
              .any()) {
        std::pair<CubitMeshSet_multiIndex::iterator, bool> p =
            cubitMeshsets.insert(base_meshset);
        if (!p.second) {
          SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                  "meshset not inserted");
        }
        if (verb > QUIET) {
          std::ostringstream ss;
          ss << "read cubit " << base_meshset << std::endl;
          // PetscSynchronizedPrintf(comm,ss.str().c_str());
          PetscPrintf(m_field.get_comm(), ss.str().c_str());
        }
      }
    }
    MoFEMFunctionReturn(0);
  }

printBcSet()

template<class CUBIT_BC_DATA_TYPE >

MoFEMErrorCode MoFEM::MeshsetsManager::printBcSet	(	CUBIT_BC_DATA_TYPE &	data,
		unsigned long int	type
	)		const

Definition at line 193 of file MeshsetsManager.hpp.

                                                          {

    MoFEMFunctionBeginHot;
    try {
      const MoFEM::Interface &m_field = cOre;
      const moab::Interface &moab = m_field.get_moab();
      for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_((*this), type, it)) {
        ierr = it->getBcDataStructure(data);
        CHKERRG(ierr);
        std::ostringstream ss;
        ss << *it << std::endl;
        ss << data << std::endl;
        Range tets, tris, edges, nodes;
        rval = moab.get_entities_by_type(it->meshset, MBTET, tets, true);
        CHKERRQ_MOAB(rval);
        rval = moab.get_entities_by_type(it->meshset, MBTRI, tris, true);
        CHKERRQ_MOAB(rval);
        rval = moab.get_entities_by_type(it->meshset, MBEDGE, edges, true);
        CHKERRQ_MOAB(rval);
        rval = moab.get_entities_by_type(it->meshset, MBVERTEX, nodes, true);
        CHKERRQ_MOAB(rval);
        ss << "name " << it->getName() << std::endl;
        ss << "msId " << it->getMeshsetId() << " nb. tets " << tets.size()
           << std::endl;
        ss << "msId " << it->getMeshsetId() << " nb. tris " << tris.size()
           << std::endl;
        ss << "msId " << it->getMeshsetId() << " nb. edges " << edges.size()
           << std::endl;
        ss << "msId " << it->getMeshsetId() << " nb. nodes " << nodes.size()
           << std::endl;
        ss << std::endl;
        PetscPrintf(m_field.get_comm(), ss.str().c_str());
      }
    } catch (MoFEMException const &e) {
      SETERRQ(PETSC_COMM_SELF, e.errorCode, e.errorMessage);
    }
    MoFEMFunctionReturnHot(0);
  }

printDisplacementSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printDisplacementSet

(

)

const

print meshsets with displacement boundary conditions data structure

Examples:

cell_forces.cpp, and elasticity.cpp.

Definition at line 124 of file MeshsetsManager.cpp.

                                                             {
    DisplacementCubitBcData mydata;
    MoFEMFunctionBeginHot;
    ierr = printBcSet(mydata, NODESET | mydata.tYpe.to_ulong());
    CHKERRG(ierr);
    MoFEMFunctionReturnHot(0);
  }

printForceSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printForceSet

(

)

const

print meshsets with force boundary conditions data structure

Examples:

cell_forces.cpp, and elasticity.cpp.

Definition at line 140 of file MeshsetsManager.cpp.

                                                      {
    ForceCubitBcData mydata;
    MoFEMFunctionBeginHot;
    ierr = printBcSet(mydata, NODESET | mydata.tYpe.to_ulong());
    CHKERRG(ierr);
    MoFEMFunctionReturnHot(0);
  }

printHeatFluxSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printHeatFluxSet

(

)

const

print meshsets with heat flux boundary conditions data structure

Definition at line 156 of file MeshsetsManager.cpp.

                                                         {
    HeatFluxCubitBcData mydata;
    MoFEMFunctionBeginHot;
    ierr = printBcSet(mydata, SIDESET | mydata.tYpe.to_ulong());
    CHKERRG(ierr);
    MoFEMFunctionReturnHot(0);
  }

printMaterialsSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printMaterialsSet

(

)

const

print meshsets with material data structure set on it

Examples:

cell_forces.cpp, and elasticity.cpp.

Definition at line 164 of file MeshsetsManager.cpp.

                                                          {
    MoFEMFunctionBeginHot;
    const Interface &m_field = cOre;
    const moab::Interface &moab = m_field.get_moab();
    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
             (*this), BLOCKSET | MAT_ELASTICSET, it)) {
      Mat_Elastic data;
      ierr = it->getAttributeDataStructure(data);
      CHKERRG(ierr);
      std::ostringstream ss;
      ss << *it << std::endl;
      ss << data;
      Range tets;
      rval = moab.get_entities_by_type(it->meshset, MBTET, tets, true);
      CHKERRQ_MOAB(rval);
      ss << "MAT_ELATIC msId " << it->getMeshsetId() << " nb. tets "
         << tets.size() << std::endl;
      ss << std::endl;
      PetscPrintf(m_field.get_comm(), ss.str().c_str());
    }

    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
             m_field, BLOCKSET | MAT_THERMALSET, it)) {
      Mat_Thermal data;
      ierr = it->getAttributeDataStructure(data);
      CHKERRG(ierr);
      std::ostringstream ss;
      ss << *it << std::endl;
      ss << data;
      PetscPrintf(m_field.get_comm(), ss.str().c_str());
    }

    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
             m_field, BLOCKSET | MAT_MOISTURESET, it)) {
      Mat_Moisture data;
      ierr = it->getAttributeDataStructure(data);
      CHKERRG(ierr);
      std::ostringstream ss;
      ss << *it << std::endl;
      ss << data;
      PetscPrintf(m_field.get_comm(), ss.str().c_str());
    }
    MoFEMFunctionReturnHot(0);
  }

printPressureSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printPressureSet

(

)

const

print meshsets with pressure boundary conditions data structure

Definition at line 132 of file MeshsetsManager.cpp.

                                                         {
    PressureCubitBcData mydata;
    MoFEMFunctionBeginHot;
    ierr = printBcSet(mydata, SIDESET | mydata.tYpe.to_ulong());
    CHKERRG(ierr);
    MoFEMFunctionReturnHot(0);
  }

printTemperatureSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printTemperatureSet

(

)

const

print meshsets with temperature boundary conditions data structure

Definition at line 148 of file MeshsetsManager.cpp.

                                                            {
    TemperatureCubitBcData mydata;
    MoFEMFunctionBeginHot;
    ierr = printBcSet(mydata, NODESET | mydata.tYpe.to_ulong());
    CHKERRG(ierr);
    MoFEMFunctionReturnHot(0);
  }

query_interface()

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

virtual

Implements MoFEM::UnknownInterface.

Definition at line 31 of file MeshsetsManager.cpp.

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

setAtributes()

MoFEMErrorCode MoFEM::MeshsetsManager::setAtributes	(	const CubitBCType	cubit_bc_type,
		const int	ms_id,
		const std::vector< double > &	attributes,
		const std::string	name = ""
	)

set attributes to cubit meshset

Parameters

cubit_bc_type	type of meshset, see CubitBC, i.e. BLOCKSET, NODESET, SIDESET
ms_id	id of meshset
attributes	attributes

Returns

error code

Definition at line 317 of file MeshsetsManager.cpp.

                                                                 {
    Interface &m_field = cOre;
    moab::Interface &moab = m_field.get_moab();
    MoFEMFunctionBeginHot;
    CubitMeshSet_multiIndex::index<
        Composite_Cubit_msId_And_MeshSetType_mi_tag>::type::iterator cit =
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>().find(
            boost::make_tuple(ms_id, cubit_bc_type.to_ulong()));
    if (cit ==
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>()
            .end()) {
      SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
               "such cubit meshset is already there", ms_id);
    }
    if (name.size() > 0) {
      bool success = cubitMeshsets.modify(
          cubitMeshsets.project<0>(cit), CubitMeshSets_change_name(moab, name));
      if (!success) {
        SETERRQ(m_field.get_comm(), MOFEM_OPERATION_UNSUCCESSFUL,
                "name to cubit meshset can not be set");
      }
    }
    bool success =
        cubitMeshsets.modify(cubitMeshsets.project<0>(cit),
                             CubitMeshSets_change_attributes(moab, attributes));
    if (!success)
      SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
              "modification unsuccessful");
    MoFEMFunctionReturnHot(0);
  }

setAtributesByDataStructure()

MoFEMErrorCode MoFEM::MeshsetsManager::setAtributesByDataStructure	(	const CubitBCType	cubit_bc_type,
		const int	ms_id,
		const GenericAttributeData &	data,
		const std::string	name = ""
	)

set (material) data structure to cubit meshset

Parameters

cubit_bc_type	type of meshset, see CubitBC, i.e. BLOCKSET, NODESET, SIDESET
ms_id	id of meshset
attributes	attributes

Returns

error code

Definition at line 350 of file MeshsetsManager.cpp.

                                                              {
    Interface &m_field = cOre;
    moab::Interface &moab = m_field.get_moab();
    MoFEMFunctionBeginHot;
    CubitMeshSet_multiIndex::index<
        Composite_Cubit_msId_And_MeshSetType_mi_tag>::type::iterator cit =
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>().find(
            boost::make_tuple(ms_id, cubit_bc_type.to_ulong()));
    if (cit ==
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>()
            .end()) {
      SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
               "such cubit meshset is already there", ms_id);
    }
    if (name.size() > 0) {
      bool success = cubitMeshsets.modify(
          cubitMeshsets.project<0>(cit), CubitMeshSets_change_name(moab, name));
      if (!success) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
                "name to cubit meshset can not be set");
      }
    }
    bool success = cubitMeshsets.modify(
        cubitMeshsets.project<0>(cit),
        CubitMeshSets_change_attributes_data_structure(moab, data));
    if (!success)
      SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
              "modification unsuccessful");
    MoFEMFunctionReturnHot(0);
  }

setBcData()

MoFEMErrorCode MoFEM::MeshsetsManager::setBcData	(	const CubitBCType	cubit_bc_type,
		const int	ms_id,
		const GenericCubitBcData &	data
	)

set boundary data structure to meshset

Parameters

cubit_bc_type	type of meshset, see CubitBC, i.e. BLOCKSET, NODESET, SIDESET
ms_id	id of meshset
data	data structure

Returns

error code

Definition at line 383 of file MeshsetsManager.cpp.

                                                                            {
    Interface &m_field = cOre;
    moab::Interface &moab = m_field.get_moab();
    MoFEMFunctionBeginHot;
    CubitMeshSet_multiIndex::index<
        Composite_Cubit_msId_And_MeshSetType_mi_tag>::type::iterator cit =
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>().find(
            boost::make_tuple(ms_id, cubit_bc_type.to_ulong()));
    if (cit ==
        cubitMeshsets.get<Composite_Cubit_msId_And_MeshSetType_mi_tag>()
            .end()) {
      SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
               "such cubit meshset is already there", ms_id);
    }
    bool success = cubitMeshsets.modify(
        cubitMeshsets.project<0>(cit),
        CubitMeshSets_change_bc_data_structure(moab, data));
    if (!success)
      SETERRQ(PETSC_COMM_SELF, MOFEM_OPERATION_UNSUCCESSFUL,
              "modification unsuccessful");
    MoFEMFunctionReturnHot(0);
  }

setMeshsetFromFile() [1/2]

MoFEMErrorCode MoFEM::MeshsetsManager::setMeshsetFromFile	(	const string	file_name,
		const bool	clean_file_options = true
	)

add blocksets reading config file

Example of config file

[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

Parameters

file_name

config file

Returns

error code

Examples:

unsaturated_transport.cpp.

setMeshsetFromFile() [2/2]

MoFEMErrorCode MoFEM::MeshsetsManager::setMeshsetFromFile

(

)

get name of config file from line command '-meshsets_config'

Returns

error code

Option is "-meshsets_config file_name.cfg"

Definition at line 955 of file MeshsetsManager.cpp.

                                                     {
    Interface &m_field = cOre;
    // moab::Interface &moab = m_field.get_moab();
    PetscBool flg_file;
    char meshset_file_name[255];
    MoFEMFunctionBeginHot;
    ierr = PetscOptionsBegin(m_field.get_comm(), "", "Set meshsets form file",
                             "none");
    CHKERRG(ierr);
    ierr = PetscOptionsString("-meshsets_config", "meshsets config  file name",
                              "", "add_cubit_meshsets.in", meshset_file_name,
                              255, &flg_file);
    CHKERRG(ierr);
    if (flg_file == PETSC_TRUE) {
      ifstream f(meshset_file_name);
      if (!f.good()) {
        SETERRQ1(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
                 "File configuring meshsets ( %s ) can not be open\n",
                 meshset_file_name);
      }
      ierr = setMeshsetFromFile(string(meshset_file_name));
      CHKERRG(ierr);
    }
    ierr = PetscOptionsEnd();
    CHKERRG(ierr);
    MoFEMFunctionReturnHot(0);
  }

updateAllMeshsetsByEntitiesChildren()

MoFEMErrorCode MoFEM::MeshsetsManager::updateAllMeshsetsByEntitiesChildren

(

const BitRefLevel &

bit

)

Definition at line 984 of file MeshsetsManager.cpp.

                                                                             {
    MoFEMFunctionBegin;
    BitRefManager *bit_mng = cOre.getInterface<BitRefManager>();
    for (_IT_CUBITMESHSETS_FOR_LOOP_((*this), iit)) {
      EntityHandle meshset = iit->getMeshset();
      ierr = bit_mng->updateMeshsetByEntitiesChildren(meshset, bit, meshset,
                                                      MBVERTEX, true);
      CHKERRG(ierr);
      ierr = bit_mng->updateMeshsetByEntitiesChildren(meshset, bit, meshset,
                                                      MBEDGE, true);
      CHKERRG(ierr);
      ierr = bit_mng->updateMeshsetByEntitiesChildren(meshset, bit, meshset,
                                                      MBTRI, true);
      CHKERRG(ierr);
      ierr = bit_mng->updateMeshsetByEntitiesChildren(meshset, bit, meshset,
                                                      MBTET, true);
      CHKERRG(ierr);
      ierr = bit_mng->updateMeshsetByEntitiesChildren(meshset, bit, meshset,
                                                      MBPRISM, true);
      CHKERRG(ierr);
    }
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

bhTag

Tag MoFEM::MeshsetsManager::bhTag

protected

Definition at line 750 of file MeshsetsManager.hpp.

bhTag_header

Tag MoFEM::MeshsetsManager::bhTag_header

protected

Definition at line 751 of file MeshsetsManager.hpp.

configFileOptionsPtr

boost::shared_ptr<boost::program_options::options_description> MoFEM::MeshsetsManager::configFileOptionsPtr

protected

config file options

Definition at line 756 of file MeshsetsManager.hpp.

cOre

MoFEM::Core& MoFEM::MeshsetsManager::cOre

Definition at line 123 of file MeshsetsManager.hpp.

cubitMeshsets

CubitMeshSet_multiIndex MoFEM::MeshsetsManager::cubitMeshsets

protected

cubit meshsets

Definition at line 754 of file MeshsetsManager.hpp.

nsTag

Tag MoFEM::MeshsetsManager::nsTag

protected

Definition at line 746 of file MeshsetsManager.hpp.

nsTag_data

Tag MoFEM::MeshsetsManager::nsTag_data

protected

Definition at line 748 of file MeshsetsManager.hpp.

ssTag

Tag MoFEM::MeshsetsManager::ssTag

protected

Definition at line 747 of file MeshsetsManager.hpp.

ssTag_data

Tag MoFEM::MeshsetsManager::ssTag_data

protected

Definition at line 749 of file MeshsetsManager.hpp.

---

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

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