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)

	~ProblemsManager ()
	Destructor. More...

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 problemThis will use one of the mesh partitioning programs available from PETSc See http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatPartitioningType.html 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)Mesh is distributed, that means that each processor keeps only own part of the mesh and shared entities. More...

MoFEMErrorCode	buildProblemOnDistributedMesh (Problem *problem_ptr, const bool square_matrix=true, int verb=VERBOSE)
	build problem data structures, assuming that mesh is distributed (collective)Mesh is distributed, that means that each processor keeps only own part of the mesh and shared entities. 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 elementsFunction which partition finite elements based on dofs partitioning. In addition it sets information about local row and cols dofs at given element on partition. 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 problemIn layout is stored information how many elements is on each processor, for more information look int petsc documentation http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/IS/PetscLayoutCreate.html#PetscLayoutCreate 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 problemRemove DOFs from problem which are on entities on the given range and given field name. On the finite element level, DOFs can be still accessed however local PETSc indices and global PETSc indices are marked with the index -1. 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.

Examples: boundary_marker.cpp, build_problems.cpp, continuity_check_on_skeleton_3d.cpp, EshelbianPlasticity.cpp, forces_and_sources_testing_edge_element.cpp, forces_and_sources_testing_flat_prism_element.cpp, forces_and_sources_testing_users_base.cpp, hcurl_check_approx_in_2d.cpp, hcurl_curl_operator.cpp, hcurl_divergence_operator_2d.cpp, mesh_insert_interface_atom.cpp, prism_elements_from_surface.cpp, prism_polynomial_approximation.cpp, quad_polynomial_approximation.cpp, reaction_diffusion_equation.cpp, Remodeling.cpp, and remove_entities_from_problem.cpp.

Definition at line 34 of file ProblemsManager.hpp.

Constructor & Destructor Documentation

◆ ProblemsManager()

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

Definition at line 74 of file ProblemsManager.cpp.

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

◆ ~ProblemsManager()

MoFEM::ProblemsManager::~ProblemsManager ( )

Destructor.

Definition at line 80 of file ProblemsManager.cpp.

80 {}

Member Function Documentation

◆ debugPartitionedProblem()

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

Definition at line 2355 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->numeredDofsRows->get<Idx_mi_tag>()),
         &(problem_ptr->numeredDofsCols->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 82 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

Examples: boundary_marker.cpp.

Definition at line 3133 of file ProblemsManager.cpp.

                                                                   {
 
   Interface &m_field = cOre;
   const Problem *problem_ptr;
   MoFEMFunctionBegin;
   CHKERR m_field.get_problem(problem_name, &problem_ptr);
   boost::shared_ptr<NumeredDofEntity_multiIndex> dofs;
   switch (rc) {
   case ROW:
     dofs = problem_ptr->getNumeredDofsRows();
     break;
   case COL:
     dofs = problem_ptr->getNumeredDofsCols();
   default:
     SETERRQ(PETSC_COMM_SELF, MOFEM_IMPOSIBLE_CASE, "Should be row or column");
   }
   marker.resize(dofs->size());
   marker.clear();
   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 2311 of file ProblemsManager.cpp.

                                                                              {
   MoFEM::Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   if (verb > 0) {
     std::ostringstream ss;
     ss << "partition_problem: rank = " << m_field.get_comm_rank()
        << " FEs row ghost dofs " << *problem_ptr << " Nb. local dof "
        << problem_ptr->getNbLocalDofsRow() << " nb global row dofs "
        << problem_ptr->getNbDofsRow() << std::endl;
     ss << "partition_problem: rank = " << m_field.get_comm_rank()
        << " FEs col ghost dofs " << *problem_ptr << " Nb. local dof "
        << problem_ptr->getNbLocalDofsCol() << " nb global col dofs "
        << problem_ptr->getNbDofsCol() << std::endl;
     PetscSynchronizedPrintf(m_field.get_comm(), ss.str().c_str());
     PetscSynchronizedFlush(m_field.get_comm(), PETSC_STDOUT);
   }
   if (verb > 1) {
     // FIXME mess if printed on more than one processors
     // std::ostringstream ss;
     std::cout << "rank = " << m_field.get_comm_rank() << " FEs row dofs "
               << *problem_ptr << " Nb. row dof " << problem_ptr->getNbDofsRow()
               << " Nb. local dof " << problem_ptr->getNbLocalDofsRow()
               << std::endl;
     NumeredDofEntity_multiIndex::iterator miit_dd_row =
         problem_ptr->numeredDofsRows->begin();
     for (; miit_dd_row != problem_ptr->numeredDofsRows->end(); miit_dd_row++) {
       std::cout << **miit_dd_row << std::endl;
     }
     std::cout << "rank = " << m_field.get_comm_rank() << " FEs col dofs "
               << *problem_ptr << " Nb. col dof " << problem_ptr->getNbDofsCol()
               << " Nb. local dof " << problem_ptr->getNbLocalDofsCol()
               << std::endl;
     NumeredDofEntity_multiIndex::iterator miit_dd_col =
         problem_ptr->numeredDofsCols->begin();
     for (; miit_dd_col != problem_ptr->numeredDofsCols->end(); miit_dd_col++) {
       std::cout << **miit_dd_col << std::endl;
     }
     // PetscSynchronizedPrintf(comm,ss.str().c_str());
     // PetscSynchronizedFlush(comm,PETSC_STDOUT);
   }
   MoFEMFunctionReturnHot(0);
 }

◆ query_interface()

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

virtual

Implements MoFEM::UnknownInterface.

Definition at line 62 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 47 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 345 of file ProblemsManager.hpp.

◆ synchroniseProblemEntities

PetscBool MoFEM::ProblemsManager::synchroniseProblemEntities

DOFs in fields, not from DOFs on elements.

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