MoFEM::Simple Struct Reference

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

#include "src/interfaces/Simple.hpp"

Inheritance diagram for MoFEM::Simple:

Collaboration diagram for MoFEM::Simple:

Public Types
typedef boost::function< MoFEMErrorCode(Interface &, const char *, const char *)>	LoadFileFunc

Public Member Functions
MoFEMErrorCode	query_interface (boost::typeindex::type_index type_index, UnknownInterface **iface) const
	Simple (const MoFEM::Core &core)
virtual	~Simple ()=default
MoFEMErrorCode	getOptions ()
	get options
MoFEMErrorCode	loadFile (const std::string options, const std::string mesh_file_name, LoadFileFunc loadFunc=defaultLoadFileFunc)
	Load mesh file.
MoFEMErrorCode	loadFile (const std::string mesh_file_name="")
	Load mesh file with parallel options if number of cores > 1.
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.
MoFEMErrorCode	addDomainBrokenField (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 broken field on domain.
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.
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.
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 data field.
MoFEMErrorCode	addMeshsetField (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 meshset field.
MoFEMErrorCode	removeDomainField (const std::string name)
	Remove field from domain.
MoFEMErrorCode	removeBoundaryField (const std::string name)
	Remove field from boundary.
MoFEMErrorCode	removeSkeletonField (const std::string name)
	Remove field from skeleton.
MoFEMErrorCode	defineFiniteElements ()
	Define finite elements.
MoFEMErrorCode	defineProblem (const PetscBool is_partitioned=PETSC_TRUE)
	define problem
MoFEMErrorCode	setFieldOrder (const std::string field_name, const int order, const Range *ents=NULL)
	Set field order.
MoFEMErrorCode	buildFields ()
	Build fields.
MoFEMErrorCode	buildFiniteElements ()
	Build finite elements.
MoFEMErrorCode	setSkeletonAdjacency (int dim=-1, std::string fe_name="")
	Set the skeleton adjacency object.
MoFEMErrorCode	buildProblem ()
	Build problem.
MoFEMErrorCode	setUp (const PetscBool is_partitioned=PETSC_TRUE)
	Setup problem.
MoFEMErrorCode	reSetUp (bool only_dm=false)
	Rebuild internal MoFEM data structures.
MoFEMErrorCode	getDM (DM *dm)
	Get DM.
SmartPetscObj< DM >	getDM ()
	Return smart DM object.
int	getDim () const
	Get the problem dimension.
void	setDim (int dim)
	Set the problem dimension.
DEPRECATED EntityHandle &	getMeshSet ()
EntityHandle &	getMeshset ()
	Get the MeshSet object.
EntityHandle &	getBoundaryMeshSet ()
	Get the BoundaryMeshSet object.
EntityHandle &	getSkeletonMeshSet ()
	Get the SkeletonMeshSet object.
BitRefLevel &	getBitRefLevel ()
	Get the BitRefLevel.
BitRefLevel &	getBitRefLevelMask ()
	Get the BitRefLevelMask.
const std::string	getDomainFEName () const
	Get the Domain FE Name.
const std::string	getBoundaryFEName () const
	Get the Boundary FE Name.
const std::string	getSkeletonFEName () const
	Get the Skeleton FE Name.
const std::string	getProblemName () const
	Get the Problem Name.
std::string &	getDomainFEName ()
	Get the Domain FE Name.
std::string &	getBoundaryFEName ()
	Get the Boundary FE Name.
std::string &	getSkeletonFEName ()
	Get the Skeleton FE Name.
std::string &	getMeshsetFEName ()
	Get the Skeleton FE Name.
auto &	getMeshsetFiniteElementEntities ()
	Range of entities added to meshset finite element.
std::string &	getProblemName ()
	Get the Problem Name.
std::vector< std::string > &	getOtherFiniteElements ()
	Get the Other Finite Elements.
MoFEMErrorCode	deleteDM ()
	Delete dm.
MoFEMErrorCode	deleteFiniteElements ()
	Delete finite elements.
bool &	getAddSkeletonFE ()
	Get the addSkeletonFE flag.
bool &	getAddBoundaryFE ()
	Get the addBoundaryFE flag.
bool &	getParentAdjacencies ()
	Get the addParentAdjacencies flag.
auto &	getBitAdjParent ()
	Get bit ref level for parent.
auto &	getBitAdjParentMask ()
	Get bit ref level mask for parent.
auto &	getBitAdjEnt ()
	Get bit ref level for entities.
auto &	getBitAdjEntMask ()
	Get bit ref level mask for entities.
MoFEMErrorCode	addFieldToEmptyFieldBlocks (const std::string row_field, const std::string col_field) const
	Add empty block to problem.
Public Member Functions inherited from MoFEM::UnknownInterface
template<class IFACE >
MoFEMErrorCode	registerInterface (bool error_if_registration_failed=true)
	Register interface.
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE *&iface) const
	Get interface reference to pointer of interface.
template<class IFACE >
MoFEMErrorCode	getInterface (IFACE **const iface) const
	Get interface pointer to pointer of interface.
template<class IFACE , typename boost::enable_if< boost::is_pointer< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get interface pointer to pointer of interface.
template<class IFACE , typename boost::enable_if< boost::is_reference< IFACE >, int >::type = 0>
IFACE	getInterface () const
	Get reference to interface.
template<class IFACE >
IFACE *	getInterface () const
	Function returning pointer to interface.
virtual	~UnknownInterface ()=default

Static Public Member Functions
static MoFEMErrorCode	defaultLoadFileFunc (Interface &m_field, const char *file_name, const char *options)
	Set function for loading mesh file.
Static Public Member Functions inherited from MoFEM::UnknownInterface
static MoFEMErrorCode	getLibVersion (Version &version)
	Get library version.
static MoFEMErrorCode	getFileVersion (moab::Interface &moab, Version &version)
	Get database major version.
static MoFEMErrorCode	setFileVersion (moab::Interface &moab, Version version=Version(MoFEM_VERSION_MAJOR, MoFEM_VERSION_MINOR, MoFEM_VERSION_BUILD))
	Get database major version.
static MoFEMErrorCode	getInterfaceVersion (Version &version)
	Get database major version.

Private Member Functions
MoFEMErrorCode	createBoundaryMeshset ()
MoFEMErrorCode	createSkeletonMeshset ()
MoFEMErrorCode	exchangeGhostCells ()
template<int DIM = -1>
MoFEMErrorCode	setSkeletonAdjacency (std::string fe_name)
template<int DIM = -1>
MoFEMErrorCode	setParentAdjacency ()
template<>
MoFEMErrorCode	setSkeletonAdjacency (std::string fe_name)
template<>
MoFEMErrorCode	setSkeletonAdjacency (std::string fe_name)
template<>
MoFEMErrorCode	setSkeletonAdjacency (std::string fe_name)
template<>
MoFEMErrorCode	setParentAdjacency ()
template<>
MoFEMErrorCode	setParentAdjacency ()
template<>
MoFEMErrorCode	setParentAdjacency ()

Private Attributes
MoFEM::Core &	cOre
BitRefLevel	bitLevel
	BitRefLevel of the problem.
BitRefLevel	bitLevelMask
	BitRefLevel of the problem.
PetscLogEvent	MOFEM_EVENT_SimpleSetUP
PetscLogEvent	MOFEM_EVENT_SimpleLoadMesh
PetscLogEvent	MOFEM_EVENT_SimpleBuildFields
PetscLogEvent	MOFEM_EVENT_SimpleBuildFiniteElements
PetscLogEvent	MOFEM_EVENT_SimpleBuildProblem
PetscLogEvent	MOFEM_EVENT_SimpleKSPSolve
EntityHandle	meshSet
	domain meshset
EntityHandle	boundaryMeshset
	meshset with boundary
EntityHandle	skeletonMeshset
	skeleton meshset with boundary
bool	addSkeletonFE
	Add skeleton FE.
bool	addBoundaryFE
	Add boundary FE.
bool	addParentAdjacencies
	If set to true parent adjacencies are build.
BitRefLevel	bitAdjParent
	bit ref level for parent
BitRefLevel	bitAdjParentMask
	bit ref level for parent parent
BitRefLevel	bitAdjEnt
	bit ref level for parent
BitRefLevel	bitAdjEntMask
	bit ref level for parent parent
std::vector< std::string >	domainFields
	domain fields
std::vector< std::string >	boundaryFields
	boundary fields
std::vector< std::string >	skeletonFields
	fields on the skeleton
std::vector< std::string >	dataFields
	Data fields.
std::vector< std::string >	meshsetFields
	meshset fields
std::vector< std::string >	noFieldFields
	NOFIELD field name.
std::vector< std::string >	noFieldDataFields
	NOFIELD field name.
std::list< std::tuple< std::string, int, Range, bool > >	fieldsOrder
std::string	nameOfProblem = "SimpleProblem"
	problem name
std::string	domainFE = "dFE"
	domain finite element
std::string	boundaryFE = "bFE"
	boundary finite element
std::string	skeletonFE = "sFE"
	skeleton finite element
std::string	meshsetFE = "mFE"
	meshset finite element
std::vector< std::string >	otherFEs
	Other finite elements.
std::vector< Range >	meshsetFiniteElementEntities
	Meshset element entities.
char	meshFileName [255]
	mesh file name
int	dIm
	dimension of problem
SmartPetscObj< DM >	dM
	Discrete manager (interface to PETSc/MoFEM functions)
boost::shared_ptr< ParentFiniteElementAdjacencyFunction< 3 > >	parentAdjFunctionDim3
boost::shared_ptr< ParentFiniteElementAdjacencyFunction< 2 > >	parentAdjFunctionDim2
boost::shared_ptr< ParentFiniteElementAdjacencyFunction< 1 > >	parentAdjFunctionDim1
boost::shared_ptr< ParentFiniteElementAdjacencyFunctionSkeleton< 2 > >	parentAdjSkeletonFunctionDim2
boost::shared_ptr< ParentFiniteElementAdjacencyFunctionSkeleton< 1 > >	parentAdjSkeletonFunctionDim1

Detailed Description

Simple interface for fast problem set-up.

Examples

mofem/atom_tests/boundary_marker.cpp, mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/field_evaluator.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_check_approx_in_2d.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_check_approx_in_3d.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/atom_tests/loop_entities.cpp, mofem/atom_tests/operators_tests.cpp, mofem/atom_tests/partition_mesh.cpp, mofem/atom_tests/scalar_check_approximation.cpp, mofem/atom_tests/schur_test_diag_mat.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/atom_tests/tensor_divergence_operator.cpp, mofem/atom_tests/test_broken_space.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tools/mesh_smoothing.cpp, mofem/tutorials/adv-0/plastic.cpp, mofem/tutorials/adv-2/thermo_elastic.cpp, mofem/tutorials/adv-3/level_set.cpp, mofem/tutorials/adv-4/dynamic_first_order_con_law.cpp, mofem/tutorials/adv-5/seepage.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/fun-0/hello_world.cpp, mofem/tutorials/fun-1/integration.cpp, mofem/tutorials/fun-2/plot_base.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/mix-1/phase.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, mofem/tutorials/scl-10/initial_diffusion.cpp, mofem/tutorials/scl-10/photon_diffusion.cpp, mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-6/heat_equation.cpp, mofem/tutorials/scl-7/wave_equation.cpp, mofem/tutorials/scl-8/radiation.cpp, mofem/tutorials/scl-9/heat_method.cpp, mofem/tutorials/vec-1/eigen_elastic.cpp, mofem/tutorials/vec-10/schur_elastic.cpp, mofem/tutorials/vec-3/nonlinear_dynamic_elastic.cpp, mofem/tutorials/vec-4/approx_sphere.cpp, mofem/tutorials/vec-4/shallow_wave.cpp, mofem/tutorials/vec-5/free_surface.cpp, mofem/tutorials/vec-6/plate.cpp, mofem/tutorials/vec-7/adjoint.cpp, mofem/users_modules/adolc-plasticity/adolc_plasticity.cpp, mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_element.cpp, and mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 27 of file Simple.hpp.

Member Typedef Documentation

LoadFileFunc

typedef boost::function<MoFEMErrorCode(Interface &, const char *, const char *)> MoFEM::Simple::LoadFileFunc

Definition at line 46 of file Simple.hpp.

Constructor & Destructor Documentation

Simple()

MoFEM::Simple::Simple

(

const MoFEM::Core &

core

)

Definition at line 162 of file Simple.cpp.

    : cOre(const_cast<Core &>(core)), bitLevel(BitRefLevel().set(0)),
      bitLevelMask(BitRefLevel().set()), meshSet(0), boundaryMeshset(0),
      skeletonMeshset(0), nameOfProblem("SimpleProblem"), domainFE("dFE"),
      boundaryFE("bFE"), skeletonFE("sFE"), dIm(-1), addSkeletonFE(false),
      addBoundaryFE(false), addParentAdjacencies(false),
      bitAdjParent(BitRefLevel().set()), bitAdjParentMask(BitRefLevel().set()),
      bitAdjEnt(BitRefLevel().set()), bitAdjEntMask(BitRefLevel().set()) {
  PetscLogEventRegister("SimpSetUp", 0, &MOFEM_EVENT_SimpleSetUP);
  PetscLogEventRegister("SimpLoadMesh", 0, &MOFEM_EVENT_SimpleLoadMesh);
  PetscLogEventRegister("SimpBuildFields", 0, &MOFEM_EVENT_SimpleBuildFields);
  PetscLogEventRegister("SimpBuildFEs", 0,
                        &MOFEM_EVENT_SimpleBuildFiniteElements);
  PetscLogEventRegister("SimpSetUp", 0, &MOFEM_EVENT_SimpleBuildProblem);
  PetscLogEventRegister("SimpKSPSolve", 0, &MOFEM_EVENT_SimpleKSPSolve);
  strcpy(meshFileName, "mesh.h5m");
}

~Simple()

virtual MoFEM::Simple::~Simple ( )

virtualdefault

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.

Adds a field to the boundary entities (entities on the mesh boundary)

Parameters

name	Name of the field
space	Function space (L2, H1, Hdiv, Hcurl)
base	Approximation base (see FieldApproximationBase)
nb_of_coefficients	Number of field coefficients per DOF
tag_type	Tag storage type (MB_TAG_DENSE is faster, MB_TAG_SPARSE more flexible)
bh	Behavior if field exists (MF_EXCL throws error, MF_ZERO ignores)
verb	Verbosity level (-1 for default)

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/simple_interface.cpp, mofem/tools/mesh_smoothing.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/mix-1/phase.cpp, and mofem/tutorials/scl-12/electrostatics.cpp.

Definition at line 355 of file Simple.cpp.

                                   {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.add_field(name, space, base, nb_of_coefficients, tag_type, bh,
                           verb);
  boundaryFields.push_back(name);
  if (space == NOFIELD)
    SETERRQ(
        PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
        "NOFIELD space for boundary field not implemented in Simple interface");
  MoFEMFunctionReturn(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 data field.

Adds a field for storing data values (typically for material properties or constants)

Parameters

name	Name of the field
space	Function space (L2, H1, Hdiv, Hcurl)
base	Approximation base (see FieldApproximationBase)
nb_of_coefficients	Number of field coefficients per DOF
tag_type	Tag storage type (MB_TAG_DENSE is faster, MB_TAG_SPARSE more flexible)
bh	Behavior if field exists (MF_EXCL throws error, MF_ZERO ignores)
verb	Verbosity level (-1 for default)

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, and mofem/tutorials/scl-12/electrostatics.cpp.

Definition at line 393 of file Simple.cpp.

                               {
 
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.add_field(name, space, base, nb_of_coefficients, tag_type, bh,
                           verb);
  if (space == NOFIELD)
    noFieldDataFields.push_back(name);
  else
    dataFields.push_back(name);
  MoFEMFunctionReturn(0);
}

addDomainBrokenField()

MoFEMErrorCode MoFEM::Simple::addDomainBrokenField	(	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 broken field on domain.

Adds a discontinuous field to the domain entities with element-wise DOF mapping

Parameters

name	Name of the field
space	Function space (L2, H1, Hdiv, Hcurl)
base	Approximation base (see FieldApproximationBase)
nb_of_coefficients	Number of field coefficients per DOF
tag_type	Tag storage type (MB_TAG_DENSE is faster, MB_TAG_SPARSE more flexible)
bh	Behavior if field exists (MF_EXCL throws error, MF_ZERO ignores)
verb	Verbosity level (-1 for default)

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/tutorials/fun-2/plot_base.cpp.

Definition at line 282 of file Simple.cpp.

                                       {
 
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  auto get_side_map_hcurl = [&]() {
    return std::vector<
 
        std::pair<EntityType,
                  std::function<MoFEMErrorCode(BaseFunction::DofsSideMap &)>
 
                  >>{
 
        {MBTRI,
         [&](BaseFunction::DofsSideMap &dofs_side_map) -> MoFEMErrorCode {
           return TriPolynomialBase::setDofsSideMap(space, DISCONTINUOUS, base,
                                                    dofs_side_map);
         }}
 
    };
  };
 
  auto get_side_map_hdiv = [&]() {
    return std::vector<
 
        std::pair<EntityType,
                  std::function<MoFEMErrorCode(BaseFunction::DofsSideMap &)>
 
                  >>{
 
        {MBTET,
         [&](BaseFunction::DofsSideMap &dofs_side_map) -> MoFEMErrorCode {
           return TetPolynomialBase::setDofsSideMap(space, DISCONTINUOUS, base,
                                                    dofs_side_map);
         }}
 
    };
  };
 
  auto get_side_map = [&](auto space) {
    switch (space) {
    case HCURL:
      return get_side_map_hcurl();
    case HDIV:
      return get_side_map_hdiv();
    default:
      MOFEM_LOG_CHANNEL("WORLD");
      MOFEM_LOG("WORLD", Sev::warning)
          << "Side dof map for broken space <" << space << "> not implemented";
    }
    return std::vector<
 
        std::pair<EntityType,
                  std::function<MoFEMErrorCode(BaseFunction::DofsSideMap &)>
 
                  >>{};
  };
 
  CHKERR m_field.add_broken_field(name, space, base, nb_of_coefficients,
                                  get_side_map(space), tag_type, bh, verb);
  domainFields.push_back(name);
  if (space == NOFIELD)
    SETERRQ(
        PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
        "NOFIELD space for boundary field not implemented in Simple interface");
  MoFEMFunctionReturn(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.

Adds a field to the domain entities (highest dimension entities in the mesh)

Parameters

name	Name of the field
space	Function space (L2, H1, Hdiv, Hcurl)
base	Approximation base (see FieldApproximationBase)
nb_of_coefficients	Number of field coefficients per DOF
tag_type	Tag storage type (MB_TAG_DENSE is faster, MB_TAG_SPARSE more flexible)
bh	Behavior if field exists (MF_EXCL throws error, MF_ZERO ignores)
verb	Verbosity level (-1 for default)

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/field_evaluator.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_check_approx_in_2d.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/atom_tests/loop_entities.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tools/mesh_smoothing.cpp, mofem/tutorials/adv-2/thermo_elastic.cpp, mofem/tutorials/adv-3/level_set.cpp, mofem/tutorials/adv-5/seepage.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/fun-2/plot_base.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/mix-1/phase.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, mofem/tutorials/scl-10/initial_diffusion.cpp, mofem/tutorials/scl-10/photon_diffusion.cpp, mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-12/electrostatics.cpp, mofem/tutorials/scl-6/heat_equation.cpp, mofem/tutorials/scl-7/wave_equation.cpp, mofem/tutorials/scl-8/radiation.cpp, mofem/tutorials/scl-9/heat_method.cpp, mofem/tutorials/vec-1/eigen_elastic.cpp, mofem/tutorials/vec-3/nonlinear_dynamic_elastic.cpp, mofem/tutorials/vec-4/approx_sphere.cpp, mofem/tutorials/vec-4/shallow_wave.cpp, mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_element.cpp, and mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 261 of file Simple.cpp.

                                 {
 
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.add_field(name, space, base, nb_of_coefficients, tag_type, bh,
                           verb);
  domainFields.push_back(name);
  if (space == NOFIELD)
    SETERRQ(
        PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
        "NOFIELD space for boundary field not implemented in Simple interface");
 
 
  MoFEMFunctionReturn(0);
}

addFieldToEmptyFieldBlocks()

MoFEMErrorCode MoFEM::Simple::addFieldToEmptyFieldBlocks	(	const std::string	row_field,
		const std::string	col_field
	)		const

Add empty block to problem.

MatrixManager assumes that all blocks, i.e. all field combinations are non-zero. This is not always the case, to optimize code and reduce memory usage user can specify which blocks are empty.

Parameters

row_field	Name of the row field
col_field	Name of the column field

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/tutorials/mix-1/phase.cpp.

Definition at line 835 of file Simple.cpp.

                                                                    {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.getInterface<ProblemsManager>()->addFieldToEmptyFieldBlocks(
      getProblemName(), row_field, col_field);
  MoFEMFunctionReturn(0);
}

addMeshsetField()

MoFEMErrorCode MoFEM::Simple::addMeshsetField	(	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 meshset field.

Add fields managed by meshset finite element for global DOFs

Parameters

name	Name of the field
space	Function space (L2, H1, Hdiv, Hcurl)
base	Approximation base (see FieldApproximationBase)
nb_of_coefficients	Number of field coefficients per DOF
tag_type	Tag storage type (MB_TAG_DENSE is faster, MB_TAG_SPARSE more flexible)
bh	Behavior if field exists (MF_EXCL throws error, MF_ZERO ignores)
verb	Verbosity level (-1 for default)

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/simple_interface.cpp, and mofem/atom_tests/simple_l2_only.cpp.

Definition at line 411 of file Simple.cpp.

                                  {
 
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.add_field(name, space, base, nb_of_coefficients, tag_type, bh,
                           verb);
  if (space != NOFIELD) {
    meshsetFields.push_back(name);
  } else {
    noFieldFields.push_back(name);
  }
  MoFEMFunctionReturn(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.

Adds a field to the skeleton entities (internal interfaces between domain elements)

Parameters

name	Name of the field
space	Function space (L2, H1, Hdiv, Hcurl)
base	Approximation base (see FieldApproximationBase)
nb_of_coefficients	Number of field coefficients per DOF
tag_type	Tag storage type (MB_TAG_DENSE is faster, MB_TAG_SPARSE more flexible)
bh	Behavior if field exists (MF_EXCL throws error, MF_ZERO ignores)
verb	Verbosity level (-1 for default)

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, and mofem/atom_tests/simple_interface.cpp.

Definition at line 373 of file Simple.cpp.

                                   {
 
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.add_field(name, space, base, nb_of_coefficients, tag_type, bh,
                           verb);
  skeletonFields.push_back(name);
  if (space == NOFIELD)
    SETERRQ(
        PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
        "NOFIELD space for boundary field not implemented in Simple interface");
 
  MoFEMFunctionReturn(0);
}

buildFields()

MoFEMErrorCode MoFEM::Simple::buildFields

(

)

Build fields.

Builds the field data structures and creates DOFs

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/tools/mesh_smoothing.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, and mofem/tutorials/scl-12/electrostatics.cpp.

Definition at line 584 of file Simple.cpp.

                                   {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  PetscLogEventBegin(MOFEM_EVENT_SimpleBuildFields, 0, 0, 0, 0);
 
  auto comm_interface_ptr = m_field.getInterface<CommInterface>();
  auto bit_ref_ptr = m_field.getInterface<BitRefManager>();
 
  // Add entities to the fields
  auto add_ents_to_field = [&](auto meshset, auto dim, auto &fields) {
    MoFEMFunctionBeginHot;
    for (auto &field : fields) {
      CHKERR m_field.add_ents_to_field_by_dim(meshset, dim, field);
      CHKERR comm_interface_ptr->synchroniseFieldEntities(field, 0);
    }
    MoFEMFunctionReturnHot(0);
  };
 
  CHKERR add_ents_to_field(meshSet, dIm, domainFields);
  CHKERR add_ents_to_field(meshSet, dIm, dataFields);
  CHKERR add_ents_to_field(boundaryMeshset, dIm - 1, boundaryFields);
  CHKERR add_ents_to_field(skeletonMeshset, dIm - 1, skeletonFields);
 
  std::set<std::string> nofield_fields;
  for (auto &f : noFieldFields)
    nofield_fields.insert(f);
  for (auto &f : noFieldDataFields)
    nofield_fields.insert(f);
 
  // Set order
 
  auto get_field_ptr = [&](auto &f) { return m_field.get_field_structure(f); };
  for (auto &t : fieldsOrder) {
    const auto f = std::get<0>(t);
    const auto order = std::get<1>(t);
 
    MOFEM_TAG_AND_LOG("WORLD", Sev::inform, "Simple")
        << "Set order to field " << f << " order " << order;
    MOFEM_LOG_CHANNEL("SYNC");
    if (std::get<3>(t)) {
      MOFEM_TAG_AND_LOG("SYNC", Sev::verbose, "Simple")
          << "To ents: " << std::endl
          << std::get<2>(t) << std::endl;
    }
    MOFEM_LOG_SEVERITY_SYNC(m_field.get_comm(), Sev::verbose);
 
    if (std::get<3>(t)) {
 
      CHKERR m_field.set_field_order(std::get<2>(t), f, order);
 
    } else {
      auto f_ptr = get_field_ptr(f);
 
      if (f_ptr->getSpace() == H1) {
        if (f_ptr->getApproxBase() == AINSWORTH_BERNSTEIN_BEZIER_BASE) {
          CHKERR m_field.set_field_order(meshSet, MBVERTEX, f, order);
        } else {
          CHKERR m_field.set_field_order(meshSet, MBVERTEX, f, 1);
        }
      }
 
      for (auto d = 1; d <= dIm; ++d) {
        for (EntityType t = CN::TypeDimensionMap[d].first;
             t <= CN::TypeDimensionMap[d].second; ++t) {
          CHKERR m_field.set_field_order(meshSet, t, f, order);
        }
      }
    }
  }
  MOFEM_LOG_CHANNEL("WORLD");
  // Build fields
  CHKERR m_field.build_fields();
  PetscLogEventEnd(MOFEM_EVENT_SimpleBuildFields, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

buildFiniteElements()

MoFEMErrorCode MoFEM::Simple::buildFiniteElements

(

)

Build finite elements.

Builds the finite element data structures and adjacencies

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/tools/mesh_smoothing.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, and mofem/tutorials/scl-12/electrostatics.cpp.

Definition at line 660 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 (addBoundaryFE || boundaryFields.size()) {
    CHKERR m_field.add_ents_to_finite_element_by_dim(boundaryMeshset, dIm - 1,
                                                     boundaryFE, true);
    CHKERR m_field.build_finite_elements(boundaryFE);
  }
  if (addSkeletonFE || skeletonFields.size()) {
    CHKERR m_field.add_ents_to_finite_element_by_dim(skeletonMeshset, dIm - 1,
                                                     skeletonFE, true);
    CHKERR m_field.build_finite_elements(skeletonFE);
  }
 
  if (noFieldFields.size() || meshsetFields.size()) {
 
    auto create_meshset = [&]() {
      EntityHandle meshset;
      auto create = [&]() {
        MoFEMFunctionBegin;
        CHKERR m_field.get_moab().create_meshset(MESHSET_SET, meshset);
        {
          ParallelComm *pcomm =
              ParallelComm::get_pcomm(&m_field.get_moab(), MYPCOMM_INDEX);
          Tag part_tag = pcomm->part_tag();
          int rank = pcomm->rank();
          CHKERR m_field.get_moab().tag_set_data(part_tag, &meshset, 1,
                                                 &rank);
          CHKERR m_field.getInterface<BitRefManager>()->setBitLevelToMeshset(
              meshset, BitRefLevel().set());
        }
        MoFEMFunctionReturn(0);
      };
      CHK_THROW_MESSAGE(create(), "Failed to create meshset");
      return meshset;
    };
 
    auto fe_meshset = create_meshset();
    CHKERR m_field.add_ents_to_finite_element_by_MESHSET(fe_meshset,
                                                         meshsetFE, false);
    // Add entities to the meshset, that will make DOFs on those entities
    // accessible on meshset finite element
    for (auto ents : meshsetFiniteElementEntities) {
      CHKERR m_field.get_moab().add_entities(fe_meshset, ents);
    }
    CHKERR m_field.build_finite_elements(meshsetFE);
 
  }
 
  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.

Builds the problem data structures and sets up DOF numbering

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/tools/mesh_smoothing.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, and mofem/tutorials/scl-12/electrostatics.cpp.

Definition at line 722 of file Simple.cpp.

                                    {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  PetscLogEventBegin(MOFEM_EVENT_SimpleBuildProblem, 0, 0, 0, 0);
  CHKERR m_field.build_adjacencies(bitLevel, bitLevelMask);
  // Set problem by the DOFs on the fields rather that by adding DOFs on the
  // elements
  m_field.getInterface<ProblemsManager>()->buildProblemFromFields = PETSC_TRUE;
  CHKERR DMSetUp_MoFEM(dM);
  m_field.getInterface<ProblemsManager>()->buildProblemFromFields = PETSC_FALSE;
  PetscLogEventEnd(MOFEM_EVENT_SimpleBuildProblem, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

createBoundaryMeshset()

MoFEMErrorCode MoFEM::Simple::createBoundaryMeshset ( )

private

Definition at line 844 of file Simple.cpp.

                                             {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  ParallelComm *pcomm =
      ParallelComm::get_pcomm(&m_field.get_moab(), MYPCOMM_INDEX);
 
  auto get_skin = [&](auto meshset) {
    // filter not owned entities, those are not on boundary
 
    Range domain_ents;
    CHKERR m_field.get_moab().get_entities_by_dimension(meshset, dIm,
                                                        domain_ents, true);
    CHKERR pcomm->filter_pstatus(domain_ents,
                                 PSTATUS_SHARED | PSTATUS_MULTISHARED,
                                 PSTATUS_NOT, -1, nullptr);
 
    Skinner skin(&m_field.get_moab());
    Range domain_skin;
    CHKERR skin.find_skin(0, domain_ents, false, domain_skin);
    CHKERR pcomm->filter_pstatus(domain_skin,
                                 PSTATUS_SHARED | PSTATUS_MULTISHARED,
                                 PSTATUS_NOT, -1, nullptr);
    return domain_skin;
  };
 
  auto create_boundary_meshset = [&](auto &&domain_skin) {
    MoFEMFunctionBeginHot;
    // create boundary meshset
    if (boundaryMeshset != 0) {
      MoFEMFunctionReturnHot(0);
    }
    CHKERR m_field.get_moab().create_meshset(MESHSET_SET, boundaryMeshset);
    CHKERR m_field.get_moab().add_entities(boundaryMeshset, domain_skin);
    for (int dd = 0; dd != dIm - 1; dd++) {
      Range adj;
      CHKERR m_field.get_moab().get_adjacencies(domain_skin, dd, false, adj,
                                                moab::Interface::UNION);
      CHKERR m_field.get_moab().add_entities(boundaryMeshset, adj);
    }
    MoFEMFunctionReturnHot(0);
  };
 
  CHKERR create_boundary_meshset(get_skin(meshSet));
 
  MoFEMFunctionReturn(0);
}

createSkeletonMeshset()

MoFEMErrorCode MoFEM::Simple::createSkeletonMeshset ( )

private

Definition at line 891 of file Simple.cpp.

                                             {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  ParallelComm *pcomm =
      ParallelComm::get_pcomm(&m_field.get_moab(), MYPCOMM_INDEX);
 
  auto create_skeleton_meshset = [&](auto meshset) {
    MoFEMFunctionBeginHot;
    // create boundary meshset
    if (skeletonMeshset != 0) {
      MoFEMFunctionReturnHot(0);
    }
    Range boundary_ents, skeleton_ents;
    CHKERR m_field.get_moab().get_entities_by_dimension(boundaryMeshset,
                                                        dIm - 1, boundary_ents);
    Range domain_ents;
    CHKERR m_field.get_moab().get_entities_by_dimension(meshset, dIm,
                                                        domain_ents, true);
    CHKERR m_field.get_moab().get_adjacencies(
        domain_ents, dIm - 1, true, skeleton_ents, moab::Interface::UNION);
    skeleton_ents = subtract(skeleton_ents, boundary_ents);
    CHKERR pcomm->filter_pstatus(skeleton_ents, PSTATUS_NOT_OWNED, PSTATUS_NOT,
                                 -1, nullptr);
    CHKERR m_field.get_moab().create_meshset(MESHSET_SET, skeletonMeshset);
    CHKERR m_field.get_moab().add_entities(skeletonMeshset, skeleton_ents);
    MoFEMFunctionReturnHot(0);
  };
 
  CHKERR create_skeleton_meshset(meshSet);
 
  MoFEMFunctionReturn(0);
}

defaultLoadFileFunc()

static MoFEMErrorCode MoFEM::Simple::defaultLoadFileFunc ( Interface &  m_field, const char *  file_name, const char *  options  )

inlinestatic

Set function for loading mesh file.

Default function for loading mesh files using MOAB interface

Parameters

m_field	MoFEM interface object
file_name	Name of the mesh file to load
options	Loader options string for MOAB

Returns

MoFEMErrorCode Error code indicating success or failure

Definition at line 57 of file Simple.hpp.

                                                                 {
    // Default behavior using the member moab.load_file
    return m_field.get_moab().load_file(file_name, 0, options);
  }

defineFiniteElements()

MoFEMErrorCode MoFEM::Simple::defineFiniteElements

(

)

Define finite elements.

Creates and defines domain, boundary, and skeleton finite elements based on added fields

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/tools/mesh_smoothing.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, and mofem/tutorials/scl-12/electrostatics.cpp.

Definition at line 471 of file Simple.cpp.

                                            {
  Interface &m_field = cOre;
  MoFEMFunctionBeginHot;
 
  auto clear_rows_and_cols = [&](auto &fe_name) {
    MoFEMFunctionBeginHot;
    auto fe_ptr = m_field.get_finite_elements();
    auto &fe_by_name = const_cast<FiniteElement_multiIndex *>(fe_ptr)
                           ->get<FiniteElement_name_mi_tag>();
    auto it_fe = fe_by_name.find(fe_name);
    if (it_fe != fe_by_name.end()) {
 
      if (!fe_by_name.modify(it_fe, FiniteElement_row_change_bit_reset()))
        SETERRQ(m_field.get_comm(), MOFEM_OPERATION_UNSUCCESSFUL,
                "modification unsuccessful");
 
      if (!fe_by_name.modify(it_fe, FiniteElement_col_change_bit_reset()))
        SETERRQ(m_field.get_comm(), MOFEM_OPERATION_UNSUCCESSFUL,
                "modification unsuccessful");
    }
    MoFEMFunctionReturnHot(0);
  };
  CHKERR clear_rows_and_cols(domainFE);
  CHKERR clear_rows_and_cols(boundaryFE);
  CHKERR clear_rows_and_cols(skeletonFE);
 
  // Define finite elements
  CHKERR m_field.add_finite_element(domainFE, MF_ZERO);
 
  auto add_fields = [&](auto &fe_name, auto &fields) {
    MoFEMFunctionBeginHot;
    for (auto &field : fields) {
      CHKERR m_field.modify_finite_element_add_field_row(fe_name, field);
      CHKERR m_field.modify_finite_element_add_field_col(fe_name, field);
      CHKERR m_field.modify_finite_element_add_field_data(fe_name, field);
    }
    MoFEMFunctionReturnHot(0);
  };
 
  auto add_data_fields = [&](auto &fe_name, auto &fields) {
    MoFEMFunctionBeginHot;
    for (auto &field : fields)
      CHKERR m_field.modify_finite_element_add_field_data(fe_name, field);
    MoFEMFunctionReturnHot(0);
  };
 
  CHKERR add_fields(domainFE, domainFields);
  CHKERR add_data_fields(domainFE, dataFields);
  CHKERR add_fields(domainFE, noFieldFields);
  CHKERR add_data_fields(domainFE, noFieldDataFields);
 
  if (addBoundaryFE || !boundaryFields.empty()) {
    CHKERR m_field.add_finite_element(boundaryFE, MF_ZERO);
    CHKERR add_fields(boundaryFE, domainFields);
    if (!boundaryFields.empty())
      CHKERR add_fields(boundaryFE, boundaryFields);
    CHKERR add_data_fields(boundaryFE, dataFields);
    CHKERR add_data_fields(boundaryFE, noFieldDataFields);
    CHKERR add_fields(boundaryFE, noFieldFields);
  }
  if (addSkeletonFE || !skeletonFields.empty()) {
    CHKERR m_field.add_finite_element(skeletonFE, MF_ZERO);
    CHKERR add_fields(skeletonFE, domainFields);
    if (!skeletonFields.empty())
      CHKERR add_fields(skeletonFE, skeletonFields);
    CHKERR add_data_fields(skeletonFE, dataFields);
    CHKERR add_data_fields(skeletonFE, noFieldDataFields);
    CHKERR add_fields(skeletonFE, noFieldFields);
  }
 
  if (noFieldFields.size() || meshsetFields.size()) {
    CHKERR m_field.add_finite_element(meshsetFE, MF_ZERO);
    CHKERR add_fields(meshsetFE, meshsetFields);
    CHKERR add_fields(meshsetFE, noFieldFields);
    CHKERR add_data_fields(meshsetFE, noFieldDataFields);
  }
 
  MoFEMFunctionReturnHot(0);
}

defineProblem()

MoFEMErrorCode MoFEM::Simple::defineProblem

(

const PetscBool

is_partitioned = PETSC_TRUE

)

define problem

Creates the problem and associates finite elements with it

Parameters

is_partitioned

Whether the mesh should be partitioned for parallel execution

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/tools/mesh_smoothing.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, and mofem/tutorials/scl-12/electrostatics.cpp.

Definition at line 551 of file Simple.cpp.

                                                                   {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  // Create dm instance
  dM = createDM(m_field.get_comm(), "DMMOFEM");
  // set dm data structure which created mofem data structures
  CHKERR DMMoFEMCreateMoFEM(dM, &m_field, nameOfProblem.c_str(), bitLevel,
                            bitLevelMask);
  CHKERR DMSetFromOptions(dM);
  CHKERR DMMoFEMAddElement(dM, domainFE);
  if (addBoundaryFE || !boundaryFields.empty()) {
    CHKERR DMMoFEMAddElement(dM, boundaryFE);
  }
  if (addSkeletonFE || !skeletonFields.empty()) {
    CHKERR DMMoFEMAddElement(dM, skeletonFE);
  }
  if (noFieldFields.size() || meshsetFields.size()) {
    CHKERR DMMoFEMAddElement(dM, meshsetFE);
  }
  CHKERR DMMoFEMAddElement(dM, otherFEs);
  CHKERR DMMoFEMSetIsPartitioned(dM, is_partitioned);
  MoFEMFunctionReturn(0);
}

deleteDM()

MoFEMErrorCode MoFEM::Simple::deleteDM

(

)

Delete dm.

Destroys the PETSc DM object and cleans up resources

Returns

MoFEMErrorCode Error code indicating success or failure

MoFEMErrorCode

Definition at line 812 of file Simple.cpp.

                                {
  MoFEMFunctionBegin;
  dM.reset();
  MoFEMFunctionReturn(0);
}

deleteFiniteElements()

MoFEMErrorCode MoFEM::Simple::deleteFiniteElements

(

)

Delete finite elements.

Removes all finite elements created by the Simple interface

Returns

MoFEMErrorCode Error code indicating success or failure

MoFEMErrorCode

Definition at line 823 of file Simple.cpp.

                                            {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  for (auto fe : {domainFE, boundaryFE, skeletonFE}) {
    if (m_field.check_finite_element(fe)) {
      CHKERR m_field.delete_finite_element(fe);
    }
  }
  MoFEMFunctionReturn(0);
}

exchangeGhostCells()

MoFEMErrorCode MoFEM::Simple::exchangeGhostCells ( )

private

get all adjacent ghosted entities

Definition at line 925 of file Simple.cpp.

                                          {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  MOFEM_LOG("WORLD", Sev::verbose) << "Exchange ghost cells";
 
  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());
 
  CHKERR pcomm->exchange_ghost_cells(getDim(), getDim() - 1, 1,
                                     3 /**get all adjacent ghosted entities */,
                                     true, true, meshSet ? &meshSet : nullptr);
 
  Range shared;
  CHKERR m_field.get_moab().get_entities_by_dimension(0, dIm, shared);
  for (auto d = dIm - 1; d >= 0; --d) {
    CHKERR m_field.get_moab().get_adjacencies(shared, d, false, shared,
                                              moab::Interface::UNION);
  }
  CHKERR pcomm->filter_pstatus(shared, PSTATUS_SHARED | PSTATUS_MULTISHARED,
                               PSTATUS_OR, -1, &shared);
  Tag part_tag = pcomm->part_tag();
  CHKERR pcomm->exchange_tags(part_tag, shared);
  CHKERR m_field.getInterface<CommInterface>()->resolveParentEntities(shared,
                                                                      VERBOSE);
 
  MoFEMFunctionReturn(0);
}

getAddBoundaryFE()

bool & MoFEM::Simple::getAddBoundaryFE ( )

inline

Get the addBoundaryFE flag.

If variable set to true, boundary element is created regardless of whether a field on boundary is added or not.

Returns

bool& Reference to flag indicating if boundary FE should be created

Examples

mofem/atom_tests/simple_l2_only.cpp, and mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp.

Definition at line 546 of file Simple.hpp.

{ return addBoundaryFE; }

getAddSkeletonFE()

bool & MoFEM::Simple::getAddSkeletonFE ( )

inline

Get the addSkeletonFE flag.

If variable set to true, skeleton element is created regardless of whether a field on skeleton is added or not.

Returns

bool& Reference to flag indicating if skeleton FE should be created

Examples

mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, and mofem/tutorials/scl-12/electrostatics.cpp.

Definition at line 536 of file Simple.hpp.

{ return addSkeletonFE; }

getBitAdjEnt()

auto & MoFEM::Simple::getBitAdjEnt ( )

inline

Get bit ref level for entities.

Returns

auto& Reference to bit refinement level for entity adjacencies

Definition at line 576 of file Simple.hpp.

{ return bitAdjEnt; }

getBitAdjEntMask()

auto & MoFEM::Simple::getBitAdjEntMask ( )

inline

Get bit ref level mask for entities.

Returns

auto& Reference to bit refinement level mask for entity adjacencies

Definition at line 583 of file Simple.hpp.

{ return bitAdjEntMask; }

getBitAdjParent()

auto & MoFEM::Simple::getBitAdjParent ( )

inline

Get bit ref level for parent.

Returns

auto& Reference to bit refinement level for parent adjacencies

Definition at line 562 of file Simple.hpp.

{ return bitAdjParent; }

getBitAdjParentMask()

auto & MoFEM::Simple::getBitAdjParentMask ( )

inline

Get bit ref level mask for parent.

Returns

auto& Reference to bit refinement level mask for parent adjacencies

Definition at line 569 of file Simple.hpp.

{ return bitAdjParentMask; }

getBitRefLevel()

BitRefLevel & MoFEM::Simple::getBitRefLevel ( )

inline

Get the BitRefLevel.

Returns

BitRefLevel& Reference to bit refinement level

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/tutorials/adv-3/level_set.cpp, mofem/tutorials/fun-2/plot_base.cpp, and mofem/tutorials/mix-0/mixed_poisson.cpp.

Definition at line 415 of file Simple.hpp.

{ return bitLevel; }

getBitRefLevelMask()

BitRefLevel & MoFEM::Simple::getBitRefLevelMask ( )

inline

Get the BitRefLevelMask.

Returns

BitRefLevel& Reference to bit refinement level mask

Examples

mofem/atom_tests/child_and_parent.cpp, and mofem/atom_tests/hanging_node_approx.cpp.

Definition at line 422 of file Simple.hpp.

{ return bitLevelMask; }

getBoundaryFEName() [1/2]

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

inline

Get the Boundary FE Name.

Returns

std::string&

Definition at line 464 of file Simple.hpp.

{ return boundaryFE; }

getBoundaryFEName() [2/2]

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

inline

Get the Boundary FE Name.

Returns

const std::string

Examples

mofem/atom_tests/simple_interface.cpp, and mofem/tools/mesh_smoothing.cpp.

Definition at line 436 of file Simple.hpp.

{ return boundaryFE; }

getBoundaryMeshSet()

EntityHandle & MoFEM::Simple::getBoundaryMeshSet ( )

inline

Get the BoundaryMeshSet object.

Returns

EntityHandle& Reference to boundary meshset handle

Examples

mofem/atom_tests/hanging_node_approx.cpp.

Definition at line 401 of file Simple.hpp.

{ return boundaryMeshset; }

getDim()

int MoFEM::Simple::getDim ( ) const

inline

Get the problem dimension.

Problem dimension is determined by highest dimension of entities on the mesh.

Returns

int Problem dimension (1D, 2D, or 3D)

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/field_evaluator.cpp, and mofem/atom_tests/hanging_node_approx.cpp.

Definition at line 373 of file Simple.hpp.

{ return dIm; }

getDM() [1/2]

SmartPetscObj< DM > MoFEM::Simple::getDM ( )

inline

Return smart DM object.

Returns a smart pointer to the DM object that handles automatic cleanup

{
   auto dm  = simple_interface->getDM();
 
   // ...
 
   // dm is automatically destroyed when out of the scope
}

Returns

SmartPetscObj<DM> Smart pointer to DM object

Definition at line 363 of file Simple.hpp.

{ return dM; }

getDM() [2/2]

MoFEMErrorCode MoFEM::Simple::getDM

(

DM *

dm

)

Get DM.

Retrieves the PETSc Discrete Manager (DM) object

Parameters

dm	Pointer to store the DM object

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/field_evaluator.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_check_approx_in_2d.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/atom_tests/loop_entities.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tools/mesh_smoothing.cpp, mofem/tutorials/adv-3/level_set.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/mix-1/phase.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-12/electrostatics.cpp, mofem/tutorials/scl-9/heat_method.cpp, mofem/tutorials/vec-4/shallow_wave.cpp, mofem/tutorials/vec-5/free_surface.cpp, mofem/tutorials/vec-7/adjoint.cpp, mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_element.cpp, and mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 800 of file Simple.cpp.

                                   {
  MoFEMFunctionBegin;
  CHKERR PetscObjectReference(getPetscObject(dM.get()));
  *dm = dM.get();
  MoFEMFunctionReturn(0);
}

getDomainFEName() [1/2]

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

inline

Get the Domain FE Name.

Returns

std::string&

Definition at line 457 of file Simple.hpp.

{ return domainFE; }

getDomainFEName() [2/2]

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

inline

Get the Domain FE Name.

Returns

const std::string

Examples

mofem/atom_tests/field_evaluator.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tools/mesh_smoothing.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, mofem/tutorials/scl-12/electrostatics.cpp, mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_element.cpp, and mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 429 of file Simple.hpp.

{ return domainFE; }

getMeshSet()

DEPRECATED EntityHandle & MoFEM::Simple::getMeshSet ( )

inline

Deprecated:

Use getMeshset

Returns

EntityHandle& Reference to domain meshset handle

Definition at line 387 of file Simple.hpp.

{ return meshSet; }

getMeshset()

EntityHandle & MoFEM::Simple::getMeshset ( )

inline

Get the MeshSet object.

Returns

EntityHandle& Reference to domain meshset handle

Definition at line 394 of file Simple.hpp.

{ return meshSet; }

getMeshsetFEName()

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

inline

Get the Skeleton FE Name.

Returns

std::string&

Examples

mofem/atom_tests/simple_interface.cpp.

Definition at line 478 of file Simple.hpp.

{ return meshsetFE; }

getMeshsetFiniteElementEntities()

auto & MoFEM::Simple::getMeshsetFiniteElementEntities ( )

inline

Range of entities added to meshset finite element.

Meshset element have not defined integration rule, if you like to integrate you have to add mesh entity to meshset finite element. In principle meshset element become standard FE.

Returns

std::vector<std::string>&

Examples

mofem/atom_tests/simple_interface.cpp, and mofem/atom_tests/simple_l2_only.cpp.

Definition at line 489 of file Simple.hpp.

                                                 {
    return meshsetFiniteElementEntities;
  }

getOptions()

MoFEMErrorCode MoFEM::Simple::getOptions

(

)

get options

Retrieves command-line options and configuration settings for the Simple interface

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/boundary_marker.cpp, mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/field_evaluator.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_check_approx_in_2d.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/atom_tests/loop_entities.cpp, mofem/atom_tests/operators_tests.cpp, mofem/atom_tests/schur_test_diag_mat.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/atom_tests/tensor_divergence_operator.cpp, mofem/atom_tests/test_broken_space.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tools/mesh_smoothing.cpp, mofem/tutorials/adv-0/plastic.cpp, mofem/tutorials/adv-2/thermo_elastic.cpp, mofem/tutorials/adv-3/level_set.cpp, mofem/tutorials/adv-4/dynamic_first_order_con_law.cpp, mofem/tutorials/adv-5/seepage.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/fun-1/integration.cpp, mofem/tutorials/fun-2/plot_base.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/mix-1/phase.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-12/electrostatics.cpp, mofem/tutorials/scl-8/radiation.cpp, mofem/tutorials/scl-9/heat_method.cpp, mofem/tutorials/vec-1/eigen_elastic.cpp, mofem/tutorials/vec-3/nonlinear_dynamic_elastic.cpp, mofem/tutorials/vec-4/approx_sphere.cpp, mofem/tutorials/vec-4/shallow_wave.cpp, mofem/tutorials/vec-6/plate.cpp, mofem/tutorials/vec-7/adjoint.cpp, mofem/users_modules/adolc-plasticity/adolc_plasticity.cpp, mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_element.cpp, and mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 180 of file Simple.cpp.

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

getOtherFiniteElements()

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

inline

Get the Other Finite Elements.

User can create finite elements using directly core interface and and them to simple problem by this function

Returns

std::vector<std::string>&

Examples

mofem/tools/mesh_smoothing.cpp.

Definition at line 508 of file Simple.hpp.

{ return otherFEs; }

getParentAdjacencies()

bool & MoFEM::Simple::getParentAdjacencies ( )

inline

Get the addParentAdjacencies flag.

If set true, parent adjacencies are added to finite elements

Returns

bool& Reference to flag indicating if parent adjacencies should be added

Examples

mofem/atom_tests/hanging_node_approx.cpp, and mofem/tutorials/vec-5/free_surface.cpp.

Definition at line 555 of file Simple.hpp.

{ return addParentAdjacencies; }

getProblemName() [1/2]

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

inline

Get the Problem Name.

Returns

std::string&

Definition at line 498 of file Simple.hpp.

{ return nameOfProblem; }

getProblemName() [2/2]

const std::string MoFEM::Simple::getProblemName ( ) const

inline

Get the Problem Name.

Returns

const std::string

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/fun-2/plot_base.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, and mofem/tutorials/scl-12/electrostatics.cpp.

Definition at line 450 of file Simple.hpp.

{ return nameOfProblem; }

getSkeletonFEName() [1/2]

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

inline

Get the Skeleton FE Name.

Returns

std::string&

Definition at line 471 of file Simple.hpp.

{ return skeletonFE; }

getSkeletonFEName() [2/2]

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

inline

Get the Skeleton FE Name.

Returns

const std::string

Examples

mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, and mofem/atom_tests/simple_interface.cpp.

Definition at line 443 of file Simple.hpp.

{ return skeletonFE; }

getSkeletonMeshSet()

EntityHandle & MoFEM::Simple::getSkeletonMeshSet ( )

inline

Get the SkeletonMeshSet object.

Returns

EntityHandle& Reference to skeleton meshset handle

Definition at line 408 of file Simple.hpp.

{ return skeletonMeshset; }

loadFile() [1/2]

MoFEMErrorCode MoFEM::Simple::loadFile

(

const std::string

mesh_file_name = ""

)

Load mesh file with parallel options if number of cores > 1.

Automatically selects appropriate loading options for parallel execution

Parameters

mesh_file_name

Name of mesh file to load (if empty, uses default or command line value)

Returns

MoFEMErrorCode Error code indicating success or failure

Definition at line 247 of file Simple.cpp.

                                                              {
  MoFEMFunctionBegin;
  Interface &m_field = cOre;
  if (m_field.get_comm_size() == 1)
    CHKERR loadFile("", meshFileName);
  else
    CHKERR loadFile("PARALLEL=READ_PART;"
                    "PARALLEL_RESOLVE_SHARED_ENTS;"
                    "PARTITION=PARALLEL_PARTITION;",
                    meshFileName);
  MoFEMFunctionReturn(0);
}

loadFile() [2/2]

MoFEMErrorCode MoFEM::Simple::loadFile	(	const std::string	options,
		const std::string	mesh_file_name,
		LoadFileFunc	loadFunc = defaultLoadFileFunc
	)

Load mesh file.

Loads mesh file with specified options and custom loading function

Parameters

options	File load options string passed to the loader
mesh_file_name	Name of mesh file to load (if empty, uses default or command line)
loadFunc	Function for loading mesh file (default: defaultLoadFileFunc)

Note

If bitRefLevel is set to any, bit ref level of loaded entities is not changed. After mesh is load, bit ref level should be set to create problem. Default setting of bit ref level is on first bit, and if is set all mesh entities on load are set to set level.

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/field_evaluator.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_check_approx_in_2d.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/atom_tests/loop_entities.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tools/mesh_smoothing.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/fun-2/plot_base.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/mix-1/phase.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-12/electrostatics.cpp, mofem/tutorials/scl-9/heat_method.cpp, mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_element.cpp, and mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 191 of file Simple.cpp.

                                                       {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  PetscLogEventBegin(MOFEM_EVENT_SimpleLoadMesh, 0, 0, 0, 0);
 
  if (!mesh_file_name.empty())
    strcpy(meshFileName, mesh_file_name.c_str());
 
  // Invoke the boost function, passing in required arguments
  CHKERR loadFunc(m_field, meshFileName, options.c_str());
  CHKERR m_field.rebuild_database();
 
  // determine problem dimension
  if (dIm == -1) {
    int nb_ents_3d;
    CHKERR m_field.get_moab().get_number_entities_by_dimension(
        meshSet, 3, nb_ents_3d, true);
    if (nb_ents_3d > 0) {
      dIm = 3;
    } else {
      int nb_ents_2d;
      CHKERR m_field.get_moab().get_number_entities_by_dimension(
          meshSet, 2, nb_ents_2d, true);
      if (nb_ents_2d > 0) {
        dIm = 2;
      } else {
        dIm = 1;
      }
    }
  }
 
  if (!boundaryMeshset)
    CHKERR createBoundaryMeshset();
  if (!skeletonMeshset)
    CHKERR createSkeletonMeshset();
  if (addSkeletonFE)
    CHKERR exchangeGhostCells();
 
  if (bitLevel.any()) {
    Range ents;
    CHKERR m_field.get_moab().get_entities_by_dimension(meshSet, dIm, ents,
                                                        true);
    CHKERR m_field.getInterface<BitRefManager>()->setBitRefLevel(ents, bitLevel,
                                                                 false);
  } else {
    MOFEM_LOG("WORLD", Sev::warning) << "BitRefLevel is none and not set";
    CHKERR m_field.getInterface<BitRefManager>()->addToDatabaseBitRefLevelByDim(
        dIm, BitRefLevel().set(), BitRefLevel().set());
  }
 
  PetscLogEventEnd(MOFEM_EVENT_SimpleLoadMesh, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

query_interface()

MoFEMErrorCode MoFEM::Simple::query_interface ( boost::typeindex::type_index  type_index, UnknownInterface **  iface  ) const

virtual

Implements MoFEM::UnknownInterface.

Definition at line 8 of file Simple.cpp.

                                                                       {
  *iface = const_cast<Simple *>(this);
  return 0;
}

removeBoundaryField()

MoFEMErrorCode MoFEM::Simple::removeBoundaryField

(

const std::string

name

)

Remove field from boundary.

Parameters

name	Name of the field to remove from boundary

Returns

MoFEMErrorCode Error code indicating success or failure

Definition at line 443 of file Simple.cpp.

                                                               {
  MoFEMFunctionBegin;
 
  auto remove_field_from_list = [&](auto &vec) {
    auto it = std::find(vec.begin(), vec.end(), name);
    if (it != vec.end())
      vec.erase(it);
  };
 
  remove_field_from_list(boundaryFields);
 
  MoFEMFunctionReturn(0);
}

removeDomainField()

MoFEMErrorCode MoFEM::Simple::removeDomainField

(

const std::string

name

)

Remove field from domain.

Parameters

name	Name of the field to remove from domain

Returns

MoFEMErrorCode Error code indicating success or failure

Definition at line 429 of file Simple.cpp.

                                                             {
  MoFEMFunctionBegin;
 
  auto remove_field_from_list = [&](auto &vec) {
    auto it = std::find(vec.begin(), vec.end(), name);
    if (it != vec.end())
      vec.erase(it);
  };
 
  remove_field_from_list(domainFields);
 
  MoFEMFunctionReturn(0);
}

removeSkeletonField()

MoFEMErrorCode MoFEM::Simple::removeSkeletonField

(

const std::string

name

)

Remove field from skeleton.

Parameters

name	Name of the field to remove from skeleton

Returns

MoFEMErrorCode Error code indicating success or failure

Definition at line 457 of file Simple.cpp.

                                                               {
  MoFEMFunctionBegin;
 
  auto remove_field_from_list = [&](auto &vec) {
    auto it = std::find(vec.begin(), vec.end(), name);
    if (it != vec.end())
      vec.erase(it);
  };
 
  remove_field_from_list(skeletonFields);
 
  MoFEMFunctionReturn(0);
}

reSetUp()

MoFEMErrorCode MoFEM::Simple::reSetUp

(

bool

only_dm = false

)

Rebuild internal MoFEM data structures.

Call this function after you add field or remove it.

Note

If you add field, or remove it, finite element and problem needs to be rebuild. However DM can remain the same.

Parameters

only_dm

If true, only rebuild DM without rebuilding fields and finite elements

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/child_and_parent.cpp.

Definition at line 762 of file Simple.cpp.

                                           {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  PetscLogEventBegin(MOFEM_EVENT_SimpleBuildProblem, 0, 0, 0, 0);
 
  if (!only_dm) {
    CHKERR defineFiniteElements();
    CHKERR buildFields();
    if (addSkeletonFE || !skeletonFields.empty()) {
      CHKERR setSkeletonAdjacency(-1, skeletonFE);
      if (addBoundaryFE || !boundaryFields.empty()) {
        CHKERR setSkeletonAdjacency(-1, boundaryFE);
      }
    }
    if (addParentAdjacencies)
      CHKERR setParentAdjacency();
    CHKERR buildFiniteElements();
  }
 
  CHKERR m_field.build_adjacencies(bitLevel, bitLevelMask);
 
  const Problem *problem_ptr;
  CHKERR DMMoFEMGetProblemPtr(dM, &problem_ptr);
  const auto problem_name = problem_ptr->getName();
  CHKERR m_field.modify_problem_ref_level_set_bit(problem_name, bitLevel);
  CHKERR m_field.modify_problem_mask_ref_level_set_bit(problem_name,
                                                       bitLevelMask);
 
  // Set problem by the DOFs on the fields rather that by adding DOFs on the
  // elements
  m_field.getInterface<ProblemsManager>()->buildProblemFromFields = PETSC_TRUE;
  CHKERR DMSetUp_MoFEM(dM);
  m_field.getInterface<ProblemsManager>()->buildProblemFromFields = PETSC_FALSE;
 
  PetscLogEventEnd(MOFEM_EVENT_SimpleBuildProblem, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

setDim()

void MoFEM::Simple::setDim ( int  dim )

inline

Set the problem dimension.

Parameters

dim	Problem dimension to set (1, 2, or 3)

Examples

mofem/tutorials/fun-2/plot_base.cpp.

Definition at line 380 of file Simple.hpp.

{ dIm = dim; };

setFieldOrder()

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

Set field order.

Sets the approximation order for a field on specified entities

Parameters

field_name	Name of the field
order	Approximation order to set
ents	Range of entities to apply order to (NULL for all entities)

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/field_evaluator.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_check_approx_in_2d.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/atom_tests/loop_entities.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tools/mesh_smoothing.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/cor-6/simple_elasticity.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/fun-2/plot_base.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/mix-1/phase.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-12/electrostatics.cpp, mofem/tutorials/scl-9/heat_method.cpp, mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_element.cpp, and mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 575 of file Simple.cpp.

                                                                         {
  MoFEMFunctionBeginHot;
  fieldsOrder.emplace_back(field_name, order,
                           ents == NULL ? Range() : Range(*ents),
                           ents == NULL ? false : true);
  MoFEMFunctionReturnHot(0);
}

setParentAdjacency() [1/4]

template<>

MoFEMErrorCode MoFEM::Simple::setParentAdjacency ( )

private

Definition at line 72 of file Simple.cpp.

                                                         {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  parentAdjFunctionDim3 =
      boost::make_shared<ParentFiniteElementAdjacencyFunction<3>>(
          bitAdjParent, bitAdjParentMask, bitAdjEnt, bitAdjEntMask);
  parentAdjFunctionDim2 =
      boost::make_shared<ParentFiniteElementAdjacencyFunction<2>>(
          bitAdjParent, bitAdjParentMask, bitAdjEnt, bitAdjEntMask);
 
  CHKERR m_field.modify_finite_element_adjacency_table(domainFE, MBTET,
                                                       *parentAdjFunctionDim3);
  CHKERR m_field.modify_finite_element_adjacency_table(domainFE, MBHEX,
                                                       *parentAdjFunctionDim3);
  if (addBoundaryFE || !boundaryFields.empty()) {
    CHKERR m_field.modify_finite_element_adjacency_table(
        boundaryFE, MBTRI, *parentAdjFunctionDim2);
    CHKERR m_field.modify_finite_element_adjacency_table(
        boundaryFE, MBQUAD, *parentAdjFunctionDim2);
  }
  if (addSkeletonFE || !skeletonFields.empty()) {
    CHKERR m_field.modify_finite_element_adjacency_table(
        skeletonFE, MBTRI, *parentAdjFunctionDim2);
    CHKERR m_field.modify_finite_element_adjacency_table(
        skeletonFE, MBQUAD, *parentAdjFunctionDim2);
  }
 
  MoFEMFunctionReturn(0);
}

setParentAdjacency() [2/4]

template<>

MoFEMErrorCode MoFEM::Simple::setParentAdjacency ( )

private

Definition at line 103 of file Simple.cpp.

                                                         {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  parentAdjFunctionDim2 =
      boost::make_shared<ParentFiniteElementAdjacencyFunction<2>>(
          bitAdjParent, bitAdjParentMask, bitAdjEnt, bitAdjEntMask);
  parentAdjFunctionDim1 =
      boost::make_shared<ParentFiniteElementAdjacencyFunction<1>>(
          bitAdjParent, bitAdjParentMask, bitAdjEnt, bitAdjEntMask);
 
  CHKERR m_field.modify_finite_element_adjacency_table(domainFE, MBTRI,
                                                       *parentAdjFunctionDim2);
  CHKERR m_field.modify_finite_element_adjacency_table(domainFE, MBQUAD,
                                                       *parentAdjFunctionDim2);
  if (addBoundaryFE || !boundaryFields.empty())
    CHKERR m_field.modify_finite_element_adjacency_table(
        boundaryFE, MBEDGE, *parentAdjFunctionDim1);
  if (addSkeletonFE || !skeletonFields.empty())
    CHKERR m_field.modify_finite_element_adjacency_table(
        skeletonFE, MBEDGE, *parentAdjFunctionDim1);
 
  MoFEMFunctionReturn(0);
}

setParentAdjacency() [3/4]

template<>

MoFEMErrorCode MoFEM::Simple::setParentAdjacency ( )

private

Definition at line 128 of file Simple.cpp.

                                                          {
  MoFEMFunctionBegin;
  switch (getDim()) {
  case 1:
    THROW_MESSAGE("Not implemented");
  case 2:
    return setParentAdjacency<2>();
  case 3:
    return setParentAdjacency<3>();
  default:
    THROW_MESSAGE("Not implemented");
  }
  MoFEMFunctionReturn(0);
}

setParentAdjacency() [4/4]

template<int DIM>

MoFEMErrorCode MoFEM::Simple::setParentAdjacency ( )

private

Definition at line 67 of file Simple.cpp.

                                                             {
  static_assert(DIM == 2 || DIM == 3, "not implemented");
  return MOFEM_NOT_IMPLEMENTED;
}

setSkeletonAdjacency() [1/5]

MoFEMErrorCode MoFEM::Simple::setSkeletonAdjacency	(	int	dim = -1,
		std::string	fe_name = ""
	)

Set the skeleton adjacency object.

Sets up adjacency relationships for skeleton finite elements

Parameters

dim	Dimension for skeleton adjacency (-1 for automatic detection)
fe_name	Name of finite element to set adjacency for (empty for default)

Returns

MoFEMErrorCode Error code indicating success or failure

Definition at line 143 of file Simple.cpp.

                                                                      {
  MoFEMFunctionBegin;
  if (dim == -1)
    dim = getDim();
 
  if (fe_name.empty())
    fe_name = skeletonFE;
 
  switch (dim) {
  case 2:
    return setSkeletonAdjacency<2>(fe_name);
  case 3:
    return setSkeletonAdjacency<3>(fe_name);
  default:
    SETERRQ(PETSC_COMM_WORLD, MOFEM_NOT_IMPLEMENTED, "Not implemented");
  }
  MoFEMFunctionReturn(0);
}

setSkeletonAdjacency() [2/5]

template<>

MoFEMErrorCode MoFEM::Simple::setSkeletonAdjacency ( std::string  fe_name )

private

Definition at line 21 of file Simple.cpp.

                                                                {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  parentAdjSkeletonFunctionDim1 =
      boost::make_shared<ParentFiniteElementAdjacencyFunctionSkeleton<1>>(
          bitAdjParent, bitAdjParentMask, bitAdjEnt, bitAdjEntMask);
  CHKERR m_field.modify_finite_element_adjacency_table(
      fe_name, MBEDGE, *parentAdjSkeletonFunctionDim1);
 
  MoFEMFunctionReturn(0);
}

setSkeletonAdjacency() [3/5]

template<>

MoFEMErrorCode MoFEM::Simple::setSkeletonAdjacency ( std::string  fe_name )

private

Definition at line 35 of file Simple.cpp.

                                                                {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  parentAdjSkeletonFunctionDim2 =
      boost::make_shared<ParentFiniteElementAdjacencyFunctionSkeleton<2>>(
          bitAdjParent, bitAdjParentMask, bitAdjEnt, bitAdjEntMask);
 
  CHKERR m_field.modify_finite_element_adjacency_table(
      fe_name, MBTRI, *parentAdjSkeletonFunctionDim2);
  CHKERR m_field.modify_finite_element_adjacency_table(
      fe_name, MBQUAD, *parentAdjSkeletonFunctionDim2);
 
  MoFEMFunctionReturn(0);
}

setSkeletonAdjacency() [4/5]

template<>

MoFEMErrorCode MoFEM::Simple::setSkeletonAdjacency ( std::string  fe_name )

private

Definition at line 52 of file Simple.cpp.

                                                                 {
  MoFEMFunctionBegin;
 
  switch (getDim()) {
  case 1:
    THROW_MESSAGE("Not implemented");
  case 2:
    return setSkeletonAdjacency<2>(fe_name);
  case 3:
    return setSkeletonAdjacency<3>(fe_name);
  default:
    THROW_MESSAGE("Not implemented");
  }
  MoFEMFunctionReturn(0);
}

setSkeletonAdjacency() [5/5]

template<int DIM>

MoFEMErrorCode MoFEM::Simple::setSkeletonAdjacency ( std::string  fe_name )

private

Definition at line 15 of file Simple.cpp.

                                                             {
  static_assert(DIM == 2 || DIM == 3, "not implemented");
  return MOFEM_NOT_IMPLEMENTED;
}

setUp()

MoFEMErrorCode MoFEM::Simple::setUp

(

const PetscBool

is_partitioned = PETSC_TRUE

)

Setup problem.

Complete setup of the problem including fields, finite elements, and problem definition

Parameters

is_partitioned

Whether the mesh should be partitioned for parallel execution

Returns

MoFEMErrorCode Error code indicating success or failure

Examples

mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/field_evaluator.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_check_approx_in_2d.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/atom_tests/loop_entities.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tutorials/clx-0/helmholtz.cpp, mofem/tutorials/cor-9/reaction_diffusion.cpp, mofem/tutorials/fun-2/plot_base.cpp, mofem/tutorials/mix-0/mixed_poisson.cpp, mofem/tutorials/mix-1/phase.cpp, mofem/tutorials/scl-0/approximaton.cpp, mofem/tutorials/scl-1/poisson_2d_homogeneous.cpp, mofem/tutorials/scl-11/poisson_2d_dis_galerkin.cpp, mofem/tutorials/scl-12/electrostatics.cpp, mofem/tutorials/scl-9/heat_method.cpp, mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_element.cpp, and mofem/users_modules/basic_finite_elements/atom_tests/testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 736 of file Simple.cpp.

                                                           {
  MoFEMFunctionBegin;
 
  PetscLogEventBegin(MOFEM_EVENT_SimpleSetUP, 0, 0, 0, 0);
 
  CHKERR defineFiniteElements();
 
  if (addSkeletonFE || !skeletonFields.empty()) {
    CHKERR setSkeletonAdjacency(-1, skeletonFE);
    if (addBoundaryFE || !boundaryFields.empty()) {
      CHKERR setSkeletonAdjacency(-1, boundaryFE);
    }
  }
 
  if (addParentAdjacencies)
    CHKERR setParentAdjacency();
 
  CHKERR defineProblem(is_partitioned);
  CHKERR buildFields();
  CHKERR buildFiniteElements();
  CHKERR buildProblem();
 
  PetscLogEventEnd(MOFEM_EVENT_SimpleSetUP, 0, 0, 0, 0);
  MoFEMFunctionReturn(0);
}

Member Data Documentation

addBoundaryFE

bool MoFEM::Simple::addBoundaryFE

private

Add boundary FE.

Definition at line 617 of file Simple.hpp.

addParentAdjacencies

bool MoFEM::Simple::addParentAdjacencies

private

If set to true parent adjacencies are build.

Definition at line 618 of file Simple.hpp.

addSkeletonFE

bool MoFEM::Simple::addSkeletonFE

private

Add skeleton FE.

Definition at line 616 of file Simple.hpp.

bitAdjEnt

BitRefLevel MoFEM::Simple::bitAdjEnt

private

bit ref level for parent

Definition at line 622 of file Simple.hpp.

bitAdjEntMask

BitRefLevel MoFEM::Simple::bitAdjEntMask

private

bit ref level for parent parent

Definition at line 623 of file Simple.hpp.

bitAdjParent

BitRefLevel MoFEM::Simple::bitAdjParent

private

bit ref level for parent

Definition at line 620 of file Simple.hpp.

bitAdjParentMask

BitRefLevel MoFEM::Simple::bitAdjParentMask

private

bit ref level for parent parent

Definition at line 621 of file Simple.hpp.

bitLevel

BitRefLevel MoFEM::Simple::bitLevel

private

BitRefLevel of the problem.

Definition at line 602 of file Simple.hpp.

bitLevelMask

BitRefLevel MoFEM::Simple::bitLevelMask

private

BitRefLevel of the problem.

Definition at line 603 of file Simple.hpp.

boundaryFE

std::string MoFEM::Simple::boundaryFE = "bFE"

private

boundary finite element

Definition at line 639 of file Simple.hpp.

boundaryFields

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

private

boundary fields

Definition at line 626 of file Simple.hpp.

boundaryMeshset

EntityHandle MoFEM::Simple::boundaryMeshset

private

meshset with boundary

Definition at line 613 of file Simple.hpp.

cOre

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

private

Definition at line 600 of file Simple.hpp.

dataFields

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

private

Data fields.

Definition at line 628 of file Simple.hpp.

dIm

int MoFEM::Simple::dIm

private

dimension of problem

Definition at line 648 of file Simple.hpp.

dM

SmartPetscObj<DM> MoFEM::Simple::dM

private

Discrete manager (interface to PETSc/MoFEM functions)

Definition at line 651 of file Simple.hpp.

domainFE

std::string MoFEM::Simple::domainFE = "dFE"

private

domain finite element

Definition at line 638 of file Simple.hpp.

domainFields

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

private

domain fields

Definition at line 625 of file Simple.hpp.

fieldsOrder

std::list<std::tuple<std::string, int, Range, bool> > MoFEM::Simple::fieldsOrder

private

fields order. 1: field name, order, range, set by range if true

Definition at line 634 of file Simple.hpp.

meshFileName

char MoFEM::Simple::meshFileName[255]

private

mesh file name

Definition at line 647 of file Simple.hpp.

meshSet

EntityHandle MoFEM::Simple::meshSet

private

domain meshset

Definition at line 612 of file Simple.hpp.

meshsetFE

std::string MoFEM::Simple::meshsetFE = "mFE"

private

meshset finite element

Definition at line 641 of file Simple.hpp.

meshsetFields

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

private

meshset fields

Definition at line 629 of file Simple.hpp.

meshsetFiniteElementEntities

std::vector<Range> MoFEM::Simple::meshsetFiniteElementEntities

private

Meshset element entities.

Definition at line 645 of file Simple.hpp.

MOFEM_EVENT_SimpleBuildFields

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleBuildFields

private

Definition at line 607 of file Simple.hpp.

MOFEM_EVENT_SimpleBuildFiniteElements

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleBuildFiniteElements

private

Definition at line 608 of file Simple.hpp.

MOFEM_EVENT_SimpleBuildProblem

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleBuildProblem

private

Definition at line 609 of file Simple.hpp.

MOFEM_EVENT_SimpleKSPSolve

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleKSPSolve

private

Definition at line 610 of file Simple.hpp.

MOFEM_EVENT_SimpleLoadMesh

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleLoadMesh

private

Definition at line 606 of file Simple.hpp.

MOFEM_EVENT_SimpleSetUP

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleSetUP

private

Definition at line 605 of file Simple.hpp.

nameOfProblem

std::string MoFEM::Simple::nameOfProblem = "SimpleProblem"

private

problem name

Definition at line 637 of file Simple.hpp.

noFieldDataFields

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

private

NOFIELD field name.

Definition at line 631 of file Simple.hpp.

noFieldFields

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

private

NOFIELD field name.

Definition at line 630 of file Simple.hpp.

otherFEs

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

private

Other finite elements.

Definition at line 643 of file Simple.hpp.

parentAdjFunctionDim1

boost::shared_ptr<ParentFiniteElementAdjacencyFunction<1> > MoFEM::Simple::parentAdjFunctionDim1

private

Definition at line 667 of file Simple.hpp.

parentAdjFunctionDim2

boost::shared_ptr<ParentFiniteElementAdjacencyFunction<2> > MoFEM::Simple::parentAdjFunctionDim2

private

Definition at line 665 of file Simple.hpp.

parentAdjFunctionDim3

boost::shared_ptr<ParentFiniteElementAdjacencyFunction<3> > MoFEM::Simple::parentAdjFunctionDim3

private

Definition at line 663 of file Simple.hpp.

parentAdjSkeletonFunctionDim1

boost::shared_ptr<ParentFiniteElementAdjacencyFunctionSkeleton<1> > MoFEM::Simple::parentAdjSkeletonFunctionDim1

private

Definition at line 672 of file Simple.hpp.

parentAdjSkeletonFunctionDim2

boost::shared_ptr<ParentFiniteElementAdjacencyFunctionSkeleton<2> > MoFEM::Simple::parentAdjSkeletonFunctionDim2

private

Definition at line 670 of file Simple.hpp.

skeletonFE

std::string MoFEM::Simple::skeletonFE = "sFE"

private

skeleton finite element

Definition at line 640 of file Simple.hpp.

skeletonFields

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

private

fields on the skeleton

Definition at line 627 of file Simple.hpp.

skeletonMeshset

EntityHandle MoFEM::Simple::skeletonMeshset

private

skeleton meshset with boundary

Definition at line 614 of file Simple.hpp.

---

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

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