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 (boost::typeindex::type_index type_index, UnknownInterface **iface) const

	Simple (const MoFEM::Core &core)

virtual	~Simple ()=default

MoFEMErrorCode	getOptions ()
	get options More...

MoFEMErrorCode	loadFile (const std::string options, const std::string mesh_file_name)
	Load mesh file. More...

MoFEMErrorCode	loadFile (const std::string mesh_file_name="")
	Load mesh file with parallel options if number of cores > 1. 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	removeDomainField (const std::string &name)
	Remove field form domain. More...

MoFEMErrorCode	removeBoundaryField (const std::string &name)
	Remove field form boundary. More...

MoFEMErrorCode	removeSkeletonField (const std::string &name)
	Remove field form skeleton. 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	setSkeletonAdjacency (int dim=-1, std::string fe_name="")
	Set the skeleton adjacency object. More...

MoFEMErrorCode	buildProblem ()
	Build problem. More...

MoFEMErrorCode	setUp (const PetscBool is_partitioned=PETSC_TRUE)
	Setup problem. More...

MoFEMErrorCode	reSetUp (bool only_dm=false)
	Rebuild internal MoFEM data structures. More...

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

SmartPetscObj< DM >	getDM ()
	Return smart DM object. More...

int	getDim () const
	Get the problem dimension. More...

void	setDim (int dim)
	Set the problem dimension. More...

EntityHandle &	getMeshSet ()
	Get the MeshSet object. More...

EntityHandle &	getBoundaryMeshSet ()
	Get the BoundaryMeshSet object. More...

EntityHandle &	getSkeletonMeshSet ()
	Get the SkeletonMeshSet object. More...

BitRefLevel &	getBitRefLevel ()
	Get the BitRefLevel. More...

BitRefLevel &	getBitRefLevelMask ()
	Get the BitRefLevel. More...

const std::string	getDomainFEName () const
	Get the Domain FE Name. More...

const std::string	getBoundaryFEName () const
	Get the Boundary FE Name. More...

const std::string	getSkeletonFEName () const
	Get the Skeleton FE Name. More...

const std::string	getProblemName () const
	Get the Problem Name. More...

std::string &	getDomainFEName ()
	Get the Domain FE Name. More...

std::string &	getBoundaryFEName ()
	Get the Boundary FE Name. More...

std::string &	getSkeletonFEName ()
	Get the Skeleton FE Name. More...

std::string &	getProblemName ()
	Get the Problem Name. More...

std::vector< std::string > &	getOtherFiniteElements ()
	Get the Other Finite Elements. More...

MoFEMErrorCode	deleteDM ()
	Delete dm. More...

MoFEMErrorCode	deleteFiniteElements ()
	Delete finite elements. More...

bool &	getAddSkeletonFE ()
	Get the addSkeletonFE. More...

bool &	getAddBoundaryFE ()
	Get the addSkeletonFE. More...

bool &	getParentAdjacencies ()
	Get the addParentAdjacencies. More...

auto &	getBitAdjParent ()
	bit ref level for parent More...

auto &	getBitAdjParentMask ()
	bit ref level for parent marent More...

auto &	getBitAdjEnt ()
	bit ref level for parent More...

auto &	getBitAdjEntMask ()
	bit ref level for parent marent More...

MoFEMErrorCode	addFieldToEmptyFieldBlocks (const std::string row_field, const std::string col_field) const
	add empty block to problem More...

Public Member Functions inherited from MoFEM::UnknownInterface
virtual MoFEMErrorCode	query_interface (boost::typeindex::type_index type_index, UnknownInterface **iface) const =0

template<class IFACE >
MoFEMErrorCode	registerInterface (bool error_if_registration_failed=true)
	Register interface. More...

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

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

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

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

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

virtual	~UnknownInterface ()=default

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 probelm. More...

BitRefLevel	bitLevelMask
	BitRefLevel of the probelm. More...

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 More...

EntityHandle	boundaryMeshset
	meshset with boundary More...

EntityHandle	skeletonMeshset
	skeleton meshset with boundary More...

bool	addSkeletonFE
	Add skeleton FE. More...

bool	addBoundaryFE
	Add boundary FE. More...

bool	addParentAdjacencies
	If set to true parent adjacencies are build. More...

BitRefLevel	bitAdjParent
	bit ref level for parent More...

BitRefLevel	bitAdjParentMask
	bit ref level for parent marent More...

BitRefLevel	bitAdjEnt
	bit ref level for parent More...

BitRefLevel	bitAdjEntMask
	bit ref level for parent marent 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::vector< std::string >	noFieldFields
	NOFIELD field name. More...

std::vector< std::string >	noFieldDataFields
	NOFIELD field name. More...

std::list< std::tuple< std::string, int, Range > >	fieldsOrder
	fields order More...

std::string	nameOfProblem
	problem name 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...

SmartPetscObj< DM >	dM
	Discrete manager (interface to PETSc/MoFEM functions) More...

boost::shared_ptr< ParentFiniteElementAdjacencyFunction< 3 > >	parentAdjFunctionDim3

boost::shared_ptr< ParentFiniteElementAdjacencyFunction< 2 > >	parentAdjFunctionDim2

boost::shared_ptr< ParentFiniteElementAdjacencyFunction< 1 > >	parentAdjFunctionDim1

Additional Inherited Members
Static Public Member Functions inherited from MoFEM::UnknownInterface
static MoFEMErrorCode	getLibVersion (Version &version)
	Get library version. More...

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

static MoFEMErrorCode	setFileVersion (moab::Interface &moab, Version version=Version(MoFEM_VERSION_MAJOR, MoFEM_VERSION_MINOR, MoFEM_VERSION_BUILD))
	Get database major version. More...

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

Detailed Description

Simple interface for fast problem set-up.

Examples: approx_sphere.cpp, boundary_marker.cpp, child_and_parent.cpp, contact.cpp, continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, eigen_elastic.cpp, field_evaluator.cpp, free_surface.cpp, hanging_node_approx.cpp, hcurl_check_approx_in_2d.cpp, hcurl_curl_operator.cpp, hdiv_divergence_operator.cpp, heat_equation.cpp, heat_method.cpp, hello_world.cpp, helmholtz.cpp, higher_derivatives.cpp, loop_entities.cpp, mesh_smoothing.cpp, mixed_poisson.cpp, nonlinear_elastic.cpp, partition_mesh.cpp, phase.cpp, photon_diffusion.cpp, plastic.cpp, plate.cpp, plot_base.cpp, poisson_2d_dis_galerkin.cpp, poisson_2d_homogeneous.cpp, reaction_diffusion.cpp, scalar_check_approximation.cpp, shallow_wave.cpp, simple_elasticity.cpp, simple_interface.cpp, simple_l2_only.cpp, test_cache_on_entities.cpp, testing_jacobian_of_hook_element.cpp, testing_jacobian_of_hook_scaled_with_density_element.cpp, thermo_elastic.cpp, and wave_equation.cpp.

Definition at line 26 of file Simple.hpp.

Constructor & Destructor Documentation

◆ Simple()

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

Definition at line 271 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("SimpleSetUp", 0, &MOFEM_EVENT_SimpleSetUP);
  PetscLogEventRegister("SimpleLoadMesh", 0, &MOFEM_EVENT_SimpleLoadMesh);
  PetscLogEventRegister("SimpleBuildFields", 0, &MOFEM_EVENT_SimpleBuildFields);
  PetscLogEventRegister("SimpleBuildFiniteElements", 0,
                        &MOFEM_EVENT_SimpleBuildFiniteElements);
  PetscLogEventRegister("SimpleSetUp", 0, &MOFEM_EVENT_SimpleBuildProblem);
  PetscLogEventRegister("SimpleKSPSolve", 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.

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: hcurl_curl_operator.cpp, hdiv_divergence_operator.cpp, mesh_smoothing.cpp, and simple_interface.cpp.

Definition at line 391 of file Simple.cpp.

                                   {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.add_field(name, space, base, nb_of_cooficients, tag_type, bh,
                           verb);
  boundaryFields.push_back(name);
  if (space == NOFIELD)
    SETERRQ(
        PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
        "NOFIELD space for boundary filed 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 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: hcurl_curl_operator.cpp, and hdiv_divergence_operator.cpp.

Definition at line 429 of file Simple.cpp.

                               {
 
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.add_field(name, space, base, nb_of_cooficients, tag_type, bh,
                           verb);
  if (space == NOFIELD)
    noFieldDataFields.push_back(name);
  else
    dataFields.push_back(name);
  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.

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: continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, field_evaluator.cpp, hcurl_check_approx_in_2d.cpp, hcurl_curl_operator.cpp, hdiv_divergence_operator.cpp, loop_entities.cpp, mesh_smoothing.cpp, reaction_diffusion.cpp, simple_elasticity.cpp, simple_interface.cpp, simple_l2_only.cpp, test_cache_on_entities.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 373 of file Simple.cpp.

                                 {
 
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.add_field(name, space, base, nb_of_cooficients, tag_type, bh,
                           verb);
  if (space == NOFIELD)
    noFieldFields.push_back(name);
  else
    domainFields.push_back(name);
  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 fields combinations are non zero. This is not always the case, to optimise code and reduce memory usage user can specifi which blocks are empty.

Parameters

row_field	row filed name
col_field	col field name

Returns: MoFEMErrorCode

Definition at line 815 of file Simple.cpp.

                                                                    {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  CHKERR m_field.getInterface<ProblemsManager>()->addFieldToEmptyFieldBlocks(
      getProblemName(), row_field, col_field);
  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.

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: continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, and simple_interface.cpp.

Definition at line 409 of file Simple.cpp.

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

◆ buildFields()

MoFEMErrorCode MoFEM::Simple::buildFields ( )

Build fields.

Returns: error code

Examples: mesh_smoothing.cpp, and simple_elasticity.cpp.

Definition at line 596 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);
  };
 
  auto make_no_field_ents = [&](auto &fields) {
    MoFEMFunctionBeginHot;
    for (auto &field : fields) {
      std::array<double, 6> coords = {0, 0, 0, 0, 0, 0};
      CHKERR m_field.create_vertices_and_add_to_field(field, coords.data(), 1);
      CHKERR comm_interface_ptr->makeFieldEntitiesMultishared(field, 0);
      CHKERR bit_ref_ptr->setFieldEntitiesBitRefLevel(field, bitLevel);
    }
    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);
 
  CHKERR make_no_field_ents(nofield_fields);
 
  // 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;
    if (!std::get<2>(t).empty()) {
      MOFEM_LOG_CHANNEL("SYNC");
      MOFEM_TAG_AND_LOG("SYNC", Sev::verbose, "Simple")
          << "To ents: " << std::endl
          << std::get<2>(t) << std::endl;
      MOFEM_LOG_SYNCHRONISE(m_field.get_comm());
    }
 
    if (std::get<2>(t).empty()) {
      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);
        }
      }
    } else {
      CHKERR m_field.set_field_order(std::get<2>(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.

Returns: error code

Examples: mesh_smoothing.cpp, and simple_elasticity.cpp.

Definition at line 683 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.empty()) {
    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.empty()) {
    CHKERR m_field.add_ents_to_finite_element_by_dim(skeletonMeshset, 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 709 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 824 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 871 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, false, 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);
}

◆ defineFiniteElements()

MoFEMErrorCode MoFEM::Simple::defineFiniteElements ( )

Define finite elements.

Returns: Error code

Examples: mesh_smoothing.cpp, and simple_elasticity.cpp.

Definition at line 492 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);
  }
  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 564 of file Simple.cpp.

                                                                   {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
  // Create dm instance
  dM = createSmartDM(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.c_str());
  if (addBoundaryFE || !boundaryFields.empty()) {
    CHKERR DMMoFEMAddElement(dM, boundaryFE.c_str());
  }
  if (addSkeletonFE || !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);
}

◆ deleteDM()

MoFEMErrorCode MoFEM::Simple::deleteDM ( )

Delete dm.

Returns: MoFEMErrorCode

Returns: MoFEMErrorCode

Definition at line 792 of file Simple.cpp.

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

◆ deleteFiniteElements()

MoFEMErrorCode MoFEM::Simple::deleteFiniteElements ( )

Delete finite elements.

Returns: MoFEMErrorCode

Returns: MoFEMErrorCode

Definition at line 803 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 905 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());
 
  Range verts;
  CHKERR m_field.get_moab().get_entities_by_type(0, MBVERTEX, verts);
  CHKERR pcomm->exchange_ghost_cells(getDim(), getDim() - 1, 1,
                                     3 /**get all adjacent ghosted entities */,
                                     true, false, meshSet ? &meshSet : nullptr);
 
  Range shared;
  CHKERR m_field.get_moab().get_entities_by_dimension(0, dIm, shared);
  for (auto d = dIm - 1; d >= 1; --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);
 
  MoFEMFunctionReturn(0);
}

◆ getAddBoundaryFE()

bool & MoFEM::Simple::getAddBoundaryFE ( )

inline

Get the addSkeletonFE.

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

Returns: true; false

Examples: simple_l2_only.cpp.

Definition at line 408 of file Simple.hpp.

408{ return addBoundaryFE; }

◆ getAddSkeletonFE()

bool & MoFEM::Simple::getAddSkeletonFE ( )

inline

Get the addSkeletonFE.

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

Returns: true; false

Examples: continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, and simple_l2_only.cpp.

Definition at line 397 of file Simple.hpp.

397{ return addSkeletonFE; }

◆ getBitAdjEnt()

auto & MoFEM::Simple::getBitAdjEnt ( )

inline

bit ref level for parent

Returns: auto&

Definition at line 439 of file Simple.hpp.

439{ return bitAdjEnt; }

◆ getBitAdjEntMask()

auto & MoFEM::Simple::getBitAdjEntMask ( )

inline

bit ref level for parent marent

Returns: auto&

Definition at line 446 of file Simple.hpp.

446{ return bitAdjEntMask; }

◆ getBitAdjParent()

auto & MoFEM::Simple::getBitAdjParent ( )

inline

bit ref level for parent

Returns: auto&

Definition at line 425 of file Simple.hpp.

425{ return bitAdjParent; }

◆ getBitAdjParentMask()

auto & MoFEM::Simple::getBitAdjParentMask ( )

inline

bit ref level for parent marent

Returns: auto&

Definition at line 432 of file Simple.hpp.

432{ return bitAdjParentMask; }

◆ getBitRefLevel()

BitRefLevel & MoFEM::Simple::getBitRefLevel ( )

inline

Get the BitRefLevel.

Returns: BitRefLevel

Definition at line 299 of file Simple.hpp.

299{ return bitLevel; }

◆ getBitRefLevelMask()

BitRefLevel & MoFEM::Simple::getBitRefLevelMask ( )

inline

Get the BitRefLevel.

Returns: BitRefLevel

Definition at line 306 of file Simple.hpp.

306{ return bitLevelMask; }

◆ getBoundaryFEName() [1/2]

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

inline

Get the Boundary FE Name.

Returns: std::string&

Definition at line 348 of file Simple.hpp.

348{ return boundaryFE; }

◆ getBoundaryFEName() [2/2]

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

inline

Get the Boundary FE Name.

Returns: const std::string

Examples: mesh_smoothing.cpp, and simple_interface.cpp.

Definition at line 320 of file Simple.hpp.

320{ return boundaryFE; }

◆ getBoundaryMeshSet()

EntityHandle & MoFEM::Simple::getBoundaryMeshSet ( )

inline

Get the BoundaryMeshSet object.

Returns: EntityHandle&

Definition at line 285 of file Simple.hpp.

285{ 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

Examples: field_evaluator.cpp.

Definition at line 264 of file Simple.hpp.

264{ return dIm; }

◆ getDM() [1/2]

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

inline

Return smart DM object.

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

Returns: SmartPetscObj<DM>

Definition at line 254 of file Simple.hpp.

254{ return dM; }

◆ getDM() [2/2]

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

Get DM.

Parameters

dm	discrete manager

Returns: error code

Examples: continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, field_evaluator.cpp, hcurl_check_approx_in_2d.cpp, loop_entities.cpp, mesh_smoothing.cpp, reaction_diffusion.cpp, simple_elasticity.cpp, simple_interface.cpp, test_cache_on_entities.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 780 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 341 of file Simple.hpp.

341{ return domainFE; }

◆ getDomainFEName() [2/2]

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

inline

Get the Domain FE Name.

Returns: const std::string

Examples: field_evaluator.cpp, mesh_smoothing.cpp, reaction_diffusion.cpp, simple_elasticity.cpp, simple_interface.cpp, simple_l2_only.cpp, test_cache_on_entities.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 313 of file Simple.hpp.

313{ return domainFE; }

◆ getMeshSet()

EntityHandle & MoFEM::Simple::getMeshSet ( )

inline

Get the MeshSet object.

Returns: EntityHandle&

Definition at line 278 of file Simple.hpp.

278{ return meshSet; }

◆ getOptions()

MoFEMErrorCode MoFEM::Simple::getOptions ( )

get options

Returns: error code

Examples: continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, field_evaluator.cpp, hcurl_check_approx_in_2d.cpp, hcurl_curl_operator.cpp, hdiv_divergence_operator.cpp, loop_entities.cpp, mesh_smoothing.cpp, reaction_diffusion.cpp, simple_elasticity.cpp, simple_interface.cpp, simple_l2_only.cpp, test_cache_on_entities.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 289 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 ( )

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: mesh_smoothing.cpp.

Definition at line 372 of file Simple.hpp.

372{ return otherFEs; }

◆ getParentAdjacencies()

bool & MoFEM::Simple::getParentAdjacencies ( )

inline

Get the addParentAdjacencies.

If set true add parent adjacencies

Returns: true; false

Definition at line 418 of file Simple.hpp.

418{ return addParentAdjacencies; }

◆ getProblemName() [1/2]

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

inline

Get the Problem Name.

Returns: std::string&

Definition at line 362 of file Simple.hpp.

362{ return nameOfProblem; }

◆ getProblemName() [2/2]

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

inline

Get the Problem Name.

Returns: const std::string

Examples: reaction_diffusion.cpp, and test_cache_on_entities.cpp.

Definition at line 334 of file Simple.hpp.

334{ return nameOfProblem; }

◆ getSkeletonFEName() [1/2]

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

inline

Get the Skeleton FE Name.

Returns: std::string&

Definition at line 355 of file Simple.hpp.

355{ return skeletonFE; }

◆ getSkeletonFEName() [2/2]

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

inline

Get the Skeleton FE Name.

Returns: const std::string

Examples: continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, and simple_interface.cpp.

Definition at line 327 of file Simple.hpp.

327{ return skeletonFE; }

◆ getSkeletonMeshSet()

EntityHandle & MoFEM::Simple::getSkeletonMeshSet ( )

inline

Get the SkeletonMeshSet object.

Returns: EntityHandle&

Definition at line 292 of file Simple.hpp.

292{ 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.

Parameters

mesh_file_name file name if not set default or set by command line is used.

Returns: error code

Definition at line 359 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
	)

Load mesh file.

Parameters

options	file load options
mesh_file_name	file name if not set default or set by command line is used.

Returns: error code

Examples: continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, field_evaluator.cpp, hcurl_check_approx_in_2d.cpp, hcurl_curl_operator.cpp, hdiv_divergence_operator.cpp, loop_entities.cpp, mesh_smoothing.cpp, reaction_diffusion.cpp, simple_elasticity.cpp, simple_interface.cpp, simple_l2_only.cpp, test_cache_on_entities.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 303 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());
 
  // This is a case of distributed mesh and algebra. In that case each processor
  // keep only part of the problem.
  CHKERR m_field.get_moab().load_file(meshFileName, 0, 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 form boundary.

Parameters

name

Returns: MoFEMErrorCode

Definition at line 462 of file Simple.cpp.

                                                                {
  Interface &m_field = cOre;
  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 form domain.

Parameters

name

Returns: MoFEMErrorCode

Definition at line 446 of file Simple.cpp.

                                                              {
  Interface &m_field = cOre;
  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(noFieldFields);
  remove_field_from_list(domainFields);
 
  MoFEMFunctionReturn(0);
}

◆ removeSkeletonField()

MoFEMErrorCode MoFEM::Simple::removeSkeletonField ( const std::string & name )

Remove field form skeleton.

Parameters

name

Returns: MoFEMErrorCode

Definition at line 477 of file Simple.cpp.

                                                                {
  Interface &m_field = cOre;
  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.

Returns: MoFEMErrorCode

Definition at line 746 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();
    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.

Returns: int

Definition at line 271 of file Simple.hpp.

271{ dIm = dim; };

◆ 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: continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, field_evaluator.cpp, hcurl_check_approx_in_2d.cpp, hcurl_curl_operator.cpp, hdiv_divergence_operator.cpp, loop_entities.cpp, mesh_smoothing.cpp, reaction_diffusion.cpp, simple_elasticity.cpp, simple_interface.cpp, simple_l2_only.cpp, test_cache_on_entities.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 588 of file Simple.cpp.

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

◆ setParentAdjacency() [1/4]

template<>

MoFEMErrorCode MoFEM::Simple::setParentAdjacency ( )

private

Definition at line 181 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 212 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 237 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 176 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.

Parameters

dim
fe_name

Returns: MoFEMErrorCode

Definition at line 252 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;
 
  auto defaultSkeletonEdge =
      [&](moab::Interface &moab, const Field &field, const EntFiniteElement &fe,
          std::vector<EntityHandle> &adjacency) -> MoFEMErrorCode {
    MoFEMFunctionBegin;
 
    CHKERR DefaultElementAdjacency::defaultEdge(moab, field, fe, adjacency);
 
    if (std::find(domainFields.begin(), domainFields.end(), field.getName()) !=
        domainFields.end()) {
 
      const EntityHandle fe_ent = fe.getEnt();
      std::vector<EntityHandle> bride_adjacency_edge;
      CHKERR moab.get_adjacencies(&fe_ent, 1, 2, false, bride_adjacency_edge);
 
      switch (field.getSpace()) {
      case H1:
        CHKERR moab.get_connectivity(&*bride_adjacency_edge.begin(),
                                     bride_adjacency_edge.size(), adjacency,
                                     true);
      case HCURL:
      case HDIV:
        CHKERR moab.get_adjacencies(&*bride_adjacency_edge.begin(),
                                    bride_adjacency_edge.size(), 1, false,
                                    adjacency, moab::Interface::UNION);
      case L2:
        adjacency.insert(adjacency.end(), bride_adjacency_edge.begin(),
                         bride_adjacency_edge.end());
        break;
      case NOFIELD:
        break;
      default:
        SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
                "this field is not implemented for TRI finite element");
      }
 
      std::sort(adjacency.begin(), adjacency.end());
      auto it = std::unique(adjacency.begin(), adjacency.end());
 
      std::vector<EntityHandle> new_adjacency(
          std::distance(adjacency.begin(), it));
      std::copy(adjacency.begin(), it, new_adjacency.begin());
 
      for (auto e : new_adjacency) {
        auto side_table = fe.getSideNumberTable();
        if (side_table.find(e) == side_table.end())
          const_cast<SideNumber_multiIndex &>(fe.getSideNumberTable())
              .insert(
                  boost::shared_ptr<SideNumber>(new SideNumber(e, -1, 0, 0)));
      }
 
      adjacency.swap(new_adjacency);
    }
 
    MoFEMFunctionReturn(0);
  };
 
  CHKERR m_field.modify_finite_element_adjacency_table(fe_name, MBEDGE,
                                                       defaultSkeletonEdge);
 
  MoFEMFunctionReturn(0);
}

◆ setSkeletonAdjacency() [3/5]

template<>

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

private

Definition at line 88 of file Simple.cpp.

                                                                {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  auto defaultSkeletonEdge =
      [&](moab::Interface &moab, const Field &field, const EntFiniteElement &fe,
          std::vector<EntityHandle> &adjacency) -> MoFEMErrorCode {
    MoFEMFunctionBegin;
 
    CHKERR DefaultElementAdjacency::defaultFace(moab, field, fe, adjacency);
 
    if (std::find(domainFields.begin(), domainFields.end(), field.getName()) !=
        domainFields.end()) {
 
      const EntityHandle fe_ent = fe.getEnt();
      std::vector<EntityHandle> bride_adjacency_edge;
      CHKERR moab.get_adjacencies(&fe_ent, 1, 2, false, bride_adjacency_edge);
 
      switch (field.getSpace()) {
      case H1:
        CHKERR moab.get_connectivity(&*bride_adjacency_edge.begin(),
                                     bride_adjacency_edge.size(), adjacency,
                                     true);
      case HCURL:
        CHKERR moab.get_adjacencies(&*bride_adjacency_edge.begin(),
                                    bride_adjacency_edge.size(), 1, false,
                                    adjacency, moab::Interface::UNION);
      case HDIV:
        CHKERR moab.get_adjacencies(&*bride_adjacency_edge.begin(),
                                    bride_adjacency_edge.size(), 2, false,
                                    adjacency, moab::Interface::UNION);
      case L2:
        adjacency.insert(adjacency.end(), bride_adjacency_edge.begin(),
                         bride_adjacency_edge.end());
        break;
      case NOFIELD:
        break;
      default:
        SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
                "this field is not implemented for TRI finite element");
      }
 
      std::sort(adjacency.begin(), adjacency.end());
      auto it = std::unique(adjacency.begin(), adjacency.end());
 
      std::vector<EntityHandle> new_adjacency(
          std::distance(adjacency.begin(), it));
      std::copy(adjacency.begin(), it, new_adjacency.begin());
 
      for (auto e : new_adjacency) {
        auto side_table = fe.getSideNumberTable();
        if (side_table.find(e) == side_table.end())
          const_cast<SideNumber_multiIndex &>(fe.getSideNumberTable())
              .insert(
                  boost::shared_ptr<SideNumber>(new SideNumber(e, -1, 0, 0)));
      }
 
      adjacency.swap(new_adjacency);
    }
 
    MoFEMFunctionReturn(0);
  };
 
  CHKERR m_field.modify_finite_element_adjacency_table(fe_name, MBTRI,
                                                       defaultSkeletonEdge);
  CHKERR m_field.modify_finite_element_adjacency_table(fe_name, MBQUAD,
                                                       defaultSkeletonEdge);
 
  MoFEMFunctionReturn(0);
}

◆ setSkeletonAdjacency() [4/5]

template<>

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

private

Definition at line 160 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.

Returns: error code

Examples: continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, field_evaluator.cpp, hcurl_check_approx_in_2d.cpp, hcurl_curl_operator.cpp, hdiv_divergence_operator.cpp, loop_entities.cpp, reaction_diffusion.cpp, simple_interface.cpp, simple_l2_only.cpp, test_cache_on_entities.cpp, testing_jacobian_of_hook_element.cpp, and testing_jacobian_of_hook_scaled_with_density_element.cpp.

Definition at line 723 of file Simple.cpp.

                                                           {
  Interface &m_field = cOre;
  MoFEMFunctionBegin;
 
  PetscLogEventBegin(MOFEM_EVENT_SimpleSetUP, 0, 0, 0, 0);
 
  CHKERR defineFiniteElements();
 
  if (addSkeletonFE || !skeletonFields.empty())
    CHKERR setSkeletonAdjacency();
 
  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 481 of file Simple.hpp.

◆ addParentAdjacencies

bool MoFEM::Simple::addParentAdjacencies

private

If set to true parent adjacencies are build.

Definition at line 482 of file Simple.hpp.

◆ addSkeletonFE

bool MoFEM::Simple::addSkeletonFE

private

Add skeleton FE.

Definition at line 480 of file Simple.hpp.

◆ bitAdjEnt

BitRefLevel MoFEM::Simple::bitAdjEnt

private

bit ref level for parent

Definition at line 486 of file Simple.hpp.

◆ bitAdjEntMask

BitRefLevel MoFEM::Simple::bitAdjEntMask

private

bit ref level for parent marent

Definition at line 487 of file Simple.hpp.

◆ bitAdjParent

BitRefLevel MoFEM::Simple::bitAdjParent

private

bit ref level for parent

Definition at line 484 of file Simple.hpp.

◆ bitAdjParentMask

BitRefLevel MoFEM::Simple::bitAdjParentMask

private

bit ref level for parent marent

Definition at line 485 of file Simple.hpp.

◆ bitLevel

BitRefLevel MoFEM::Simple::bitLevel

private

BitRefLevel of the probelm.

Definition at line 466 of file Simple.hpp.

◆ bitLevelMask

BitRefLevel MoFEM::Simple::bitLevelMask

private

BitRefLevel of the probelm.

Definition at line 467 of file Simple.hpp.

◆ boundaryFE

std::string MoFEM::Simple::boundaryFE

private

boundary finite element

Definition at line 501 of file Simple.hpp.

◆ boundaryFields

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

private

boundary fields

Definition at line 490 of file Simple.hpp.

◆ boundaryMeshset

EntityHandle MoFEM::Simple::boundaryMeshset

private

meshset with boundary

Definition at line 477 of file Simple.hpp.

◆ cOre

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

private

Definition at line 464 of file Simple.hpp.

◆ dataFields

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

private

Data fields.

Definition at line 492 of file Simple.hpp.

◆ dIm

int MoFEM::Simple::dIm

private

dimension of problem

Definition at line 507 of file Simple.hpp.

◆ dM

SmartPetscObj<DM> MoFEM::Simple::dM

private

Discrete manager (interface to PETSc/MoFEM functions)

Definition at line 510 of file Simple.hpp.

◆ domainFE

std::string MoFEM::Simple::domainFE

private

domain finite element

Definition at line 500 of file Simple.hpp.

◆ domainFields

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

private

domain fields

Definition at line 489 of file Simple.hpp.

◆ fieldsOrder

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

private

fields order

Definition at line 497 of file Simple.hpp.

◆ meshFileName

char MoFEM::Simple::meshFileName[255]

private

mesh file name

Definition at line 506 of file Simple.hpp.

◆ meshSet

EntityHandle MoFEM::Simple::meshSet

private

domain meshset

Definition at line 476 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleBuildFields

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleBuildFields

private

Definition at line 471 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleBuildFiniteElements

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleBuildFiniteElements

private

Definition at line 472 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleBuildProblem

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleBuildProblem

private

Definition at line 473 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleKSPSolve

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleKSPSolve

private

Definition at line 474 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleLoadMesh

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleLoadMesh

private

Definition at line 470 of file Simple.hpp.

◆ MOFEM_EVENT_SimpleSetUP

PetscLogEvent MoFEM::Simple::MOFEM_EVENT_SimpleSetUP

private

Definition at line 469 of file Simple.hpp.

◆ nameOfProblem

std::string MoFEM::Simple::nameOfProblem

private

problem name

Definition at line 499 of file Simple.hpp.

◆ noFieldDataFields

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

private

NOFIELD field name.

Definition at line 494 of file Simple.hpp.

◆ noFieldFields

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

private

NOFIELD field name.

Definition at line 493 of file Simple.hpp.

◆ otherFEs

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

private

Other finite elements.

Definition at line 504 of file Simple.hpp.

◆ parentAdjFunctionDim1

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

private

Definition at line 526 of file Simple.hpp.

◆ parentAdjFunctionDim2

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

private

Definition at line 524 of file Simple.hpp.

◆ parentAdjFunctionDim3

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

private

Definition at line 522 of file Simple.hpp.

◆ skeletonFE

std::string MoFEM::Simple::skeletonFE

private

skeleton finite element

Definition at line 502 of file Simple.hpp.

◆ skeletonFields

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

private

fields on the skeleton

Definition at line 491 of file Simple.hpp.

◆ skeletonMeshset

EntityHandle MoFEM::Simple::skeletonMeshset

private

skeleton meshset with boundary

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

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ Simple()

◆ ~Simple()

Member Function Documentation

◆ addBoundaryField()

◆ addDataField()

◆ addDomainField()

◆ addFieldToEmptyFieldBlocks()

◆ addSkeletonField()

◆ buildFields()

◆ buildFiniteElements()

◆ buildProblem()

◆ createBoundaryMeshset()

◆ createSkeletonMeshset()

◆ defineFiniteElements()

◆ defineProblem()

◆ deleteDM()

◆ deleteFiniteElements()

◆ exchangeGhostCells()

◆ getAddBoundaryFE()

◆ getAddSkeletonFE()

◆ getBitAdjEnt()

◆ getBitAdjEntMask()

◆ getBitAdjParent()

◆ getBitAdjParentMask()

◆ getBitRefLevel()

◆ getBitRefLevelMask()

◆ getBoundaryFEName() [1/2]

◆ getBoundaryFEName() [2/2]

◆ getBoundaryMeshSet()

◆ getDim()

◆ getDM() [1/2]

◆ getDM() [2/2]

◆ getDomainFEName() [1/2]

◆ getDomainFEName() [2/2]

◆ getMeshSet()

◆ getOptions()

◆ getOtherFiniteElements()

◆ getParentAdjacencies()

◆ getProblemName() [1/2]

◆ getProblemName() [2/2]

◆ getSkeletonFEName() [1/2]

◆ getSkeletonFEName() [2/2]

◆ getSkeletonMeshSet()

◆ loadFile() [1/2]

◆ loadFile() [2/2]

◆ query_interface()

◆ removeBoundaryField()

◆ removeDomainField()

◆ removeSkeletonField()

◆ reSetUp()

◆ setDim()

◆ setFieldOrder()

◆ setParentAdjacency() [1/4]

◆ setParentAdjacency() [2/4]

◆ setParentAdjacency() [3/4]

◆ setParentAdjacency() [4/4]

◆ setSkeletonAdjacency() [1/5]

◆ setSkeletonAdjacency() [2/5]

◆ setSkeletonAdjacency() [3/5]

◆ setSkeletonAdjacency() [4/5]

◆ setSkeletonAdjacency() [5/5]

◆ setUp()

Member Data Documentation

◆ addBoundaryFE

◆ addParentAdjacencies

◆ addSkeletonFE

◆ bitAdjEnt

◆ bitAdjEntMask

◆ bitAdjParent

◆ bitAdjParentMask

◆ bitLevel

◆ bitLevelMask

◆ boundaryFE

◆ boundaryFields