MoFEM::Problem Struct Reference

keeps basic data about problem More...

#include <src/multi_indices/ProblemsMultiIndices.hpp>

Collaboration diagram for MoFEM::Problem:

Classes
struct	SubProblemData
	Subproblem problem data. More...

Public Types
typedef multi_index_container< boost::weak_ptr< std::vector< NumeredDofEntity > >, indexed_by< sequenced<> > >	SequenceDofContainer
using	EmptyFieldBlocks = std::pair< BitFieldId, BitFieldId >

Public Member Functions
auto &	getNumeredRowDofsPtr () const
	get access to numeredRowDofsPtr storing DOFs on rows More...
auto &	getNumeredColDofsPtr () const
	get access to numeredColDofsPtr storing DOFs on cols More...
const auto &	getNumeredFiniteElementsPtr () const
	get access to reference for multi-index storing finite elements More...
boost::shared_ptr< SubProblemData > &	getSubData () const
	Get main problem of sub-problem is. More...
auto &	getComposedProblemsData () const
	Het composed problems data structure. More...
MoFEMErrorCode	getDofByNameEntAndEntDofIdx (const int field_bit_number, const EntityHandle ent, const int ent_dof_idx, const RowColData row_or_col, boost::shared_ptr< NumeredDofEntity > &dof_ptr) const
	get DOFs from problem More...
auto	getNumeredFEsBegin (std::string fe_name) const
auto	getNumeredFEsEnd (std::string fe_name) const
NumeredDofEntity_multiIndex::iterator	getNumeredRowDofsBegin () const
NumeredDofEntity_multiIndex::iterator	getNumeredRowDofsEnd () const
NumeredDofEntity_multiIndex::iterator	getNumeredColDofsBegin () const
NumeredDofEntity_multiIndex::iterator	getNumeredColDofsEnd () const
auto	getNumeredRowDofsByLocIdxBegin (const DofIdx locidx) const
auto	getNumeredRowDofsByLocIdxEnd (const DofIdx locidx) const
auto	getNumeredColDofsByLocIdxBegin (const DofIdx locidx) const
auto	getNumeredColDofsByLocIdxEnd (const DofIdx locidx) const
auto	getNumeredRowDofsByEntBegin (const EntityHandle ent) const
auto	getNumeredRowDofsByEntEnd (const EntityHandle ent) const
auto	getNumeredColDofsByEntBegin (const EntityHandle ent) const
auto	getNumeredColDofsByEntEnd (const EntityHandle ent) const
	Problem (Interface &moab, const EntityHandle meshset)
virtual	~Problem ()=default
BitProblemId	getId () const
auto	getName () const
DofIdx	getNbDofsRow () const
DofIdx	getNbDofsCol () const
DofIdx	getNbLocalDofsRow () const
DofIdx	getNbLocalDofsCol () const
DofIdx	getNbGhostDofsRow () const
DofIdx	getNbGhostDofsCol () const
BitRefLevel	getBitRefLevel () const
BitRefLevel	getBitRefLevelMask () const
MoFEMErrorCode	getRowDofsByPetscGlobalDofIdx (DofIdx idx, const NumeredDofEntity **dof_ptr, MoFEMTypes bh=MF_EXIST) const
	Get the Row Dofs By Petsc Global Dof Idx object. More...
MoFEMErrorCode	getColDofsByPetscGlobalDofIdx (DofIdx idx, const NumeredDofEntity **dof_ptr, MoFEMTypes bh=MF_EXIST) const
	Get the Col Dofs By Petsc Global Dof Idx object. More...
boost::weak_ptr< NumeredDofEntity >	getRowDofsByPetscGlobalDofIdx (DofIdx idx) const
	Get the Row Dofs By Petsc Global Dof Idx object. More...
boost::weak_ptr< NumeredDofEntity >	getColDofsByPetscGlobalDofIdx (DofIdx idx) const
	Get the Col Dofs By Petsc Global Dof Idx object. More...
BitFEId	getBitFEId () const
MoFEMErrorCode	getNumberOfElementsByNameAndPart (MPI_Comm comm, const std::string name, PetscLayout *layout) const
	Get number of finite elements by name on processors. More...
MoFEMErrorCode	getNumberOfElementsByPart (MPI_Comm comm, PetscLayout *layout) const
	Get number of finite elements on processors. More...
auto &	getRowDofsSequence () const
	Get reference to sequence data numbered dof container. More...
auto &	getColDofsSequence () const
	Get reference to sequence data numbered dof container. More...
EmptyFieldBlocks &	getEmptyFieldBlocks () const
	Get the empty field blocks. More...
EmptyFieldBlocks &	addFieldToEmptyFieldBlocks (const EmptyFieldBlocks add_fields) const
	Add fields to the empty field blocks. More...

Public Attributes
EntityHandle	meshset
BitProblemId *	tagId
	Unique problem ID. More...
const char *	tagName
	Problem name. More...
int	tagNameSize
	Size of problem name. More...
BitFEId *	tagBitFEId
	IDs of finite elements in problem. More...
BitRefLevel *	tagBitRefLevel
	BitRef level of finite elements in problem. More...
BitRefLevel *	tagBitRefLevelMask
	BItRefMask of elements in problem. More...
DofIdx	nbDofsRow
	Global number of DOFs in row. More...
DofIdx	nbDofsCol
	Global number of DOFs in col. More...
DofIdx	nbLocDofsRow
	Local number of DOFs in row. More...
DofIdx	nbLocDofsCol
	Local number of DOFs in colIs. More...
DofIdx	nbGhostDofsRow
	Number of ghost DOFs in row. More...
DofIdx	nbGhostDofsCol
	Number of ghost DOFs in col. More...
boost::shared_ptr< NumeredDofEntity_multiIndex >	numeredRowDofsPtr
	store DOFs on rows for this problem More...
boost::shared_ptr< NumeredDofEntity_multiIndex >	numeredColDofsPtr
	store DOFs on columns for this problem More...
boost::shared_ptr< NumeredEntFiniteElement_multiIndex >	numeredFiniteElementsPtr
	store finite elements More...
boost::shared_ptr< SubProblemData >	subProblemData
boost::shared_ptr< ComposedProblemsData >	composedProblemsData

Private Attributes
boost::shared_ptr< SequenceDofContainer >	sequenceRowDofContainer
boost::shared_ptr< SequenceDofContainer >	sequenceColDofContainer
EmptyFieldBlocks	emptyFieldBlocks

Friends
std::ostream &	operator<< (std::ostream &os, const Problem &e)

Detailed Description

keeps basic data about problem

This is low level structure with information about problem, what elements compose problem and what DOFs are on rows and columns.

Todo:

fix names following name convention

Examples

EshelbianPlasticity.cpp, MagneticElement.hpp, Remodeling.cpp, cell_forces.cpp, field_evaluator.cpp, loop_entities.cpp, lorentz_force.cpp, nonlinear_dynamics.cpp, partition_mesh.cpp, poisson_2d_dis_galerkin.cpp, and test_cache_on_entities.cpp.

Definition at line 41 of file ProblemsMultiIndices.hpp.

Member Typedef Documentation

EmptyFieldBlocks

using MoFEM::Problem::EmptyFieldBlocks = std::pair<BitFieldId, BitFieldId>

Definition at line 527 of file ProblemsMultiIndices.hpp.

SequenceDofContainer

typedef multi_index_container<boost::weak_ptr<std::vector<NumeredDofEntity> >, indexed_by<sequenced<> > > MoFEM::Problem::SequenceDofContainer

Definition at line 493 of file ProblemsMultiIndices.hpp.

Constructor & Destructor Documentation

Problem()

MoFEM::Problem::Problem	(	Interface &	moab,
		const EntityHandle	meshset
	)

Definition at line 18 of file ProblemsMultiIndices.cpp.

    : meshset(meshset), numeredRowDofsPtr(new NumeredDofEntity_multiIndex()),
      numeredColDofsPtr(new NumeredDofEntity_multiIndex()),
      numeredFiniteElementsPtr(new NumeredEntFiniteElement_multiIndex()),
      sequenceRowDofContainer(new SequenceDofContainer()),
      sequenceColDofContainer(new SequenceDofContainer()) {
  ErrorCode rval;
  Tag th_ProblemId;
  rval = moab.tag_get_handle("_ProblemId", th_ProblemId);
  MOAB_THROW(rval);
  rval = moab.tag_get_by_ptr(th_ProblemId, &meshset, 1, (const void **)&tagId);
  MOAB_THROW(rval);
  Tag th_ProblemName;
  rval = moab.tag_get_handle("_ProblemName", th_ProblemName);
  MOAB_THROW(rval);
  rval = moab.tag_get_by_ptr(th_ProblemName, &meshset, 1,
                             (const void **)&tagName, &tagNameSize);
  MOAB_THROW(rval);
  Tag th_ProblemFEId;
  rval = moab.tag_get_handle("_ProblemFEId", th_ProblemFEId);
  MOAB_THROW(rval);
  rval = moab.tag_get_by_ptr(th_ProblemFEId, &meshset, 1,
                             (const void **)&tagBitFEId);
  MOAB_THROW(rval);
  Tag th_RefBitLevel;
  rval = moab.tag_get_handle("_RefBitLevel", th_RefBitLevel);
  MOAB_THROW(rval);
  rval = moab.tag_get_by_ptr(th_RefBitLevel, &meshset, 1,
                             (const void **)&tagBitRefLevel);
  MOAB_THROW(rval);
  Tag th_RefBitLevel_Mask;
  rval = moab.tag_get_handle("_RefBitLevelMask", th_RefBitLevel_Mask);
  MOAB_THROW(rval);
  rval = moab.tag_get_by_ptr(th_RefBitLevel_Mask, &meshset, 1,
                             (const void **)&tagBitRefLevelMask);
  MOAB_THROW(rval);
}

~Problem()

virtual MoFEM::Problem::~Problem ( )

virtualdefault

Member Function Documentation

addFieldToEmptyFieldBlocks()

EmptyFieldBlocks & MoFEM::Problem::addFieldToEmptyFieldBlocks ( const EmptyFieldBlocks  add_fields ) const

inline

Add fields to the empty field blocks.

Emtpy field blocks is a pair contains IDs of the fields for which matrix has zero entries.

Parameters

add_fields

Returns

EmptyFieldBlocks&

Definition at line 551 of file ProblemsMultiIndices.hpp.

                                                                      {
    emptyFieldBlocks.first |= add_fields.first;
    emptyFieldBlocks.second |= add_fields.second;
    return emptyFieldBlocks;
  }

getBitFEId()

BitFEId MoFEM::Problem::getBitFEId

(

)

const

Definition at line 64 of file ProblemsMultiIndices.cpp.

{ return *tagBitFEId; }

getBitRefLevel()

BitRefLevel MoFEM::Problem::getBitRefLevel ( ) const

inline

Definition at line 386 of file ProblemsMultiIndices.hpp.

{ return *tagBitRefLevel; }

getBitRefLevelMask()

BitRefLevel MoFEM::Problem::getBitRefLevelMask ( ) const

inline

Definition at line 387 of file ProblemsMultiIndices.hpp.

{ return *tagBitRefLevelMask; }

getColDofsByPetscGlobalDofIdx() [1/2]

boost::weak_ptr< NumeredDofEntity > MoFEM::Problem::getColDofsByPetscGlobalDofIdx

(

DofIdx

idx

)

const

Get the Col Dofs By Petsc Global Dof Idx object.

Parameters

idx

Returns

boost::weak_ptr<NumeredDofEntity>

Definition at line 78 of file ProblemsMultiIndices.cpp.

                                                       {
  MoFEMFunctionBeginHot;
  boost::weak_ptr<NumeredDofEntity> dof_weak_ptr;
  NumeredDofEntity_multiIndex::index<PetscGlobalIdx_mi_tag>::type::iterator dit;
  dit = numeredColDofsPtr->get<PetscGlobalIdx_mi_tag>().find(idx);
  if (dit != numeredColDofsPtr->get<PetscGlobalIdx_mi_tag>().end())
    dof_weak_ptr = *dit;
  return dof_weak_ptr;
}

getColDofsByPetscGlobalDofIdx() [2/2]

MoFEMErrorCode MoFEM::Problem::getColDofsByPetscGlobalDofIdx	(	DofIdx	idx,
		const NumeredDofEntity **	dof_ptr,
		MoFEMTypes	bh = MF_EXIST
	)		const

Get the Col Dofs By Petsc Global Dof Idx object.

Parameters

idx
dof_ptr
bh

Returns

MoFEMErrorCode

Definition at line 102 of file ProblemsMultiIndices.cpp.

                                                                       {
  MoFEMFunctionBeginHot;
  auto weak_dof_ptr = getColDofsByPetscGlobalDofIdx(idx);
  if (auto shared_dof_ptr = weak_dof_ptr.lock()) {
    *dof_ptr = shared_dof_ptr.get();
  } else {
    if (bh == MF_EXIST)
      SETERRQ1(PETSC_COMM_SELF, MOFEM_NOT_FOUND, "row dof <%d> not found", idx);
    *dof_ptr = nullptr;
  }
  MoFEMFunctionReturnHot(0);
}

getColDofsSequence()

auto & MoFEM::Problem::getColDofsSequence ( ) const

inline

Get reference to sequence data numbered dof container.

In sequence data container data are physically stored. The purpose of this is to allocate NumeredDofEntity data in bulk, having only one allocation instead each time entity is inserted. That makes code efficient.

The vector in sequence is destroyed if last entity inside that vector is destroyed. All MoFEM::NumeredDofEntity have aliased shared_ptr which points to the vector.

Returns

MoFEM::Problem::SequenceDofContainer

Definition at line 525 of file ProblemsMultiIndices.hpp.

{ return sequenceColDofContainer; }

getComposedProblemsData()

auto & MoFEM::Problem::getComposedProblemsData ( ) const

inline

Het composed problems data structure.

Definition at line 110 of file ProblemsMultiIndices.hpp.

{ return composedProblemsData; }

getDofByNameEntAndEntDofIdx()

MoFEMErrorCode MoFEM::Problem::getDofByNameEntAndEntDofIdx	(	const int	field_bit_number,
		const EntityHandle	ent,
		const int	ent_dof_idx,
		const RowColData	row_or_col,
		boost::shared_ptr< NumeredDofEntity > &	dof_ptr
	)		const

get DOFs from problem

Note that ent_dof_idx is not coefficient number, is local number of DOFs on the entity. The coefficient number and local index of DOFs or entity are the same on vertices and H1 approximation.

Parameters

field_bit_number	field name field_bit_number = (use m_field.get_field_bit_number(field_name);
ent	entity handle
ent_dof_idx	index of DOFs on entity
row_or_col	ROW or COL
dof_ptr	shared pointer to DOFs if found

Returns

error code

Examples

cell_forces.cpp.

Definition at line 151 of file ProblemsMultiIndices.cpp.

                                                      {
  MoFEMFunctionBegin;
  decltype(numeredRowDofsPtr) numered_dofs;
  switch (row_or_col) {
  case ROW:
    if (!numeredRowDofsPtr) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_INVALID_DATA,
              "Row numbered index in problem not allocated");
    }
    numered_dofs = numeredRowDofsPtr;
    break;
  case COL:
    if (!numeredColDofsPtr) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_INVALID_DATA,
              "Col numbered index in problem not allocated");
    }
    numered_dofs = numeredColDofsPtr;
    break;
  default:
    SETERRQ(PETSC_COMM_SELF, MOFEM_INVALID_DATA,
            "Only ROW and COL is possible for 3rd argument");
  }
 
  auto it =
      numered_dofs->get<Unique_mi_tag>().find(DofEntity::getUniqueIdCalculate(
          ent_dof_idx,
          FieldEntity::getLocalUniqueIdCalculate(field_bit_number, ent)));
  if (it != numered_dofs->get<Unique_mi_tag>().end()) {
    dof_ptr = *it;
  } else {
    dof_ptr = boost::shared_ptr<NumeredDofEntity>();
  }
  MoFEMFunctionReturn(0);
}

getEmptyFieldBlocks()

EmptyFieldBlocks & MoFEM::Problem::getEmptyFieldBlocks ( ) const

inline

Get the empty field blocks.

Emtpy field blocks is a pair contains IDs of the fields for which matrix has zero entries.

Returns

EmptyFieldBlocks&

Definition at line 537 of file ProblemsMultiIndices.hpp.

                                                       {
    return emptyFieldBlocks;
  }

getId()

BitProblemId MoFEM::Problem::getId ( ) const

inline

Definition at line 373 of file ProblemsMultiIndices.hpp.

{ return *((BitProblemId *)tagId); }

getName()

auto MoFEM::Problem::getName ( ) const

inline

Examples

EshelbianPlasticity.cpp, MagneticElement.hpp, Remodeling.cpp, field_evaluator.cpp, and lorentz_force.cpp.

Definition at line 375 of file ProblemsMultiIndices.hpp.

                              {
    return std::string((char *)tagName, tagNameSize);
  }

getNbDofsCol()

DofIdx MoFEM::Problem::getNbDofsCol ( ) const

inline

Definition at line 380 of file ProblemsMultiIndices.hpp.

{ return nbDofsCol; }

getNbDofsRow()

DofIdx MoFEM::Problem::getNbDofsRow ( ) const

inline

Examples

nonlinear_dynamics.cpp, and poisson_2d_dis_galerkin.cpp.

Definition at line 379 of file ProblemsMultiIndices.hpp.

{ return nbDofsRow; }

getNbGhostDofsCol()

DofIdx MoFEM::Problem::getNbGhostDofsCol ( ) const

inline

Definition at line 384 of file ProblemsMultiIndices.hpp.

{ return nbGhostDofsCol; }

getNbGhostDofsRow()

DofIdx MoFEM::Problem::getNbGhostDofsRow ( ) const

inline

Definition at line 383 of file ProblemsMultiIndices.hpp.

{ return nbGhostDofsRow; }

getNbLocalDofsCol()

DofIdx MoFEM::Problem::getNbLocalDofsCol ( ) const

inline

Examples

nonlinear_dynamics.cpp.

Definition at line 382 of file ProblemsMultiIndices.hpp.

{ return nbLocDofsCol; }

getNbLocalDofsRow()

DofIdx MoFEM::Problem::getNbLocalDofsRow ( ) const

inline

Examples

nonlinear_dynamics.cpp.

Definition at line 381 of file ProblemsMultiIndices.hpp.

{ return nbLocDofsRow; }

getNumberOfElementsByNameAndPart()

MoFEMErrorCode MoFEM::Problem::getNumberOfElementsByNameAndPart	(	MPI_Comm	comm,
		const std::string	name,
		PetscLayout *	layout
	)		const

Get number of finite elements by name on processors.

Size retuned IS is equal to size of processors.

What PetscLayout see, http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/IS/index.html

Example of usage for layout

PetscInt rstart, rend;
ierr = PetscLayoutGetRange(layout, &rstart, &rend); CHKERRG(ierr);
int global_size;
ierr = PetscLayoutGetSize(layout,&global_size); CHKERRG(ierr);

Parameters

comm	Communicator
name	Finite element name
layout	Get number of elements on each processor

Returns

error code

Definition at line 117 of file ProblemsMultiIndices.cpp.

                                                                     {
  MoFEMFunctionBegin;
  int size, rank;
  MPI_Comm_size(comm, &size);
  MPI_Comm_rank(comm, &rank);
  CHKERR PetscLayoutCreate(comm, layout);
  CHKERR PetscLayoutSetBlockSize(*layout, 1);
  const auto &fe_by_name_and_part =
      numeredFiniteElementsPtr->get<Composite_Name_And_Part_mi_tag>();
  int nb_elems;
  nb_elems = fe_by_name_and_part.count(boost::make_tuple(name, rank));
  CHKERR PetscLayoutSetLocalSize(*layout, nb_elems);
  CHKERR PetscLayoutSetUp(*layout);
  MoFEMFunctionReturn(0);
}

getNumberOfElementsByPart()

MoFEMErrorCode MoFEM::Problem::getNumberOfElementsByPart	(	MPI_Comm	comm,
		PetscLayout *	layout
	)		const

Get number of finite elements on processors.

Size retuned IS is equal to size of processors.

What PetscLayout see, http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/IS/index.html

Example of usage for layout

PetscInt rstart, rend;
ierr = PetscLayoutGetRange(layout, &rstart, &rend); CHKERRG(ierr);
int global_size;
ierr = PetscLayoutGetSize(layout,&global_size); CHKERRG(ierr);

Parameters

comm	Communicator
layout	Get number of elements on each processor

Returns

error code

Definition at line 134 of file ProblemsMultiIndices.cpp.

                                                                             {
  MoFEMFunctionBegin;
  int size, rank;
  MPI_Comm_size(comm, &size);
  MPI_Comm_rank(comm, &rank);
  CHKERR PetscLayoutCreate(comm, layout);
  CHKERR PetscLayoutSetBlockSize(*layout, 1);
  typedef NumeredEntFiniteElement_multiIndex::index<Part_mi_tag>::type FeByPart;
  const FeByPart &fe_by_part = numeredFiniteElementsPtr->get<Part_mi_tag>();
  int nb_elems;
  nb_elems = fe_by_part.count(rank);
  CHKERR PetscLayoutSetLocalSize(*layout, nb_elems);
  CHKERR PetscLayoutSetUp(*layout);
  MoFEMFunctionReturn(0);
}

getNumeredColDofsBegin()

NumeredDofEntity_multiIndex::iterator MoFEM::Problem::getNumeredColDofsBegin ( ) const

inline

get begin iterator for numeredColDofsPtr (insted you can use #_IT_NUMEREDDOFMOFEMENTITY_COL_FOR_LOOP_ for loops)

Definition at line 204 of file ProblemsMultiIndices.hpp.

                                                                     {
    return numeredColDofsPtr->begin();
  }

getNumeredColDofsByEntBegin()

auto MoFEM::Problem::getNumeredColDofsByEntBegin ( const EntityHandle  ent ) const

inline

get begin iterator for numeredColDofsPtr (insted you can use _IT_NUMEREDDOF_COL_BY_ENT_FOR_LOOP_ for loops)

Definition at line 319 of file ProblemsMultiIndices.hpp.

                                                                 {
    return numeredColDofsPtr->get<Ent_mi_tag>().lower_bound(ent);
  }

getNumeredColDofsByEntEnd()

auto MoFEM::Problem::getNumeredColDofsByEntEnd ( const EntityHandle  ent ) const

inline

get end iterator for numeredColDofsPtr (insted you can use _IT_NUMEREDDOF_COL_BY_ENT_FOR_LOOP_ for loops)

Definition at line 325 of file ProblemsMultiIndices.hpp.

                                                               {
    return numeredColDofsPtr->get<Ent_mi_tag>().upper_bound(ent);
  }

getNumeredColDofsByLocIdxBegin()

auto MoFEM::Problem::getNumeredColDofsByLocIdxBegin ( const DofIdx  locidx ) const

inline

get begin iterator for numeredColDofsPtr (insted you can use _IT_NUMEREDDOF_COL_FOR_LOOP_ for loops)

Definition at line 263 of file ProblemsMultiIndices.hpp.

                                                                 {
    return numeredColDofsPtr->get<PetscLocalIdx_mi_tag>().lower_bound(locidx);
  }

getNumeredColDofsByLocIdxEnd()

auto MoFEM::Problem::getNumeredColDofsByLocIdxEnd ( const DofIdx  locidx ) const

inline

get end iterator for numeredColDofsPtr (insted you can use _IT_NUMEREDDOF_COL_FOR_LOOP_ for loops)

Definition at line 269 of file ProblemsMultiIndices.hpp.

                                                               {
    return numeredColDofsPtr->get<PetscLocalIdx_mi_tag>().upper_bound(locidx);
  }

getNumeredColDofsEnd()

NumeredDofEntity_multiIndex::iterator MoFEM::Problem::getNumeredColDofsEnd ( ) const

inline

get end iterator for numeredColDofsPtr (insted you can use #_IT_NUMEREDDOFMOFEMENTITY_COL_FOR_LOOP_ for loops)

Definition at line 210 of file ProblemsMultiIndices.hpp.

                                                                   {
    return numeredColDofsPtr->end();
  }

getNumeredColDofsPtr()

auto & MoFEM::Problem::getNumeredColDofsPtr ( ) const

inline

get access to numeredColDofsPtr storing DOFs on cols

Definition at line 74 of file ProblemsMultiIndices.hpp.

{ return numeredColDofsPtr; }

getNumeredFEsBegin()

auto MoFEM::Problem::getNumeredFEsBegin ( std::string  fe_name ) const

inline

Definition at line 146 of file ProblemsMultiIndices.hpp.

                                                   {
    return numeredFiniteElementsPtr->get<FiniteElement_name_mi_tag>()
        .lower_bound(fe_name);
  }

getNumeredFEsEnd()

auto MoFEM::Problem::getNumeredFEsEnd ( std::string  fe_name ) const

inline

Definition at line 151 of file ProblemsMultiIndices.hpp.

                                                 {
    return numeredFiniteElementsPtr->get<FiniteElement_name_mi_tag>()
        .upper_bound(fe_name);
  }

getNumeredFiniteElementsPtr()

const auto & MoFEM::Problem::getNumeredFiniteElementsPtr ( ) const

inline

get access to reference for multi-index storing finite elements

Definition at line 79 of file ProblemsMultiIndices.hpp.

                                                         {
    return numeredFiniteElementsPtr;
  }

getNumeredRowDofsBegin()

NumeredDofEntity_multiIndex::iterator MoFEM::Problem::getNumeredRowDofsBegin ( ) const

inline

get begin iterator for numeredRowDofsPtr (insted you can use #_IT_NUMEREDDOFMOFEMENTITY_ROW_FOR_LOOP_ for loops)

Definition at line 192 of file ProblemsMultiIndices.hpp.

                                                                     {
    return numeredRowDofsPtr->begin();
  }

getNumeredRowDofsByEntBegin()

auto MoFEM::Problem::getNumeredRowDofsByEntBegin ( const EntityHandle  ent ) const

inline

get begin iterator for numeredRowDofsPtr (insted you can use _IT_NUMEREDDOF_ROW_BY_ENT_FOR_LOOP_ for loops)

Definition at line 307 of file ProblemsMultiIndices.hpp.

                                                                 {
    return numeredRowDofsPtr->get<Ent_mi_tag>().lower_bound(ent);
  }

getNumeredRowDofsByEntEnd()

auto MoFEM::Problem::getNumeredRowDofsByEntEnd ( const EntityHandle  ent ) const

inline

get end iterator for numeredRowDofsPtr (insted you can use _IT_NUMEREDDOF_ROW_BY_ENT_FOR_LOOP_ for loops)

Definition at line 313 of file ProblemsMultiIndices.hpp.

                                                               {
    return numeredRowDofsPtr->get<Ent_mi_tag>().upper_bound(ent);
  }

getNumeredRowDofsByLocIdxBegin()

auto MoFEM::Problem::getNumeredRowDofsByLocIdxBegin ( const DofIdx  locidx ) const

inline

get begin iterator for numeredRowDofsPtr (insted you can use _IT_NUMEREDDOF_ROW_FOR_LOOP_ for loops)

Definition at line 251 of file ProblemsMultiIndices.hpp.

                                                                 {
    return numeredRowDofsPtr->get<PetscLocalIdx_mi_tag>().lower_bound(locidx);
  }

getNumeredRowDofsByLocIdxEnd()

auto MoFEM::Problem::getNumeredRowDofsByLocIdxEnd ( const DofIdx  locidx ) const

inline

get end iterator for numeredRowDofsPtr (insted you can use _IT_NUMEREDDOF_ROW_FOR_LOOP_ for loops)

Definition at line 257 of file ProblemsMultiIndices.hpp.

                                                               {
    return numeredRowDofsPtr->get<PetscLocalIdx_mi_tag>().upper_bound(locidx);
  }

getNumeredRowDofsEnd()

NumeredDofEntity_multiIndex::iterator MoFEM::Problem::getNumeredRowDofsEnd ( ) const

inline

get end iterator for numeredRowDofsPtr (insted you can use #_IT_NUMEREDDOFMOFEMENTITY_ROW_FOR_LOOP_ for loops)

Definition at line 198 of file ProblemsMultiIndices.hpp.

                                                                   {
    return numeredRowDofsPtr->end();
  }

getNumeredRowDofsPtr()

auto & MoFEM::Problem::getNumeredRowDofsPtr ( ) const

inline

get access to numeredRowDofsPtr storing DOFs on rows

Examples

eigen_elastic.cpp, and test_cache_on_entities.cpp.

Definition at line 69 of file ProblemsMultiIndices.hpp.

{ return numeredRowDofsPtr; }

getRowDofsByPetscGlobalDofIdx() [1/2]

boost::weak_ptr< NumeredDofEntity > MoFEM::Problem::getRowDofsByPetscGlobalDofIdx

(

DofIdx

idx

)

const

Get the Row Dofs By Petsc Global Dof Idx object.

Parameters

idx

Returns

boost::weak_ptr<NumeredDofEntity>

Definition at line 67 of file ProblemsMultiIndices.cpp.

                                                       {
  MoFEMFunctionBeginHot;
  boost::weak_ptr<NumeredDofEntity> dof_weak_ptr;
  NumeredDofEntity_multiIndex::index<PetscGlobalIdx_mi_tag>::type::iterator dit;
  dit = numeredRowDofsPtr->get<PetscGlobalIdx_mi_tag>().find(idx);
  if (dit != numeredRowDofsPtr->get<PetscGlobalIdx_mi_tag>().end())
    dof_weak_ptr = *dit;
  return dof_weak_ptr;
}

getRowDofsByPetscGlobalDofIdx() [2/2]

MoFEMErrorCode MoFEM::Problem::getRowDofsByPetscGlobalDofIdx	(	DofIdx	idx,
		const NumeredDofEntity **	dof_ptr,
		MoFEMTypes	bh = MF_EXIST
	)		const

Get the Row Dofs By Petsc Global Dof Idx object.

Parameters

idx
dof_ptr
bh

Returns

MoFEMErrorCode

Definition at line 88 of file ProblemsMultiIndices.cpp.

                                                                       {
  MoFEMFunctionBeginHot;
  auto weak_dof_ptr = getRowDofsByPetscGlobalDofIdx(idx);
  if (auto shared_dof_ptr = weak_dof_ptr.lock()) {
    *dof_ptr = shared_dof_ptr.get();
  } else {
    if (bh == MF_EXIST)
      SETERRQ1(PETSC_COMM_SELF, MOFEM_NOT_FOUND, "row dof <%d> not found", idx);
    *dof_ptr = nullptr;
  }
  MoFEMFunctionReturnHot(0);
}

getRowDofsSequence()

auto & MoFEM::Problem::getRowDofsSequence ( ) const

inline

Get reference to sequence data numbered dof container.

In sequence data container data are physically stored. The purpose of this is to allocate NumeredDofEntity data in bulk, having only one allocation instead each time entity is inserted. That makes code efficient.

The vector in sequence is destroyed if last entity inside that vector is destroyed. All MoFEM::NumeredDofEntity have aliased shared_ptr which points to the vector.

Returns

MoFEM::Problem::SequenceDofContainer

Definition at line 509 of file ProblemsMultiIndices.hpp.

{ return sequenceRowDofContainer; }

getSubData()

boost::shared_ptr< SubProblemData > & MoFEM::Problem::getSubData ( ) const

inline

Get main problem of sub-problem is.

Returns

sub problem data structure

Examples

cell_forces.cpp.

Definition at line 98 of file ProblemsMultiIndices.hpp.

                                                             {
    return subProblemData;
  }

Friends And Related Function Documentation

operator<<

std::ostream & operator<< ( std::ostream &  os, const Problem &  e  )

friend

Definition at line 56 of file ProblemsMultiIndices.cpp.

                                                         {
  os << "Problem id " << e.getId() << " name " << e.getName() << endl;
  os << "FiniteElement id " << e.getBitFEId() << endl;
  os << "BitRefLevel " << e.getBitRefLevel() << endl;
  os << "BitRefLevelMask " << e.getBitRefLevelMask();
  return os;
}

Member Data Documentation

composedProblemsData

boost::shared_ptr<ComposedProblemsData> MoFEM::Problem::composedProblemsData

mutable

Pointer to data structure from which this problem is composed

Definition at line 105 of file ProblemsMultiIndices.hpp.

emptyFieldBlocks

EmptyFieldBlocks MoFEM::Problem::emptyFieldBlocks

mutableprivate

Definition at line 562 of file ProblemsMultiIndices.hpp.

meshset

EntityHandle MoFEM::Problem::meshset

Problem meshset (on tag of this meshset all data related to problem are stored)

Definition at line 43 of file ProblemsMultiIndices.hpp.

nbDofsCol

DofIdx MoFEM::Problem::nbDofsCol

mutable

Global number of DOFs in col.

Definition at line 53 of file ProblemsMultiIndices.hpp.

nbDofsRow

DofIdx MoFEM::Problem::nbDofsRow

mutable

Global number of DOFs in row.

Examples

partition_mesh.cpp.

Definition at line 52 of file ProblemsMultiIndices.hpp.

nbGhostDofsCol

DofIdx MoFEM::Problem::nbGhostDofsCol

mutable

Number of ghost DOFs in col.

Definition at line 57 of file ProblemsMultiIndices.hpp.

nbGhostDofsRow

DofIdx MoFEM::Problem::nbGhostDofsRow

mutable

Number of ghost DOFs in row.

Definition at line 56 of file ProblemsMultiIndices.hpp.

nbLocDofsCol

DofIdx MoFEM::Problem::nbLocDofsCol

mutable

Local number of DOFs in colIs.

Definition at line 55 of file ProblemsMultiIndices.hpp.

nbLocDofsRow

DofIdx MoFEM::Problem::nbLocDofsRow

mutable

Local number of DOFs in row.

Definition at line 54 of file ProblemsMultiIndices.hpp.

numeredColDofsPtr

boost::shared_ptr<NumeredDofEntity_multiIndex> MoFEM::Problem::numeredColDofsPtr

mutable

store DOFs on columns for this problem

Definition at line 62 of file ProblemsMultiIndices.hpp.

numeredFiniteElementsPtr

boost::shared_ptr<NumeredEntFiniteElement_multiIndex> MoFEM::Problem::numeredFiniteElementsPtr

mutable

store finite elements

Definition at line 64 of file ProblemsMultiIndices.hpp.

numeredRowDofsPtr

boost::shared_ptr<NumeredDofEntity_multiIndex> MoFEM::Problem::numeredRowDofsPtr

mutable

store DOFs on rows for this problem

Examples

cell_forces.cpp.

Definition at line 60 of file ProblemsMultiIndices.hpp.

sequenceColDofContainer

boost::shared_ptr<SequenceDofContainer> MoFEM::Problem::sequenceColDofContainer

mutableprivate

Definition at line 560 of file ProblemsMultiIndices.hpp.

sequenceRowDofContainer

boost::shared_ptr<SequenceDofContainer> MoFEM::Problem::sequenceRowDofContainer

mutableprivate

Definition at line 559 of file ProblemsMultiIndices.hpp.

subProblemData

boost::shared_ptr<SubProblemData> MoFEM::Problem::subProblemData

mutable

Pointer to data structure. This pointer has allocated data only for sub problems.

Examples

EshelbianPlasticity.cpp.

Definition at line 92 of file ProblemsMultiIndices.hpp.

tagBitFEId

BitFEId* MoFEM::Problem::tagBitFEId

IDs of finite elements in problem.

Definition at line 48 of file ProblemsMultiIndices.hpp.

tagBitRefLevel

BitRefLevel* MoFEM::Problem::tagBitRefLevel

BitRef level of finite elements in problem.

Definition at line 49 of file ProblemsMultiIndices.hpp.

tagBitRefLevelMask

BitRefLevel* MoFEM::Problem::tagBitRefLevelMask

BItRefMask of elements in problem.

Definition at line 50 of file ProblemsMultiIndices.hpp.

tagId

BitProblemId* MoFEM::Problem::tagId

Unique problem ID.

Definition at line 45 of file ProblemsMultiIndices.hpp.

tagName

const char* MoFEM::Problem::tagName

Problem name.

Definition at line 46 of file ProblemsMultiIndices.hpp.

tagNameSize

int MoFEM::Problem::tagNameSize

Size of problem name.

Definition at line 47 of file ProblemsMultiIndices.hpp.

---

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

• src/multi_indices/ProblemsMultiIndices.hpp
• src/multi_indices/impl/ProblemsMultiIndices.cpp