Problem manager is used to build and partition problems. More...

#include <src/interfaces/ProblemsManager.hpp>

Inheritance diagram for MoFEM::ProblemsManager:

Collaboration diagram for MoFEM::ProblemsManager:

Public Member Functions
MoFEMErrorCode	query_interface (const MOFEMuuid &uuid, UnknownInterface **iface) const

	ProblemsManager (const MoFEM::Core &core)

virtual	~ProblemsManager ()=default

MoFEMErrorCode	getOptions ()

MoFEMErrorCode	partitionMesh (const Range &ents, const int dim, const int adj_dim, const int n_parts, Tag th_vertex_weights=nullptr, Tag th_edge_weights=nullptr, Tag *th_part_weights=nullptr, int verb=VERBOSE, const bool debug=false)
	Set partition tag to each finite element in the problem. More...

MoFEMErrorCode	buildProblem (const std::string name, const bool square_matrix, int verb=VERBOSE)
	build problem data structures More...

MoFEMErrorCode	buildProblem (Problem *problem_ptr, const bool square_matrix, int verb=VERBOSE)
	build problem data structures More...

MoFEMErrorCode	buildProblemOnDistributedMesh (const std::string name, const bool square_matrix, int verb=VERBOSE)
	build problem data structures, assuming that mesh is distributed (collective) More...

MoFEMErrorCode	buildProblemOnDistributedMesh (Problem *problem_ptr, const bool square_matrix=true, int verb=VERBOSE)
	build problem data structures, assuming that mesh is distributed (collective) More...

MoFEMErrorCode	buildSubProblem (const std::string out_name, const std::vector< std::string > &fields_row, const std::vector< std::string > &fields_col, const std::string main_problem, const bool square_matrix=true, const map< std::string, std::pair< EntityType, EntityType >> entityMapRow=nullptr, const map< std::string, std::pair< EntityType, EntityType >> entityMapCol=nullptr, int verb=VERBOSE)
	build sub problem More...

MoFEMErrorCode	buildCompsedProblem (const std::string out_name, const std::vector< std::string > add_row_problems, const std::vector< std::string > add_col_problems, const bool square_matrix=true, int verb=1)
	build composite problem More...

MoFEMErrorCode	inheritPartition (const std::string name, const std::string problem_for_rows, bool copy_rows, const std::string problem_for_cols, bool copy_cols, int verb=VERBOSE)
	build indexing and partition problem inheriting indexing and partitioning from two other problems More...

MoFEMErrorCode	partitionSimpleProblem (const std::string name, int verb=VERBOSE)
	partition problem dofs More...

MoFEMErrorCode	partitionProblem (const std::string name, int verb=VERBOSE)
	partition problem dofs (collective) More...

MoFEMErrorCode	printPartitionedProblem (const Problem *problem_ptr, int verb=VERBOSE)

MoFEMErrorCode	debugPartitionedProblem (const Problem *problem_ptr, int verb=VERBOSE)

MoFEMErrorCode	partitionFiniteElements (const std::string name, bool part_from_moab=false, int low_proc=-1, int hi_proc=-1, int verb=VERBOSE)
	partition finite elements More...

MoFEMErrorCode	partitionGhostDofs (const std::string name, int verb=VERBOSE)
	determine ghost nodes More...

MoFEMErrorCode	partitionGhostDofsOnDistributedMesh (const std::string name, int verb=VERBOSE)
	determine ghost nodes on distributed meshes More...

MoFEMErrorCode	getFEMeshset (const std::string prb_name, const std::string fe_name, EntityHandle *meshset) const
	create add entities of finite element in the problem More...

MoFEMErrorCode	getProblemElementsLayout (const std::string name, const std::string fe_name, PetscLayout *layout) const
	Get layout of elements in the problem. More...

MoFEMErrorCode	removeDofsOnEntities (const std::string problem_name, const std::string field_name, const Range ents, const int lo_coeff=0, const int hi_coeff=MAX_DOFS_ON_ENTITY, int verb=VERBOSE, const bool debug=false)
	Remove DOFs from problem. More...

MoFEMErrorCode	markDofs (const std::string problem_name, RowColData rc, const Range ents, std::vector< bool > &marker)
	Create vector with marked indices. More...

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

template<class IFACE , bool VERIFY = false>
MoFEMErrorCode	getInterface (const MOFEMuuid &uuid, IFACE *&iface) const
	Get interface by uuid and return reference to pointer of interface. More...

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

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

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

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

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

virtual	~UnknownInterface ()=default

virtual MoFEMErrorCode	getLibVersion (Version &version) const
	Get library version. More...

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

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

template<>
MoFEMErrorCode	getInterface (const MOFEMuuid &uuid, UnknownInterface *&iface) const

Public Attributes
MoFEM::Core &	cOre

PetscBool	buildProblemFromFields

PetscBool	synchroniseProblemEntities
	DOFs in fields, not from DOFs on elements. More...

Private Attributes
PetscLogEvent	MOFEM_EVENT_ProblemsManager

Additional Inherited Members
Protected Member Functions inherited from MoFEM::UnknownInterface
boost::typeindex::type_index	getClassIdx (const MOFEMuuid &uid) const
	Get type name for interface Id. More...

MOFEMuuid	getUId (const boost::typeindex::type_index &class_idx) const
	Get interface Id for class name. More...

Detailed Description

Problem manager is used to build and partition problems.

Definition at line 34 of file ProblemsManager.hpp.

Constructor & Destructor Documentation

◆ ProblemsManager()

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

Definition at line 82 of file ProblemsManager.cpp.

     : cOre(const_cast<MoFEM::Core &>(core)),
       buildProblemFromFields(PETSC_FALSE),
       synchroniseProblemEntities(PETSC_FALSE) {
   PetscLogEventRegister("ProblemsManager", 0, &MOFEM_EVENT_ProblemsManager);
 }

◆ ~ProblemsManager()

virtual MoFEM::ProblemsManager::~ProblemsManager ( )

virtualdefault

Member Function Documentation

◆ debugPartitionedProblem()

MoFEMErrorCode MoFEM::ProblemsManager::debugPartitionedProblem	(	const Problem *	problem_ptr,
		int	verb = `VERBOSE`
	)

Definition at line 2354 of file ProblemsManager.cpp.

                                                                              {
   MoFEM::Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   if (debug > 0) {
  
     typedef NumeredDofEntity_multiIndex::index<Idx_mi_tag>::type
         NumeredDofEntitysByIdx;
     NumeredDofEntitysByIdx::iterator dit, hi_dit;
     const NumeredDofEntitysByIdx *numered_dofs_ptr[] = {
         &(problem_ptr->numeredRowDofsPtr->get<Idx_mi_tag>()),
         &(problem_ptr->numeredColDofsPtr->get<Idx_mi_tag>())};
  
     int *nbdof_ptr[] = {&problem_ptr->nbDofsRow, &problem_ptr->nbDofsCol};
     int *local_nbdof_ptr[] = {&problem_ptr->nbLocDofsRow,
                               &problem_ptr->nbLocDofsCol};
  
     for (int ss = 0; ss < 2; ss++) {
  
       dit = numered_dofs_ptr[ss]->begin();
       hi_dit = numered_dofs_ptr[ss]->end();
       for (; dit != hi_dit; dit++) {
         if ((*dit)->getPart() == (unsigned int)m_field.get_comm_rank()) {
           if ((*dit)->getPetscLocalDofIdx() < 0) {
             std::ostringstream zz;
             zz << "rank " << m_field.get_comm_rank() << " " << **dit;
             SETERRQ2(PETSC_COMM_SELF, MOFEM_IMPOSIBLE_CASE,
                      "local dof index for %d (0-row, 1-col) not set, i.e. has "
                      "negative value\n %s",
                      ss, zz.str().c_str());
           }
           if ((*dit)->getPetscLocalDofIdx() >= *local_nbdof_ptr[ss]) {
             std::ostringstream zz;
             zz << "rank " << m_field.get_comm_rank() << " " << **dit;
             SETERRQ2(PETSC_COMM_SELF, MOFEM_IMPOSIBLE_CASE,
                      "local dofs for %d (0-row, 1-col) out of range\n %s", ss,
                      zz.str().c_str());
           }
         } else {
           if ((*dit)->getPetscGlobalDofIdx() < 0) {
             std::ostringstream zz;
             zz << "rank " << m_field.get_comm_rank() << " "
                << dit->get()->getBitRefLevel() << " " << **dit;
             SETERRQ2(PETSC_COMM_SELF, MOFEM_IMPOSIBLE_CASE,
                      "global dof index for %d (0-row, 1-col) row not set, i.e. "
                      "has negative value\n %s",
                      ss, zz.str().c_str());
           }
           if ((*dit)->getPetscGlobalDofIdx() >= *nbdof_ptr[ss]) {
             std::ostringstream zz;
             zz << "rank " << m_field.get_comm_rank() << " nb_dofs "
                << *nbdof_ptr[ss] << " " << **dit;
             SETERRQ2(PETSC_COMM_SELF, MOFEM_IMPOSIBLE_CASE,
                      "global dofs for %d (0-row, 1-col) out of range\n %s", ss,
                      zz.str().c_str());
           }
         }
       }
     }
   }
   MoFEMFunctionReturnHot(0);
 }

◆ getOptions()

MoFEMErrorCode MoFEM::ProblemsManager::getOptions ( )

Definition at line 89 of file ProblemsManager.cpp.

                                            {
   MoFEM::Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   CHKERR PetscOptionsBegin(m_field.get_comm(), "", "Problem manager", "none");
   {
     CHKERR PetscOptionsBool(
         "-problem_build_from_fields",
         "Add DOFs to problem directly from fields not through DOFs on elements",
         "", buildProblemFromFields, &buildProblemFromFields, NULL);
   }
   ierr = PetscOptionsEnd();
   CHKERRG(ierr);
   MoFEMFunctionReturnHot(0);
 }

◆ markDofs()

MoFEMErrorCode MoFEM::ProblemsManager::markDofs	(	const std::string	problem_name,
		RowColData	rc,
		const Range	ents,
		std::vector< bool > &	marker
	)

Create vector with marked indices.

Vector with local DOFs marked by entities

Parameters

problem_name
row
ents
marker

Returns: MoFEMErrorCode

Definition at line 3078 of file ProblemsManager.cpp.

                                                                   {
  
   Interface &m_field = cOre;
   const Problem *problem_ptr;
   ProblemManagerFunctionBegin;
   CHKERR m_field.get_problem(problem_name, &problem_ptr);
   boost::shared_ptr<NumeredDofEntity_multiIndex> dofs;
   switch (rc) {
   case ROW:
     dofs = problem_ptr->getNumeredRowDofsPtr();
     break;
   case COL:
     dofs = problem_ptr->getNumeredColDofsPtr();
   default:
     SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSIBLE_CASE, "Should be row or column");
   }
   marker.resize(dofs->size());
   std::fill(marker.begin(), marker.end(), false);
   for (auto p = ents.pair_begin(); p != ents.pair_end(); ++p) {
     auto lo = dofs->get<Ent_mi_tag>().lower_bound(p->first);
     auto hi = dofs->get<Ent_mi_tag>().upper_bound(p->second);
     for (; lo != hi; ++lo)
       marker[(*lo)->getPetscLocalDofIdx()] = true;
   }
   MoFEMFunctionReturn(0);
 }

◆ printPartitionedProblem()

MoFEMErrorCode MoFEM::ProblemsManager::printPartitionedProblem	(	const Problem *	problem_ptr,
		int	verb = `VERBOSE`
	)

Definition at line 2335 of file ProblemsManager.cpp.

                                                                              {
   MoFEM::Interface &m_field = cOre;
   ProblemManagerFunctionBegin;
  
   if (verb > QUIET) {
  
     MOFEM_LOG("SYNC", Sev::inform)
         << problem_ptr->getName() << " Nb. local dof "
         << problem_ptr->getNbLocalDofsRow() << " by "
         << problem_ptr->getNbLocalDofsCol() << " nb global dofs "
         << problem_ptr->getNbDofsRow() << " by " << problem_ptr->getNbDofsCol();
  
     MOFEM_LOG_SYNCHRONISE(m_field.get_comm())
   }
  
   MoFEMFunctionReturn(0);
 }

◆ query_interface()

MoFEMErrorCode MoFEM::ProblemsManager::query_interface	(	const MOFEMuuid &	uuid,
		UnknownInterface **	iface
	)		const

virtual

Implements MoFEM::UnknownInterface.

Definition at line 70 of file ProblemsManager.cpp.

                                                                  {
   MoFEMFunctionBeginHot;
   *iface = NULL;
   if (uuid == IDD_MOFEMProblemsManager) {
     *iface = const_cast<ProblemsManager *>(this);
     MoFEMFunctionReturnHot(0);
   }
   SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "unknown interface");
   MoFEMFunctionReturnHot(0);
 }

Member Data Documentation

◆ buildProblemFromFields

PetscBool MoFEM::ProblemsManager::buildProblemFromFields

If set to true, problem is build from

Definition at line 43 of file ProblemsManager.hpp.

◆ cOre

MoFEM::Core& MoFEM::ProblemsManager::cOre

Definition at line 39 of file ProblemsManager.hpp.

◆ MOFEM_EVENT_ProblemsManager

PetscLogEvent MoFEM::ProblemsManager::MOFEM_EVENT_ProblemsManager

private

Definition at line 341 of file ProblemsManager.hpp.

◆ synchroniseProblemEntities

PetscBool MoFEM::ProblemsManager::synchroniseProblemEntities

DOFs in fields, not from DOFs on elements.

Definition at line 46 of file ProblemsManager.hpp.

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

src/interfaces/ProblemsManager.hpp
src/interfaces/impl/ProblemsManager.cpp

Public Member Functions

Public Attributes

Private Attributes

Additional Inherited Members