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...
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, remove_entities_from_problem.cpp, remove_entities_from_problem_not_partitioned.cpp, and test_cache_on_entities.cpp.

Definition at line 40 of file ProblemsMultiIndices.hpp.

Member Typedef Documentation

EmptyFieldBlocks

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

Definition at line 502 of file ProblemsMultiIndices.hpp.

SequenceDofContainer

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

Definition at line 468 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 526 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 361 of file ProblemsMultiIndices.hpp.

{ return *tagBitRefLevel; }

getBitRefLevelMask()

BitRefLevel MoFEM::Problem::getBitRefLevelMask ( ) const

inline

Definition at line 362 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 500 of file ProblemsMultiIndices.hpp.

{ return sequenceColDofContainer; }

getComposedProblemsData()

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

inline

Het composed problems data structure.

Definition at line 109 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 512 of file ProblemsMultiIndices.hpp.

                                                       {
    return emptyFieldBlocks;
  }

getId()

BitProblemId MoFEM::Problem::getId ( ) const

inline

Definition at line 348 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 350 of file ProblemsMultiIndices.hpp.

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

getNbDofsCol()

DofIdx MoFEM::Problem::getNbDofsCol ( ) const

inline

Definition at line 355 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 354 of file ProblemsMultiIndices.hpp.

{ return nbDofsRow; }

getNbGhostDofsCol()

DofIdx MoFEM::Problem::getNbGhostDofsCol ( ) const

inline

Definition at line 359 of file ProblemsMultiIndices.hpp.

{ return nbGhostDofsCol; }

getNbGhostDofsRow()

DofIdx MoFEM::Problem::getNbGhostDofsRow ( ) const

inline

Definition at line 358 of file ProblemsMultiIndices.hpp.

{ return nbGhostDofsRow; }

getNbLocalDofsCol()

DofIdx MoFEM::Problem::getNbLocalDofsCol ( ) const

inline

Examples

nonlinear_dynamics.cpp.

Definition at line 357 of file ProblemsMultiIndices.hpp.

{ return nbLocDofsCol; }

getNbLocalDofsRow()

DofIdx MoFEM::Problem::getNbLocalDofsRow ( ) const

inline

Examples

nonlinear_dynamics.cpp.

Definition at line 356 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 179 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 294 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 300 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 238 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 244 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 185 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 73 of file ProblemsMultiIndices.hpp.

{ return numeredColDofsPtr; }

getNumeredFiniteElementsPtr()

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

inline

get access to reference for multi-index storing finite elements

Definition at line 78 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 167 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 282 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 288 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 226 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 232 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 173 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 68 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 484 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, remove_entities_from_problem.cpp, and remove_entities_from_problem_not_partitioned.cpp.

Definition at line 97 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 104 of file ProblemsMultiIndices.hpp.

emptyFieldBlocks

EmptyFieldBlocks MoFEM::Problem::emptyFieldBlocks

mutableprivate

Definition at line 537 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 42 of file ProblemsMultiIndices.hpp.

nbDofsCol

DofIdx MoFEM::Problem::nbDofsCol

mutable

Global number of DOFs in col.

Definition at line 52 of file ProblemsMultiIndices.hpp.

nbDofsRow

DofIdx MoFEM::Problem::nbDofsRow

mutable

Global number of DOFs in row.

Examples

partition_mesh.cpp.

Definition at line 51 of file ProblemsMultiIndices.hpp.

nbGhostDofsCol

DofIdx MoFEM::Problem::nbGhostDofsCol

mutable

Number of ghost DOFs in col.

Definition at line 56 of file ProblemsMultiIndices.hpp.

nbGhostDofsRow

DofIdx MoFEM::Problem::nbGhostDofsRow

mutable

Number of ghost DOFs in row.

Definition at line 55 of file ProblemsMultiIndices.hpp.

nbLocDofsCol

DofIdx MoFEM::Problem::nbLocDofsCol

mutable

Local number of DOFs in colIs.

Definition at line 54 of file ProblemsMultiIndices.hpp.

nbLocDofsRow

DofIdx MoFEM::Problem::nbLocDofsRow

mutable

Local number of DOFs in row.

Definition at line 53 of file ProblemsMultiIndices.hpp.

numeredColDofsPtr

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

mutable

store DOFs on columns for this problem

Definition at line 61 of file ProblemsMultiIndices.hpp.

numeredFiniteElementsPtr

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

mutable

store finite elements

Definition at line 63 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 59 of file ProblemsMultiIndices.hpp.

sequenceColDofContainer

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

mutableprivate

Definition at line 535 of file ProblemsMultiIndices.hpp.

sequenceRowDofContainer

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

mutableprivate

Definition at line 534 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 91 of file ProblemsMultiIndices.hpp.

tagBitFEId

BitFEId* MoFEM::Problem::tagBitFEId

IDs of finite elements in problem.

Definition at line 47 of file ProblemsMultiIndices.hpp.

tagBitRefLevel

BitRefLevel* MoFEM::Problem::tagBitRefLevel

BitRef level of finite elements in problem.

Definition at line 48 of file ProblemsMultiIndices.hpp.

tagBitRefLevelMask

BitRefLevel* MoFEM::Problem::tagBitRefLevelMask

BItRefMask of elements in problem.

Definition at line 49 of file ProblemsMultiIndices.hpp.

tagId

BitProblemId* MoFEM::Problem::tagId

Unique problem ID.

Definition at line 44 of file ProblemsMultiIndices.hpp.

tagName

const char* MoFEM::Problem::tagName

Problem name.

Definition at line 45 of file ProblemsMultiIndices.hpp.

tagNameSize

int MoFEM::Problem::tagNameSize

Size of problem name.

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