MoFEM::BcManager Struct Reference

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...
MoFEMErrorCode	addBlockDOFsToMPCs (const std::string problem_name, const std::string field_name, bool get_low_dim_ents=false, bool block_name_field_prefix=false, bool is_distributed_mesh=false)
MoFEMErrorCode	pushMarkSideDofs (const std::string problem_name, const std::string block_name, const std::string field_name, int bridge_dim, int lo, int hi)
	Mark side DOFs. More...
MoFEMErrorCode	removeSideDOFs (const std::string problem_name, const std::string block_name, const std::string field_name, int bridge_dim, int lo, int hi, bool is_distributed_mesh=true)
	Remove side DOFs. More...
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, bool is_distributed_mesh=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...
template<typename T >
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents=true, bool block_name_field_prefix=false, bool is_distributed_mesh=true)
	Mark block DOFs. More...
template<typename T >
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents=true, bool block_name_field_prefix=false)
	Mark block DOFs. More...
template<typename T >
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string block_name, const std::string field_name, bool get_low_dim_ents=true, bool is_distributed_mesh=true)
	Mark block DOFs. More...
template<typename T >
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string block_name, const std::string field_name, 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...
Range	getMergedBlocksRange (std::vector< std::regex > bc_regex_vec)
	Merge block ranges. More...
auto	getMergedBlocksRange (std::vector< string > bc_names)
	Merge block ranges. 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< std::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...
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string block_name, const std::string field_name, bool get_low_dim_ents, bool is_distributed_mesh)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, const std::string block_name, bool get_low_dim_ents)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string block_name, const std::string field_name, bool get_low_dim_ents, bool is_distributed_mesh)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string block_name, const std::string field_name, bool get_low_dim_ents)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	pushMarkDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
template<>
MoFEMErrorCode	removeBlockDOFsOnEntities (const std::string problem_name, const std::string field_name, bool get_low_dim_ents, bool block_name_field_prefix, bool is_distributed_mesh)
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 reference 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

Static Public Member Functions
static std::pair< std::string, std::string >	extractStringFromBlockId (const std::string block_id, const std::string prb_name)
	Extract block name and block name form block id. More...
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...

Private Attributes
MoFEM::Core &	cOre
BcMapByBlockName	bcMapByBlockName

Detailed Description

Simple interface for fast problem set-up.

Examples

adolc_plasticity.cpp, dynamic_first_order_con_law.cpp, EshelbianPlasticity.cpp, free_surface.cpp, heat_equation.cpp, photon_diffusion.cpp, plastic.cpp, plot_base.cpp, poisson_2d_homogeneous.cpp, schur_test_diag_mat.cpp, seepage.cpp, shallow_wave.cpp, tensor_divergence_operator.cpp, test_broken_space.cpp, thermo_elastic.cpp, and wave_equation.cpp.

Definition at line 25 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 225 of file BcManager.hpp.

BcMarkerPtr

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

Definition at line 227 of file BcManager.hpp.

Constructor & Destructor Documentation

BcManager()

MoFEM::BcManager::BcManager

(

const MoFEM::Core &

core

)

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

addBlockDOFsToMPCs()

MoFEMErrorCode MoFEM::BcManager::addBlockDOFsToMPCs	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents = false,
		bool	block_name_field_prefix = false,
		bool	is_distributed_mesh = false
	)

Definition at line 185 of file BcManager.cpp.

                                                                       {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  if (block_name_field_prefix)
    MOFEM_LOG("BcMngWorld", Sev::warning)
        << "Argument block_name_field_prefix=true has no effect";
  if (is_distributed_mesh)
    MOFEM_LOG("BcMngWorld", Sev::warning)
        << "Argument is_distributed_mesh=true has no effect";
  if (get_low_dim_ents)
    MOFEM_LOG("BcMngWorld", Sev::warning)
        << "Argument get_low_dim_ents=true has no effect";
 
  auto get_dim = [&](const Range &ents) {
    for (auto d : {3, 2, 1})
      if (ents.num_of_dimension(d))
        return d;
    return 0;
  };
 
  auto iterate_mpc_meshsets = [&]() {
    MoFEMFunctionBegin;
 
    auto mpc_meshset_ptr =
        m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(
            std::regex((boost::format("%s(.*)") % "MPC_(.*)").str()));
 
    for (auto m : mpc_meshset_ptr) {
 
      if (std::regex_match(m->getName(),
                           std::regex("(.*)COUPLING_LINKS(.*)"))) {
 
        auto bc = boost::make_shared<BCs>();
        bc->mpcPtr = boost::make_shared<MPCsType>();
        bc->mpcPtr->mpcType = MPC::COUPLING;
 
        std::string const corresponding_master_ms =
            std::regex_replace(m->getName(), std::regex("LINKS"), "MASTER");
 
        Range links_ents;
        CHKERR m_field.get_moab().get_entities_by_handle(m->getMeshset(),
                                                         links_ents, true);
 
        Range master_nodes;
        if (m_field.getInterface<MeshsetsManager>()->checkMeshset(
                corresponding_master_ms)) {
          const CubitMeshSets *l;
          CHKERR m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(
              corresponding_master_ms, &l);
          bc->mpcPtr->isReprocitical = false;
 
          CHKERR m_field.get_moab().get_entities_by_handle(l->getMeshset(),
                                                           master_nodes, true);
          // if (master_nodes.subset_by_dimension(0).size() <
          // links_ents.subset_by_dimension(1))
          {
            auto low_dim_ents = BcManagerImplTools::get_adj_ents(
                m_field.get_moab(), master_nodes);
            low_dim_ents = low_dim_ents.subset_by_dimension(0);
            master_nodes.swap(low_dim_ents);
          }
 
          MOFEM_LOG("BcMngWorld", Sev::verbose)
              << "Found block MASTER LINKS block: " << l->getName()
              << " Entities size: " << master_nodes.size();
 
        } else {
          MOFEM_LOG("BcMngWorld", Sev::warning)
              << "MASTER LINKS block not found: " << corresponding_master_ms
              << " setting reprocitical constraint. ("
              << bc->mpcPtr->isReprocitical << ").";
        }
 
        // if (std::regex_match(bc_id, std::regex("(.*)TIE(.*)")))
        //   bc->mpcPtr->mpcType = MPC::TIE;
        // if (std::regex_match(bc_id, std::regex("(.*)RIGID_BODY(.*)")))
        //   bc->mpcPtr->mpcType = MPC::RIGID_BODY;
 
        // std::cout << "links_ents range is: " << links_ents << "\n";
        for (auto &link : links_ents.subset_by_dimension(1)) {
          Range verts;
          CHKERR m_field.get_moab().get_connectivity(&link, 1, verts, true);
          // std::cout << "verts range is: " << verts << "\n";
          if (bc->mpcPtr->isReprocitical) {
            bc->mpcPtr->mpcMasterEnts.insert(verts[0]);
            bc->mpcPtr->mpcSlaveEnts.insert(verts[1]);
          } else {
            for (auto &m_node : verts)
              if (master_nodes.find(m_node) != master_nodes.end()) {
                // std::cout << "found ent: " << m_node << "\n";
                bc->mpcPtr->mpcMasterEnts.insert(m_node);
                bc->mpcPtr->mpcSlaveEnts.merge(
                    subtract(verts, Range(m_node, m_node)));
                break;
              }
          }
        }
 
        MOFEM_LOG("BcMngWorld", Sev::verbose)
            << "Found block MPC LINKS block: " << m->getName()
            << " Entities size (edges): " << links_ents.size()
            << " Entities size (nodes): "
            << bc->mpcPtr->mpcMasterEnts.size() +
                   bc->mpcPtr->mpcSlaveEnts.size()
            << " (" << bc->mpcPtr->mpcMasterEnts.size() << " "
            << bc->mpcPtr->mpcSlaveEnts.size() << ")";
 
        MOFEM_LOG("BcMngSync", Sev::noisy) << *bc->mpcPtr;
        MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::noisy);
 
        // CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(
        //     bc->mpcPtr->mpcMasterEnts);
        // CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(
        //     bc->mpcPtr->mpcSlaveEnts);
        vector<double> mAttributes;
        CHKERR m->getAttributes(mAttributes);
 
        auto setFlags = [&](const auto &flags) {
          auto &d = bc->mpcPtr->data;
          if (flags.empty()) {
            d.flag1 = d.flag2 = d.flag3 = d.flag4 = d.flag5 = d.flag6 = true;
            return;
          }
          for (size_t i = 0; i < std::min(flags.size(), size_t(6)); ++i)
            (&d.flag1)[i] = flags[i] > 0.0;
        };
 
        setFlags(mAttributes);
 
        MOFEM_LOG("BcMngWorld", Sev::verbose)
            << "Found block " << m->getName() << " number of entities "
            << bc->mpcPtr->mpcMasterEnts.size() << " number of attributes "
            << mAttributes.size() << " highest dim of entities "
            << get_dim(bc->mpcPtr->mpcMasterEnts);
 
        // NOTE: we are not using markers at the moment
        // for (int i = 0; i != mAttributes.size(); ++i) {
        //   if (mAttributes[i] > 0.0)
        //     CHKERR mark_fix_dofs(bc->mpcPtr->bcMasterMarkers, i, i);
        // }
 
        // if (mAttributes.empty())
        //   CHKERR mark_fix_dofs(bc->mpcPtr->bcMasterMarkers, 0,
        //                        MAX_DOFS_ON_ENTITY);
 
        const std::string bc_id =
            problem_name + "_" + field_name + "_" + m->getName();
        bcMapByBlockName[bc_id] = bc;
      }
    }
    MoFEMFunctionReturn(0);
  };
 
  CHKERR iterate_mpc_meshsets();
 
  MoFEMFunctionReturn(0);
}

checkBlock() [1/2]

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

inline

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 317 of file BcManager.hpp.

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

checkBlock() [2/2]

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

inline

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 305 of file BcManager.hpp.

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

extractStringFromBlockId()

std::pair< std::string, std::string > MoFEM::BcManager::extractStringFromBlockId ( const std::string  block_id, const std::string  prb_name  )

static

Extract block name and block name form block id.

Parameters

block_id
prb_name

Returns

std::pair<std::string, std::string>

Definition at line 1381 of file BcManager.cpp.

                                                              {
 
  // Assumes that field name is consist with letters and numbers.
  // No special characters.
  auto field_rgx_str =
      (boost::format("%s_([a-zA-Z0-9]*)_(.*)") % prb_name).str();
  std::regex field_rgx(field_rgx_str);
  std::smatch match_field_name;
  std::string field_name;
  std::string block_name;
 
  if (std::regex_search(block_id, match_field_name, field_rgx)) {
    field_name = match_field_name[1];
    block_name = match_field_name[2];
  } else {
    CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                      "Field name and block name can not be resolved");
  }
 
  return std::make_pair(field_name, block_name);
}

getBcMapByBlockName()

BcMapByBlockName& MoFEM::BcManager::getBcMapByBlockName ( )

inline

Get the bc map.

Returns

auto

Examples

photon_diffusion.cpp.

Definition at line 243 of file BcManager.hpp.

{ return bcMapByBlockName; }

getBcStructure()

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

inline

Get the bc structure object.

Parameters

block_name

Returns

auto

Definition at line 234 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 405 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("BcMngWorld", Sev::verbose)
          << "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 445 of file BcManager.cpp.

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

getMergedBlocksMarker() [1/3]

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

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 374 of file BcManager.cpp.

                                                                 {
  BcManager::BcMarkerPtr boundary_marker_ptr;
  if (bcMapByBlockName.size()) {
    for (auto b : bcMapByBlockName) {
      for (auto &reg_name : bc_regex_vec) {
        if (std::regex_match(b.first, reg_name)) {
          if (!boundary_marker_ptr)
            boundary_marker_ptr =
                boost::make_shared<std::vector<char unsigned>>();
          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 )

inline

Get the Merged Boundary Marker object.

Parameters

bc_names

vector of boundary names

Returns

boundaryMarker

Definition at line 281 of file BcManager.hpp.

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

getMergedBlocksRange() [1/2]

Range MoFEM::BcManager::getMergedBlocksRange

(

std::vector< std::regex >

bc_regex_vec

)

Merge block ranges.

Parameters

bc_regex_vec

Returns

Range

Definition at line 359 of file BcManager.cpp.

                                                                      {
  Range ents;
  if (bcMapByBlockName.size()) {
    for (auto b : bcMapByBlockName) {
      for (auto &reg_name : bc_regex_vec) {
        if (std::regex_match(b.first, reg_name)) {
          ents.merge(b.second->bcEnts);
        }
      }
    }
  }
  return ents;
}

getMergedBlocksRange() [2/2]

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

inline

Merge block ranges.

Parameters

bc_names

Returns

auto

Definition at line 259 of file BcManager.hpp.

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

getOptions()

MoFEMErrorCode MoFEM::BcManager::getOptions

(

)

get options

Returns

error code

Definition at line 65 of file BcManager.cpp.

                                     {
  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 349 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() [1/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	block_name,
		const std::string	field_name,
		bool	get_low_dim_ents
	)

Definition at line 1248 of file BcManager.cpp.

                                                       {
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities(problem_name, block_name, field_name, 0,
                                MAX_DOFS_ON_ENTITY, get_low_dim_ents);
 
  auto regex_str =
      (boost::format("%s_%s_%s(.*)") % problem_name % field_name % block_name)
          .str();
 
  for (auto &m : bcMapByBlockName) {
 
    auto &bc_id = m.first;
 
    if (std::regex_match(bc_id, std::regex(regex_str))) {
 
      auto &bc = m.second;
      bc->tempBcPtr = boost::make_shared<TemperatureCubitBcData>();
      bc->tempBcPtr->data.flag1 = 1;
      if (bc->bcAttributes.empty()) {
        bc->tempBcPtr->data.value1 = 0;
        MOFEM_LOG("BcMngWorld", Sev::warning)
            << "Expected one attribute on block but have "
            << bc->bcAttributes.size();
      } else if (bc->bcAttributes.size() >= 1) {
        bc->tempBcPtr->data.value1 = bc->bcAttributes[0];
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [2/25]

template<typename T >

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

Mark block DOFs.

Parameters

problem_name
field_name
block_name
get_low_dim_ents	get lower dimension entities

Returns

MoFEMErrorCode

pushMarkDOFsOnEntities() [3/25]

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 field name

Returns

MoFEMErrorCode

Definition at line 111 of file BcManager.cpp.

                                                                        {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  // if(problem_name.size())
  //   MOFEM_LOG("BcMngWorld", Sev::warning)
  //       << "Argument problem_name has no effect";
 
  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);
    };
 
    auto iterate_meshsets = [&](auto &&meshset_vec_ptr) {
      MoFEMFunctionBegin;
      for (auto m : meshset_vec_ptr) {
        auto bc = boost::make_shared<BCs>();
        CHKERR m_field.get_moab().get_entities_by_handle(m->getMeshset(),
                                                         bc->bcEnts, true);
        CHKERR m->getAttributes(bc->bcAttributes);
 
        if (problem_name.size())
          CHKERR mark_fix_dofs(bc->bcMarkers, lo, hi);
        MOFEM_LOG("BcMngWorld", Sev::verbose)
            << "Found block " << m->getName() << " number of attributes "
            << bc->bcAttributes.size();
 
        if (get_low_dim_ents) {
          auto low_dim_ents =
              BcManagerImplTools::get_adj_ents(m_field.get_moab(), bc->bcEnts);
          bc->bcEnts.swap(low_dim_ents);
        }
 
        CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(
            bc->bcEnts);
        if (problem_name.size())
          CHKERR prb_mng->markDofs(problem_name, ROW, ProblemsManager::AND,
                                   bc->bcEnts, bc->bcMarkers);
 
        const std::string bc_id =
            problem_name + "_" + field_name + "_" + m->getName();
        bcMapByBlockName[bc_id] = bc;
      }
      MoFEMFunctionReturn(0);
    };
 
    CHKERR iterate_meshsets(
 
        m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(std::regex(
 
            (boost::format("%s(.*)") % block_name).str()
 
                ))
 
    );
 
    MoFEMFunctionReturn(0);
  };
 
  CHKERR fix_disp();
 
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [4/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [5/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [6/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [7/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [8/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [9/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [10/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [11/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [12/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [13/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

pushMarkDOFsOnEntities() [14/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 753 of file BcManager.cpp.

                                                         {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  // if(problem_name.size()==0)
  //   MOFEM_LOG("BcMngWorld", Sev::warning)
  //       << "Argument problem_name has no effect";
 
  if (block_name_field_prefix)
    MOFEM_LOG("BcMngWorld", Sev::warning)
        << "Argument block_name_field_prefix=true has no effect";
 
  auto fix_disp = [&]() {
    MoFEMFunctionBegin;
 
    auto iterate_meshsets = [&](auto &&meshset_vec_ptr) {
      MoFEMFunctionBegin;
      for (auto m : meshset_vec_ptr) {
        auto bc = boost::make_shared<BCs>();
        CHKERR m_field.get_moab().get_entities_by_handle(m->getMeshset(),
                                                         bc->bcEnts, true);
        bc->dispBcPtr = boost::make_shared<DisplacementCubitBcData>();
        CHKERR m->getBcDataStructure(*(bc->dispBcPtr));
 
        MOFEM_LOG("BcMngWorld", Sev::verbose)
            << "Found block DISPLACEMENTSET id = " << m->getMeshsetId();
        MOFEM_LOG("BcMngWorld", Sev::verbose) << *bc->dispBcPtr;
 
        MOFEM_LOG("BcMngSync", Sev::noisy)
            << "Found block DISPLACEMENTSET id = " << m->getMeshsetId()
            << " nb. of entities " << bc->bcEnts.size()
            << " highest dim of entities "
            << BcManagerImplTools::get_dim(bc->bcEnts);
        MOFEM_LOG("BcMngSync", Sev::noisy) << *bc->dispBcPtr;
        MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::noisy);
 
        if (problem_name.size()) {
 
          if (bc->dispBcPtr->data.flag1)
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 0, 0,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          if (bc->dispBcPtr->data.flag2)
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 1, 1,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          if (bc->dispBcPtr->data.flag3)
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 2, 2,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          if (bc->dispBcPtr->data.flag4) {
 
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 1, 1,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 2, 2,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          }
          if (bc->dispBcPtr->data.flag5) {
 
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 0, 0,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 2, 2,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          }
          if (bc->dispBcPtr->data.flag6) {
 
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 0, 0,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 1, 1,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          }
        }
 
        if (get_low_dim_ents) {
          auto low_dim_ents =
              BcManagerImplTools::get_adj_ents(m_field.get_moab(), bc->bcEnts);
          bc->bcEnts.swap(low_dim_ents);
        }
 
        CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(
            bc->bcEnts);
        if (problem_name.size())
          CHKERR prb_mng->markDofs(problem_name, ROW, ProblemsManager::AND,
                                   bc->bcEnts, bc->bcMarkers);
 
        const std::string bc_id =
            problem_name + "_" + field_name + "_DISPLACEMENTSET" +
            boost::lexical_cast<std::string>(m->getMeshsetId());
        bcMapByBlockName[bc_id] = bc;
      }
      MoFEMFunctionReturn(0);
    };
 
    CHKERR iterate_meshsets(
 
        m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(
            NODESET | DISPLACEMENTSET)
 
    );
 
    MoFEMFunctionReturn(0);
  };
 
  CHKERR fix_disp();
 
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [15/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 871 of file BcManager.cpp.

                                                         {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  if (block_name_field_prefix)
    MOFEM_LOG("BcMngWorld", Sev::warning)
        << "Argument block_name_field_prefix=true has no effect";
 
  auto fix_temp = [&]() {
    MoFEMFunctionBegin;
 
    auto iterate_meshsets = [&](auto &&meshset_vec_ptr) {
      MoFEMFunctionBegin;
      for (auto m : meshset_vec_ptr) {
        auto bc = boost::make_shared<BCs>();
        CHKERR m_field.get_moab().get_entities_by_handle(m->getMeshset(),
                                                         bc->bcEnts, true);
        bc->tempBcPtr = boost::make_shared<TemperatureCubitBcData>();
        CHKERR m->getBcDataStructure(*(bc->tempBcPtr));
 
        MOFEM_LOG("BcMngWorld", Sev::verbose)
            << "Found block TEMPERATURESET id = " << m->getMeshsetId();
        MOFEM_LOG("BcMngWorld", Sev::verbose) << *bc->tempBcPtr;
 
        CHKERR prb_mng->modifyMarkDofs(
            problem_name, ROW, field_name, 0, MAX_DOFS_ON_ENTITY,
            ProblemsManager::MarkOP::OR, 1, bc->bcMarkers);
 
        if (get_low_dim_ents) {
          auto low_dim_ents =
              BcManagerImplTools::get_adj_ents(m_field.get_moab(), bc->bcEnts);
          bc->bcEnts.swap(low_dim_ents);
        }
 
        CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(
            bc->bcEnts);
        CHKERR prb_mng->markDofs(problem_name, ROW, ProblemsManager::AND,
                                 bc->bcEnts, bc->bcMarkers);
 
        const std::string bc_id =
            problem_name + "_" + field_name + "_TEMPERATURESET" +
            boost::lexical_cast<std::string>(m->getMeshsetId());
        bcMapByBlockName[bc_id] = bc;
      }
      MoFEMFunctionReturn(0);
    };
 
    CHKERR iterate_meshsets(
 
        m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(
            NODESET | TEMPERATURESET)
 
    );
 
    MoFEMFunctionReturn(0);
  };
 
  CHKERR fix_temp();
 
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [16/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 937 of file BcManager.cpp.

                                                         {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  if (block_name_field_prefix)
    MOFEM_LOG("BcMngWorld", Sev::warning)
        << "Argument block_name_field_prefix=true has no effect";
 
  auto fix_disp = [&]() {
    MoFEMFunctionBegin;
 
    auto iterate_meshsets = [&](auto &&meshset_vec_ptr) {
      MoFEMFunctionBegin;
      for (auto m : meshset_vec_ptr) {
        auto bc = boost::make_shared<BCs>();
        CHKERR m_field.get_moab().get_entities_by_handle(m->getMeshset(),
                                                         bc->bcEnts, true);
        bc->heatFluxBcPtr = boost::make_shared<HeatFluxCubitBcData>();
        CHKERR m->getBcDataStructure(*(bc->heatFluxBcPtr));
 
        CHKERR prb_mng->modifyMarkDofs(
            problem_name, ROW, field_name, 0, MAX_DOFS_ON_ENTITY,
            ProblemsManager::MarkOP::OR, 1, bc->bcMarkers);
 
        MOFEM_LOG("BcMngWorld", Sev::inform)
            << "Found block HEATFLUX id = " << m->getMeshsetId();
        MOFEM_LOG("BcMngWorld", Sev::inform) << *bc->heatFluxBcPtr;
 
        if (get_low_dim_ents) {
          auto low_dim_ents =
              BcManagerImplTools::get_adj_ents(m_field.get_moab(), bc->bcEnts);
          bc->bcEnts.swap(low_dim_ents);
        }
 
        CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(
            bc->bcEnts);
        CHKERR prb_mng->markDofs(problem_name, ROW, ProblemsManager::AND,
                                 bc->bcEnts, bc->bcMarkers);
 
        const std::string bc_id =
            problem_name + "_" + field_name + "_HEATFLUXSET" +
            boost::lexical_cast<std::string>(m->getMeshsetId());
        bcMapByBlockName[bc_id] = bc;
      }
      MoFEMFunctionReturn(0);
    };
 
    CHKERR iterate_meshsets(
 
        m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(SIDESET |
                                                                    HEATFLUXSET)
 
    );
 
    MoFEMFunctionReturn(0);
  };
 
  CHKERR fix_disp();
 
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [17/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 1004 of file BcManager.cpp.

                                                         {
  MoFEMFunctionBegin;
 
  auto mark_dofs = [&](const string block_name, const int &idx_0,
                       const int &idx_1) {
    MoFEMFunctionBeginHot;
    if (block_name_field_prefix) {
      const string field_block = field_name + "_" + block_name;
      CHKERR pushMarkDOFsOnEntities(problem_name, field_block, field_name,
                                    idx_0, idx_1, get_low_dim_ents);
    } else {
 
      CHKERR pushMarkDOFsOnEntities(problem_name, block_name, field_name, idx_0,
                                    idx_1, get_low_dim_ents);
    }
    MoFEMFunctionReturnHot(0);
  };
 
  // displacement
  CHKERR mark_dofs("FIX_X", 0, 0);
  CHKERR mark_dofs("FIX_Y", 1, 1);
  CHKERR mark_dofs("FIX_Z", 2, 2);
  CHKERR mark_dofs("FIX_ALL", 0, MAX_DOFS_ON_ENTITY);
 
  // rotation
  CHKERR mark_dofs("ROTATE_X", 1, 1);
  CHKERR mark_dofs("ROTATE_X", 2, 2);
  CHKERR mark_dofs("ROTATE_Y", 0, 0);
  CHKERR mark_dofs("ROTATE_Y", 2, 2);
  CHKERR mark_dofs("ROTATE_Z", 0, 0);
  CHKERR mark_dofs("ROTATE_Z", 1, 1);
  CHKERR mark_dofs("ROTATE_ALL", 0, MAX_DOFS_ON_ENTITY);
 
  std::string regex_str;
  if (block_name_field_prefix) {
    regex_str = (boost::format("%s_%s_%s_(.*)") % problem_name % field_name %
                 field_name)
                    .str();
  } else {
    regex_str = (boost::format("%s_%s_(.*)") % problem_name % field_name).str();
  }
 
  for (auto &m : bcMapByBlockName) {
    auto &bc_id = m.first;
    if (std::regex_match(bc_id, std::regex(regex_str))) {
      auto &bc = m.second;
      if (std::regex_match(bc_id, std::regex("(.*)_FIX_X(.*)"))) {
        bc->dispBcPtr = boost::make_shared<DisplacementCubitBcData>();
        bc->dispBcPtr->data.flag1 = 1;
        if (bc->bcAttributes.empty()) {
          bc->dispBcPtr->data.value1 = 0;
          MOFEM_LOG("BcMngWorld", Sev::warning)
              << "Expected one attribute on block but have "
              << bc->bcAttributes.size();
        } else if (bc->bcAttributes.size() >= 1) {
          bc->dispBcPtr->data.value1 = bc->bcAttributes[0];
        }
        MOFEM_LOG("BcMngWorld", Sev::inform) << "Add X " << bc_id;
        MOFEM_LOG("BcMngWorld", Sev::inform) << *bc->dispBcPtr;
      } else if (std::regex_match(bc_id, std::regex("(.*)_FIX_Y(.*)"))) {
        bc->dispBcPtr = boost::make_shared<DisplacementCubitBcData>();
        bc->dispBcPtr->data.flag2 = 1;
        if (bc->bcAttributes.empty()) {
          bc->dispBcPtr->data.value2 = 0;
          MOFEM_LOG("BcMngWorld", Sev::warning)
              << "Expected one attribute on block but have "
              << bc->bcAttributes.size();
        } else if (bc->bcAttributes.size() == 1) {
          bc->dispBcPtr->data.value2 = bc->bcAttributes[0];
        } else if (bc->bcAttributes.size() >= 2) {
          bc->dispBcPtr->data.value2 = bc->bcAttributes[1];
        }
        MOFEM_LOG("BcMngWorld", Sev::inform) << "Add Y " << bc_id;
        MOFEM_LOG("BcMngWorld", Sev::inform) << *(bc->dispBcPtr);
      } else if (std::regex_match(bc_id, std::regex("(.*)_FIX_Z(.*)"))) {
        bc->dispBcPtr = boost::make_shared<DisplacementCubitBcData>();
        bc->dispBcPtr->data.flag3 = 1;
        if (bc->bcAttributes.empty()) {
          bc->dispBcPtr->data.value3 = 0;
          MOFEM_LOG("BcMngWorld", Sev::warning)
              << "Expected one attribute on block but have "
              << bc->bcAttributes.size();
        } else if (bc->bcAttributes.size() == 1) {
          bc->dispBcPtr->data.value3 = bc->bcAttributes[0];
        } else if (bc->bcAttributes.size() == 3) {
          bc->dispBcPtr->data.value3 = bc->bcAttributes[2];
        }
        MOFEM_LOG("BcMngWorld", Sev::inform) << "Add Z " << bc_id;
        MOFEM_LOG("BcMngWorld", Sev::inform) << *(bc->dispBcPtr);
      } else if (std::regex_match(bc_id, std::regex("(.*)_FIX_ALL(.*)"))) {
        bc->dispBcPtr = boost::make_shared<DisplacementCubitBcData>();
        bc->dispBcPtr->data.flag1 = 1;
        bc->dispBcPtr->data.flag2 = 1;
        bc->dispBcPtr->data.flag3 = 1;
        if (bc->bcAttributes.size() >= 1) {
          bc->dispBcPtr->data.value1 = bc->bcAttributes[0];
        }
        if (bc->bcAttributes.size() >= 2) {
          bc->dispBcPtr->data.value2 = bc->bcAttributes[1];
        }
        if (bc->bcAttributes.size() >= 3) {
          bc->dispBcPtr->data.value3 = bc->bcAttributes[2];
        }
        MOFEM_LOG("BcMngWorld", Sev::inform) << "Add ALL " << bc_id;
        MOFEM_LOG("BcMngWorld", Sev::inform) << *(bc->dispBcPtr);
      } else if (std::regex_match(bc_id, std::regex("(.*)_ROTATE_X(.*)"))) {
        bc->dispBcPtr =
            boost::make_shared<DisplacementCubitBcDataWithRotation>();
        bc->dispBcPtr->data.flag4 = 1;
        bc->dispBcPtr->data.flag5 = 0;
        bc->dispBcPtr->data.flag6 = 0;
        // for the ROTATE_X block the angles can be specified with either one or
        // three attributes, e.g. 1, coords or 1,0,0,coords
        if (bc->bcAttributes.empty()) {
          bc->dispBcPtr->data.value4 = 0;
          MOFEM_LOG("BcMngWorld", Sev::warning)
              << "Expected one attribute on block on block (angle (1 or 3), "
                 "center coords(3) but have "
              << bc->bcAttributes.size();
        } else if (bc->bcAttributes.size() >= 1) {
          bc->dispBcPtr->data.value4 = bc->bcAttributes[0];
        }
        MOFEM_LOG("BcMngWorld", Sev::inform) << "Add X " << bc_id;
        MOFEM_LOG("BcMngWorld", Sev::inform) << *bc->dispBcPtr;
        if (bc->bcAttributes.size() == 4 || bc->bcAttributes.size() == 6) {
          if (auto ext_disp_bc =
                  dynamic_cast<DisplacementCubitBcDataWithRotation *>(
                      bc->dispBcPtr.get())) {
            auto &o = ext_disp_bc->rotOffset;
            for (int a = 0; a != 3; ++a)
              o[a] = bc->bcAttributes[bc->bcAttributes.size() - 3 + a];
            MOFEM_LOG("BcMngWorld", Sev::inform)
                << "Add Rotate X Center: " << o[0] << " " << o[1] << " "
                << o[2];
          }
        }
      } else if (std::regex_match(bc_id, std::regex("(.*)_ROTATE_Y(.*)"))) {
        bc->dispBcPtr =
            boost::make_shared<DisplacementCubitBcDataWithRotation>();
        bc->dispBcPtr->data.flag4 = 0;
        bc->dispBcPtr->data.flag5 = 1;
        bc->dispBcPtr->data.flag6 = 0;
        // for the ROTATE_Y block the angles can be specified with either one or
        // three attributes, e.g. 1, coords or 0,1,0,coords
        if (bc->bcAttributes.empty()) {
          bc->dispBcPtr->data.value5 = 0;
          MOFEM_LOG("BcMngWorld", Sev::warning)
              << "Expected one attribute on block on block (angle (1 or 3), "
                 "center coords(3) but have "
              << bc->bcAttributes.size();
        } else if (bc->bcAttributes.size() == 1 ||
                   bc->bcAttributes.size() == 4) {
          bc->dispBcPtr->data.value5 = bc->bcAttributes[0];
        } else if (bc->bcAttributes.size() == 6) {
          bc->dispBcPtr->data.value5 = bc->bcAttributes[1];
        }
        MOFEM_LOG("BcMngWorld", Sev::inform) << "Add Y " << bc_id;
        MOFEM_LOG("BcMngWorld", Sev::inform) << *(bc->dispBcPtr);
        if (bc->bcAttributes.size() == 4 || bc->bcAttributes.size() == 6) {
          if (auto ext_disp_bc =
                  dynamic_cast<DisplacementCubitBcDataWithRotation *>(
                      bc->dispBcPtr.get())) {
            auto &o = ext_disp_bc->rotOffset;
            for (int a = 0; a != 3; ++a)
              o[a] = bc->bcAttributes[bc->bcAttributes.size() - 3 + a];
            MOFEM_LOG("BcMngWorld", Sev::inform)
                << "Add Rotate Y Center: " << o[0] << " " << o[1] << " "
                << o[2];
          }
        }
      } else if (std::regex_match(bc_id, std::regex("(.*)_ROTATE_Z(.*)"))) {
        bc->dispBcPtr =
            boost::make_shared<DisplacementCubitBcDataWithRotation>();
        bc->dispBcPtr->data.flag4 = 0;
        bc->dispBcPtr->data.flag5 = 0;
        bc->dispBcPtr->data.flag6 = 1;
        // for the ROTATE_Z block the angles can be specified with either one or
        // three attributes, e.g. 1, coords or 0,0,1,coords
        if (bc->bcAttributes.empty()) {
          bc->dispBcPtr->data.value6 = 0;
          MOFEM_LOG("BcMngWorld", Sev::warning)
              << "Expected one attribute on block (angle (1 or 3), center "
                 "coords(3) but have "
              << bc->bcAttributes.size();
        } else if (bc->bcAttributes.size() == 1 ||
                   bc->bcAttributes.size() == 4) {
          bc->dispBcPtr->data.value6 = bc->bcAttributes[0];
        } else if (bc->bcAttributes.size() == 3 ||
                   bc->bcAttributes.size() == 6) {
          bc->dispBcPtr->data.value6 = bc->bcAttributes[2];
        }
        MOFEM_LOG("BcMngWorld", Sev::inform) << "Add Z " << bc_id;
        MOFEM_LOG("BcMngWorld", Sev::inform) << *(bc->dispBcPtr);
        if (bc->bcAttributes.size() == 4 || bc->bcAttributes.size() == 6) {
          if (auto ext_disp_bc =
                  dynamic_cast<DisplacementCubitBcDataWithRotation *>(
                      bc->dispBcPtr.get())) {
            auto &o = ext_disp_bc->rotOffset;
            for (int a = 0; a != 3; ++a)
              o[a] = bc->bcAttributes[bc->bcAttributes.size() - 3 + a];
            MOFEM_LOG("BcMngWorld", Sev::inform)
                << "Add Rotate Z Center: " << o[0] << " " << o[1] << " "
                << o[2];
          }
        }
      } else if (std::regex_match(bc_id, std::regex("(.*)_ROTATE_ALL(.*)"))) {
        bc->dispBcPtr =
            boost::make_shared<DisplacementCubitBcDataWithRotation>();
        bc->dispBcPtr->data.flag4 = 1;
        bc->dispBcPtr->data.flag5 = 1;
        bc->dispBcPtr->data.flag6 = 1;
        if (bc->bcAttributes.size() >= 1) {
          bc->dispBcPtr->data.value4 = bc->bcAttributes[0];
        }
        if (bc->bcAttributes.size() >= 2) {
          bc->dispBcPtr->data.value5 = bc->bcAttributes[1];
        }
        if (bc->bcAttributes.size() >= 3) {
          bc->dispBcPtr->data.value6 = bc->bcAttributes[2];
        }
        MOFEM_LOG("BcMngWorld", Sev::inform) << "Add ALL " << bc_id;
        MOFEM_LOG("BcMngWorld", Sev::inform) << *(bc->dispBcPtr);
        if (bc->bcAttributes.size() > 3) {
          if (auto ext_disp_bc =
                  dynamic_cast<DisplacementCubitBcDataWithRotation *>(
                      bc->dispBcPtr.get())) {
            auto &o = ext_disp_bc->rotOffset;
            for (int a = 0; a != 3; ++a)
              o[a] = bc->bcAttributes[3 + a];
            MOFEM_LOG("BcMngWorld", Sev::inform)
                << "Add Rotate ALL Center: " << o[0] << " " << o[1] << " "
                << o[2];
          }
        }
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [18/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 1284 of file BcManager.cpp.

                                                         {
  MoFEMFunctionBegin;
  // that marks DOFs and create data when are set by cubit nodesets.
  CHKERR pushMarkDOFsOnEntities<BcMeshsetType<DISPLACEMENTSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
  // that marks DOFs and create data when are set by blocsket.
  CHKERR pushMarkDOFsOnEntities<BcDisplacementMeshsetType<BLOCKSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [19/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 1316 of file BcManager.cpp.

                                                         {
  MoFEMFunctionBegin;
  CHKERR pushMarkDOFsOnEntities<BcMeshsetType<TEMPERATURESET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
 
  auto get_block_name = [&]() {
    if (block_name_field_prefix)
      return (boost::format("%s_FIX_SCALAR") % field_name).str();
    else
      return field_name;
  };
 
  CHKERR pushMarkDOFsOnEntities<BcScalarMeshsetType<BLOCKSET>>(
      problem_name, get_block_name(), field_name, get_low_dim_ents);
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [20/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 1359 of file BcManager.cpp.

                                                         {
  MoFEMFunctionBegin;
  CHKERR pushMarkDOFsOnEntities<BcMeshsetType<HEATFLUXSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [21/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 1405 of file BcManager.cpp.

                                                         {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  if (problem_name.empty())
    MOFEM_LOG("BcMngWorld", Sev::warning) << "Argument problem_name is empty";
 
  if (block_name_field_prefix)
    MOFEM_LOG("BcMngWorld", Sev::warning)
        << "Argument block_name_field_prefix=true has no effect";
 
  auto fix_force = [&]() {
    MoFEMFunctionBegin;
 
    auto iterate_meshsets = [&](auto &&meshset_vec_ptr) {
      MoFEMFunctionBegin;
      for (auto m : meshset_vec_ptr) {
        auto bc = boost::make_shared<BCs>();
        CHKERR m_field.get_moab().get_entities_by_handle(m->getMeshset(),
                                                         bc->bcEnts, true);
        bc->forceBcPtr = boost::make_shared<ForceCubitBcData>();
        CHKERR m->getBcDataStructure(*(bc->forceBcPtr));
 
        MOFEM_LOG("BcMngWorld", Sev::verbose)
            << "Found block FORCESET id = " << m->getMeshsetId();
        MOFEM_LOG("BcMngWorld", Sev::verbose) << *bc->forceBcPtr;
 
        MOFEM_LOG("BcMngSync", Sev::noisy)
            << "Found block FORCESET id = " << m->getMeshsetId()
            << " nb. of entities " << bc->bcEnts.size()
            << " highest dim of entities "
            << BcManagerImplTools::get_dim(bc->bcEnts);
        MOFEM_LOG("BcMngSync", Sev::noisy) << *bc->forceBcPtr;
        MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::noisy);
 
        if (problem_name.size()) {
 
          if (bc->forceBcPtr->data.value2 > 0)
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 0, 0,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          if (bc->forceBcPtr->data.value3 > 0)
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 1, 1,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          if (bc->forceBcPtr->data.value4 > 0)
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 2, 2,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
 
          if (bc->forceBcPtr->data.value5 > 0) {
 
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 1, 1,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 2, 2,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          }
          if (bc->forceBcPtr->data.value5) {
 
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 0, 0,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 2, 2,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          }
          if (bc->forceBcPtr->data.value6) {
 
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 0, 0,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
            CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 1, 1,
                                           ProblemsManager::MarkOP::OR, 1,
                                           bc->bcMarkers);
          }
        }
 
        if (get_low_dim_ents) {
          auto low_dim_ents =
              BcManagerImplTools::get_adj_ents(m_field.get_moab(), bc->bcEnts);
          bc->bcEnts.swap(low_dim_ents);
        }
 
        CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(
            bc->bcEnts);
        if (problem_name.size())
          CHKERR prb_mng->markDofs(problem_name, ROW, ProblemsManager::AND,
                                   bc->bcEnts, bc->bcMarkers);
 
        const std::string bc_id =
            problem_name + "_" + field_name + "_FORCESET" +
            boost::lexical_cast<std::string>(m->getMeshsetId());
        bcMapByBlockName[bc_id] = bc;
      }
      MoFEMFunctionReturn(0);
    };
 
    CHKERR iterate_meshsets(
 
        m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(NODESET |
                                                                    FORCESET)
 
    );
 
    MoFEMFunctionReturn(0);
  };
 
  CHKERR fix_force();
 
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [22/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 1523 of file BcManager.cpp.

                                                         {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  if (problem_name.size() == 0)
    MOFEM_LOG("BcMngWorld", Sev::warning) << "Argument problem_name is empty";
 
  auto get_force_block = [&](auto block_name) {
    MoFEMFunctionBegin;
 
    for (auto m :
         m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(std::regex(
 
             (boost::format("%s(.*)") % block_name).str()
 
                 ))
 
    ) {
 
      const auto block_name = m->getName();
 
      MOFEM_LOG("BcMngWorld", Sev::inform)
          << "Found force block " << block_name;
 
      auto bc = boost::make_shared<BCs>();
      CHKERR m_field.get_moab().get_entities_by_handle(m->getMeshset(),
                                                       bc->bcEnts, true);
 
      CHKERR m->getAttributes(bc->bcAttributes);
      if (bc->bcAttributes.size() != 3) {
        SETERRQ(m_field.get_comm(), MOFEM_DATA_INCONSISTENCY,
                "Expect three block attributes for force block");
      }
 
      bc->forceBcPtr = boost::make_shared<ForceCubitBcData>();
      // For details look at ForceCubitBcData in
      // mofem/src/multi_indices/BCData.hpp
      bc->forceBcPtr->data.value1 = 1;
      bc->forceBcPtr->data.value3 = bc->bcAttributes[0];
      bc->forceBcPtr->data.value4 = bc->bcAttributes[1];
      bc->forceBcPtr->data.value5 = bc->bcAttributes[2];
 
      MOFEM_LOG("BcMngWorld", Sev::inform) << *bc->forceBcPtr;
      MOFEM_LOG("BcMngSync", Sev::noisy)
          << "Found block FORCESET id = " << m->getMeshsetId()
          << " nb. of entities " << bc->bcEnts.size()
          << " highest dim of entities "
          << BcManagerImplTools::get_dim(bc->bcEnts);
      MOFEM_LOG("BcMngSync", Sev::noisy) << *bc->forceBcPtr;
      MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::noisy);
 
      if (problem_name.size()) {
 
        if (bc->forceBcPtr->data.value2 > 0)
          CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 0, 0,
                                         ProblemsManager::MarkOP::OR, 1,
                                         bc->bcMarkers);
        if (bc->forceBcPtr->data.value3 > 0)
          CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 1, 1,
                                         ProblemsManager::MarkOP::OR, 1,
                                         bc->bcMarkers);
        if (bc->forceBcPtr->data.value4 > 0)
          CHKERR prb_mng->modifyMarkDofs(problem_name, ROW, field_name, 2, 2,
                                         ProblemsManager::MarkOP::OR, 1,
                                         bc->bcMarkers);
      }
 
      if (get_low_dim_ents) {
        auto low_dim_ents =
            BcManagerImplTools::get_adj_ents(m_field.get_moab(), bc->bcEnts);
        bc->bcEnts.swap(low_dim_ents);
      }
 
      CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(
          bc->bcEnts);
      if (problem_name.size())
        CHKERR prb_mng->markDofs(problem_name, ROW, ProblemsManager::AND,
                                 bc->bcEnts, bc->bcMarkers);
 
      const std::string bc_id =
          problem_name + "_" + field_name + "_" + block_name;
      bcMapByBlockName[bc_id] = bc;
    }
    MoFEMFunctionReturn(0);
  };
 
  CHKERR get_force_block("FORCE");
 
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [23/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix
	)

Definition at line 1618 of file BcManager.cpp.

                                                         {
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities<BcMeshsetType<FORCESET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
 
  CHKERR pushMarkDOFsOnEntities<BcForceMeshsetType<BLOCKSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
 
  MoFEMFunctionReturn(0);
}

pushMarkDOFsOnEntities() [24/25]

template<typename T >

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents = true,
		bool	block_name_field_prefix = false
	)

Mark block DOFs.

Parameters

problem_name
field_name
get_low_dim_ents	get lower dimension entities
block_name_field_prefix

Returns

MoFEMErrorCode

pushMarkDOFsOnEntities() [25/25]

template<>

MoFEMErrorCode MoFEM::BcManager::pushMarkDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		const std::string	block_name,
		bool	get_low_dim_ents
	)

pushMarkSideDofs()

MoFEMErrorCode MoFEM::BcManager::pushMarkSideDofs	(	const std::string	problem_name,
		const std::string	block_name,
		const std::string	field_name,
		int	bridge_dim,
		int	lo,
		int	hi
	)

Mark side DOFs.

Parameters

problem_name
block_name
field_name
lo
hi

Returns

MoFEMErrorCode

Definition at line 1779 of file BcManager.cpp.

  {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  if (problem_name.empty())
    MoFEMFunctionReturnHot(0);
 
  auto iterate_meshsets = [&](auto &&meshset_vec_ptr) {
    auto prb_mng = m_field.getInterface<ProblemsManager>();
    MoFEMFunctionBegin;
    for (auto m : meshset_vec_ptr) {
      auto bc = boost::make_shared<BCs>();
      CHKERR m_field.get_moab().get_entities_by_handle(m->getMeshset(),
                                                       bc->bcEnts, true);
      CHKERR m->getAttributes(bc->bcAttributes);
 
      bc->dofsViewPtr = boost::make_shared<BCs::DofsView>();
 
      CHKERR m_field.getInterface<CommInterface>()->synchroniseEntities(
          bc->bcEnts);
      CHKERR prb_mng->getSideDofsOnBrokenSpaceEntities(
          *(bc->dofsViewPtr), problem_name, ROW, field_name, bc->bcEnts,
          bridge_dim, lo, hi);
      CHKERR prb_mng->markDofs(problem_name, ROW, *(bc->dofsViewPtr),
                               ProblemsManager::OR, bc->bcMarkers);
 
      MOFEM_LOG("BcMngWorld", Sev::inform)
          << "Found block " << m->getName() << " number of attributes "
          << bc->bcAttributes.size() << " number of entities "
          << bc->bcEnts.size();
 
      const std::string bc_id =
          problem_name + "_" + field_name + "_" + m->getName();
      bcMapByBlockName[bc_id] = bc;
    }
    MoFEMFunctionReturn(0);
  };
 
  CHKERR iterate_meshsets(
 
      m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(std::regex(
 
          (boost::format("%s(.*)") % block_name).str()
 
              ))
 
  );
 
  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 34 of file BcManager.cpp.

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

removeBlockDOFsOnEntities() [1/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	block_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [2/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	block_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	is_distributed_mesh
	)

Definition at line 697 of file BcManager.cpp.

                              {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities<BcScalarMeshsetType<BLOCKSET>>(
      problem_name, block_name, field_name, get_low_dim_ents);
 
  Range ents_to_remove;
 
  for (auto m :
 
       m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(
           BLOCKSET | UNKNOWNNAME)) {
 
    std::string bc_id = problem_name + "_" + field_name + "_" + m->getName();
 
    auto str = boost::format("%s_%s_%s(.*)")
 
               % problem_name % field_name % block_name;
 
    if (std::regex_match(bc_id, std::regex(str.str()))) {
 
      auto bc = bcMapByBlockName.at(bc_id);
 
      if (auto temp_bc = bc->tempBcPtr) {
        if (temp_bc->data.flag1) {
          ents_to_remove.merge(bc->bcEnts);
        }
      } else {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "BC type not implemented");
      }
    }
  }
 
  auto remove_dofs_on_ents = [&](const Range &ents, const int lo,
                                 const int hi) {
    if (is_distributed_mesh)
      return prb_mng->removeDofsOnEntities(problem_name, field_name, ents, lo,
                                           hi);
    else
      return prb_mng->removeDofsOnEntitiesNotDistributed(
          problem_name, field_name, ents, lo, hi);
  };
 
  CHKERR remove_dofs_on_ents(ents_to_remove, 0, MAX_DOFS_ON_ENTITY);
 
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [3/25]

template<typename T >

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	block_name,
		const std::string	field_name,
		bool	get_low_dim_ents = true,
		bool	is_distributed_mesh = true
	)

Mark block DOFs.

Template Parameters

BCSET

Parameters

problem_name
field_name
block_name
get_low_dim_ents
is_distributed_mesh

Returns

MoFEMErrorCode

removeBlockDOFsOnEntities() [4/25]

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,
		bool	is_distributed_mesh = 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
is_distributed_mesh	distributed mesh

Returns

MoFEMErrorCode

Definition at line 73 of file BcManager.cpp.

                              {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities(problem_name, block_name, field_name, lo, hi,
                                get_low_dim_ents);
 
  auto remove_dofs_on_ents = [&](const Range &ents, const int lo,
                                 const int hi) {
    if (is_distributed_mesh)
      return prb_mng->removeDofsOnEntities(problem_name, field_name, ents, lo,
                                           hi);
    else
      return prb_mng->removeDofsOnEntitiesNotDistributed(
          problem_name, field_name, ents, lo, hi);
  };
 
  for (auto m :
       m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(std::regex(
 
           (boost::format("%s(.*)") % block_name).str()
 
               ))
 
  ) {
    const std::string bc_id =
        problem_name + "_" + field_name + "_" + m->getName();
    CHKERR remove_dofs_on_ents(bcMapByBlockName.at(bc_id)->bcEnts, lo, hi);
    bcMapByBlockName.at(bc_id)->bcMarkers = std::vector<unsigned char>();
  }
 
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [5/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [6/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [7/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [8/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [9/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [10/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [11/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 455 of file BcManager.cpp.

                              {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities<BcMeshsetType<DISPLACEMENTSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
 
  std::array<Range, 3> ents_to_remove;
 
  for (auto m :
 
       m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(
           NODESET | DISPLACEMENTSET)) {
 
    const std::string bc_id =
        problem_name + "_" + field_name + "_DISPLACEMENTSET" +
        boost::lexical_cast<std::string>(m->getMeshsetId());
 
    auto bc = bcMapByBlockName.at(bc_id);
 
    if (bc->dispBcPtr) {
      if (bc->dispBcPtr->data.flag1) {
        ents_to_remove[0].merge(bc->bcEnts);
      }
      if (bc->dispBcPtr->data.flag2) {
        ents_to_remove[1].merge(bc->bcEnts);
      }
      if (bc->dispBcPtr->data.flag3) {
        ents_to_remove[2].merge(bc->bcEnts);
      }
      if (bc->dispBcPtr->data.flag4) {
        ents_to_remove[1].merge(bc->bcEnts);
        ents_to_remove[2].merge(bc->bcEnts);
      }
      if (bc->dispBcPtr->data.flag5) {
        ents_to_remove[0].merge(bc->bcEnts);
        ents_to_remove[2].merge(bc->bcEnts);
      }
      if (bc->dispBcPtr->data.flag6) {
        ents_to_remove[0].merge(bc->bcEnts);
        ents_to_remove[1].merge(bc->bcEnts);
      }
    }
    bc->bcMarkers = std::vector<unsigned char>();
  }
 
  auto remove_dofs_on_ents = [&](const Range &ents, const int lo,
                                 const int hi) {
    if (is_distributed_mesh)
      return prb_mng->removeDofsOnEntities(problem_name, field_name, ents, lo,
                                           hi);
    else
      return prb_mng->removeDofsOnEntitiesNotDistributed(
          problem_name, field_name, ents, lo, hi);
  };
 
  CHKERR remove_dofs_on_ents(ents_to_remove[0], 0, 0);
  CHKERR remove_dofs_on_ents(ents_to_remove[1], 1, 1);
  CHKERR remove_dofs_on_ents(ents_to_remove[2], 2, 2);
 
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [12/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [13/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [14/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [15/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

removeBlockDOFsOnEntities() [16/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 524 of file BcManager.cpp.

                              {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities<BcMeshsetType<TEMPERATURESET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
 
  Range ents_to_remove;
 
  for (auto m :
 
       m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(
           NODESET | TEMPERATURESET)
 
  ) {
    const std::string bc_id =
        problem_name + "_" + field_name + "_TEMPERATURESET" +
        boost::lexical_cast<std::string>(m->getMeshsetId());
    auto bc = bcMapByBlockName.at(bc_id);
    ents_to_remove.merge(bc->bcEnts);
    bc->bcMarkers = std::vector<unsigned char>();
  }
 
  auto remove_dofs_on_ents = [&](const Range &ents, const int lo,
                                 const int hi) {
    if (is_distributed_mesh)
      return prb_mng->removeDofsOnEntities(problem_name, field_name, ents, lo,
                                           hi);
    else
      return prb_mng->removeDofsOnEntitiesNotDistributed(
          problem_name, field_name, ents, lo, hi);
  };
 
  CHKERR remove_dofs_on_ents(ents_to_remove, 0, MAX_DOFS_ON_ENTITY);
 
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [17/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 567 of file BcManager.cpp.

                              {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities<BcMeshsetType<HEATFLUXSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
 
  Range ents_to_remove;
 
  for (auto m :
 
       m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(NODESET |
                                                                   HEATFLUXSET)
 
  ) {
    const std::string bc_id =
        problem_name + "_" + field_name + "_HEATFLUXSET" +
        boost::lexical_cast<std::string>(m->getMeshsetId());
    auto bc = bcMapByBlockName.at(bc_id);
    ents_to_remove.merge(bc->bcEnts);
    bc->bcMarkers = std::vector<unsigned char>();
  }
 
  auto remove_dofs_on_ents = [&](const Range &ents, const int lo,
                                 const int hi) {
    if (is_distributed_mesh)
      return prb_mng->removeDofsOnEntities(problem_name, field_name, ents, lo,
                                           hi);
    else
      return prb_mng->removeDofsOnEntitiesNotDistributed(
          problem_name, field_name, ents, lo, hi);
  };
 
  CHKERR remove_dofs_on_ents(ents_to_remove, 0, MAX_DOFS_ON_ENTITY);
 
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [18/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 611 of file BcManager.cpp.

                              {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities<BcDisplacementMeshsetType<BLOCKSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix);
 
  std::array<Range, 3> ents_to_remove;
 
  for (auto m :
 
       m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(
           BLOCKSET | UNKNOWNNAME)) {
 
    const auto block_name = m->getName();
 
    std::string bc_id = problem_name + "_" + field_name + "_" + block_name;
    std::string regex_str;
    if (block_name_field_prefix) {
      regex_str = (boost::format("%s_%s_%s_((FIX_(ALL|X|Y|Z))|("
                                 "DISPLACEMENT|ROTATE))(.*)") %
                   problem_name % field_name % field_name)
                      .str();
    } else {
      regex_str = (boost::format("%s_%s_((FIX_(ALL|X|Y|Z))|("
                                 "DISPLACEMENT|ROTATE))(.*)") %
                   problem_name % field_name)
                      .str();
    }
 
    if (std::regex_match(bc_id, std::regex(regex_str))) {
 
      auto bc = bcMapByBlockName.at(bc_id);
 
      if (auto disp_bc = bc->dispBcPtr) {
        if (disp_bc->data.flag1) {
          ents_to_remove[0].merge(bc->bcEnts);
        }
        if (disp_bc->data.flag2) {
          ents_to_remove[1].merge(bc->bcEnts);
        }
        if (disp_bc->data.flag3) {
          ents_to_remove[2].merge(bc->bcEnts);
        }
        if (disp_bc->data.flag4) {
          ents_to_remove[1].merge(bc->bcEnts);
          ents_to_remove[2].merge(bc->bcEnts);
        }
        if (disp_bc->data.flag5) {
          ents_to_remove[0].merge(bc->bcEnts);
          ents_to_remove[2].merge(bc->bcEnts);
        }
        if (disp_bc->data.flag6) {
          ents_to_remove[0].merge(bc->bcEnts);
          ents_to_remove[1].merge(bc->bcEnts);
        }
      } else {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "BC type not implemented");
      }
    }
  }
 
  auto remove_dofs_on_ents = [&](const Range &ents, const int lo,
                                 const int hi) {
    if (is_distributed_mesh)
      return prb_mng->removeDofsOnEntities(problem_name, field_name, ents, lo,
                                           hi);
    else
      return prb_mng->removeDofsOnEntitiesNotDistributed(
          problem_name, field_name, ents, lo, hi);
  };
 
  CHKERR remove_dofs_on_ents(ents_to_remove[0], 0, 0);
  CHKERR remove_dofs_on_ents(ents_to_remove[1], 1, 1);
  CHKERR remove_dofs_on_ents(ents_to_remove[2], 2, 2);
 
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [19/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 1298 of file BcManager.cpp.

                              {
  MoFEMFunctionBegin;
  // that remove DOFs when are set by cubit nodesets.
  CHKERR removeBlockDOFsOnEntities<BcMeshsetType<DISPLACEMENTSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix,
      is_distributed_mesh);
  // that remove DOFs when are by blocksets
  CHKERR removeBlockDOFsOnEntities<BcDisplacementMeshsetType<BLOCKSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix,
      is_distributed_mesh);
  // add more ways to remove bcs when appropiate
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [20/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 1336 of file BcManager.cpp.

                              {
  MoFEMFunctionBegin;
  CHKERR removeBlockDOFsOnEntities<BcMeshsetType<TEMPERATURESET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix,
      is_distributed_mesh);
 
  auto get_block_name = [&]() {
    if (block_name_field_prefix)
      return (boost::format("%s_FIX_SCALAR") % field_name).str();
    else
      return field_name;
  };
 
  CHKERR removeBlockDOFsOnEntities<BcScalarMeshsetType<BLOCKSET>>(
      problem_name, get_block_name(), field_name, get_low_dim_ents,
      is_distributed_mesh);
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [21/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 1369 of file BcManager.cpp.

                              {
  MoFEMFunctionBegin;
  CHKERR removeBlockDOFsOnEntities<BcMeshsetType<HEATFLUXSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix,
      is_distributed_mesh);
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [22/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 1633 of file BcManager.cpp.

                              {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities<BcMeshsetType<FORCESET>>(
      problem_name, field_name, get_low_dim_ents);
 
  std::array<Range, 3> ents_to_remove;
 
  for (auto m :
 
       m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(NODESET |
                                                                   FORCESET)) {
 
    const auto block_name = m->getName();
    std::string bc_id = problem_name + "_" + field_name + "_" + block_name;
 
    auto str = boost::format("%s_%s_%s(.*)")
 
               % problem_name % field_name % block_name;
 
    if (std::regex_match(bc_id, std::regex(str.str()))) {
 
      auto bc = bcMapByBlockName.at(bc_id);
 
      if (auto force_bc = bc->forceBcPtr) {
        if (force_bc->data.value3 > 0) {
          ents_to_remove[0].merge(bc->bcEnts);
        }
        if (force_bc->data.value4 > 0) {
          ents_to_remove[1].merge(bc->bcEnts);
        }
        if (force_bc->data.value5 > 0) {
          ents_to_remove[2].merge(bc->bcEnts);
        }
      } else {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "BC type not implemented");
      }
    }
  }
 
  auto remove_dofs_on_ents = [&](const Range &ents, const int lo,
                                 const int hi) {
    if (is_distributed_mesh)
      return prb_mng->removeDofsOnEntities(problem_name, field_name, ents, lo,
                                           hi);
    else
      return prb_mng->removeDofsOnEntitiesNotDistributed(
          problem_name, field_name, ents, lo, hi);
  };
 
  CHKERR remove_dofs_on_ents(ents_to_remove[0], 0, 0);
  CHKERR remove_dofs_on_ents(ents_to_remove[1], 1, 1);
  CHKERR remove_dofs_on_ents(ents_to_remove[2], 2, 2);
 
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [23/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 1698 of file BcManager.cpp.

                              {
  Interface &m_field = cOre;
  auto prb_mng = m_field.getInterface<ProblemsManager>();
  MoFEMFunctionBegin;
 
  CHKERR pushMarkDOFsOnEntities<BcForceMeshsetType<BLOCKSET>>(
      problem_name, field_name, get_low_dim_ents);
 
  std::array<Range, 3> ents_to_remove;
 
  for (auto m :
 
       m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(
           BLOCKSET | UNKNOWNNAME)) {
 
    const auto block_name = m->getName();
    std::string bc_id = problem_name + "_" + field_name + "_" + block_name;
 
    auto str = boost::format("%s_%s_%s(.*)")
 
               % problem_name % field_name % block_name;
 
    if (std::regex_match(bc_id, std::regex(str.str()))) {
 
      auto bc = bcMapByBlockName.at(bc_id);
 
      if (auto force_bc = bc->forceBcPtr) {
        if (force_bc->data.value3 > 0) {
          ents_to_remove[0].merge(bc->bcEnts);
        }
        if (force_bc->data.value4 > 0) {
          ents_to_remove[1].merge(bc->bcEnts);
        }
        if (force_bc->data.value5 > 0) {
          ents_to_remove[2].merge(bc->bcEnts);
        }
      } else {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "BC type not implemented");
      }
    }
  }
 
  auto remove_dofs_on_ents = [&](const Range &ents, const int lo,
                                 const int hi) {
    if (is_distributed_mesh)
      return prb_mng->removeDofsOnEntities(problem_name, field_name, ents, lo,
                                           hi);
    else
      return prb_mng->removeDofsOnEntitiesNotDistributed(
          problem_name, field_name, ents, lo, hi);
  };
 
  CHKERR remove_dofs_on_ents(ents_to_remove[0], 0, 0);
  CHKERR remove_dofs_on_ents(ents_to_remove[1], 1, 1);
  CHKERR remove_dofs_on_ents(ents_to_remove[2], 2, 2);
 
  MoFEMFunctionReturn(0);
}

removeBlockDOFsOnEntities() [24/25]

template<>

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents,
		bool	block_name_field_prefix,
		bool	is_distributed_mesh
	)

Definition at line 1762 of file BcManager.cpp.

                              {
  MoFEMFunctionBegin;
 
  CHKERR removeBlockDOFsOnEntities<BcMeshsetType<FORCESET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix,
      is_distributed_mesh);
 
  CHKERR removeBlockDOFsOnEntities<BcForceMeshsetType<BLOCKSET>>(
      problem_name, field_name, get_low_dim_ents, block_name_field_prefix,
      is_distributed_mesh);
 
  MoFEMFunctionReturn(0);
};

removeBlockDOFsOnEntities() [25/25]

template<typename T >

MoFEMErrorCode MoFEM::BcManager::removeBlockDOFsOnEntities	(	const std::string	problem_name,
		const std::string	field_name,
		bool	get_low_dim_ents = true,
		bool	block_name_field_prefix = false,
		bool	is_distributed_mesh = true
	)

Mark block DOFs.

Template Parameters

BCSET

Parameters

problem_name
field_name
get_low_dim_ents
is_distributed_mesh
block_name_field_prefix

Returns

MoFEMErrorCode

removeSideDOFs()

MoFEMErrorCode MoFEM::BcManager::removeSideDOFs	(	const std::string	problem_name,
		const std::string	block_name,
		const std::string	field_name,
		int	bridge_dim,
		int	lo,
		int	hi,
		bool	is_distributed_mesh = true
	)

Remove side DOFs.

Parameters

problem_name
block_name
field_name
lo
hi
is_distributed_mesh

Returns

MoFEMErrorCode

Definition at line 1835 of file BcManager.cpp.

                                                                   {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  CHKERR pushMarkSideDofs(problem_name, block_name, field_name, bridge_dim, lo,
                          hi);
 
  auto iterate_meshsets = [&](auto &&meshset_vec_ptr) {
    MoFEMFunctionBegin;
    auto prb_mng = m_field.getInterface<ProblemsManager>();
    for (auto m : meshset_vec_ptr) {
      const std::string bc_id =
          problem_name + "_" + field_name + "_" + m->getName();
      auto &bc = bcMapByBlockName.at(bc_id);
      CHKERR prb_mng->removeDofs(problem_name, ROW, *(bc->dofsViewPtr), lo, hi);
      CHKERR prb_mng->removeDofs(problem_name, COL, *(bc->dofsViewPtr), lo, hi);
    }
    MoFEMFunctionReturn(0);
  };
 
  if (is_distributed_mesh) {
 
    CHKERR iterate_meshsets(
 
        m_field.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(std::regex(
 
            (boost::format("%s(.*)") % block_name).str()
 
                ))
 
    );
  } else {
    SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Not implemented");
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

bcMapByBlockName

BcMapByBlockName MoFEM::BcManager::bcMapByBlockName

private

Definition at line 370 of file BcManager.hpp.

cOre

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

private

Definition at line 368 of file BcManager.hpp.

---

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

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