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 (boost::typeindex::type_index type_index, UnknownInterface **iface) const

	MeshsetsManager (const MoFEM::Core &core)

virtual	~MeshsetsManager ()=default

MoFEMErrorCode	getTags (int verb=-1)
	get tags handlers used on meshsets

Tag	get_nsTag () const
	get tag handle used to store "id" of NODESET

Tag	get_ssTag () const
	get tag handle used to store "id" of SIDESET

Tag	get_nsTag_data () const
	get tag handle used to store boundary data on NODESET

Tag	get_ssTag_data () const
	get tag handle used to store boundary data on SIDESET

Tag	get_bhTag () const
	get tag handle used to store "id" of BLOCKSET

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

MeshsetsManager *	get_meshsets_manager_ptr ()
	return pointer to meshset manager

const MeshsetsManager *	get_meshsets_manager_ptr () const
	return pointer to meshset manager

MoFEMErrorCode	clearMap ()
	clear multi-index container

MoFEMErrorCode	initialiseDatabaseFromMesh (int verb=DEFAULT_VERBOSITY)

MoFEMErrorCode	readMeshsets (int verb=DEFAULT_VERBOSITY)
	Read meshsets.

MoFEMErrorCode	broadcastMeshsets (int verb=DEFAULT_VERBOSITY)
	Broadcast meshsets.

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

MoFEMErrorCode	printPressureSet () const
	print meshsets with pressure boundary conditions data structure

MoFEMErrorCode	printForceSet () const
	print meshsets with force boundary conditions data structure

MoFEMErrorCode	printTemperatureSet () const
	print meshsets with temperature boundary conditions data structure

MoFEMErrorCode	printHeatFluxSet () const
	print meshsets with heat flux boundary conditions data structure

MoFEMErrorCode	printMaterialsSet () const
	print meshsets with material data structure set on it

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)

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)

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)

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)

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)

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)

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)

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)

bool	checkMeshset (const int ms_id, const CubitBCType cubit_bc_type) const
	check for CUBIT Id and CUBIT type

bool	checkMeshset (const string name, int *const number_of_meshsets_ptr=NULL) const
	check if meshset of given name exist

MoFEMErrorCode	addMeshset (const CubitBCType cubit_bc_type, const int ms_id, const std::string name="")
	add cubit meshset

MoFEMErrorCode	addEntitiesToMeshset (const CubitBCType cubit_bc_type, const int ms_id, const Range &ents)
	add entities to cubit meshset

MoFEMErrorCode	addEntitiesToMeshset (const CubitBCType cubit_bc_type, const int ms_id, const EntityHandle *ents, const int nb_ents)
	add entities to cubit meshset

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

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

MoFEMErrorCode	setBcData (const CubitBCType cubit_bc_type, const int ms_id, const GenericCubitBcData &data)
	set boundary data structure to meshset

MoFEMErrorCode	deleteMeshset (const CubitBCType cubit_bc_type, const int ms_id, const MoFEMTypes bh=MF_EXIST)
	delete cubit meshset

MoFEMErrorCode	getCubitMeshsetPtr (const int ms_id, const CubitBCType cubit_bc_type, const CubitMeshSets **cubit_meshset_ptr) const
	get cubit meshset

const CubitMeshSets *	getCubitMeshsetPtr (const int ms_id, const CubitBCType cubit_bc_type) const
	get cubit meshset

MoFEMErrorCode	getCubitMeshsetPtr (const CubitBCType cubit_bc_type, std::vector< const CubitMeshSets * > &vec_ptr) const
	Get vector of pointer by bc type.

std::vector< const CubitMeshSets * >	getCubitMeshsetPtr (const CubitBCType cubit_bc_type) const
	Get vector of pointer by bc type.

MoFEMErrorCode	getCubitMeshsetPtr (const string name, const CubitMeshSets **cubit_meshset_ptr) const
	get cubit meshset

MoFEMErrorCode	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 * >	getCubitMeshsetPtr (const std::regex reg_exp_name) const
	Get vector of pointer to blocksets with name satisfying regular expression.

MoFEMErrorCode	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	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	getMeshset (const int ms_id, const unsigned int cubit_bc_type, EntityHandle &meshset) const
	get meshset from CUBIT Id and CUBIT type

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)
	Check if meshset constains entities.

MoFEMErrorCode	getMeshsetsByType (const unsigned int cubit_bc_type, Range &meshsets) const
	get all CUBIT meshsets by CUBIT type

MoFEMErrorCode	setMeshsetFromFile (const string file_name, const bool clean_file_options=true)
	add blocksets reading config file

MoFEMErrorCode	setMeshsetFromFile ()
	get name of config file from line command '-meshsets_config'

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
	save cubit meshset entities on the moab mesh

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
	save cubit meshset entities on the moab mesh

boost::shared_ptr< boost::program_options::options_description > &	getConfigFileOptionsPtr ()
	Get config file options, use with care.

MoFEMErrorCode	updateAllMeshsetsByEntitiesChildren (const BitRefLevel &bit)
	Update all blolsets, sidesets and node sets.

Public Member Functions inherited from MoFEM::UnknownInterface
template<class IFACE >
MoFEMErrorCode	registerInterface (bool error_if_registration_failed=true)
	Register interface.

template<class IFACE >
MoFEMErrorCode	getInterface (IFACE *&iface) const
	Get interface reference to pointer of interface.

template<class IFACE >
MoFEMErrorCode	getInterface (IFACE **const iface) const
	Get interface pointer to pointer of interface.

template<class IFACE , typename boost::enable_if< boost::is_pointer< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get interface pointer to pointer of interface.

template<class IFACE , typename boost::enable_if< boost::is_reference< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get reference to interface.

template<class IFACE >
IFACE *	getInterface () const
	Function returning pointer to interface.

virtual	~UnknownInterface ()=default

Static Public Member Functions
static MoFEMErrorCode	getTags (moab::Interface &moab, Tag &nsTag, Tag &ssTag, Tag &nsTag_data, Tag &ssTag_data, Tag &bhTag, Tag &bhTag_header, int verb=-1)

static void	sortMeshsets (std::vector< const CubitMeshSets * > &vec_ptr)

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

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

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

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

Public Attributes
MoFEM::Core &	cOre

Static Public Attributes
static bool	broadcastMeshsetsOn = true

Protected Attributes
Tag	nsTag

Tag	ssTag

Tag	nsTag_data

Tag	ssTag_data

Tag	bhTag

Tag	bhTag_header

CubitMeshSet_multiIndex	cubitMeshsets
	cubit meshsets

boost::shared_ptr< boost::program_options::options_description >	configFileOptionsPtr
	config file options

Detailed Description

Interface for managing meshsets containing materials and boundary conditions.

Definition at line 104 of file MeshsetsManager.hpp.

Constructor & Destructor Documentation

◆ MeshsetsManager()

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

Definition at line 34 of file MeshsetsManager.cpp.

    : cOre(const_cast<Core &>(core)) {
 
  if (!LogManager::checkIfChannelExist("MeshsetMngWorld")) {
    auto core_log = logging::core::get();
 
    core_log->add_sink(
        LogManager::createSink(LogManager::getStrmWorld(), "MeshsetMngWorld"));
    core_log->add_sink(
        LogManager::createSink(LogManager::getStrmSync(), "MeshsetMngSync"));
    core_log->add_sink(
        LogManager::createSink(LogManager::getStrmSelf(), "MeshsetMngSelf"));
 
    LogManager::setLog("MeshsetMngWorld");
    LogManager::setLog("MeshsetMngSync");
    LogManager::setLog("MeshsetMngSelf");
 
    MOFEM_LOG_TAG("MeshsetMngWorld", "MeshsetMng");
    MOFEM_LOG_TAG("MeshsetMngSync", "MeshsetMng");
    MOFEM_LOG_TAG("MeshsetMngSelf", "MeshsetMng");
  }
 
  MOFEM_LOG("MeshsetMngWorld", Sev::noisy) << "Mashset manager created";
}

◆ ~MeshsetsManager()

virtual MoFEM::MeshsetsManager::~MeshsetsManager ( )

virtualdefault

Member Function Documentation

◆ addEntitiesToMeshset() [1/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 449 of file MeshsetsManager.cpp.

                                                         {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
  auto 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()) {
    SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "Cannot find Cubit meshset with id: %d", ms_id);
  }
  EntityHandle meshset = cit->getMeshset();
  CHKERR moab.add_entities(meshset, ents, nb_ents);
 
  MoFEMFunctionReturn(0);
}

◆ addEntitiesToMeshset() [2/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 429 of file MeshsetsManager.cpp.

                                                                          {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
  auto 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()) {
    SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "Cannot find Cubit meshset with id: %d", ms_id);
  }
  EntityHandle meshset = cit->getMeshset();
  CHKERR moab.add_entities(meshset, ents);
 
  MoFEMFunctionReturn(0);
}

◆ broadcastMeshsets()

MoFEMErrorCode MoFEM::MeshsetsManager::broadcastMeshsets ( int verb = DEFAULT_VERBOSITY )

Broadcast meshsets.

Parameters

verb

Returns: MoFEMErrorCode

Definition at line 129 of file MeshsetsManager.cpp.

                                                          {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
 
  ParallelComm *pcomm = ParallelComm::get_pcomm(&moab, MYPCOMM_INDEX);
  if (pcomm == NULL)
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
            "MOAB communicator not set");
 
  const double coords[] = {0, 0, 0};
 
  auto set_tags_dummy_node = [&](const EntityHandle dummy_node,
                                 const EntityHandle meshset) {
    MoFEMFunctionBegin;
    std::vector<Tag> tag_handles;
    CHKERR moab.tag_get_tags_on_entity(meshset, tag_handles);
    for (auto th : tag_handles) {
      void *data[1];
      int tag_size;
      CHKERR moab.tag_get_by_ptr(th, &meshset, 1, (const void **)data,
                                 &tag_size);
      if (tag_size > 0)
        CHKERR moab.tag_set_by_ptr(th, &dummy_node, 1, data, &tag_size);
    }
    MoFEMFunctionReturn(0);
  };
 
  for (auto t : {NODESET, SIDESET, BLOCKSET}) {
    std::vector<const CubitMeshSets *> vec_ptr;
    for (auto &m : cubitMeshsets) {
      if ((m.getBcType() & CubitBCType(t)).any()) {
        vec_ptr.push_back(&m);
      }
    }
    sortMeshsets(vec_ptr);
 
    std::vector<EntityHandle> vec_meshsets;
    vec_meshsets.reserve(vec_ptr.size());
    for (auto m_ptr : vec_ptr) {
      vec_meshsets.push_back(m_ptr->getMeshset());
    }
 
    switch (t) {
    case NODESET:
      MOFEM_LOG("MeshsetMngSync", Sev::verbose)
          << "broadcast NODESET " << vec_meshsets.size();
      break;
    case SIDESET:
      MOFEM_LOG("MeshsetMngSync", Sev::verbose)
          << "broadcast SIDESET " << vec_meshsets.size();
      break;
    case BLOCKSET:
      MOFEM_LOG("MeshsetMngSync", Sev::verbose)
          << "broadcast BLOCKSET " << vec_meshsets.size();
      break;
    default:
      MOFEM_LOG("MeshsetMngSync", Sev::verbose)
          << "broadcast Unknown " << vec_meshsets.size();
    }
 
    for (int from_proc = 0; from_proc < pcomm->size(); ++from_proc) {
 
      Range r_dummy_nodes;
 
      if (from_proc == pcomm->rank()) {
        std::vector<EntityHandle> dummy_nodes(vec_meshsets.size(), 0);
        int i = 0;
        for (auto m : vec_meshsets) {
          CHKERR moab.create_vertex(coords, dummy_nodes[i]);
          CHKERR set_tags_dummy_node(dummy_nodes[i], m);
          ++i;
        }
        r_dummy_nodes.insert_list(dummy_nodes.begin(), dummy_nodes.end());
      }
 
      CHKERR pcomm->broadcast_entities(from_proc, r_dummy_nodes, false, true);
 
      if (from_proc != pcomm->rank()) {
 
        for (auto dummy_node : r_dummy_nodes) {
          CubitMeshSets broadcast_block(moab, dummy_node);
          auto id = broadcast_block.getMeshsetId();
          if (id != -1) {
            auto mit =
                cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>()
                    .find(boost::make_tuple(
                        id, broadcast_block.getMaskedBcTypeULong()));
            if (mit ==
                cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>()
                    .end()) {
              EntityHandle m;
              CHKERR moab.create_meshset(MESHSET_SET, m);
              CHKERR set_tags_dummy_node(m, dummy_node);
              auto hint = cubitMeshsets.end();
              /*auto p = */ cubitMeshsets.emplace_hint(hint, moab, m);
            }
          } else {
            MOFEM_LOG("MeshsetMngSync", Sev::warning)
                << "broadcasted vertex " << dummy_node << " has negative id";
          }
          
        }
 
      } else {
        MOFEM_LOG("MeshsetMngSync", Sev::verbose)
            << "broadcast send from " << from_proc;
      }
 
      CHKERR moab.delete_entities(r_dummy_nodes);
    }
 
    MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::verbose);
  }
 
  MoFEMFunctionReturn(0);
}

◆ checkIfMeshsetContainsEntities()

bool MoFEM::MeshsetsManager::checkIfMeshsetContainsEntities	(	const int	ms_id,
		const unsigned int	cubit_bc_type,
		const EntityHandle *	entities,
		int	num_entities,
		const int	operation_type = moab::Interface::INTERSECT )

Check if meshset constains entities.

Parameters

ms_id
cubit_bc_type
entities
num_entities
operation_type

Returns: true; false

Definition at line 737 of file MeshsetsManager.cpp.

                                                                              {
  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()) {
    Interface &m_field = cOre;
    return m_field.get_moab().contains_entities(miit->meshset, entities,
                                                num_entities, operation_type);
  } else
    return false;
}

◆ 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

Definition at line 383 of file MeshsetsManager.cpp.

                                                                            {
  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) {
    return false;
  }
  if (number_of_meshsets_ptr) {
    *number_of_meshsets_ptr = std::distance(miit, hi_miit);
  }
  return true;
}

◆ clearMap()

MoFEMErrorCode MoFEM::MeshsetsManager::clearMap ( )

clear multi-index container

Returns: error code

Definition at line 59 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 564 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.to_ulong()));
  if (miit ==
      cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().end()) {
    if (bh & MF_EXIST) {
      SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
               "meshset not found id = %d", ms_id);
    } else {
      MoFEMFunctionReturnHot(0);
    }
  }
  EntityHandle meshset = miit->getMeshset();
  cubitMeshsets.get<Composite_Cubit_msId_And_MeshsetType_mi_tag>().erase(miit);
  CHKERR moab.delete_entities(&meshset, 1);
 
  MoFEMFunctionReturn(0);
}

◆ get_bhTag()

Tag MoFEM::MeshsetsManager::get_bhTag ( ) const

inline

get tag handle used to store "id" of BLOCKSET

Definition at line 153 of file MeshsetsManager.hpp.

153{ return bhTag; }

MoFEM::MeshsetsManager::bhTag

Tag bhTag

Definition MeshsetsManager.hpp:895

◆ get_bhTag_header()

Tag MoFEM::MeshsetsManager::get_bhTag_header ( ) const

inline

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

Definition at line 158 of file MeshsetsManager.hpp.

158{ return bhTag_header; }

MoFEM::MeshsetsManager::bhTag_header

Tag bhTag_header

Definition MeshsetsManager.hpp:896

◆ get_meshsets_manager_ptr() [1/2]

MeshsetsManager * MoFEM::MeshsetsManager::get_meshsets_manager_ptr ( )

inline

return pointer to meshset manager

Definition at line 163 of file MeshsetsManager.hpp.

163{ return this; }

◆ get_meshsets_manager_ptr() [2/2]

const MeshsetsManager * MoFEM::MeshsetsManager::get_meshsets_manager_ptr ( ) const

inline

return pointer to meshset manager

Definition at line 168 of file MeshsetsManager.hpp.

168{ return this; }

◆ get_nsTag()

Tag MoFEM::MeshsetsManager::get_nsTag ( ) const

inline

get tag handle used to store "id" of NODESET

Definition at line 133 of file MeshsetsManager.hpp.

133{ return nsTag; }

MoFEM::MeshsetsManager::nsTag

Tag nsTag

Definition MeshsetsManager.hpp:891

◆ get_nsTag_data()

Tag MoFEM::MeshsetsManager::get_nsTag_data ( ) const

inline

get tag handle used to store boundary data on NODESET

Definition at line 143 of file MeshsetsManager.hpp.

143{ return nsTag_data; }

MoFEM::MeshsetsManager::nsTag_data

Tag nsTag_data

Definition MeshsetsManager.hpp:893

◆ get_ssTag()

Tag MoFEM::MeshsetsManager::get_ssTag ( ) const

inline

get tag handle used to store "id" of SIDESET

Definition at line 138 of file MeshsetsManager.hpp.

138{ return ssTag; }

MoFEM::MeshsetsManager::ssTag

Tag ssTag

Definition MeshsetsManager.hpp:892

◆ get_ssTag_data()

Tag MoFEM::MeshsetsManager::get_ssTag_data ( ) const

inline

get tag handle used to store boundary data on SIDESET

Definition at line 148 of file MeshsetsManager.hpp.

148{ return ssTag_data; }

MoFEM::MeshsetsManager::ssTag_data

Tag ssTag_data

Definition MeshsetsManager.hpp:894

◆ getConfigFileOptionsPtr()

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

inline

Get config file options, use with care.

Returns: error code

Definition at line 873 of file MeshsetsManager.hpp.

                            {
    return configFileOptionsPtr;
  }

◆ getMeshsetsMultindex()

CubitMeshSet_multiIndex & MoFEM::MeshsetsManager::getMeshsetsMultindex ( )

inline

Examples: meshset_to_vtk.cpp.

Definition at line 233 of file MeshsetsManager.hpp.

                                                         {
    return cubitMeshsets;
  }

◆ getTags() [1/2]

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

                                                {
  MoFEMFunctionBegin;
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  CHKERR getTags(moab, nsTag, ssTag, nsTag_data, ssTag_data, bhTag,
                 bhTag_header, verb);
  MoFEMFunctionReturn(0);
}

◆ getTags() [2/2]

MoFEMErrorCode MoFEM::MeshsetsManager::getTags	(	moab::Interface &	moab,
		Tag &	nsTag,
		Tag &	ssTag,
		Tag &	nsTag_data,
		Tag &	ssTag_data,
		Tag &	bhTag,
		Tag &	bhTag_header,
		int	verb = -1 )

static

Definition at line 256 of file MeshsetsManager.cpp.

                                                                     {
  MoFEMFunctionBegin;
 
  int default_val = -1;
  CHKERR moab.tag_get_handle(DIRICHLET_SET_TAG_NAME, 1, MB_TYPE_INTEGER, nsTag,
                             MB_TAG_SPARSE | MB_TAG_CREAT, &default_val);
 
  CHKERR moab.tag_get_handle(NEUMANN_SET_TAG_NAME, 1, MB_TYPE_INTEGER, ssTag,
                             MB_TAG_SPARSE | MB_TAG_CREAT, &default_val);
 
  const int def_bc_data_len = 0;
  std::string tag_name = std::string(DIRICHLET_SET_TAG_NAME) + "__BC_DATA";
  CHKERR 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);
 
  tag_name = std::string(NEUMANN_SET_TAG_NAME) + "__BC_DATA";
  CHKERR 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);
 
  CHKERR moab.tag_get_handle(MATERIAL_SET_TAG_NAME, 1, MB_TYPE_INTEGER, bhTag,
                             MB_TAG_SPARSE | MB_TAG_CREAT, &default_val);
 
  std::vector<unsigned int> def_uint_zero(3, 0);
  CHKERR 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]);
 
  Tag block_attribs;
  int def_Block_Attributes_length = 0;
  CHKERR 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);
 
  Tag entity_name_tag;
  CHKERR moab.tag_get_handle(NAME_TAG_NAME, NAME_TAG_SIZE, MB_TYPE_OPAQUE,
                             entity_name_tag, MB_TAG_SPARSE | MB_TAG_CREAT);
 
  MoFEMFunctionReturn(0);
}

◆ initialiseDatabaseFromMesh()

MoFEMErrorCode MoFEM::MeshsetsManager::initialiseDatabaseFromMesh ( int verb = DEFAULT_VERBOSITY )

\brier initialize container form data on mesh

Returns: [description]

Definition at line 65 of file MeshsetsManager.cpp.

                                                                   {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR readMeshsets(verb);
  if (broadcastMeshsetsOn)
    CHKERR broadcastMeshsets(verb);
 
  std::vector<const CubitMeshSets *> vec_ptr;
  for (auto &m : cubitMeshsets) {
    vec_ptr.push_back(&m);
  }
  sortMeshsets(vec_ptr);
 
  for (auto m_ptr : vec_ptr) {
    MOFEM_LOG("MeshsetMngWorld", Sev::inform)
        << "meshset in database " << *m_ptr;
  }
  // Verbose synchronised print
  for (auto m_ptr : vec_ptr) {
    MOFEM_LOG("MeshsetMngSync", Sev::verbose)
        << "meshset in database " << *m_ptr;
  }
  MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::verbose);
 
  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 906 of file MeshsetsManager.hpp.

                                                                         {
  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);
}

◆ printDisplacementSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printDisplacementSet ( ) const

print meshsets with displacement boundary conditions data structure

Examples: elasticity.cpp, and elasticity_mixed_formulation.cpp.

Definition at line 301 of file MeshsetsManager.cpp.

                                                           {
  DisplacementCubitBcData mydata;
  MoFEMFunctionBegin;
  CHKERR printBcSet(mydata, NODESET | mydata.tYpe.to_ulong());
  MoFEMFunctionReturn(0);
}

◆ printForceSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printForceSet ( ) const

print meshsets with force boundary conditions data structure

Examples: elasticity.cpp, and elasticity_mixed_formulation.cpp.

Definition at line 315 of file MeshsetsManager.cpp.

                                                    {
  ForceCubitBcData mydata;
  MoFEMFunctionBegin;
  CHKERR printBcSet(mydata, NODESET | mydata.tYpe.to_ulong());
  MoFEMFunctionReturn(0);
}

◆ printHeatFluxSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printHeatFluxSet ( ) const

print meshsets with heat flux boundary conditions data structure

Definition at line 329 of file MeshsetsManager.cpp.

                                                       {
  HeatFluxCubitBcData mydata;
  MoFEMFunctionBegin;
  CHKERR printBcSet(mydata, SIDESET | mydata.tYpe.to_ulong());
  MoFEMFunctionReturn(0);
}

◆ printMaterialsSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printMaterialsSet ( ) const

print meshsets with material data structure set on it

Examples: elasticity.cpp, and elasticity_mixed_formulation.cpp.

Definition at line 336 of file MeshsetsManager.cpp.

                                                        {
  MoFEMFunctionBegin;
  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;
    CHKERR it->getAttributeDataStructure(data);
    MOFEM_LOG("MeshsetMngWorld", Sev::inform) << *it;
    MOFEM_LOG("MeshsetMngWorld", Sev::inform) << data;
    Range tets;
    CHKERR moab.get_entities_by_dimension(it->meshset, 3, tets, true);
    MOFEM_LOG("MeshsetMngWorld", Sev::inform)
        << "MAT_ELATIC msId " << it->getMeshsetId() << " nb. volumes "
        << tets.size();
  }
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
           m_field, BLOCKSET | MAT_THERMALSET, it)) {
    Mat_Thermal data;
    CHKERR it->getAttributeDataStructure(data);
    MOFEM_LOG("MeshsetMngWorld", Sev::inform) << *it;
    MOFEM_LOG("MeshsetMngWorld", Sev::inform) << data;
  }
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
           m_field, BLOCKSET | MAT_MOISTURESET, it)) {
    Mat_Moisture data;
    CHKERR it->getAttributeDataStructure(data);
    MOFEM_LOG("MeshsetMngWorld", Sev::inform) << *it;
    MOFEM_LOG("MeshsetMngWorld", Sev::inform) << data;
  }
  MoFEMFunctionReturn(0);
}

◆ printPressureSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printPressureSet ( ) const

print meshsets with pressure boundary conditions data structure

Definition at line 308 of file MeshsetsManager.cpp.

                                                       {
  PressureCubitBcData mydata;
  MoFEMFunctionBegin;
  CHKERR printBcSet(mydata, SIDESET | mydata.tYpe.to_ulong());
  MoFEMFunctionReturn(0);
}

◆ printTemperatureSet()

MoFEMErrorCode MoFEM::MeshsetsManager::printTemperatureSet ( ) const

print meshsets with temperature boundary conditions data structure

Definition at line 322 of file MeshsetsManager.cpp.

                                                          {
  TemperatureCubitBcData mydata;
  MoFEMFunctionBegin;
  CHKERR printBcSet(mydata, NODESET | mydata.tYpe.to_ulong());
  MoFEMFunctionReturn(0);
}

◆ query_interface()

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

virtual

Implements MoFEM::UnknownInterface.

Definition at line 28 of file MeshsetsManager.cpp.

                                                                 {
  *iface = const_cast<MeshsetsManager *>(this);
  return 0;
}

◆ readMeshsets()

MoFEMErrorCode MoFEM::MeshsetsManager::readMeshsets ( int verb = DEFAULT_VERBOSITY )

Read meshsets.

Parameters

verb

Returns: MoFEMErrorCode

Definition at line 92 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 (auto m : meshsets) {
    // check if meshset is cubit meshset
    CubitMeshSets block(moab, m);
    if ((block.cubitBcType & CubitBCType(NODESET | SIDESET | BLOCKSET)).any()) {
      auto p = cubitMeshsets.insert(block);
      if (!p.second) {
        // blockset/nodeset/sideset set exist, could be created on other
        // processor.
        Range ents;
        CHKERR m_field.get_moab().get_entities_by_handle(m, ents, true);
        CHKERR m_field.get_moab().add_entities(p.first->getMeshset(), ents);
        CHKERR m_field.get_moab().delete_entities(&m, 1);
      }
    }
  }
 
  std::vector<const CubitMeshSets *> vec_ptr;
  for (auto &m : cubitMeshsets) {
    vec_ptr.push_back(&m);
  }
  sortMeshsets(vec_ptr);
 
  for (auto m_ptr : vec_ptr) {
    MOFEM_LOG("MeshsetMngSync", Sev::verbose) << "read " << *m_ptr;
  }
  MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::verbose);
 
  MoFEMFunctionReturn(0);
}

◆ saveMeshsetToFile() [1/2]

MoFEMErrorCode MoFEM::MeshsetsManager::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

save cubit meshset entities on the moab mesh

Parameters

ms_id	id of the cubit meshset
cubit_bc_type	type of a cubit mesheset (NODESET SIDESET BLOCKSET)
dim	dimension of the entities
file_name	optional name for the file
file_type	optional file type for moab (VTK MOAB)
options	optional parameters for moab writer (PARALLEL=WRITE_PART)

Returns: MoFEMErrorCode

Definition at line 1430 of file MeshsetsManager.cpp.

                                                              {
 
  MoFEMFunctionBegin;
  MoFEM::Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  Range entities;
  CHKERR getEntitiesByDimension(ms_id, cubit_bc_type, dim, entities, recursive);
  EntityHandle meshset;
  CHKERR moab.create_meshset(MESHSET_SET, meshset);
  CHKERR moab.add_entities(meshset, entities);
  CHKERR moab.write_file(file_name.c_str(), file_type.c_str(), options.c_str(),
                         &meshset, 1);
  CHKERR moab.delete_entities(&meshset, 1);
  MoFEMFunctionReturn(0);
}

◆ saveMeshsetToFile() [2/2]

MoFEMErrorCode MoFEM::MeshsetsManager::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

save cubit meshset entities on the moab mesh

Parameters

ms_id	id of the cubit meshset (NODESET SIDESET BLOCKSET)
cubit_bc_type	type of a cubit mesheset
file_name	optional name for the file
file_type	optional file type for moab (VTK MOAB)
options	optional parameters for moab writer (PARALLEL=WRITE_PART)

Returns: MoFEMErrorCode

Examples: meshset_to_vtk.cpp.

Definition at line 1415 of file MeshsetsManager.cpp.

                                   {
 
  MoFEMFunctionBegin;
  MoFEM::Interface &m_field = cOre;
  const CubitMeshSets *cubit_meshset_ptr;
  CHKERR getCubitMeshsetPtr(ms_id, cubit_bc_type, &cubit_meshset_ptr);
  EntityHandle meshset = cubit_meshset_ptr->getMeshset();
  CHKERR m_field.get_moab().write_file(file_name.c_str(), file_type.c_str(),
                                       options.c_str(), &meshset, 1);
  MoFEMFunctionReturn(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
name	set name to blockset

Returns: error code

Examples: add_cubit_meshsets.cpp.

Definition at line 470 of file MeshsetsManager.cpp.

                                                    {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
  auto 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()) {
    SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "Cannot find Cubit meshset with id: %d", 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");
 
  std::ostringstream ss;
  ss << "Block " << cit->getName();
  ss << " Add attr: ";
  for (unsigned int ii = 0; ii != attributes.size(); ii++) {
    ss << attributes[ii] << " ";
  }
  MOFEM_LOG("MeshsetMngSelf", Sev::noisy) << ss.str();
 
  MoFEMFunctionReturn(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

Examples: add_cubit_meshsets.cpp.

Definition at line 510 of file MeshsetsManager.cpp.

                                                            {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
  auto 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()) {
    SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
             "Cannot find Cubit meshset with id: %d", 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");
  MoFEMFunctionReturn(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

Examples: add_cubit_meshsets.cpp.

Definition at line 541 of file MeshsetsManager.cpp.

                                                                          {
  Interface &m_field = cOre;
  moab::Interface &moab = m_field.get_moab();
  MoFEMFunctionBegin;
  auto 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()) {
    SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
            "Cubit meshset with id is already there id = %d", 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");
  MoFEMFunctionReturn(0);
}

◆ setMeshsetFromFile() [1/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 1386 of file MeshsetsManager.cpp.

                                                   {
  Interface &m_field = cOre;
  PetscBool flg_file;
  char meshset_file_name[255] = "config_file.cfg";
  MoFEMFunctionBegin;
  PetscOptionsBegin(m_field.get_comm(), "", "Set meshsets form file",
                           "none");
  CHKERR PetscOptionsString("-meshsets_config", "meshsets config  file name",
                            "", "add_cubit_meshsets.in", meshset_file_name, 255,
                            &flg_file);
  if (flg_file == PETSC_TRUE) {
    ifstream f(meshset_file_name);
    if (!f.good()) {
      SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
               "File configuring meshsets ( %s ) can not be open\n",
               meshset_file_name);
    }
    CHKERR setMeshsetFromFile(string(meshset_file_name));
    MOFEM_LOG("MeshsetMngWorld", Sev::inform)
        << "Printing all meshsets (old and added from meshsets "
           "configurational file)";
    for (auto cit = getBegin(); cit != getEnd(); cit++) {
      MOFEM_LOG("MeshsetMngWorld", Sev::inform) << *cit;
    }
  }
  PetscOptionsEnd();
  MoFEMFunctionReturn(0);
}

◆ setMeshsetFromFile() [2/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
 
 
[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

Parameters

file_name config file

Returns: error code

Examples: add_cubit_meshsets.cpp, meshset_to_vtk.cpp, and unsaturated_transport.cpp.

Definition at line 802 of file MeshsetsManager.cpp.

                                                                            {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  std::ifstream ini_file(file_name.c_str(), std::ifstream::in);
  po::variables_map vm;
  if (clean_file_options) {
    configFileOptionsPtr =
        boost::shared_ptr<boost::program_options::options_description>(
            new po::options_description());
  }
 
  auto add_block_attributes = [&](auto prefix, auto &block_lists, auto &it) {
    // remove spacec from blockset name
    configFileOptionsPtr->add_options()(
        (prefix + ".number_of_attributes").c_str(),
        po::value<int>(&block_lists[it->getMeshsetId()].numberOfAttributes)
            ->default_value(-1),
        "Number of blockset attribute");
    for (int ii = 1; ii <= 10; ii++) {
      std::string surfix = ".user" + boost::lexical_cast<std::string>(ii);
      configFileOptionsPtr->add_options()(
          (prefix + surfix).c_str(),
          po::value<double>(&block_lists[it->getMeshsetId()].aTtr[ii - 1])
              ->default_value(0.0),
          "Add block attribute");
    }
  };
 
  // Add blocks
  map<int, BlockData> block_lists;
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
    block_lists[it->getMeshsetId()].cubitMeshset = it->getMeshset();
    std::string prefix =
        "block_" + boost::lexical_cast<std::string>(it->getMeshsetId());
    configFileOptionsPtr->add_options()(
        (prefix + ".add").c_str(),
        po::value<string>(&block_lists[it->getMeshsetId()].addType)
            ->default_value("UNKNOWNSET"),
        "Add block set")(
        (prefix + ".id").c_str(),
        po::value<int>(&block_lists[it->getMeshsetId()].iD)->default_value(-1),
        "Id of meshset")(
        (prefix + ".name").c_str(),
        po::value<string>(&block_lists[it->getMeshsetId()].nAme)
            ->default_value(""),
        "Name of the meshset");
 
    // Block attributes
    add_block_attributes(prefix, block_lists, it);
 
    // Mat elastic
    {
      // double Young;          ///< Young's modulus
      // double Poisson;        ///< Poisson's ratio
      // double ThermalExpansion;   ///< Thermal expansion
      configFileOptionsPtr->add_options()(
          (prefix + ".young").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matElastic.data.Young)
              ->default_value(-1),
          "Young modulus")(
          (prefix + ".poisson").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matElastic.data.Poisson)
              ->default_value(-2),
          "Poisson ratio")(
          (prefix + ".thermalexpansion").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matElastic.data.ThermalExpansion)
              ->default_value(-1),
          "Thermal expansion");
      // TODO Add users parameters
    }
    // Mat Trans Iso
    {
      // Young's modulus in xy plane (Ep)
      // Young's modulus in z-direction (Ez)
      // Poisson's ratio in xy plane (vp)
      // Poisson's ratio in z-direction (vpz)
      // Shear modulus in z-direction (Gzp)
      configFileOptionsPtr->add_options()(
          (prefix + ".Youngp").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matTransIso.data.Youngp)
              ->default_value(-1),
          "Youngp")(
          (prefix + ".Youngz").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matTransIso.data.Youngz)
              ->default_value(-1),
          "Youngz")(
          (prefix + ".Poissonp").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matTransIso.data.Poissonp)
              ->default_value(0),
          "Poissonp")(
          (prefix + ".Poissonpz").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matTransIso.data.Poissonpz)
              ->default_value(0),
          "Poissonpz")(
          (prefix + ".Shearzp").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matTransIso.data.Shearzp)
              ->default_value(-1),
          "Shearzp");
      // TODO Add users parameters
    }
    // Mat thermal
    {
      // double Conductivity; ///< Thermal conductivity
      // double HeatCapacity; ///< Heat Capacity
      configFileOptionsPtr->add_options()(
          (prefix + ".conductivity").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matThermal.data.Conductivity)
              ->default_value(-1),
          "Conductivity")(
          (prefix + ".capacity").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matThermal.data.HeatCapacity)
              ->default_value(-1),
          "Capacity");
      // TODO Add users parameters
    }
    // Mat interface
    {
      // double alpha; ///< Elastic modulus multiplier
      // double beta;  ///< Damage Coupling multiplier between normal and
      // shear (g=sqrt(gn^2 + beta(gt1^2 + gt2^2))) double ft;    ///< Maximum
      // stress of crack double Gf;    ///< Fracture Energy
      configFileOptionsPtr->add_options()(
          (prefix + ".interface_alpha").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].matInterf.data.alpha)
              ->default_value(-1),
          "alpha")((prefix + ".interface_beta").c_str(),
                   po::value<double>(
                       &block_lists[it->getMeshsetId()].matInterf.data.beta)
                       ->default_value(-1),
                   "beta")(
          (prefix + ".interface_ft").c_str(),
          po::value<double>(&block_lists[it->getMeshsetId()].matInterf.data.ft)
              ->default_value(-1),
          "ft")(
          (prefix + ".interface_Gf").c_str(),
          po::value<double>(&block_lists[it->getMeshsetId()].matInterf.data.Gf)
              ->default_value(-1),
          "Gf");
      // TODO Add users parameters
    }
 
    // Displacement bc
    {
      // char flag1; //< Flag for X-Translation (0: N/A, 1: specified)
      // char flag2; //< Flag for Y-Translation (0: N/A, 1: specified)
      // char flag3; //< Flag for Z-Translation (0: N/A, 1: specified)
      // char flag4; //< Flag for X-Rotation (0: N/A, 1: specified)
      // char flag5; //< Flag for Y-Rotation (0: N/A, 1: specified)
      // char flag6; //< Flag for Z-Rotation (0: N/A, 1: specified)
      // double value1; //< Value for X-Translation
      // double value2; //< Value for Y-Translation
      // double value3; //< Value for Z-Translation
      // double value4; //< Value for X-Rotation
      // double value5; //< Value for Y-Rotation
      // double value6; //< Value for Z-Rotation
      configFileOptionsPtr->add_options()(
          (prefix + ".disp_flag1").c_str(),
          po::value<char>(&block_lists[it->getMeshsetId()].dispBc.data.flag1)
              ->default_value(0),
          "flag1")(
          (prefix + ".disp_flag2").c_str(),
          po::value<char>(&block_lists[it->getMeshsetId()].dispBc.data.flag2)
              ->default_value(0),
          "flag2")(
          (prefix + ".disp_flag3").c_str(),
          po::value<char>(&block_lists[it->getMeshsetId()].dispBc.data.flag3)
              ->default_value(0),
          "flag3")(
          (prefix + ".disp_flag4").c_str(),
          po::value<char>(&block_lists[it->getMeshsetId()].dispBc.data.flag4)
              ->default_value(0),
          "flag4")(
          (prefix + ".disp_flag5").c_str(),
          po::value<char>(&block_lists[it->getMeshsetId()].dispBc.data.flag5)
              ->default_value(0),
          "flag5")(
          (prefix + ".disp_flag6").c_str(),
          po::value<char>(&block_lists[it->getMeshsetId()].dispBc.data.flag6)
              ->default_value(0),
          "flag6")(
          (prefix + ".disp_ux").c_str(),
          po::value<double>(&block_lists[it->getMeshsetId()].dispBc.data.value1)
              ->default_value(0),
          "value1")(
          (prefix + ".disp_uy").c_str(),
          po::value<double>(&block_lists[it->getMeshsetId()].dispBc.data.value2)
              ->default_value(0),
          "value2")(
          (prefix + ".disp_uz").c_str(),
          po::value<double>(&block_lists[it->getMeshsetId()].dispBc.data.value3)
              ->default_value(0),
          "value3")(
          (prefix + ".disp_rx").c_str(),
          po::value<double>(&block_lists[it->getMeshsetId()].dispBc.data.value4)
              ->default_value(0),
          "value4")(
          (prefix + ".disp_ry").c_str(),
          po::value<double>(&block_lists[it->getMeshsetId()].dispBc.data.value5)
              ->default_value(0),
          "value5")(
          (prefix + ".disp_rz").c_str(),
          po::value<double>(&block_lists[it->getMeshsetId()].dispBc.data.value6)
              ->default_value(0),
          "value6");
    }
    // Force BC data
    {
      // char zero[3]; //< 3 zeros
      // double value1; //< Force magnitude
      // double value2; //< Moment magnitude
      // double value3; //< X-component of force direction vector
      // double value4; //< Y-component of force direction vector
      // double value5; //< Z-component of force direction vector
      // double value6; //< X-component of moment direction vector
      // double value7; //< Y-component of moment direction vector
      // double value8; //< Z-component of moment direction vector
      // char zero2; // 0
      configFileOptionsPtr->add_options()(
          (prefix + ".force_magnitude").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].forceBc.data.value1)
              ->default_value(0),
          "value1")((prefix + ".moment_magnitude").c_str(),
                    po::value<double>(
                        &block_lists[it->getMeshsetId()].forceBc.data.value2)
                        ->default_value(0),
                    "value2")(
          (prefix + ".force_fx").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].forceBc.data.value3)
              ->default_value(0),
          "value3")((prefix + ".force_fy").c_str(),
                    po::value<double>(
                        &block_lists[it->getMeshsetId()].forceBc.data.value4)
                        ->default_value(0),
                    "value4")(
          (prefix + ".force_fz").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].forceBc.data.value5)
              ->default_value(0),
          "value5")((prefix + ".moment_mx").c_str(),
                    po::value<double>(
                        &block_lists[it->getMeshsetId()].forceBc.data.value6)
                        ->default_value(0),
                    "value6")(
          (prefix + ".moment_my").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].forceBc.data.value7)
              ->default_value(0),
          "value7")((prefix + ".moment_mz").c_str(),
                    po::value<double>(
                        &block_lists[it->getMeshsetId()].forceBc.data.value8)
                        ->default_value(0),
                    "value8");
    }
    {
      // char name[11]; //< 11 characters for "Temperature"
      // char pre1; //< This is always zero
      // char pre2; //< 0: temperature is not applied on thin shells
      // (default); 1: temperature is applied on thin shells char flag1; //<
      // 0: N/A, 1: temperature value applied (not on thin shells) char flag2;
      // //< 0: N/A, 1: temperature applied on thin shell middle char flag3;
      // //< 0: N/A, 1: thin shell temperature gradient specified char flag4;
      // //< 0: N/A, 1: top thin shell temperature char flag5; //< 0: N/A, 1:
      // bottom thin shell temperature char flag6; //< This is always zero
      // double value1; //< Temperature (default case - no thin shells)
      // double value2; //< Temperature for middle of thin shells
      // double value3; //< Temperature gradient for thin shells
      // double value4; //< Temperature for top of thin shells
      // double value5; //< Temperature for bottom of thin shells
      // double value6; //< This is always zero, i.e. ignore
      configFileOptionsPtr->add_options()(
          (prefix + ".temperature_flag1").c_str(),
          po::value<char>(
              &block_lists[it->getMeshsetId()].temperatureBc.data.flag1)
              ->default_value(0),
          "flag1")(
          (prefix + ".temperature_t").c_str(),
          po::value<double>(
              &block_lists[it->getMeshsetId()].temperatureBc.data.value1)
              ->default_value(0),
          "value1");
      // TODO: Add more cases, see above
    }
    // Sideset
    {
      // char name[8];   //< 8 characters for "Pressure"
      // char zero;      //< This is always zero
      // char flag2;     //< 0: Pressure is interpreted as pure pressure 1:
      // pressure is interpreted as total force double value1;  //< Pressure
      // value
      configFileOptionsPtr->add_options()(
          (prefix + ".pressure_flag2").c_str(),
          po::value<char>(
              &block_lists[it->getMeshsetId()].pressureBc.data.flag2)
              ->default_value(0),
          "flag2")((prefix + ".pressure_magnitude").c_str(),
                   po::value<double>(
                       &block_lists[it->getMeshsetId()].pressureBc.data.value1)
                       ->default_value(0),
                   "value1");
    }
    {
      // char name[8]; //< 8 characters for "HeatFlux" (no space)
      // char pre1; //< This is always zero
      // char pre2; //< 0: heat flux is not applied on thin shells (default);
      // 1: heat flux is applied on thin shells char flag1; //< 0: N/A, 1:
      // normal heat flux case (i.e. single value, case without thin shells)
      // char flag2; //< 0: N/A, 1: Thin shell top heat flux specified
      // char flag3; //< 0: N/A, 1: Thin shell bottom heat flux specidied
      // double value1; //< Heat flux value for default case (no thin shells)
      // double value2; //< Heat flux (thin shell top)
      // double value3; //< Heat flux (thin shell bottom)
      configFileOptionsPtr->add_options()(
          (prefix + ".heatflux_flag1").c_str(),
          po::value<char>(
              &block_lists[it->getMeshsetId()].heatFluxBc.data.flag1)
              ->default_value(0),
          "flag1")((prefix + ".heatflux_magnitude").c_str(),
                   po::value<double>(
                       &block_lists[it->getMeshsetId()].heatFluxBc.data.value1)
                       ->default_value(0),
                   "value1");
    }
    // Interface set
    {
      configFileOptionsPtr->add_options()(
          (prefix + ".interface_type").c_str(),
          po::value<char>(&block_lists[it->getMeshsetId()].cfgBc.data.type)
              ->default_value(0),
          "type");
    }
  }
 
  map<int, BlockData> block_set_attributes;
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
    block_set_attributes[it->getMeshsetId()].cubitMeshset = it->getMeshset();
    block_set_attributes[it->getMeshsetId()].iD = it->getMeshsetId();
    block_set_attributes[it->getMeshsetId()].bcType = BLOCKSET;
    std::string block_name = it->getName();
    block_name.erase(
        std::remove_if(block_name.begin(), block_name.end(), ::isspace),
        block_name.end());
    block_set_attributes[it->getMeshsetId()].nAme = block_name;
    // Only blocks which have name
    if (block_name.compare("NoNameSet") != 0) {
      std::string prefix = "SET_ATTR_" + block_name;
      // Block attributes
      add_block_attributes(prefix, block_set_attributes, it);
    }
  }
 
  po::parsed_options parsed =
      parse_config_file(ini_file, *configFileOptionsPtr, true);
  store(parsed, vm);
  po::notify(vm);
 
  // Set type from name
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
 
    CubitBCType bc_type;
    unsigned jj = 0;
    while (1 << jj != LASTSET_BC) {
      if (string(CubitBCNames[jj + 1]) ==
          block_lists[it->getMeshsetId()].addType) {
        bc_type = 1 << jj;
      }
      ++jj;
    }
    if (bc_type.none()) {
      block_lists[it->getMeshsetId()].bcType = UNKNOWNSET;
      // Skip the bockset nothing is defined for it
      continue;
    }
 
    if (bc_type.to_ulong() == BLOCKSET)
      block_lists[it->getMeshsetId()].bcType = BLOCKSET;
    else if (bc_type.to_ulong() == NODESET)
      block_lists[it->getMeshsetId()].bcType = NODESET;
    else if (bc_type.to_ulong() == SIDESET)
      block_lists[it->getMeshsetId()].bcType = SIDESET;
    else {
      SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
               "Not yet implemented type %s\n",
               block_lists[it->getMeshsetId()].addType.c_str());
    }
    if (block_lists[it->getMeshsetId()].iD == -1) {
      SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
               "Unset iD number %d\n", block_lists[it->getMeshsetId()].iD);
    }
  }
 
  std::vector<std::string> additional_parameters;
  additional_parameters =
      collect_unrecognized(parsed.options, po::include_positional);
  for (std::vector<std::string>::iterator vit = additional_parameters.begin();
       vit != additional_parameters.end(); vit++) {
    MOFEM_LOG_C("MeshsetMngSelf", Sev::warning, "Unrecognized option %s",
                vit->c_str());
  }
  for (map<int, BlockData>::iterator mit = block_lists.begin();
       mit != block_lists.end(); mit++) {
    CubitMeshSet_multiIndex::iterator cubit_meshset_it =
        cubitMeshsets.find(mit->second.cubitMeshset);
    if (cubit_meshset_it == cubitMeshsets.end()) {
      SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
              "Data inconsistency\n");
    }
    switch (mit->second.bcType) {
    case UNKNOWNSET:
      break;
    case BLOCKSET: {
      if ((CubitBCType(mit->second.bcType) & cubit_meshset_it->getBcType())
              .any() &&
          mit->second.iD == cubit_meshset_it->getMeshsetId()) {
        // Meshset is the same, only modification
      } else {
        CHKERR addMeshset(mit->second.bcType, mit->second.iD, mit->second.nAme);
        EntityHandle meshset = cubit_meshset_it->getMeshset();
        CHKERR addEntitiesToMeshset(mit->second.bcType, mit->second.iD,
                                    &meshset, 1);
      }
      // Add attributes
      CHKERR setAtributes(mit->second.bcType, mit->second.iD, mit->second.aTtr);
      // Add material elastic data if value are physical (i.e. Young > 0,
      // Poisson in (-1.0.5) and ThermalExpansion>0)
      if (mit->second.matElastic.data.Young != -1) {
        CHKERR setAtributesByDataStructure(mit->second.bcType, mit->second.iD,
                                           mit->second.matElastic);
      }
      if (mit->second.matTransIso.data.Youngp != -1) {
        CHKERR setAtributesByDataStructure(mit->second.bcType, mit->second.iD,
                                           mit->second.matTransIso);
      }
      if (mit->second.matThermal.data.Conductivity != -1) {
        CHKERR setAtributesByDataStructure(mit->second.bcType, mit->second.iD,
                                           mit->second.matThermal);
      }
      if (mit->second.matInterf.data.ft != -1) {
        CHKERR setAtributesByDataStructure(mit->second.bcType, mit->second.iD,
                                           mit->second.matInterf);
      }
    } break;
    case NODESET: {
      if ((CubitBCType(mit->second.bcType) & cubit_meshset_it->getBcType())
              .any() &&
          mit->second.iD == cubit_meshset_it->getMeshsetId()) {
        // Meshset is the same, only modification
      } else {
        CHKERR addMeshset(mit->second.bcType, mit->second.iD);
        EntityHandle meshset = cubit_meshset_it->getMeshset();
        CHKERR addEntitiesToMeshset(mit->second.bcType, mit->second.iD,
                                    &meshset, 1);
      }
      // Add displacement bc
      if (mit->second.dispBc.data.flag1 || mit->second.dispBc.data.flag2 ||
          mit->second.dispBc.data.flag3 || mit->second.dispBc.data.flag4 ||
          mit->second.dispBc.data.flag5 || mit->second.dispBc.data.flag6) {
        if (mit->second.dispBc.data.flag1 == '0')
          mit->second.dispBc.data.flag1 = 0;
        if (mit->second.dispBc.data.flag1 == 'N')
          mit->second.dispBc.data.flag1 = 0;
        if (mit->second.dispBc.data.flag1)
          mit->second.dispBc.data.flag1 = 1;
        if (mit->second.dispBc.data.flag2 == '0')
          mit->second.dispBc.data.flag2 = 0;
        if (mit->second.dispBc.data.flag2 == 'N')
          mit->second.dispBc.data.flag2 = 0;
        if (mit->second.dispBc.data.flag2)
          mit->second.dispBc.data.flag2 = 1;
        if (mit->second.dispBc.data.flag3 == '0')
          mit->second.dispBc.data.flag3 = 0;
        if (mit->second.dispBc.data.flag3 == 'N')
          mit->second.dispBc.data.flag3 = 0;
        if (mit->second.dispBc.data.flag3)
          mit->second.dispBc.data.flag3 = 1;
        if (mit->second.dispBc.data.flag4 == '0')
          mit->second.dispBc.data.flag4 = 0;
        if (mit->second.dispBc.data.flag4 == 'N')
          mit->second.dispBc.data.flag4 = 0;
        if (mit->second.dispBc.data.flag4)
          mit->second.dispBc.data.flag4 = 1;
        if (mit->second.dispBc.data.flag5 == '0')
          mit->second.dispBc.data.flag5 = 0;
        if (mit->second.dispBc.data.flag5 == 'N')
          mit->second.dispBc.data.flag5 = 0;
        if (mit->second.dispBc.data.flag5)
          mit->second.dispBc.data.flag5 = 1;
        if (mit->second.dispBc.data.flag6 == '0')
          mit->second.dispBc.data.flag6 = 0;
        if (mit->second.dispBc.data.flag6 == 'N')
          mit->second.dispBc.data.flag6 = 0;
        if (mit->second.dispBc.data.flag6)
          mit->second.dispBc.data.flag6 = 1;
        CHKERR setBcData(mit->second.bcType, mit->second.iD,
                         mit->second.dispBc);
      }
      if (mit->second.forceBc.data.value1 != 0 ||
          mit->second.forceBc.data.value2 != 0) {
        CHKERR setBcData(mit->second.bcType, mit->second.iD,
                         mit->second.forceBc);
      }
      // Add temperature boundary condition
      if (mit->second.temperatureBc.data.flag1) {
        if (mit->second.temperatureBc.data.flag1 == '0')
          mit->second.temperatureBc.data.flag1 = 0;
        if (mit->second.temperatureBc.data.flag1 == 'N')
          mit->second.temperatureBc.data.flag1 = 0;
        if (mit->second.temperatureBc.data.flag1)
          mit->second.temperatureBc.data.flag1 = 1;
        CHKERR setBcData(mit->second.bcType, mit->second.iD,
                         mit->second.temperatureBc);
      }
    } break;
    case SIDESET: {
      if ((CubitBCType(mit->second.bcType) & cubit_meshset_it->getBcType())
              .any() &&
          mit->second.iD == cubit_meshset_it->getMeshsetId()) {
        // Meshset is the same, only modification
      } else {
        CHKERR addMeshset(mit->second.bcType, mit->second.iD);
        EntityHandle meshset = cubit_meshset_it->getMeshset();
        CHKERR addEntitiesToMeshset(mit->second.bcType, mit->second.iD,
                                    &meshset, 1);
      }
      // Add pressure
      if (mit->second.pressureBc.data.value1 != 0) {
        if (mit->second.pressureBc.data.flag2 == '0')
          mit->second.pressureBc.data.flag2 = 0;
        if (mit->second.pressureBc.data.flag2 == 'N')
          mit->second.pressureBc.data.flag2 = 0;
        if (mit->second.pressureBc.data.flag2)
          mit->second.pressureBc.data.flag2 = 1;
        CHKERR setBcData(mit->second.bcType, mit->second.iD,
                         mit->second.pressureBc);
      }
      // Add heat flux
      if (mit->second.heatFluxBc.data.value1 != 0) {
        if (mit->second.heatFluxBc.data.flag1 == '0')
          mit->second.heatFluxBc.data.flag1 = 0;
        if (mit->second.heatFluxBc.data.flag1 == 'N')
          mit->second.heatFluxBc.data.flag1 = 0;
        if (mit->second.heatFluxBc.data.flag1)
          mit->second.heatFluxBc.data.flag1 = 1;
        CHKERR setBcData(mit->second.bcType, mit->second.iD,
                         mit->second.heatFluxBc);
      }
      // Add Interface
      if (mit->second.cfgBc.data.type != 0) {
        CHKERR setBcData(mit->second.bcType, mit->second.iD, mit->second.cfgBc);
      }
    } break;
    default:
      SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
              "Not yet implemented type\n");
    }
  }
 
  for (auto set_attr : block_set_attributes) {
    // Add attributes
    if (set_attr.second.numberOfAttributes > 0) {
      MOFEM_LOG("MeshsetMngSelf", Sev::verbose)
          << "Set attributes to blockset " << set_attr.second.nAme;
      set_attr.second.aTtr.resize(set_attr.second.numberOfAttributes);
      CHKERR setAtributes(set_attr.second.bcType, set_attr.second.iD,
                          set_attr.second.aTtr);
    }
  }
 
  MoFEMFunctionReturn(0);
}

◆ sortMeshsets()

void MoFEM::MeshsetsManager::sortMeshsets ( std::vector< const CubitMeshSets * > & vec_ptr )

static

Definition at line 12 of file MeshsetsManager.cpp.

                                               {
  std::sort(vec_ptr.begin(), vec_ptr.end(),
            [](const CubitMeshSets *a, const CubitMeshSets *b) {
              if (a->getBcType().to_ullong() < b->getBcType().to_ullong())
                return true;
              if (a->getBcType().to_ullong() > b->getBcType().to_ullong())
                return false;
              else
                return a->getMeshsetId() < b->getMeshsetId();
            });
}

◆ updateAllMeshsetsByEntitiesChildren()

MoFEMErrorCode MoFEM::MeshsetsManager::updateAllMeshsetsByEntitiesChildren ( const BitRefLevel & bit )

Update all blolsets, sidesets and node sets.

Parameters

bit

Returns: MoFEMErrorCode

Definition at line 1450 of file MeshsetsManager.cpp.

                                                                           {
  MoFEMFunctionBegin;
  BitRefManager *bit_mng = cOre.getInterface<BitRefManager>();
  for (_IT_CUBITMESHSETS_FOR_LOOP_((*this), iit)) {
    EntityHandle meshset = iit->getMeshset();
    for (EntityType t = MBVERTEX; t != MBENTITYSET; ++t)
      CHKERR bit_mng->updateMeshsetByEntitiesChildren(meshset, bit, meshset, t,
                                                      true);
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

◆ bhTag

Tag MoFEM::MeshsetsManager::bhTag

protected

Definition at line 895 of file MeshsetsManager.hpp.

◆ bhTag_header

Tag MoFEM::MeshsetsManager::bhTag_header

protected

Definition at line 896 of file MeshsetsManager.hpp.

◆ broadcastMeshsetsOn

bool MoFEM::MeshsetsManager::broadcastMeshsetsOn = true

static

if true meshsets are synchronised between processors

Definition at line 885 of file MeshsetsManager.hpp.

◆ configFileOptionsPtr

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

protected

config file options

Definition at line 901 of file MeshsetsManager.hpp.

◆ cOre

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

Definition at line 109 of file MeshsetsManager.hpp.

◆ cubitMeshsets

CubitMeshSet_multiIndex MoFEM::MeshsetsManager::cubitMeshsets

protected

cubit meshsets

Definition at line 899 of file MeshsetsManager.hpp.

◆ nsTag

Tag MoFEM::MeshsetsManager::nsTag

protected

Definition at line 891 of file MeshsetsManager.hpp.

◆ nsTag_data

Tag MoFEM::MeshsetsManager::nsTag_data

protected

Definition at line 893 of file MeshsetsManager.hpp.

◆ ssTag

Tag MoFEM::MeshsetsManager::ssTag

protected

Definition at line 892 of file MeshsetsManager.hpp.

◆ ssTag_data

Tag MoFEM::MeshsetsManager::ssTag_data

protected

Definition at line 894 of file MeshsetsManager.hpp.

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

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

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ MeshsetsManager()

◆ ~MeshsetsManager()

Member Function Documentation

◆ addEntitiesToMeshset() [1/2]

◆ addEntitiesToMeshset() [2/2]

◆ broadcastMeshsets()

◆ checkIfMeshsetContainsEntities()

◆ checkMeshset()

◆ clearMap()

◆ deleteMeshset()

◆ get_bhTag()

◆ get_bhTag_header()

◆ get_meshsets_manager_ptr() [1/2]

◆ get_meshsets_manager_ptr() [2/2]

◆ get_nsTag()

◆ get_nsTag_data()

◆ get_ssTag()

◆ get_ssTag_data()

◆ getConfigFileOptionsPtr()

◆ getMeshsetsMultindex()

◆ getTags() [1/2]

◆ getTags() [2/2]

◆ initialiseDatabaseFromMesh()

◆ printBcSet()

◆ printDisplacementSet()

◆ printForceSet()

◆ printHeatFluxSet()

◆ printMaterialsSet()

◆ printPressureSet()

◆ printTemperatureSet()

◆ query_interface()

◆ readMeshsets()

◆ saveMeshsetToFile() [1/2]

◆ saveMeshsetToFile() [2/2]

◆ setAtributes()

◆ setAtributesByDataStructure()

◆ setBcData()

◆ setMeshsetFromFile() [1/2]

◆ setMeshsetFromFile() [2/2]

◆ sortMeshsets()

◆ updateAllMeshsetsByEntitiesChildren()

Member Data Documentation

◆ bhTag

◆ bhTag_header

◆ broadcastMeshsetsOn

◆ configFileOptionsPtr

◆ cOre

◆ cubitMeshsets

◆ nsTag

◆ nsTag_data

◆ ssTag

◆ ssTag_data