Simple interface for fast problem set-up. More...

#include <src/interfaces/BcManager.hpp>

Inheritance diagram for MoFEM::BcManager:

Collaboration diagram for MoFEM::BcManager:

Classes
struct	BCs
	Data structure storing bc markers and atributes. More...

Public Types
using	BcMapByBlockName = std::map< string, boost::shared_ptr< BCs > >

using	BcMarkerPtr = boost::shared_ptr< std::vector< char unsigned > >

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

	BcManager (const MoFEM::Core &core)

virtual	~BcManager ()=default

MoFEMErrorCode	getOptions ()
	get options More...

Range	getAdjEnts (Range ents)

MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string block_name, const std::string field_name, int lo, int hi, bool get_low_dim_ents=true)
	Remove DOFs from problem. More...

MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string block_name, const std::string field_name, int lo, int hi, bool get_low_dim_ents=true)
	Mark block dofs. More...

boost::shared_ptr< BCs >	popMarkDOFsOnEntities (const std::string block_name)
	Get bc data and remove element. More...

auto	getBcStructure (const std::string bc_id)
	Get the bc structure object. More...

BcMapByBlockName &	getBcMapByBlockName ()
	Get the bc map. More...

BcMarkerPtr	getMergedBlocksMarker (std::vector< std::regex > bc_regex_vec)
	Get the Merged Boundary Marker object. More...

auto	getMergedBlocksMarker (std::vector< string > bc_names)
	Get the Merged Boundary Marker object. More...

BcMarkerPtr	getMergedBlocksMarker (const std::vector< BcMarkerPtr > &boundary_markers_ptr_vec)
	Get the Merged Blocks Marker object. More...

auto	checkBlock (const std::pair< string, boost::shared_ptr< BCs >> &bc, std::regex reg)
	check if given boundary condition name is in the map bc element More...

auto	checkBlock (const std::pair< string, boost::shared_ptr< BCs >> &bc, std::string name)
	check if given boundary condition name is in the map bc element More...

SmartPetscObj< IS >	getBlockIS (const std::string block_prefix, const std::string block_name, const std::string field_name, const std::string problem_name, int lo, int hi, SmartPetscObj< IS > is_expand=SmartPetscObj< IS >())
	Get block is. More...

SmartPetscObj< IS >	getBlockIS (const std::string problem_name, const std::string block_name, const std::string field_name, int lo, int hi, SmartPetscObj< IS > is_expand=SmartPetscObj< IS >())
	Get block is. More...

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

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

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

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

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

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

virtual	~UnknownInterface ()=default

Private Attributes
MoFEM::Core &	cOre

BcMapByBlockName	bcMapByBlockName

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

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

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

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

Detailed Description

Simple interface for fast problem set-up.

Definition at line 27 of file BcManager.hpp.

Member Typedef Documentation

◆ BcMapByBlockName

using MoFEM::BcManager::BcMapByBlockName = std::map<string, boost::shared_ptr<BCs> >

Todo:: Add markers for standard BCs from cubit on Nodests and Sidesets used bu cubit for displacements, forces, etc. Also add function to add blockset by id and type.

Definition at line 113 of file BcManager.hpp.

◆ BcMarkerPtr

using MoFEM::BcManager::BcMarkerPtr = boost::shared_ptr<std::vector<char unsigned> >

Definition at line 115 of file BcManager.hpp.

Constructor & Destructor Documentation

◆ BcManager()

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

Definition at line 30 of file BcManager.cpp.

                                      : cOre(const_cast<Core &>(core)) {
  
   if (!LogManager::checkIfChannelExist("BcMngWorld")) {
     auto core_log = logging::core::get();
  
     core_log->add_sink(
         LogManager::createSink(LogManager::getStrmWorld(), "BcMngWorld"));
     core_log->add_sink(
         LogManager::createSink(LogManager::getStrmSync(), "BcMngSync"));
     core_log->add_sink(
         LogManager::createSink(LogManager::getStrmSelf(), "BcMngSelf"));
  
     LogManager::setLog("BcMngWorld");
     LogManager::setLog("BcMngSync");
     LogManager::setLog("BcMngSelf");
  
     MOFEM_LOG_TAG("BcMngWorld", "BcMng");
     MOFEM_LOG_TAG("BcMngSync", "BcMng");
     MOFEM_LOG_TAG("BcMngSelf", "BcMng");
   }
  
   MOFEM_LOG("BcMngWorld", Sev::noisy) << "BC manager created";
 }

◆ ~BcManager()

virtual MoFEM::BcManager::~BcManager ( )

virtualdefault

Member Function Documentation

◆ checkBlock() [1/2]

auto MoFEM::BcManager::checkBlock	(	const std::pair< string, boost::shared_ptr< BCs >> &	bc,
		std::regex	reg
	)

check if given boundary condition name is in the map bc element

Parameters

bc	element of the map
reg	bc regex

Returns: auto

Definition at line 169 of file BcManager.hpp.

                                        {
     return std::regex_match(bc.first, reg);
   }

◆ checkBlock() [2/2]

auto MoFEM::BcManager::checkBlock	(	const std::pair< string, boost::shared_ptr< BCs >> &	bc,
		std::string	name
	)

check if given boundary condition name is in the map bc element

Parameters

bc	element of the map
name	bc name

Returns: auto

Definition at line 180 of file BcManager.hpp.

                                          {
     string full_name = string("(.*)_") + name + string("(.*)");
     return checkBlock(bc, std::regex(full_name));
   }

◆ getAdjEnts()

Range MoFEM::BcManager::getAdjEnts ( Range ents )

◆ getBcMapByBlockName()

BcMapByBlockName& MoFEM::BcManager::getBcMapByBlockName ( )

Get the bc map.

Returns: auto

Definition at line 131 of file BcManager.hpp.

131 { return bcMapByBlockName; }

MoFEM::BcManager::bcMapByBlockName

BcMapByBlockName bcMapByBlockName

Definition: BcManager.hpp:222

◆ getBcStructure()

auto MoFEM::BcManager::getBcStructure ( const std::string bc_id )

Get the bc structure object.

Parameters

block_name

Returns: auto

Definition at line 122 of file BcManager.hpp.

                                                     {
     return bcMapByBlockName.at(bc_id);
   }

◆ getBlockIS() [1/2]

SmartPetscObj< IS > MoFEM::BcManager::getBlockIS	(	const std::string	block_prefix,
		const std::string	block_name,
		const std::string	field_name,
		const std::string	problem_name,
		int	lo,
		int	hi,
		SmartPetscObj< IS >	is_expand = `SmartPetscObj<IS>()`
	)

Get block is.

Parameters

block_prefix	for hashmap
block_name	for hash map
field_name	for hash map and IS
problem_name	for IS
lo
hi
is_expand	is to extend

Returns: SmartPetscObj<IS>

Definition at line 238 of file BcManager.cpp.

                                                                              {
   Interface &m_field = cOre;
  
   const std::string bc_id =
       block_prefix + "_" + field_name + "_" + block_name + "(.*)";
  
   Range bc_ents;
   for (auto bc : getBcMapByBlockName()) {
     if (std::regex_match(bc.first, std::regex(bc_id))) {
       bc_ents.merge(*(bc.second->getBcEntsPtr()));
       MOFEM_LOG("BcMngSelf", Sev::noisy)
           << "Get entities from block and add to IS. Block name " << bc.first;
     }
   }
  
   SmartPetscObj<IS> is_bc;
   auto get_is = [&]() {
     MoFEMFunctionBegin;
     CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(bc_ents);
     CHKERR m_field.getInterface<ISManager>()->isCreateProblemFieldAndRank(
         problem_name, ROW, field_name, lo, hi, is_bc, &bc_ents);
     if (is_expand) {
       IS is_tmp;
       CHKERR ISExpand(is_bc, is_expand, &is_tmp);
       is_bc = SmartPetscObj<IS>(is_tmp);
     }
     CHKERR ISSort(is_bc);
     MoFEMFunctionReturn(0);
   };
  
   if (get_is())
     CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY, "IS is not created");
  
   return is_bc;
 }

◆ getBlockIS() [2/2]

SmartPetscObj< IS > MoFEM::BcManager::getBlockIS	(	const std::string	problem_name,
		const std::string	block_name,
		const std::string	field_name,
		int	lo,
		int	hi,
		SmartPetscObj< IS >	is_expand = `SmartPetscObj<IS>()`
	)

Get block is.

Parameters

problem_name
block_name
field_name
lo
hi
is_expand	is to extend

Returns: SmartPetscObj<IS>

Definition at line 278 of file BcManager.cpp.

                                                                              {
   return getBlockIS(problem_name, block_name, field_name, problem_name, lo, hi,
                     is_expand);
 }

◆ getMergedBlocksMarker() [1/3]

BcManager::BcMarkerPtr MoFEM::BcManager::getMergedBlocksMarker ( const std::vector< BcMarkerPtr > & boundary_markers_ptr_vec )

Get the Merged Blocks Marker object.

Parameters

boundary_markers_ptr_vec vector of boundary markers to merge

Returns: BcMarkerPtr

Definition at line 227 of file BcManager.cpp.

                                                                      {
   auto boundary_marker_ptr = boost::make_shared<std::vector<char unsigned>>();
   for (auto &bcm : boundary_markers_ptr_vec) {
     boundary_marker_ptr->resize(bcm->size(), 0);
     for (int i = 0; i != bcm->size(); ++i)
       (*boundary_marker_ptr)[i] |= (*bcm)[i];
   }
   return boundary_marker_ptr;
 }

◆ getMergedBlocksMarker() [2/3]

BcManager::BcMarkerPtr MoFEM::BcManager::getMergedBlocksMarker ( std::vector< std::regex > bc_regex_vec )

Get the Merged Boundary Marker object.

Parameters

bc_regex_vec boundary name regex vector

Returns: boundaryMarker

Definition at line 209 of file BcManager.cpp.

                                                                  {
   BcManager::BcMarkerPtr boundary_marker_ptr;
   if (bcMapByBlockName.size()) {
     boundary_marker_ptr = boost::make_shared<std::vector<char unsigned>>();
     for (auto b : bcMapByBlockName) {
       for (auto &reg_name : bc_regex_vec) {
         if (std::regex_match(b.first, reg_name)) {
           boundary_marker_ptr->resize(b.second->bcMarkers.size(), 0);
           for (int i = 0; i != b.second->bcMarkers.size(); ++i) {
             (*boundary_marker_ptr)[i] |= b.second->bcMarkers[i];
           }
         }
       }
     }
   }
   return boundary_marker_ptr;
 }

◆ getMergedBlocksMarker() [3/3]

auto MoFEM::BcManager::getMergedBlocksMarker ( std::vector< string > bc_names )

Get the Merged Boundary Marker object.

Parameters

bc_names vector of boundary names

Returns: boundaryMarker

Definition at line 146 of file BcManager.hpp.

                                                                 {
     std::vector<std::regex> reg_vec(bc_names.size());
     for (int i = 0; i != bc_names.size(); ++i) {
       string full_name = string("(.*)_") + bc_names[i] + string("(.*)");
       reg_vec[i] = std::regex(full_name);
     }
     return getMergedBlocksMarker(reg_vec);
   }

◆ getOptions()

MoFEMErrorCode MoFEM::BcManager::getOptions ( )

get options

Returns: error code

Definition at line 54 of file BcManager.cpp.

                                      {
   PetscBool flg = PETSC_TRUE;
   MoFEMFunctionBeginHot;
   ierr = PetscOptionsBegin(PETSC_COMM_WORLD, "", "BcManager options", "none");
   ierr = PetscOptionsEnd();
   CHKERRG(ierr);
   MoFEMFunctionReturnHot(0);
 }

◆ popMarkDOFsOnEntities()

boost::shared_ptr< BcManager::BCs > MoFEM::BcManager::popMarkDOFsOnEntities ( const std::string block_name )

Get bc data and remove element.

Parameters

block_name

Returns: boost::shared_ptr<BCs>

Definition at line 198 of file BcManager.cpp.

                                                            {
   auto bc_it = bcMapByBlockName.find(block_name);
   if (bc_it != bcMapByBlockName.end()) {
     auto bc = bc_it->second;
     bcMapByBlockName.erase(bc_it);
     return bc;
   }
   return boost::shared_ptr<BCs>();
 }

◆ pushMarkDOFsOnEntities()

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	block_name,
		const std::string	field_name,
		int	lo,
		int	hi,
		bool	get_low_dim_ents = `true`
	)

Mark block dofs.

Parameters

problem_name
block_name
field_name
lo	lowest coefficient
hi	highest coefficient
get_low_dim_ents	get lower dimension entities

Returns: MoFEMErrorCode

Definition at line 121 of file BcManager.cpp.

                                                                         {
   Interface &m_field = cOre;
   auto prb_mng = m_field.getInterface<ProblemsManager>();
   MoFEMFunctionBegin;
  
   auto get_dim = [&](const Range &ents) {
     for (auto d : {3, 2, 1})
       if (ents.num_of_dimension(d))
         return d;
     return 0;
   };
  
   auto get_adj_ents = [&](const Range &ents) {
     Range verts;
     CHKERR m_field.get_moab().get_connectivity(ents, verts, true);
     const auto dim = get_dim(ents);
     for (size_t d = 1; d < dim; ++d)
       CHKERR m_field.get_moab().get_adjacencies(ents, d, false, verts,
                                                 moab::Interface::UNION);
     verts.merge(ents);
     CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(verts);
     return verts;
   };
  
   auto fix_disp = [&]() {
     MoFEMFunctionBegin;
  
     auto mark_fix_dofs = [&](std::vector<unsigned char> &marked_field_dofs,
                              const auto lo, const auto hi) {
       return prb_mng->modifyMarkDofs(problem_name, ROW, field_name, lo, hi,
                                      ProblemsManager::MarkOP::OR, 1,
                                      marked_field_dofs);
     };
  
     for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
       if (it->getName().compare(0, block_name.length(), block_name) == 0) {
  
         const std::string bc_id =
             problem_name + "_" + field_name + "_" + it->getName();
  
         auto bc = boost::make_shared<BCs>();
         CHKERR m_field.get_moab().get_entities_by_handle(it->meshset,
                                                          bc->bcEnts, true);
         CHKERR it->getAttributes(bc->bcAttributes);
  
         MOFEM_LOG("BcMngWorld", Sev::verbose)
             << "Found block " << block_name << " number of entities "
             << bc->bcEnts.size() << " number of attributes "
             << bc->bcAttributes.size() << " highest dim of entities "
             << get_dim(bc->bcEnts);
  
         CHKERR mark_fix_dofs(bc->bcMarkers, lo, hi);
         if (get_low_dim_ents) {
           auto low_dim_ents = get_adj_ents(bc->bcEnts);
           CHKERR prb_mng->markDofs(problem_name, ROW, ProblemsManager::AND,
                                    low_dim_ents, bc->bcMarkers);
           bc->bcEnts.swap(low_dim_ents);
         } else
           CHKERR prb_mng->markDofs(problem_name, ROW, ProblemsManager::AND,
                                    bc->bcEnts, bc->bcMarkers);
  
         bcMapByBlockName[bc_id] = bc;
       }
     }
     MoFEMFunctionReturn(0);
   };
  
   CHKERR fix_disp();
  
   MoFEMFunctionReturn(0);
 }

◆ query_interface()

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

virtual

Implements MoFEM::UnknownInterface.

Definition at line 23 of file BcManager.cpp.

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

◆ removeBlockDOFsOnEntities()

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	block_name,
		const std::string	field_name,
		int	lo,
		int	hi,
		bool	get_low_dim_ents = `true`
	)

Remove DOFs from problem.

Parameters

problem_name
block_name
field_name
lo	lowest coefficient
hi	highest coefficient
get_low_dim_ents	get lower dimension entities

Returns: MoFEMErrorCode

Definition at line 63 of file BcManager.cpp.

                                                                        {
   Interface &m_field = cOre;
   auto prb_mng = m_field.getInterface<ProblemsManager>();
   MoFEMFunctionBegin;
  
   auto get_dim = [&](const Range &ents) {
     for (auto d : {3, 2, 1})
       if (ents.num_of_dimension(d))
         return d;
     return 0;
   };
  
   auto get_adj_ents = [&](const Range &ents) {
     Range verts;
     CHKERR m_field.get_moab().get_connectivity(ents, verts, true);
     const auto dim = get_dim(ents);
     for (size_t d = 1; d < dim; ++d)
       CHKERR m_field.get_moab().get_adjacencies(ents, d, false, verts,
                                                 moab::Interface::UNION);
     verts.merge(ents);
     CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(verts);
     return verts;
   };
  
   auto remove_dofs_on_ents = [&](const Range &&ents, const int lo,
                                  const int hi) {
     MoFEMFunctionBegin;
     CHKERR prb_mng->removeDofsOnEntities(problem_name, field_name, ents, lo,
                                          hi);
     MoFEMFunctionReturn(0);
   };
  
   auto get_block_ents = [&](const std::string block_name) {
     Range remove_ents;
     for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
       if (it->getName().compare(0, block_name.length(), block_name) == 0) {
         CHKERR m_field.get_moab().get_entities_by_handle(it->meshset,
                                                          remove_ents, true);
         MOFEM_LOG("BcMngWorld", Sev::verbose)
             << "Found block to remove " << block_name << " number of entities "
             << remove_ents.size() << " highest dim of entities "
             << get_dim(remove_ents);
       }
     }
     return remove_ents;
   };
  
   if (get_lod_dim_ents)
     CHKERR remove_dofs_on_ents(get_adj_ents(get_block_ents(block_name)), lo,
                                hi);
   else
     CHKERR remove_dofs_on_ents(get_block_ents(block_name), lo, hi);
  
   MoFEMFunctionReturn(0);
 }

Member Data Documentation

◆ bcMapByBlockName

BcMapByBlockName MoFEM::BcManager::bcMapByBlockName

private

Definition at line 222 of file BcManager.hpp.

◆ cOre

MoFEM::Core& MoFEM::BcManager::cOre

private

Definition at line 220 of file BcManager.hpp.

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

src/interfaces/BcManager.hpp
src/interfaces/impl/BcManager.cpp

Classes

Public Types

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ BcMapByBlockName

◆ BcMarkerPtr

Constructor & Destructor Documentation

◆ BcManager()

◆ ~BcManager()

Member Function Documentation

◆ checkBlock() [1/2]

◆ checkBlock() [2/2]

◆ getAdjEnts()

◆ getBcMapByBlockName()

◆ getBcStructure()

◆ getBlockIS() [1/2]

◆ getBlockIS() [2/2]

◆ getMergedBlocksMarker() [1/3]

◆ getMergedBlocksMarker() [2/3]

◆ getMergedBlocksMarker() [3/3]

◆ getOptions()

◆ popMarkDOFsOnEntities()

◆ pushMarkDOFsOnEntities()

◆ query_interface()

◆ removeBlockDOFsOnEntities()

Member Data Documentation

◆ bcMapByBlockName

◆ cOre