Simple interface for fast problem set-up. More...

#include <src/interfaces/Simple.hpp>

Inheritance diagram for MoFEM::Simple:

Collaboration diagram for MoFEM::Simple:

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

	Simple (const MoFEM::Core &core)

	~Simple ()
	Destructor. More...

MoFEMErrorCode	getOptions ()
	get options More...

MoFEMErrorCode	loadFile ()
	Load mesh file. More...

MoFEMErrorCode	addDomainField (const std::string &name, const FieldSpace space, const FieldApproximationBase base, const FieldCoefficientsNumber nb_of_coefficients, const TagType tag_type=MB_TAG_SPARSE, const enum MoFEMTypes bh=MF_ZERO, int verb=-1)
	Add field on domain. More...

MoFEMErrorCode	addBoundaryField (const std::string &name, const FieldSpace space, const FieldApproximationBase base, const FieldCoefficientsNumber nb_of_coefficients, const TagType tag_type=MB_TAG_SPARSE, const enum MoFEMTypes bh=MF_ZERO, int verb=-1)
	Add field on boundary. More...

MoFEMErrorCode	addSkeletonField (const std::string &name, const FieldSpace space, const FieldApproximationBase base, const FieldCoefficientsNumber nb_of_coefficients, const TagType tag_type=MB_TAG_SPARSE, const enum MoFEMTypes bh=MF_ZERO, int verb=-1)
	Add field on skeleton. More...

MoFEMErrorCode	addDataField (const std::string &name, const FieldSpace space, const FieldApproximationBase base, const FieldCoefficientsNumber nb_of_coefficients, const TagType tag_type=MB_TAG_SPARSE, const enum MoFEMTypes bh=MF_ZERO, int verb=-1)
	Add field on domain. More...

MoFEMErrorCode	defineFiniteElements ()
	Define finite elements. More...

MoFEMErrorCode	defineProblem (const PetscBool is_partitioned=PETSC_TRUE)
	define problem More...

MoFEMErrorCode	setFieldOrder (const std::string field_name, const int order, const Range *ents=NULL)
	Set field order. More...

MoFEMErrorCode	buildFields ()
	Build fields. More...

MoFEMErrorCode	buildFiniteElements ()
	Build finite elements. More...

MoFEMErrorCode	buildProblem ()
	Build problem. More...

MoFEMErrorCode	setUp ()
	Setup problem. More...

MoFEMErrorCode	getDM (DM *dm)
	Get DM. More...

int	getDim () const

const std::string &	getDomainFEName () const

const std::string &	getBoundaryFEName () const

const std::string &	getSkeletonFEName () const

std::string &	getDomainFEName ()

std::string &	getBoundaryFEName ()

std::string &	getSkeletonFEName ()

std::vector< std::string > &	getOtherFiniteElements ()

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 ()

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

Private Attributes
const BitRefLevel	bitLevel

PetscLogEvent	MOFEM_EVENT_SimpleLoadMesh

PetscLogEvent	MOFEM_EVENT_SimpleBuildFields

PetscLogEvent	MOFEM_EVENT_SimpleBuildFiniteElements

PetscLogEvent	MOFEM_EVENT_SimpleBuildProblem

PetscLogEvent	MOFEM_EVENT_SimpleKSPSolve

EntityHandle	meshSet
	domain meshset More...

EntityHandle	boundaryMeshset
	meshset with boundary More...

std::string	nameOfProblem
	problem name More...

std::vector< std::string >	domainFields
	domain fields More...

std::vector< std::string >	boundaryFields
	boundary fields More...

std::vector< std::string >	skeletonFields
	fields on the skeleton More...

std::vector< std::string >	dataFields
	Data fields. More...

std::map< std::string, std::pair< int, Range > >	fieldsOrder
	fields order More...

std::string	domainFE
	domain finite element More...

std::string	boundaryFE
	boundary finite element More...

std::string	skeletonFE
	skeleton finite element More...

std::vector< std::string >	otherFEs
	Other finite elements. More...

char	meshFileName [255]
	mesh file name More...

int	dIm
	dimension of problem More...

DM	dM

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

Simple interface for fast problem set-up.

Examples:: hello_world.cpp, mesh_smoothing.cpp, simple_elasticity.cpp, simple_interface.cpp, and testing_jacobian_of_hook_element.cpp.

Definition at line 33 of file Simple.hpp.

Constructor & Destructor Documentation

◆ Simple()

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

Definition at line 40 of file Simple.cpp.

     : cOre(const_cast<Core &>(core)), bitLevel(BitRefLevel().set(0)),
       meshSet(0), boundaryMeshset(0), domainFE("dFE"), boundaryFE("bFE"),
       skeletonFE("sFE"), dIm(-1), dM(PETSC_NULL) {
   PetscLogEventRegister("LoadMesh", 0, &MOFEM_EVENT_SimpleLoadMesh);
   PetscLogEventRegister("buildFields", 0, &MOFEM_EVENT_SimpleBuildFields);
   PetscLogEventRegister("buildFiniteElements", 0,
                         &MOFEM_EVENT_SimpleBuildFiniteElements);
   PetscLogEventRegister("SimpleSetUp", 0, &MOFEM_EVENT_SimpleBuildProblem);
   PetscLogEventRegister("SimpleKSPSolve", 0, &MOFEM_EVENT_SimpleKSPSolve);
   strcpy(meshFileName, "mesh.h5m");
 }

◆ ~Simple()

MoFEM::Simple::~Simple ( )

Destructor.

Definition at line 52 of file Simple.cpp.

                 {
 
   if (dM != PETSC_NULL) {
     ierr = DMDestroy(&dM);
     CHKERRABORT(PETSC_COMM_WORLD, ierr);
   }
 }

Member Function Documentation

◆ addBoundaryField()

MoFEMErrorCode MoFEM::Simple::addBoundaryField	(	const std::string &	name,
		const FieldSpace	space,
		const FieldApproximationBase	base,
		const FieldCoefficientsNumber	nb_of_coefficients,
		const TagType	tag_type = `MB_TAG_SPARSE`,
		const enum MoFEMTypes	bh = `MF_ZERO`,
		int	verb = `-1`
	)

Add field on boundary.

Parameters

name	name of the filed
space	space (L2,H1,Hdiv,Hcurl)
base	approximation base, see FieldApproximationBase
nb_of_coefficients	number of field coefficients
tag_type	type of the tag MB_TAG_DENSE or MB_TAG_SPARSE (DENSE is faster and uses less memory, SPARSE is more flexible if you define field on subdomains)
bh	if MF_EXCL throws error if field exits, MF_ZERO no error if field exist
verb	verbosity level

Returns: error code

Examples:: hello_world.cpp, mesh_smoothing.cpp, and simple_interface.cpp.

Definition at line 139 of file Simple.cpp.

                                    {
 
   Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   ierr = m_field.add_field(name, space, base, nb_of_cooficients, tag_type, bh,
                            verb);
   CHKERRG(ierr);
   boundaryFields.push_back(name);
   MoFEMFunctionReturnHot(0);
 }

◆ addDataField()

MoFEMErrorCode MoFEM::Simple::addDataField	(	const std::string &	name,
		const FieldSpace	space,
		const FieldApproximationBase	base,
		const FieldCoefficientsNumber	nb_of_coefficients,
		const TagType	tag_type = `MB_TAG_SPARSE`,
		const enum MoFEMTypes	bh = `MF_ZERO`,
		int	verb = `-1`
	)

Add field on domain.

Parameters

name	name of the filed
space	space (L2,H1,Hdiv,Hcurl)
base	approximation base, see FieldApproximationBase
nb_of_coefficients	number of field coefficients
tag_type	type of the tag MB_TAG_DENSE or MB_TAG_SPARSE (DENSE is faster and uses less memory, SPARSE is more flexible if you define field on subdomains)
bh	if MF_EXCL throws error if field exits, MF_ZERO no error if field exist
verb	verbosity level

Returns: error code

Definition at line 171 of file Simple.cpp.

                                {
 
   Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   ierr = m_field.add_field(name, space, base, nb_of_cooficients, tag_type, bh,
                            verb);
   CHKERRG(ierr);
   dataFields.push_back(name);
   MoFEMFunctionReturnHot(0);
 }

◆ addDomainField()

MoFEMErrorCode MoFEM::Simple::addDomainField	(	const std::string &	name,
		const FieldSpace	space,
		const FieldApproximationBase	base,
		const FieldCoefficientsNumber	nb_of_coefficients,
		const TagType	tag_type = `MB_TAG_SPARSE`,
		const enum MoFEMTypes	bh = `MF_ZERO`,
		int	verb = `-1`
	)

Add field on domain.

Parameters

name	name of the filed
space	space (L2,H1,Hdiv,Hcurl)
base	approximation base, see FieldApproximationBase
nb_of_coefficients	number of field coefficients
tag_type	type of the tag MB_TAG_DENSE or MB_TAG_SPARSE (DENSE is faster and uses less memory, SPARSE is more flexible if you define field on subdomains)
bh	if MF_EXCL throws error if field exits, MF_ZERO no error if field exist
verb	verbosity level

Returns: error code

Examples:: hello_world.cpp, mesh_smoothing.cpp, simple_elasticity.cpp, simple_interface.cpp, and testing_jacobian_of_hook_element.cpp.

Definition at line 123 of file Simple.cpp.

                                  {
 
   Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   ierr = m_field.add_field(name, space, base, nb_of_cooficients, tag_type, bh,
                            verb);
   CHKERRG(ierr);
   domainFields.push_back(name);
   MoFEMFunctionReturnHot(0);
 }

◆ addSkeletonField()

MoFEMErrorCode MoFEM::Simple::addSkeletonField	(	const std::string &	name,
		const FieldSpace	space,
		const FieldApproximationBase	base,
		const FieldCoefficientsNumber	nb_of_coefficients,
		const TagType	tag_type = `MB_TAG_SPARSE`,
		const enum MoFEMTypes	bh = `MF_ZERO`,
		int	verb = `-1`
	)

Add field on skeleton.

Parameters

name	name of the filed
space	space (L2,H1,Hdiv,Hcurl)
base	approximation base, see FieldApproximationBase
nb_of_coefficients	number of field coefficients
tag_type	type of the tag MB_TAG_DENSE or MB_TAG_SPARSE (DENSE is faster and uses less memory, SPARSE is more flexible if you define field on subdomains)
bh	if MF_EXCL throws error if field exits, MF_ZERO no error if field exist
verb	verbosity level

Returns: error code

Examples:: hello_world.cpp.

Definition at line 155 of file Simple.cpp.

                                    {
 
   Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   ierr = m_field.add_field(name, space, base, nb_of_cooficients, tag_type, bh,
                            verb);
   CHKERRG(ierr);
   skeletonFields.push_back(name);
   MoFEMFunctionReturnHot(0);
 }

◆ buildFields()

MoFEMErrorCode MoFEM::Simple::buildFields ( )

Build fields.

Returns: error code

Examples:: mesh_smoothing.cpp, and simple_elasticity.cpp.

Definition at line 320 of file Simple.cpp.

                                    {
   Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   PetscLogEventBegin(MOFEM_EVENT_SimpleBuildFields, 0, 0, 0, 0);
   // take skin
   {
     Range domain_ents;
     rval = m_field.get_moab().get_entities_by_dimension(meshSet, dIm,
                                                         domain_ents, true);
     CHKERRQ_MOAB(rval);
     Skinner skin(&m_field.get_moab());
     Range domain_skin;
     rval = skin.find_skin(0, domain_ents, false, domain_skin);
     CHKERRQ_MOAB(rval);
     // filter not owned entities, those are not on boundary
     ParallelComm *pcomm =
         ParallelComm::get_pcomm(&m_field.get_moab(), MYPCOMM_INDEX);
     Range proc_domain_skin;
     rval =
         pcomm->filter_pstatus(domain_skin, PSTATUS_SHARED | PSTATUS_MULTISHARED,
                               PSTATUS_NOT, -1, &proc_domain_skin);
     CHKERRQ_MOAB(rval);
     // cerr << proc_domain_skin << endl;
     // create boundary meshset
     if (boundaryMeshset != 0) {
       rval = m_field.get_moab().delete_entities(&boundaryMeshset, 1);
       CHKERRQ_MOAB(rval);
     }
     rval = m_field.get_moab().create_meshset(MESHSET_SET, boundaryMeshset);
     CHKERRQ_MOAB(rval);
     rval = m_field.get_moab().add_entities(boundaryMeshset, proc_domain_skin);
     CHKERRQ_MOAB(rval);
     for (int dd = 0; dd != dIm - 1; dd++) {
       Range adj;
       rval = m_field.get_moab().get_adjacencies(proc_domain_skin, dd, false,
                                                 adj, moab::Interface::UNION);
       CHKERRQ_MOAB(rval);
       rval = m_field.get_moab().add_entities(boundaryMeshset, adj);
       CHKERRQ_MOAB(rval);
     }
   }
   // Add entities to the fields
   for (unsigned int ff = 0; ff != domainFields.size(); ff++) {
     ierr = m_field.add_ents_to_field_by_dim(meshSet, dIm, domainFields[ff]);
     CHKERRG(ierr);
     ierr = m_field.synchronise_field_entities(domainFields[ff], 0);
     CHKERRG(ierr);
   }
   for (unsigned int ff = 0; ff != dataFields.size(); ff++) {
     ierr = m_field.add_ents_to_field_by_dim(meshSet, dIm, dataFields[ff]);
     CHKERRG(ierr);
     ierr = m_field.synchronise_field_entities(dataFields[ff], 0);
     CHKERRG(ierr);
   }
   for (unsigned int ff = 0; ff != boundaryFields.size(); ff++) {
     ierr = m_field.add_ents_to_field_by_dim(boundaryMeshset, dIm - 1,
                                             boundaryFields[ff]);
     CHKERRG(ierr);
     ierr = m_field.synchronise_field_entities(boundaryFields[ff], 0);
     CHKERRG(ierr);
   }
   for (unsigned int ff = 0; ff != skeletonFields.size(); ff++) {
     ierr =
         m_field.add_ents_to_field_by_dim(meshSet, dIm - 1, skeletonFields[ff]);
     CHKERRG(ierr);
     ierr = m_field.synchronise_field_entities(skeletonFields[ff], 0);
     CHKERRG(ierr);
   }
   // Set order
   for (unsigned int ff = 0; ff != domainFields.size(); ff++) {
     if (fieldsOrder.find(domainFields[ff]) == fieldsOrder.end()) {
       SETERRQ1(PETSC_COMM_WORLD, MOFEM_INVALID_DATA,
                "Order for field not set %s", domainFields[ff].c_str());
     }
     int dds = 0;
     const Field *field = m_field.get_field_structure(domainFields[ff]);
     switch (field->getSpace()) {
     case L2:
       dds = dIm;
       break;
     case HDIV:
       dds = 2;
       break;
     case HCURL:
       dds = 1;
       break;
     case H1:
       dds = 1;
       break;
     default:
       SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
               "Huston we have a problem");
     }
     if (field->getSpace() == H1) {
       ierr = m_field.set_field_order(meshSet, MBVERTEX, domainFields[ff], 1);
       CHKERRG(ierr);
     }
     for (int dd = dds; dd <= dIm; dd++) {
       Range ents;
       ierr =
           m_field.get_field_entities_by_dimension(domainFields[ff], dd, ents);
       CHKERRG(ierr);
       if (!fieldsOrder.at(domainFields[ff]).second.empty()) {
         ents = intersect(ents, fieldsOrder.at(domainFields[ff]).second);
       }
       ierr = m_field.set_field_order(ents, domainFields[ff],
                                      fieldsOrder.at(domainFields[ff]).first);
       CHKERRG(ierr);
     }
   }
   // Set order to data fields
   for (unsigned int ff = 0; ff != dataFields.size(); ff++) {
     if (fieldsOrder.find(dataFields[ff]) == fieldsOrder.end()) {
       SETERRQ1(PETSC_COMM_WORLD, MOFEM_INVALID_DATA,
                "Order for field not set %s", dataFields[ff].c_str());
     }
     int dds = 0;
     const Field *field = m_field.get_field_structure(dataFields[ff]);
     switch (field->getSpace()) {
     case L2:
       dds = dIm;
       break;
     case HDIV:
       dds = 2;
       break;
     case HCURL:
       dds = 1;
       break;
     case H1:
       dds = 1;
       break;
     default:
       SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
               "Huston we have a problem");
     }
     if (field->getSpace() == H1) {
       ierr = m_field.set_field_order(meshSet, MBVERTEX, dataFields[ff], 1);
       CHKERRG(ierr);
     }
     for (int dd = dds; dd <= dIm; dd++) {
       Range ents;
       ierr = m_field.get_field_entities_by_dimension(dataFields[ff], dd, ents);
       CHKERRG(ierr);
       ierr = m_field.set_field_order(ents, dataFields[ff],
                                      fieldsOrder.at(dataFields[ff]).first);
       CHKERRG(ierr);
     }
   }
   // Set order to boundary
   for (unsigned int ff = 0; ff != boundaryFields.size(); ff++) {
     if (fieldsOrder.find(boundaryFields[ff]) == fieldsOrder.end()) {
       SETERRQ1(PETSC_COMM_WORLD, MOFEM_INVALID_DATA,
                "Order for field not set %s", boundaryFields[ff].c_str());
     }
     int dds = 0;
     const Field *field = m_field.get_field_structure(boundaryFields[ff]);
     switch (field->getSpace()) {
     case L2:
       dds = dIm - 1;
       break;
     case HDIV:
       dds = 2;
       break;
     case HCURL:
       dds = 1;
       break;
     case H1:
       dds = 1;
       break;
     default:
       SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
               "Huston we have a problem");
     }
     if (field->getSpace() == H1) {
       ierr = m_field.set_field_order(meshSet, MBVERTEX, boundaryFields[ff], 1);
       CHKERRG(ierr);
     }
     for (int dd = dds; dd <= dIm - 1; dd++) {
       Range ents;
       ierr =
           m_field.get_field_entities_by_dimension(boundaryFields[ff], dd, ents);
       CHKERRG(ierr);
       ierr = m_field.set_field_order(ents, boundaryFields[ff],
                                      fieldsOrder.at(boundaryFields[ff]).first);
       CHKERRG(ierr);
     }
   }
   // Set order to skeleton
   for (unsigned int ff = 0; ff != skeletonFields.size(); ff++) {
     if (fieldsOrder.find(skeletonFields[ff]) == fieldsOrder.end()) {
       SETERRQ1(PETSC_COMM_WORLD, MOFEM_INVALID_DATA,
                "Order for field not set %s", skeletonFields[ff].c_str());
     }
     int dds = 0;
     const Field *field = m_field.get_field_structure(skeletonFields[ff]);
     switch (field->getSpace()) {
     case L2:
       dds = dIm - 1;
       break;
     case HDIV:
       dds = 2;
       break;
     case HCURL:
       dds = 1;
       break;
     case H1:
       dds = 1;
       break;
     default:
       SETERRQ(PETSC_COMM_WORLD, MOFEM_DATA_INCONSISTENCY,
               "Huston we have a problem");
     }
     if (field->getSpace() == H1) {
       ierr = m_field.set_field_order(meshSet, MBVERTEX, skeletonFields[ff], 1);
       CHKERRG(ierr);
     }
     for (int dd = dds; dd <= dIm - 1; dd++) {
       Range ents;
       ierr =
           m_field.get_field_entities_by_dimension(skeletonFields[ff], dd, ents);
       CHKERRG(ierr);
       ierr = m_field.set_field_order(ents, skeletonFields[ff],
                                      fieldsOrder.at(skeletonFields[ff]).first);
       CHKERRG(ierr);
     }
   }
   // Build fields
   ierr = m_field.build_fields();
   CHKERRG(ierr);
   PetscLogEventEnd(MOFEM_EVENT_SimpleBuildFields, 0, 0, 0, 0);
   MoFEMFunctionReturnHot(0);
 }

◆ buildFiniteElements()

MoFEMErrorCode MoFEM::Simple::buildFiniteElements ( )

Build finite elements.

Returns: error code

Examples:: mesh_smoothing.cpp, and simple_elasticity.cpp.

Definition at line 553 of file Simple.cpp.

                                            {
   Interface &m_field = cOre;
   MoFEMFunctionBegin;
   PetscLogEventBegin(MOFEM_EVENT_SimpleBuildFiniteElements, 0, 0, 0, 0);
   // Add finite elements
   CHKERR m_field.add_ents_to_finite_element_by_dim(meshSet, dIm, domainFE,
                                                    true);
   CHKERR m_field.build_finite_elements(domainFE);
   if (!boundaryFields.empty()) {
     CHKERR m_field.add_ents_to_finite_element_by_dim(boundaryMeshset, dIm - 1,
                                                      boundaryFE, true);
     CHKERR m_field.build_finite_elements(boundaryFE);
   }
   if (!skeletonFields.empty()) {
     CHKERR m_field.add_ents_to_finite_element_by_dim(meshSet, dIm - 1,
                                                      skeletonFE, true);
     CHKERR m_field.build_finite_elements(skeletonFE);
   }
   for (std::vector<std::string>::iterator fit = otherFEs.begin();
        fit != otherFEs.end(); ++fit) {
     CHKERR m_field.build_finite_elements(*fit);
   }
   PetscLogEventEnd(MOFEM_EVENT_SimpleBuildFiniteElements, 0, 0, 0, 0);
   MoFEMFunctionReturn(0);
 }

◆ buildProblem()

MoFEMErrorCode MoFEM::Simple::buildProblem ( )

Build problem.

Returns: error code

Examples:: mesh_smoothing.cpp, and simple_elasticity.cpp.

Definition at line 579 of file Simple.cpp.

                                     {
   Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   PetscLogEventBegin(MOFEM_EVENT_SimpleBuildProblem, 0, 0, 0, 0);
   ierr = m_field.build_adjacencies(bitLevel);
   CHKERRG(ierr);
   // Set problem by the DOFs on the fields rather that by adding DOFs on the
   // elements
   ierr = m_field.getInterface<ProblemsManager>()->buildProblemFromFields =
       PETSC_TRUE;
   ierr = DMSetUp(dM);
   CHKERRG(ierr);
   PetscLogEventEnd(MOFEM_EVENT_SimpleBuildProblem, 0, 0, 0, 0);
   MoFEMFunctionReturnHot(0);
 }

◆ defineFiniteElements()

MoFEMErrorCode MoFEM::Simple::defineFiniteElements ( )

Define finite elements.

Returns: Error code

Examples:: mesh_smoothing.cpp, and simple_elasticity.cpp.

Definition at line 186 of file Simple.cpp.

                                             {
   Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   // Define finite elements
   ierr = m_field.add_finite_element(domainFE);
   CHKERRG(ierr);
   for (unsigned int ff = 0; ff != domainFields.size(); ff++) {
     ierr =
         m_field.modify_finite_element_add_field_row(domainFE, domainFields[ff]);
     CHKERRG(ierr);
     ierr =
         m_field.modify_finite_element_add_field_col(domainFE, domainFields[ff]);
     CHKERRG(ierr);
     ierr = m_field.modify_finite_element_add_field_data(domainFE,
                                                         domainFields[ff]);
     CHKERRG(ierr);
   }
   for (unsigned int ff = 0; ff != dataFields.size(); ff++) {
     ierr =
         m_field.modify_finite_element_add_field_data(domainFE, dataFields[ff]);
     CHKERRG(ierr);
   }
   if (!boundaryFields.empty()) {
     ierr = m_field.add_finite_element(boundaryFE);
     CHKERRG(ierr);
     for (unsigned int ff = 0; ff != domainFields.size(); ff++) {
       ierr = m_field.modify_finite_element_add_field_row(boundaryFE,
                                                          domainFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_col(boundaryFE,
                                                          domainFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_data(boundaryFE,
                                                           domainFields[ff]);
       CHKERRG(ierr);
     }
     for (unsigned int ff = 0; ff != boundaryFields.size(); ff++) {
       ierr = m_field.modify_finite_element_add_field_row(boundaryFE,
                                                          boundaryFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_col(boundaryFE,
                                                          boundaryFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_data(boundaryFE,
                                                           boundaryFields[ff]);
       CHKERRG(ierr);
     }
     for (unsigned int ff = 0; ff != skeletonFields.size(); ff++) {
       ierr = m_field.modify_finite_element_add_field_row(boundaryFE,
                                                          skeletonFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_col(boundaryFE,
                                                          skeletonFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_data(boundaryFE,
                                                           skeletonFields[ff]);
       CHKERRG(ierr);
     }
   }
   if (!skeletonFields.empty()) {
     ierr = m_field.add_finite_element(skeletonFE);
     CHKERRG(ierr);
     for (unsigned int ff = 0; ff != domainFields.size(); ff++) {
       ierr = m_field.modify_finite_element_add_field_row(skeletonFE,
                                                          domainFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_col(skeletonFE,
                                                          domainFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_data(skeletonFE,
                                                           domainFields[ff]);
       CHKERRG(ierr);
     }
     for (unsigned int ff = 0; ff != boundaryFields.size(); ff++) {
       ierr = m_field.modify_finite_element_add_field_row(skeletonFE,
                                                          boundaryFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_col(skeletonFE,
                                                          boundaryFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_data(skeletonFE,
                                                           boundaryFields[ff]);
       CHKERRG(ierr);
     }
     for (unsigned int ff = 0; ff != skeletonFields.size(); ff++) {
       ierr = m_field.modify_finite_element_add_field_row(skeletonFE,
                                                          skeletonFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_col(skeletonFE,
                                                          skeletonFields[ff]);
       CHKERRG(ierr);
       ierr = m_field.modify_finite_element_add_field_data(skeletonFE,
                                                           skeletonFields[ff]);
       CHKERRG(ierr);
     }
   }
   MoFEMFunctionReturnHot(0);
 }

◆ defineProblem()

MoFEMErrorCode MoFEM::Simple::defineProblem ( const PetscBool is_partitioned = PETSC_TRUE )

define problem

Returns: error code

Examples:: mesh_smoothing.cpp, and simple_elasticity.cpp.

Definition at line 285 of file Simple.cpp.

                                                                    {
   Interface &m_field = cOre;
   MoFEMFunctionBegin;
   if (dM != PETSC_NULL) {
     CHKERR DMDestroy(&dM);
   }
   // Create dm instance
   CHKERR DMCreate(m_field.get_comm(), &dM);
   CHKERR DMSetType(dM, "DMMOFEM");
   // set dm data structure which created mofem data structures
   CHKERR DMMoFEMCreateMoFEM(dM, &m_field, "SimpleProblem", bitLevel);
   CHKERR DMSetFromOptions(dM);
   CHKERR DMMoFEMAddElement(dM, domainFE.c_str());
   if (!boundaryFields.empty()) {
     CHKERR DMMoFEMAddElement(dM, boundaryFE.c_str());
   }
   if (!skeletonFields.empty()) {
     CHKERR DMMoFEMAddElement(dM, skeletonFE.c_str());
   }
   for (std::vector<std::string>::iterator fit = otherFEs.begin();
        fit != otherFEs.end(); ++fit) {
     CHKERR DMMoFEMAddElement(dM, fit->c_str());
   }
   CHKERR DMMoFEMSetIsPartitioned(dM, is_partitioned);
   MoFEMFunctionReturn(0);
 }

◆ getBoundaryFEName() [1/2]

const std::string& MoFEM::Simple::getBoundaryFEName ( ) const

Examples:: hello_world.cpp, mesh_smoothing.cpp, and simple_interface.cpp.

Definition at line 198 of file Simple.hpp.

198 { return boundaryFE; }

MoFEM::Simple::boundaryFE

std::string boundaryFE

boundary finite element

Definition: Simple.hpp:230

◆ getBoundaryFEName() [2/2]

std::string& MoFEM::Simple::getBoundaryFEName ( )

Definition at line 202 of file Simple.hpp.

202 { return boundaryFE; }

MoFEM::Simple::boundaryFE

std::string boundaryFE

boundary finite element

Definition: Simple.hpp:230

◆ getDim()

int MoFEM::Simple::getDim ( ) const

Examples:: simple_elasticity.cpp.

Definition at line 196 of file Simple.hpp.

196 { return dIm; }

MoFEM::Simple::dIm

int dIm

dimension of problem

Definition: Simple.hpp:236

◆ getDM()

MoFEMErrorCode MoFEM::Simple::getDM ( DM * dm )

Get DM.

Parameters

dm	discrete manager

Returns: error code

Examples:: hello_world.cpp, mesh_smoothing.cpp, simple_elasticity.cpp, simple_interface.cpp, and testing_jacobian_of_hook_element.cpp.

Definition at line 605 of file Simple.cpp.

                                    {
   MoFEMFunctionBeginHot;
   PetscObjectReference((PetscObject)dM);
   *dm = dM;
   MoFEMFunctionReturnHot(0);
 }

◆ getDomainFEName() [1/2]

const std::string& MoFEM::Simple::getDomainFEName ( ) const

Examples:: hello_world.cpp, mesh_smoothing.cpp, simple_elasticity.cpp, simple_interface.cpp, and testing_jacobian_of_hook_element.cpp.

Definition at line 197 of file Simple.hpp.

197 { return domainFE; }

MoFEM::Simple::domainFE

std::string domainFE

domain finite element

Definition: Simple.hpp:229

◆ getDomainFEName() [2/2]

std::string& MoFEM::Simple::getDomainFEName ( )

Definition at line 201 of file Simple.hpp.

201 { return domainFE; }

MoFEM::Simple::domainFE

std::string domainFE

domain finite element

Definition: Simple.hpp:229

◆ getOptions()

MoFEMErrorCode MoFEM::Simple::getOptions ( )

get options

Returns: error code

Examples:: hello_world.cpp, mesh_smoothing.cpp, simple_elasticity.cpp, simple_interface.cpp, and testing_jacobian_of_hook_element.cpp.

Definition at line 60 of file Simple.cpp.

                                   {
 
   PetscBool flg = PETSC_TRUE;
   MoFEMFunctionBeginHot;
   ierr = PetscOptionsBegin(PETSC_COMM_WORLD, "", "Simple interface options",
                            "none");
   CHKERRG(ierr);
   ierr = PetscOptionsString("-file_name", "file name", "", "mesh.h5m",
                             meshFileName, 255, &flg);
   CHKERRG(ierr);
   ierr = PetscOptionsEnd();
   CHKERRG(ierr);
   MoFEMFunctionReturnHot(0);
 }

◆ getOtherFiniteElements()

std::vector<std::string>& MoFEM::Simple::getOtherFiniteElements ( )

Examples:: mesh_smoothing.cpp.

Definition at line 205 of file Simple.hpp.

                                                           {
       return otherFEs;
     }

◆ getSkeletonFEName() [1/2]

const std::string& MoFEM::Simple::getSkeletonFEName ( ) const

Examples:: hello_world.cpp.

Definition at line 199 of file Simple.hpp.

199 { return skeletonFE; }

MoFEM::Simple::skeletonFE

std::string skeletonFE

skeleton finite element

Definition: Simple.hpp:231

◆ getSkeletonFEName() [2/2]

std::string& MoFEM::Simple::getSkeletonFEName ( )

Definition at line 203 of file Simple.hpp.

203 { return skeletonFE; }

MoFEM::Simple::skeletonFE

std::string skeletonFE

skeleton finite element

Definition: Simple.hpp:231

◆ loadFile()

MoFEMErrorCode MoFEM::Simple::loadFile ( )

Load mesh file.

Parameters

field_name file name

Returns: error code

Examples:: hello_world.cpp, mesh_smoothing.cpp, simple_elasticity.cpp, simple_interface.cpp, and testing_jacobian_of_hook_element.cpp.

Definition at line 75 of file Simple.cpp.

                                 {
   Interface &m_field = cOre;
   MoFEMFunctionBeginHot;
   PetscLogEventBegin(MOFEM_EVENT_SimpleLoadMesh, 0, 0, 0, 0);
   // This is a case of distributed mesh and algebra. In that case each processor
   // keep only part of the problem.
   const char *option;
   option = "PARALLEL=READ_PART;"
            "PARALLEL_RESOLVE_SHARED_ENTS;"
            "PARTITION=PARALLEL_PARTITION;";
   rval = m_field.get_moab().load_file(meshFileName, 0, option);
   CHKERRQ_MOAB(rval);
   ierr = m_field.rebuild_database();
   CHKERRG(ierr);
   // determine problem dimension
   if (dIm == -1) {
     int nb_ents_3d;
     rval = m_field.get_moab().get_number_entities_by_dimension(
         meshSet, 3, nb_ents_3d, true);
     CHKERRQ_MOAB(rval);
     if (nb_ents_3d > 0) {
       dIm = 3;
     } else {
       int nb_ents_2d;
       rval = m_field.get_moab().get_number_entities_by_dimension(
           meshSet, 2, nb_ents_2d, true);
       CHKERRQ_MOAB(rval);
       if (nb_ents_2d > 0) {
         dIm = 2;
       } else {
         dIm = 1;
       }
     }
   }
   Range ents;
   ierr = m_field.get_moab().get_entities_by_dimension(meshSet, dIm, ents, true);
   ierr = m_field.getInterface<BitRefManager>()->setBitRefLevel(ents, bitLevel,
                                                                false);
   CHKERRG(ierr);
   ParallelComm *pcomm =
       ParallelComm::get_pcomm(&m_field.get_moab(), MYPCOMM_INDEX);
   if (pcomm == NULL)
     pcomm = new ParallelComm(&m_field.get_moab(), m_field.get_comm());
   PetscLogEventEnd(MOFEM_EVENT_SimpleLoadMesh, 0, 0, 0, 0);
   MoFEMFunctionReturnHot(0);
 }

◆ query_interface()

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

virtual

Implements MoFEM::UnknownInterface.

Definition at line 28 of file Simple.cpp.

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

◆ setFieldOrder()

MoFEMErrorCode MoFEM::Simple::setFieldOrder	(	const std::string	field_name,
		const int	order,
		const Range *	ents = `NULL`
	)

Set field order.

Parameters

std::field_name	field name
order	order
range	of entities to which order is set (If null it sat to all entities)

Returns: error code

Examples:: hello_world.cpp, mesh_smoothing.cpp, simple_elasticity.cpp, simple_interface.cpp, and testing_jacobian_of_hook_element.cpp.

Definition at line 312 of file Simple.cpp.

                                                                          {
   MoFEMFunctionBeginHot;
   fieldsOrder[field_name] =
       std::pair<int, Range>(order, ents == NULL ? Range() : Range(*ents));
   MoFEMFunctionReturnHot(0);
 }

◆ setUp()

MoFEMErrorCode MoFEM::Simple::setUp ( )

Setup problem.

Returns: error code

Examples:: hello_world.cpp, simple_interface.cpp, and testing_jacobian_of_hook_element.cpp.

Definition at line 595 of file Simple.cpp.

                              {
   MoFEMFunctionBegin;
   CHKERR defineFiniteElements();
   CHKERR defineProblem();
   CHKERR buildFields();
   CHKERR buildFiniteElements();
   CHKERR buildProblem();
   MoFEMFunctionReturn(0);
 }

Member Data Documentation

◆ bitLevel

const BitRefLevel MoFEM::Simple::bitLevel

private

Definition at line 211 of file Simple.hpp.

◆ boundaryFE

std::string MoFEM::Simple::boundaryFE

private

boundary finite element

Definition at line 230 of file Simple.hpp.

◆ boundaryFields

std::vector<std::string> MoFEM::Simple::boundaryFields

private

boundary fields

Definition at line 223 of file Simple.hpp.

◆ boundaryMeshset

EntityHandle MoFEM::Simple::boundaryMeshset

private

meshset with boundary

Definition at line 220 of file Simple.hpp.

◆ cOre

MoFEM::Core& MoFEM::Simple::cOre

Definition at line 37 of file Simple.hpp.

◆ dataFields

std::vector<std::string> MoFEM::Simple::dataFields

private

Data fields.

Definition at line 225 of file Simple.hpp.

◆ dIm

int MoFEM::Simple::dIm

private

dimension of problem

Definition at line 236 of file Simple.hpp.

◆ dM

DM MoFEM::Simple::dM

private

Definition at line 238 of file Simple.hpp.

◆ domainFE

std::string MoFEM::Simple::domainFE

private

domain finite element

Definition at line 229 of file Simple.hpp.

◆ domainFields

std::vector<std::string> MoFEM::Simple::domainFields

private

domain fields

Definition at line 222 of file Simple.hpp.

◆ fieldsOrder

std::map<std::string,std::pair<int,Range> > MoFEM::Simple::fieldsOrder

private

fields order

Definition at line 227 of file Simple.hpp.

◆ meshFileName

char MoFEM::Simple::meshFileName[255]

private

mesh file name

Definition at line 235 of file Simple.hpp.

◆ meshSet

EntityHandle MoFEM::Simple::meshSet

private

domain meshset

Definition at line 219 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleBuildFields

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleBuildFields

private

Definition at line 214 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleBuildFiniteElements

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleBuildFiniteElements

private

Definition at line 215 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleBuildProblem

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleBuildProblem

private

Definition at line 216 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleKSPSolve

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleKSPSolve

private

Definition at line 217 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleLoadMesh

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleLoadMesh

private

Definition at line 213 of file Simple.hpp.

◆ nameOfProblem

std::string MoFEM::Simple::nameOfProblem

private

problem name

Definition at line 221 of file Simple.hpp.

◆ otherFEs

std::vector<std::string> MoFEM::Simple::otherFEs

private

Other finite elements.

Definition at line 233 of file Simple.hpp.

◆ skeletonFE

std::string MoFEM::Simple::skeletonFE

private

skeleton finite element

Definition at line 231 of file Simple.hpp.

◆ skeletonFields

std::vector<std::string> MoFEM::Simple::skeletonFields

private

fields on the skeleton

Definition at line 224 of file Simple.hpp.

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

src/interfaces/Simple.hpp
src/interfaces/impl/Simple.cpp

Public Member Functions

Public Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Simple()

◆ ~Simple()

Member Function Documentation

◆ addBoundaryField()

◆ addDataField()

◆ addDomainField()

◆ addSkeletonField()

◆ buildFields()

◆ buildFiniteElements()

◆ buildProblem()

◆ defineFiniteElements()

◆ defineProblem()

◆ getBoundaryFEName() [1/2]

◆ getBoundaryFEName() [2/2]

◆ getDim()

◆ getDM()

◆ getDomainFEName() [1/2]

◆ getDomainFEName() [2/2]

◆ getOptions()

◆ getOtherFiniteElements()

◆ getSkeletonFEName() [1/2]

◆ getSkeletonFEName() [2/2]

◆ loadFile()

◆ query_interface()

◆ setFieldOrder()

◆ setUp()

Member Data Documentation

◆ bitLevel

◆ boundaryFE

◆ boundaryFields

◆ boundaryMeshset

◆ cOre

◆ dataFields

◆ dIm

◆ dM

◆ domainFE

◆ domainFields

◆ fieldsOrder

◆ meshFileName

◆ meshSet

◆ MOFEM_EVENT_SimpleBuildFields

◆ MOFEM_EVENT_SimpleBuildFiniteElements

◆ MOFEM_EVENT_SimpleBuildProblem

◆ MOFEM_EVENT_SimpleKSPSolve

◆ MOFEM_EVENT_SimpleLoadMesh

◆ nameOfProblem

◆ otherFEs

◆ skeletonFE

◆ skeletonFields