structure to get information form mofem into DataForcesAndSourcesCore More...

#include <src/finite_elements/ForcesAndSourcesCore.hpp>

Inheritance diagram for MoFEM::ForcesAndSourcesCore:

Collaboration diagram for MoFEM::ForcesAndSourcesCore:

Classes
struct	UserDataOperator
	Data operator to do calculations at integration points.Is inherited and implemented by user to do calculations. It can be used in many different ways but typically is used to integrate matrices (f.e. stiffness matrix) and the right hand vector (f.e. force vector). More...

Public Types
typedef boost::function< int(int order_row, int order_col, int order_data)>	RuleHookFun

Public Types inherited from MoFEM::KspMethod
enum	KSPContext { CTX_SETFUNCTION, CTX_OPERATORS, CTX_KSPNONE }
	pass information about context of KSP/DM for with finite element is computed More...

Public Types inherited from MoFEM::SnesMethod
enum	SNESContext { CTX_SNESSETFUNCTION, CTX_SNESSETJACOBIAN, CTX_SNESNONE }

Public Types inherited from MoFEM::TSMethod
enum	TSContext { CTX_TSSETRHSFUNCTION, CTX_TSSETRHSJACOBIAN, CTX_TSSETIFUNCTION, CTX_TSSETIJACOBIAN, CTX_TSTSMONITORSET, CTX_TSNONE }

Public Member Functions
	ForcesAndSourcesCore (Interface &m_field)

virtual	~ForcesAndSourcesCore ()

MoFEMErrorCode	getNumberOfNodes (int &num_nodes) const

int	getMaxDataOrder () const
	Get max order of approximation for data fields. More...

int	getMaxRowOrder () const
	Get max order of approximation for field in rows. More...

int	getMaxColOrder () const
	Get max order of approximation for field in columns. More...

MoFEMErrorCode	getSense (EntityType type, boost::ptr_vector< DataForcesAndSourcesCore::EntData > &data) const
	get sense (orientation) of entity More...

MoFEMErrorCode	getDataOrder (const EntityType type, const FieldSpace space, boost::ptr_vector< DataForcesAndSourcesCore::EntData > &data) const
	get maximal approximation order of approximation on the entity More...

template<EntityType type>
MoFEMErrorCode	getEntitySense (DataForcesAndSourcesCore &data) const

template<EntityType type>
MoFEMErrorCode	getEntityFieldDataOrder (DataForcesAndSourcesCore &data, const FieldSpace space) const

MoFEMErrorCode	getNodesIndices (const boost::string_ref field_name, FENumeredDofEntity_multiIndex &dofs, VectorInt &nodes_indices, VectorInt &local_nodes_indices) const
	get node indices More...

MoFEMErrorCode	getRowNodesIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get row node indices from FENumeredDofEntity_multiIndex More...

MoFEMErrorCode	getColNodesIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get col node indices from FENumeredDofEntity_multiIndex More...

MoFEMErrorCode	getEntityIndices (DataForcesAndSourcesCore &data, const std::string &field_name, FENumeredDofEntity_multiIndex &dofs, const EntityType type_lo=MBVERTEX, const EntityType type_hi=MBPOLYHEDRON) const

MoFEMErrorCode	getEntityRowIndices (DataForcesAndSourcesCore &data, const std::string &field_name, const EntityType type_lo=MBVERTEX, const EntityType type_hi=MBPOLYHEDRON) const

MoFEMErrorCode	getEntityColIndices (DataForcesAndSourcesCore &data, const std::string &field_name, const EntityType type_lo=MBVERTEX, const EntityType type_hi=MBPOLYHEDRON) const

MoFEMErrorCode	getNoFieldIndices (const std::string &field_name, FENumeredDofEntity_multiIndex &dofs, VectorInt &nodes_indices) const
	get NoField indices More...

MoFEMErrorCode	getNoFieldRowIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get col NoField indices More...

MoFEMErrorCode	getNoFieldColIndices (DataForcesAndSourcesCore &data, const std::string &field_name) const
	get col NoField indices More...

MoFEMErrorCode	getNodesFieldData (const boost::string_ref field_name, FEDofEntity_multiIndex &dofs, VectorDouble &nodes_data, VectorDofs &nodes_dofs, FieldSpace &space, FieldApproximationBase &base) const
	Get field data on nodes. More...

MoFEMErrorCode	getNoFieldFieldData (const boost::string_ref field_name, FEDofEntity_multiIndex &dofs, VectorDouble &ent_field_data, VectorDofs &ent_field_dofs) const

MoFEMErrorCode	getNoFieldFieldData (DataForcesAndSourcesCore &data, const boost::string_ref field_name) const

MoFEMErrorCode	getNodesFieldData (DataForcesAndSourcesCore &data, const std::string &field_name) const
	Get data on nodes. More...

MoFEMErrorCode	getEntityFieldData (DataForcesAndSourcesCore &data, const std::string &field_name, const EntityType type_lo=MBVERTEX, const EntityType type_hi=MBPOLYHEDRON) const

MoFEMErrorCode	getFaceTriNodes (DataForcesAndSourcesCore &data) const
	Get nodes on triangles. More...

MoFEMErrorCode	getSpacesAndBaseOnEntities (DataForcesAndSourcesCore &data) const
	Get field approximation space and base on entities. More...

MoFEMErrorCode	getProblemNodesIndices (const std::string &field_name, const NumeredDofEntity_multiIndex &dofs, VectorInt &nodes_indices) const
	get indices of nodal indices which are declared for problem but not this particular element More...

MoFEMErrorCode	getProblemTypeIndices (const std::string &field_name, const NumeredDofEntity_multiIndex &dofs, EntityType type, int side_number, VectorInt &indices) const
	get indices by type (generic function) which are declared for problem but not this particular element More...

MoFEMErrorCode	getProblemNodesRowIndices (const std::string &field_name, VectorInt &nodes_indices) const

MoFEMErrorCode	getProblemTypeRowIndices (const std::string &field_name, EntityType type, int side_number, VectorInt &indices) const

MoFEMErrorCode	getProblemNodesColIndices (const std::string &field_name, VectorInt &nodes_indices) const

MoFEMErrorCode	getProblemTypeColIndices (const std::string &field_name, EntityType type, int side_number, VectorInt &indices) const

virtual int	getRule (int order)
	set integration rule for finite element More...

virtual int	getRule (int order_row, int order_col, int order_data)
	another variant of getRule More...

virtual MoFEMErrorCode	setGaussPts (int order)
	It will be removed in the future use other variant. More...

virtual MoFEMErrorCode	setGaussPts (int order_row, int order_col, int order_data)
	set user specific integration rule More...

boost::ptr_vector< UserDataOperator > &	getOpPtrVector ()
	Use to push back operator for row operator. More...

auto &	getElementPolynomialBase ()
	Get the Entity Polynomial Base object. More...

auto &	getUserPolynomialBase ()
	Get the User Polynomial Base object. More...

virtual MoFEMErrorCode	preProcess ()
	function is run at the beginning of loop More...

virtual MoFEMErrorCode	operator() ()
	function is run for every finite element More...

virtual MoFEMErrorCode	postProcess ()
	function is run at the end of loop More...

MoFEMErrorCode	calculateBaseFunctionsOnElement (const FieldApproximationBase b)
	Calculate base functions. More...

MoFEMErrorCode	calculateBaseFunctionsOnElement ()
	Calculate base functions. More...

MoFEMErrorCode	createDataOnElement ()
	Create a entity data on element object. More...

MoFEMErrorCode	loopOverOperators ()
	Iterate user data operators. More...

Public Member Functions inherited from MoFEM::FEMethod
MoFEMErrorCode	query_interface (const MOFEMuuid &uuid, UnknownInterface **iface) const

	FEMethod ()

template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const std::string &field_name, const EntityType type, const int side_number) const

template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const std::string &field_name, const EntityType type, const int side_number) const

template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const std::string &field_name, const EntityType type) const

template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const std::string &field_name, const EntityType type) const

template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const std::string &field_name) const

template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const std::string &field_name) const

template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const EntityHandle ent) const

template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const EntityHandle ent) const

template<class MULTIINDEX >
MULTIINDEX::iterator	get_begin (const MULTIINDEX &index, const std::string &field_name, const EntityHandle ent) const

template<class MULTIINDEX >
MULTIINDEX::iterator	get_end (const MULTIINDEX &index, const std::string &field_name, const EntityHandle ent) const

Public Member Functions inherited from MoFEM::BasicMethod
	BasicMethod ()

virtual	~BasicMethod ()

int	getNinTheLoop () const
	get number of evaluated element in the loop More...

int	getLoopSize () const
	get loop size More...

MoFEMErrorCode	copyBasicMethod (const BasicMethod &basic)

Public Member Functions inherited from MoFEM::KspMethod
	KspMethod ()

virtual	~KspMethod ()

MoFEMErrorCode	setKspCtx (const KSPContext &ctx)
	set operator type More...

MoFEMErrorCode	setKsp (KSP ksp)
	set solver More...

MoFEMErrorCode	copyKsp (const KspMethod &ksp)
	copy data form another method More...

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

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

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

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

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

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

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

virtual	~UnknownInterface ()

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

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

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

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

Public Member Functions inherited from MoFEM::SnesMethod
	SnesMethod ()

virtual	~SnesMethod ()

MoFEMErrorCode	setSnesCtx (const SNESContext &ctx)
	Set SNES context. More...

MoFEMErrorCode	setSnes (SNES snes)
	Set SNES instance. More...

MoFEMErrorCode	copySnes (const SnesMethod &snes)
	Copy snes data. More...

Public Member Functions inherited from MoFEM::TSMethod
	TSMethod ()

virtual	~TSMethod ()

MoFEMErrorCode	setTsCtx (const TSContext &ctx)
	Set Ts context. More...

MoFEMErrorCode	copyTs (const TSMethod &ts)
	Copy TS solver data. More...

MoFEMErrorCode	setTs (TS _ts)
	Set TS solver. More...

Public Attributes
Interface &	mField

const boost::shared_ptr< DataForcesAndSourcesCore >	dataOnElement [LASTSPACE]
	Entity data on element entity rows fields. More...

const boost::shared_ptr< DataForcesAndSourcesCore >	derivedDataOnElement [LASTSPACE]
	Entity data on element entity columns fields. More...

DataForcesAndSourcesCore &	dataNoField

DataForcesAndSourcesCore &	dataH1

DataForcesAndSourcesCore &	dataHcurl

DataForcesAndSourcesCore &	dataHdiv

DataForcesAndSourcesCore &	dataL2

RuleHookFun	getRuleHook
	Hook to get rule. More...

boost::ptr_vector< UserDataOperator >	opPtrVector
	Vector of finite element users data operators. More...

MatrixDouble	gaussPts
	Matrix of integration points. More...

Public Attributes inherited from MoFEM::FEMethod
std::string	feName

boost::shared_ptr< const NumeredEntFiniteElement >	numeredEntFiniteElementPtr

boost::shared_ptr< const FENumeredDofEntity_multiIndex >	rowPtr

boost::shared_ptr< const FENumeredDofEntity_multiIndex >	colPtr

boost::shared_ptr< const FEDofEntity_multiIndex >	dataPtr

Public Attributes inherited from MoFEM::BasicMethod
int	nInTheLoop

int	loopSize

int	rAnk

int	sIze

const RefEntity_multiIndex *	refinedEntitiesPtr

const RefElement_multiIndex *	refinedFiniteElementsPtr

const Problem *	problemPtr

const Field_multiIndex *	fieldsPtr

const FieldEntity_multiIndex *	entitiesPtr

const DofEntity_multiIndex *	dofsPtr

const FiniteElement_multiIndex *	finiteElementsPtr

const EntFiniteElement_multiIndex *	finiteElementsEntitiesPtr

const FieldEntityEntFiniteElementAdjacencyMap_multiIndex *	adjacenciesPtr

boost::function< MoFEMErrorCode()>	preProcessHook

boost::function< MoFEMErrorCode()>	postProcessHook

boost::function< MoFEMErrorCode()>	operatorHook

Public Attributes inherited from MoFEM::KspMethod
KSPContext	ksp_ctx
	Context. More...

KSP	ksp
	KSP solver. More...

Vec	ksp_f
	the right hand side vector More...

Mat	ksp_A
	matrix More...

Mat	ksp_B
	preconditioner matrix More...

Public Attributes inherited from MoFEM::SnesMethod
SNESContext	snes_ctx

SNES	snes

Vec	snes_x

Vec	snes_f

Mat	snes_A

Mat	snes_B

Public Attributes inherited from MoFEM::TSMethod
TSContext	ts_ctx

TS	ts

Vec	ts_u

Vec	ts_u_t

Vec	ts_F

Mat	ts_A

Mat	ts_B

PetscInt	ts_step

PetscReal	ts_a

PetscReal	ts_t

Private Attributes
EntityType	lastEvaluatedElementEntityType
	Last evaluated type of element entity. More...

boost::shared_ptr< BaseFunction >	elementPolynomialBasePtr
	Pointer to entity polynomial base. More...

boost::shared_ptr< BaseFunction >	userPolynomialBasePtr
	Pointer to user polynomail base. More...

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

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

Detailed Description

structure to get information form mofem into DataForcesAndSourcesCore

Examples:: hello_world.cpp.

Definition at line 34 of file ForcesAndSourcesCore.hpp.

Member Typedef Documentation

◆ RuleHookFun

typedef boost::function<int(int order_row, int order_col, int order_data)> MoFEM::ForcesAndSourcesCore::RuleHookFun

Definition at line 278 of file ForcesAndSourcesCore.hpp.

Constructor & Destructor Documentation

◆ ForcesAndSourcesCore()

MoFEM::ForcesAndSourcesCore::ForcesAndSourcesCore ( Interface & m_field )

Definition at line 72 of file ForcesAndSourcesCore.cpp.

     :
 
       mField(m_field),
       dataOnElement{
 
           nullptr,
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DataForcesAndSourcesCore(MBENTITYSET)), // NOFIELD
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DataForcesAndSourcesCore(MBENTITYSET)), // H1
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DataForcesAndSourcesCore(MBENTITYSET)), // HCURL
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DataForcesAndSourcesCore(MBENTITYSET)), // HDIV
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DataForcesAndSourcesCore(MBENTITYSET)) // L2
 
       },
       derivedDataOnElement{
 
           nullptr,
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DerivedDataForcesAndSourcesCore(dataOnElement[NOFIELD])),
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DerivedDataForcesAndSourcesCore(dataOnElement[H1])),
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DerivedDataForcesAndSourcesCore(dataOnElement[HCURL])),
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DerivedDataForcesAndSourcesCore(dataOnElement[HDIV])),
           boost::shared_ptr<DataForcesAndSourcesCore>(
               new DerivedDataForcesAndSourcesCore(dataOnElement[L2]))
 
       },
       dataNoField(*dataOnElement[NOFIELD].get()),
       dataH1(*dataOnElement[H1].get()), dataHcurl(*dataOnElement[HCURL].get()),
       dataHdiv(*dataOnElement[HDIV].get()), dataL2(*dataOnElement[L2].get()),
       getRuleHook(0), lastEvaluatedElementEntityType(MBMAXTYPE) {}

◆ ~ForcesAndSourcesCore()

virtual MoFEM::ForcesAndSourcesCore::~ForcesAndSourcesCore ( )

virtual

Definition at line 63 of file ForcesAndSourcesCore.hpp.

63 {}

Member Function Documentation

◆ calculateBaseFunctionsOnElement() [1/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::calculateBaseFunctionsOnElement ( const FieldApproximationBase b )

Calculate base functions.

Returns: Error code

Definition at line 840 of file ForcesAndSourcesCore.cpp.

                                     {
   MoFEMFunctionBegin;
   if (dataOnElement[H1]->bAse.test(b)) {
     switch (static_cast<FieldApproximationBase>(b)) {
     case NOBASE:
       break;
     case AINSWORTH_LEGENDRE_BASE:
     case AINSWORTH_LOBATTO_BASE:
     case DEMKOWICZ_JACOBI_BASE:
       if (!getElementPolynomialBase()) {
         SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                 "Functions genrating approximation base not defined");
       }
       for (int space = H1; space != LASTSPACE; ++space) {
         if (dataOnElement[H1]->sPace.test(space) &&
             dataOnElement[H1]->bAse.test(b) &&
             dataOnElement[H1]->basesOnSpaces[space].test(b)) {
           CHKERR getElementPolynomialBase()->getValue(
               gaussPts,
               boost::shared_ptr<BaseFunctionCtx>(new EntPolynomialBaseCtx(
                   *dataOnElement[space], static_cast<FieldSpace>(space),
                   static_cast<FieldApproximationBase>(b), NOBASE)));
         }
       }
       break;
     case USER_BASE:
       if (!getUserPolynomialBase()) {
         SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                 "Functions genrating user approximation base not defined");
       }
       for (int space = H1; space != LASTSPACE; ++space)
         if (dataOnElement[H1]->sPace.test(space) &&
             dataOnElement[H1]->bAse.test(b) &&
             dataOnElement[H1]->basesOnSpaces[space].test(b)) {
           CHKERR getUserPolynomialBase()->getValue(
               gaussPts,
               boost::shared_ptr<BaseFunctionCtx>(new EntPolynomialBaseCtx(
                   *dataOnElement[space], static_cast<FieldSpace>(space),
                   static_cast<FieldApproximationBase>(b), NOBASE)));
         }
       break;
     default:
       SETERRQ1(mField.get_comm(), MOFEM_DATA_INCONSISTENCY,
                "Base <%s> not yet implemented",
                ApproximationBaseNames[static_cast<FieldApproximationBase>(b)]);
     }
   }
   MoFEMFunctionReturn(0);
 }

◆ calculateBaseFunctionsOnElement() [2/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::calculateBaseFunctionsOnElement ( )

Calculate base functions.

Returns: Error code

Definition at line 891 of file ForcesAndSourcesCore.cpp.

                                                                      {
   MoFEMFunctionBegin;
   /// Use the some node base. Node base is usually used for construction other
   /// bases.
   for (int space = HCURL; space != LASTSPACE; ++space) {
     dataOnElement[space]->dataOnEntities[MBVERTEX][0].getNSharedPtr(NOBASE) =
         dataOnElement[H1]->dataOnEntities[MBVERTEX][0].getNSharedPtr(NOBASE);
   }
   for (int b = AINSWORTH_LEGENDRE_BASE; b != LASTBASE; b++) {
     CHKERR calculateBaseFunctionsOnElement(
         static_cast<FieldApproximationBase>(b));
   }
   MoFEMFunctionReturn(0);
 }

◆ createDataOnElement()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::createDataOnElement ( )

Create a entity data on element object.

Returns: MoFEMErrorCode

Definition at line 906 of file ForcesAndSourcesCore.cpp.

                                                          {
   MoFEMFunctionBegin;
 
   const EntityType type = numeredEntFiniteElementPtr->getEntType();
   if (type == lastEvaluatedElementEntityType)
     MoFEMFunctionReturnHot(0);
 
   // Data on elements for proper spaces
   for (int space = H1; space != LASTSPACE; ++space) {
     dataOnElement[space]->setElementType(type);
     derivedDataOnElement[space]->setElementType(type);
   }
 
   lastEvaluatedElementEntityType = type;
 
   MoFEMFunctionReturn(0);
 }

◆ getColNodesIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getColNodesIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get col node indices from FENumeredDofEntity_multiIndex

Definition at line 313 of file ForcesAndSourcesCore.cpp.

                                                                             {
   return getNodesIndices(field_name,
                          const_cast<FENumeredDofEntity_multiIndex &>(
                              numeredEntFiniteElementPtr->getColsDofs()),
                          data.dataOnEntities[MBVERTEX][0].getIndices(),
                          data.dataOnEntities[MBVERTEX][0].getLocalIndices());
 }

◆ getDataOrder()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getDataOrder	(	const EntityType	type,
		const FieldSpace	space,
		boost::ptr_vector< DataForcesAndSourcesCore::EntData > &	data
	)		const

get maximal approximation order of approximation on the entity

Definition at line 200 of file ForcesAndSourcesCore.cpp.

                                                                   {
   MoFEMFunctionBegin;
   SideNumber_multiIndex &side_table = const_cast<SideNumber_multiIndex &>(
       numeredEntFiniteElementPtr->getSideNumberTable());
   if (PetscUnlikely(data.size() < side_table.get<2>().count(type))) {
     // prims has 9 edges, some of edges for "flat" prism are not active
     SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "data structure too small to keep data %d < %d", data.size(),
              side_table.get<2>().count(type));
   }
   for (unsigned int side = 0; side != data.size(); ++side) {
     data[side].getDataOrder() = 0;
   }
   auto &data_dofs = const_cast<FEDofEntity_multiIndex::index<
       Composite_EntType_and_Space_mi_tag>::type &>(
       numeredEntFiniteElementPtr->getDataDofs()
           .get<Composite_EntType_and_Space_mi_tag>());
   auto tuple = boost::make_tuple(type, space);
   auto dit = data_dofs.lower_bound(tuple);
   if (dit != data_dofs.end()) {
     auto hi_dit = data_dofs.upper_bound(tuple);
     for (; dit != hi_dit; dit++) {
       auto &dof = **dit;
       if (dof.getEntDofIdx() == 0) {
         ApproximationOrder ent_order = dof.getMaxOrder();
         int side_number = dof.sideNumberPtr->side_number;
         if (PetscUnlikely(side_number < 0)) {
           SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                   "side number can not be negative");
         }
         auto &dat = data[side_number];
         dat.getDataOrder() =
             dat.getDataOrder() > ent_order ? dat.getDataOrder() : ent_order;
         if (dof.sideNumberPtr->brother_side_number != -1) {
           if (PetscUnlikely(
                   data.size() <
                   (unsigned int)(*dit)->sideNumberPtr->brother_side_number)) {
             SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                     "data size not big enough to keep brother data");
           }
           data[dof.sideNumberPtr->brother_side_number].getDataOrder() =
               dat.getDataOrder();
         }
       }
     }
   }
   MoFEMFunctionReturn(0);
 }

◆ getElementPolynomialBase()

auto& MoFEM::ForcesAndSourcesCore::getElementPolynomialBase ( )

Get the Entity Polynomial Base object.

Returns: boost::shared_ptr<BaseFunction>&&

Definition at line 626 of file ForcesAndSourcesCore.hpp.

626 { return elementPolynomialBasePtr; }

MoFEM::ForcesAndSourcesCore::elementPolynomialBasePtr

boost::shared_ptr< BaseFunction > elementPolynomialBasePtr

Pointer to entity polynomial base.

Definition: ForcesAndSourcesCore.hpp:710

◆ getEntityColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEntityColIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name,
		const EntityType	type_lo = `MBVERTEX`,
		const EntityType	type_hi = `MBPOLYHEDRON`
	)		const

Definition at line 146 of file ForcesAndSourcesCore.hpp.

                                                                      {
     return getEntityIndices(data, field_name,
                             const_cast<FENumeredDofEntity_multiIndex &>(
                                 numeredEntFiniteElementPtr->getColsDofs()),
                             type_lo, type_hi);
   }

◆ getEntityFieldData()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEntityFieldData	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name,
		const EntityType	type_lo = `MBVERTEX`,
		const EntityType	type_hi = `MBPOLYHEDRON`
	)		const

Definition at line 608 of file ForcesAndSourcesCore.cpp.

                                                               {
   MoFEMFunctionBegin;
 
   for (EntityType t = type_lo; t != type_hi; ++t) {
     for (auto &dat : data.dataOnEntities[t]) {
       dat.getDataOrder() = 0;
       dat.getBase() = NOBASE;
       dat.getSpace() = NOSPACE;
       dat.getFieldData().resize(0, false);
       dat.getFieldDofs().resize(0, false);
     }
   }
 
   auto &dofs = const_cast<FEDofEntity_multiIndex &>(
       numeredEntFiniteElementPtr->getDataDofs());
   auto &dofs_by_type = dofs.get<Composite_Name_And_Type_mi_tag>();
   auto dit = dofs_by_type.lower_bound(boost::make_tuple(field_name, type_lo));
   if (dit == dofs_by_type.end())
     MoFEMFunctionReturnHot(0);
   auto hi_dit =
       dofs_by_type.lower_bound(boost::make_tuple(field_name, type_hi));
   std::vector<boost::weak_ptr<FEDofEntity>> brother_dofs_vec;
   for (; dit != hi_dit;) {
 
     auto &dof = **dit;
     const int nb_dofs_on_ent = dof.getNbDofsOnEnt();
     if (nb_dofs_on_ent) {
 
       const EntityType type = dof.getEntType();
       const int side = dof.sideNumberPtr->side_number;
       auto &dat = data.dataOnEntities[type][side];
 
       auto &ent_field_dofs = dat.getFieldDofs();
       auto &ent_field_data = dat.getFieldData();
       const int brother_side = dof.sideNumberPtr->brother_side_number;
       if (brother_side != -1) 
         brother_dofs_vec.emplace_back(*dit);
 
       if (ent_field_data.empty()) {
 
         dat.getBase() = dof.getApproxBase();
         dat.getSpace() = dof.getSpace();
         const int ent_order = dof.getMaxOrder();
         dat.getDataOrder() =
             dat.getDataOrder() > ent_order ? dat.getDataOrder() : ent_order;
         const auto dof_ent_field_data = dof.getEntFieldData();
         ent_field_data.resize(nb_dofs_on_ent, false);
         noalias(ent_field_data) = dof.getEntFieldData();
         ent_field_dofs.resize(nb_dofs_on_ent, false);
         for (int ii = 0; ii != nb_dofs_on_ent; ++ii) {
           ent_field_dofs[ii] = *dit;
           ++dit;
         }
       }
     }
   }
 
   for(auto &dof_ptr : brother_dofs_vec) {
     const EntityType type = dof_ptr.lock()->getEntType();
     const int side = dof_ptr.lock()->sideNumberPtr->side_number;
     const int brother_side = dof_ptr.lock()->sideNumberPtr->brother_side_number;
     auto &dat = data.dataOnEntities[type][side];
     auto &dat_brother = data.dataOnEntities[type][brother_side];
     dat_brother.getBase() = dat.getBase();
     dat_brother.getSpace() = dat.getSpace();
     dat_brother.getDataOrder() = dat.getDataOrder();
     dat_brother.getFieldData() = dat.getFieldData();
     dat_brother.getFieldDofs() = dat.getFieldDofs();
   }
 
   MoFEMFunctionReturn(0);
 }

◆ getEntityFieldDataOrder()

template<EntityType type>

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEntityFieldDataOrder	(	DataForcesAndSourcesCore &	data,
		const FieldSpace	space
	)		const

Definition at line 108 of file ForcesAndSourcesCore.hpp.

                                                                   {
     return getDataOrder(type, space, data.dataOnEntities[type]);
   }

◆ getEntityIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEntityIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name,
		FENumeredDofEntity_multiIndex &	dofs,
		const EntityType	type_lo = `MBVERTEX`,
		const EntityType	type_hi = `MBPOLYHEDRON`
	)		const

Definition at line 322 of file ForcesAndSourcesCore.cpp.

                                     {
   MoFEMFunctionBegin;
 
   for (EntityType t = type_lo; t != type_hi; ++t) {
     for (auto &dat : data.dataOnEntities[t]) {
       dat.getIndices().resize(0, false);
       dat.getLocalIndices().resize(0, false);
     }
   }
 
   auto &dofs_by_type = dofs.get<Composite_Name_And_Type_mi_tag>();
   auto dit = dofs_by_type.lower_bound(boost::make_tuple(field_name, type_lo));
   if (dit == dofs_by_type.end())
     MoFEMFunctionReturnHot(0);
   auto hi_dit =
       dofs_by_type.lower_bound(boost::make_tuple(field_name, type_hi));
   for (; dit != hi_dit; ++dit) {
     auto &dof = **dit;
     const EntityType type = dof.getEntType();
     const int side = dof.sideNumberPtr->side_number;
     auto &dat = data.dataOnEntities[type][side];
 
     const int nb_dofs_on_ent = dof.getNbDofsOnEnt();
     if (nb_dofs_on_ent) {
       const int brother_side = dof.sideNumberPtr->brother_side_number;
       auto &ent_field_indices = dat.getIndices();
       auto &ent_field_local_indices = dat.getLocalIndices();
       if (ent_field_indices.empty()) {
         ent_field_indices.resize(nb_dofs_on_ent, false);
         ent_field_local_indices.resize(nb_dofs_on_ent, false);
         std::fill(ent_field_indices.data().begin(),
                   ent_field_indices.data().end(), -1);
         std::fill(ent_field_local_indices.data().begin(),
                   ent_field_local_indices.data().end(), -1);
       }
       const int idx = dof.getEntDofIdx();
       ent_field_indices[idx] = dof.getPetscGlobalDofIdx();
       ent_field_local_indices[idx] = dof.getPetscLocalDofIdx();
       if (brother_side != -1) {
         auto &dat_brother = data.dataOnEntities[type][brother_side];
         dat_brother.getIndices().resize(nb_dofs_on_ent, false);
         dat_brother.getLocalIndices().resize(nb_dofs_on_ent, false); 
         dat_brother.getIndices()[idx] = dat.getIndices()[idx];
         dat_brother.getLocalIndices()[idx] = dat.getLocalIndices()[idx];
       }
     }
   }
 
   MoFEMFunctionReturn(0);
 }

◆ getEntityRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEntityRowIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name,
		const EntityType	type_lo = `MBVERTEX`,
		const EntityType	type_hi = `MBPOLYHEDRON`
	)		const

Definition at line 135 of file ForcesAndSourcesCore.hpp.

                                                                      {
     return getEntityIndices(data, field_name,
                             const_cast<FENumeredDofEntity_multiIndex &>(
                                 numeredEntFiniteElementPtr->getRowsDofs()),
                             type_lo, type_hi);
   }

◆ getEntitySense()

template<EntityType type>

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getEntitySense ( DataForcesAndSourcesCore & data ) const

Definition at line 103 of file ForcesAndSourcesCore.hpp.

                                                                       {
     return getSense(type, data.dataOnEntities[type]);
   }

◆ getFaceTriNodes()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getFaceTriNodes ( DataForcesAndSourcesCore & data ) const

Get nodes on triangles.

Definition at line 718 of file ForcesAndSourcesCore.cpp.

                                                                           {
   MoFEMFunctionBegin;
   // PetscAttachDebugger();
   data.facesNodes.resize(4, 3, false);
   auto &side_table = const_cast<SideNumber_multiIndex &>(
       numeredEntFiniteElementPtr->getSideNumberTable());
   auto siit = side_table.get<1>().lower_bound(boost::make_tuple(MBTRI, 0));
   auto hi_siit = side_table.get<1>().upper_bound(boost::make_tuple(MBTRI, 4));
   if (std::distance(siit, hi_siit) != 4) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "Should be 4 triangles on tet, side_table not initialized");
   }
   const int canonical_face_sense_p1[4][3] = {
       {0, 1, 3},
       {1, 2, 3},
       {0, 3, 2} /**/,
       {0, 2, 1} /**/}; // second index is offset (positive sense)
   const int canonical_face_sense_m1[4][3] = {
       {0, 3, 1},
       {1, 3, 2},
       {0, 2, 3},
       {0, 1, 2}}; // second index is offset (negative sense
   for (; siit != hi_siit; siit++) {
     const boost::shared_ptr<SideNumber> side = *siit;
     int face_conn[3] = {-1, -1, -1};
     if (side->offset == 0) {
       face_conn[0] = side->sense == 1
                          ? canonical_face_sense_p1[(int)side->side_number][0]
                          : canonical_face_sense_m1[(int)side->side_number][0];
       face_conn[1] = side->sense == 1
                          ? canonical_face_sense_p1[(int)side->side_number][1]
                          : canonical_face_sense_m1[(int)side->side_number][1];
       face_conn[2] = side->sense == 1
                          ? canonical_face_sense_p1[(int)side->side_number][2]
                          : canonical_face_sense_m1[(int)side->side_number][2];
     }
     if (side->offset == 1) {
       face_conn[0] =
           side->sense == 1
               ? canonical_face_sense_p1[(int)side->side_number][1]
               : canonical_face_sense_m1[(int)side->side_number][2] /**/;
       face_conn[1] = side->sense == 1
                          ? canonical_face_sense_p1[(int)side->side_number][2]
                          : canonical_face_sense_m1[(int)side->side_number][0];
       face_conn[2] = side->sense == 1
                          ? canonical_face_sense_p1[(int)side->side_number][0]
                          : canonical_face_sense_m1[(int)side->side_number][1];
     }
     if (side->offset == 2) {
       face_conn[0] =
           side->sense == 1
               ? canonical_face_sense_p1[(int)side->side_number][2]
               : canonical_face_sense_m1[(int)side->side_number][1] /**/;
       face_conn[1] = side->sense == 1
                          ? canonical_face_sense_p1[(int)side->side_number][0]
                          : canonical_face_sense_m1[(int)side->side_number][2];
       face_conn[2] = side->sense == 1
                          ? canonical_face_sense_p1[(int)side->side_number][1]
                          : canonical_face_sense_m1[(int)side->side_number][0];
     }
     for (int nn = 0; nn < 3; nn++)
       data.facesNodes(side->side_number, nn) = face_conn[nn];
     {
       const EntityHandle *conn_tet;
       int num_nodes_tet;
       EntityHandle ent = numeredEntFiniteElementPtr->getEnt();
       CHKERR mField.get_moab().get_connectivity(ent, conn_tet, num_nodes_tet,
                                                 true);
       if (num_nodes_tet != 4)
         SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                 "data inconsistency");
       int num_nodes_face;
       const EntityHandle *conn_face;
       CHKERR mField.get_moab().get_connectivity(side->ent, conn_face,
                                                 num_nodes_face, true);
       if (num_nodes_face != 3)
         SETERRQ(PETSC_COMM_SELF, 1, "data inconsistency");
       if (conn_face[0] != conn_tet[data.facesNodes(side->side_number, 0)])
         SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                 "data inconsistency");
       if (conn_face[1] != conn_tet[data.facesNodes(side->side_number, 1)])
         SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                 "data inconsistency");
       if (conn_face[2] != conn_tet[data.facesNodes(side->side_number, 2)])
         SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                 "data inconsistency");
     }
   }
   MoFEMFunctionReturn(0);
 }

◆ getMaxColOrder()

int MoFEM::ForcesAndSourcesCore::getMaxColOrder ( ) const

Get max order of approximation for field in columns.

Definition at line 196 of file ForcesAndSourcesCore.cpp.

                                                {
   return getMaxOrder(numeredEntFiniteElementPtr->getColsDofs());
 }

◆ getMaxDataOrder()

int MoFEM::ForcesAndSourcesCore::getMaxDataOrder ( ) const

Get max order of approximation for data fields.

Method getMaxDataOrder () return maximal order on entities, for all data on the element. So for example if finite element is triangle, and triangle base function have order 4 and on edges base function have order 2, this function return 4.

If finite element has for example 2 or more approximated fields, for example Pressure (order 3) and displacement field (order 5), this function returns 5.

Definition at line 188 of file ForcesAndSourcesCore.cpp.

                                                 {
   return getMaxOrder(numeredEntFiniteElementPtr->getDataDofs());
 }

◆ getMaxRowOrder()

int MoFEM::ForcesAndSourcesCore::getMaxRowOrder ( ) const

Get max order of approximation for field in rows.

Definition at line 192 of file ForcesAndSourcesCore.cpp.

                                                {
   return getMaxOrder(numeredEntFiniteElementPtr->getRowsDofs());
 }

◆ getNodesFieldData() [1/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNodesFieldData	(	const boost::string_ref	field_name,
		FEDofEntity_multiIndex &	dofs,
		VectorDouble &	nodes_data,
		VectorDofs &	nodes_dofs,
		FieldSpace &	space,
		FieldApproximationBase &	base
	)		const

Get field data on nodes.

Parameters

field_name	Name of field
dofs	Dofs (element) multi index
nodes_data	Returned DOFs values
nodes_dofs	Vector of pointers to DOFs data structure
space	Get space on nodes (Only H! is valid)
base	Get base on nodes

Returns: Error code

Definition at line 528 of file ForcesAndSourcesCore.cpp.

                                         {
   MoFEMFunctionBegin;
   auto &dofs_by_name_and_type = dofs.get<Composite_Name_And_Type_mi_tag>();
   auto tuple = boost::make_tuple(field_name, MBVERTEX);
   auto dit = dofs_by_name_and_type.lower_bound(tuple);
   if(dit == dofs_by_name_and_type.end())
     SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "No nodal dofs on element");
   auto hi_dit = dofs.get<Composite_Name_And_Type_mi_tag>().upper_bound(tuple);
 
   if (dit != hi_dit) {
     auto &first_dof = **dit;
     space = first_dof.getSpace();
     base = first_dof.getApproxBase();
     int num_nodes;
     CHKERR getNumberOfNodes(num_nodes);
     const int nb_dof_idx = first_dof.getNbOfCoeffs();
     const int max_nb_dofs = nb_dof_idx * num_nodes;
     nodes_data.resize(max_nb_dofs, false);
     nodes_dofs.resize(max_nb_dofs, false);
     nodes_data.clear();
 
     std::vector<boost::weak_ptr<FEDofEntity>> brother_dofs_vec;
     for (; dit != hi_dit;) {
       const auto &dof_ptr = *dit;
       const auto &dof = *dof_ptr;
       const auto &sn = *dof.sideNumberPtr;
       const int side_number = sn.side_number;
       const int brother_side_number = sn.brother_side_number;
       if(brother_side_number != -1)
         brother_dofs_vec.emplace_back(dof_ptr);
 
       int pos = side_number * nb_dof_idx;
       auto ent_filed_data_vec = dof.getEntFieldData();
       for (int ii = 0; ii != nb_dof_idx; ++ii) {
         nodes_data[pos] = ent_filed_data_vec[ii];
         nodes_dofs[pos] = *dit;
         ++pos;
         ++dit;
       }
     }
 
     for(auto &dof_ptr : brother_dofs_vec) {
       const auto &dof = *dof_ptr.lock();
       const auto &sn = *dof.sideNumberPtr;
       const int side_number = sn.side_number;
       const int brother_side_number = sn.brother_side_number;
       int pos = side_number * nb_dof_idx;
       int brother_pos = brother_side_number * nb_dof_idx;
       for (int ii = 0; ii != nb_dof_idx; ++ii) {
         nodes_data[brother_pos] = nodes_data[pos];
         nodes_dofs[brother_pos] = nodes_dofs[pos];        
         ++pos;
         ++brother_pos;
       }
     }
 
   } else {
     nodes_data.resize(0, false);
     nodes_dofs.resize(0, false);
   }
 
   MoFEMFunctionReturn(0);
 }

◆ getNodesFieldData() [2/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNodesFieldData	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

Get data on nodes.

Parameters

data	Data structure
field_name	Field name

Returns: Error code

Definition at line 597 of file ForcesAndSourcesCore.cpp.

                                                                            {
   return getNodesFieldData(field_name,
                            const_cast<FEDofEntity_multiIndex &>(
                                numeredEntFiniteElementPtr->getDataDofs()),
                            data.dataOnEntities[MBVERTEX][0].getFieldData(),
                            data.dataOnEntities[MBVERTEX][0].getFieldDofs(),
                            data.dataOnEntities[MBVERTEX][0].getSpace(),
                            data.dataOnEntities[MBVERTEX][0].getBase());
 }

◆ getNodesIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNodesIndices	(	const boost::string_ref	field_name,
		FENumeredDofEntity_multiIndex &	dofs,
		VectorInt &	nodes_indices,
		VectorInt &	local_nodes_indices
	)		const

get node indices

Definition at line 253 of file ForcesAndSourcesCore.cpp.

                                                                     {
   MoFEMFunctionBegin;
   auto &dofs_by_type = dofs.get<Composite_Name_And_Type_mi_tag>();
   auto tuple = boost::make_tuple(field_name, MBVERTEX);
   auto dit = dofs_by_type.lower_bound(tuple);
   auto hi_dit = dofs_by_type.upper_bound(tuple);
 
   if (dit != hi_dit) {
 
     int num_nodes;
     CHKERR getNumberOfNodes(num_nodes);
     int max_nb_dofs = 0;
     const int nb_dofs_on_vert = (*dit)->getNbOfCoeffs();
     max_nb_dofs = nb_dofs_on_vert * num_nodes;
     nodes_indices.resize(max_nb_dofs, false);
     local_nodes_indices.resize(max_nb_dofs, false);
     if (std::distance(dit, hi_dit) != max_nb_dofs) {
       std::fill(nodes_indices.begin(), nodes_indices.end(), -1);
       std::fill(local_nodes_indices.begin(), local_nodes_indices.end(), -1);
     }
 
     for (; dit != hi_dit; dit++) {
       auto &dof = **dit;
       const int idx = dof.getPetscGlobalDofIdx();
       const int local_idx = dof.getPetscLocalDofIdx();
       const int side_number = dof.sideNumberPtr->side_number;
       const int pos = side_number * nb_dofs_on_vert + dof.getDofCoeffIdx();
       nodes_indices[pos] = idx;
       local_nodes_indices[pos] = local_idx;
       const int brother_side_number =
           (*dit)->sideNumberPtr->brother_side_number;
       if (brother_side_number != -1) {
         const int elem_idx =
             brother_side_number * nb_dofs_on_vert + (*dit)->getDofCoeffIdx();
         nodes_indices[elem_idx] = idx;
         local_nodes_indices[elem_idx] = local_idx;
       }
     }
 
   } else {
     nodes_indices.resize(0, false);
     local_nodes_indices.resize(0, false);
   }
 
   MoFEMFunctionReturn(0);
 }

◆ getNoFieldColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldColIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get col NoField indices

Definition at line 404 of file ForcesAndSourcesCore.cpp.

                                                                        {
   MoFEMFunctionBegin;
   if (data.dataOnEntities[MBENTITYSET].size() == 0) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
   }
   CHKERR getNoFieldIndices(field_name,
                            const_cast<FENumeredDofEntity_multiIndex &>(
                                numeredEntFiniteElementPtr->getColsDofs()),
                            data.dataOnEntities[MBENTITYSET][0].getIndices());
   MoFEMFunctionReturn(0);
 }

◆ getNoFieldFieldData() [1/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldFieldData	(	const boost::string_ref	field_name,
		FEDofEntity_multiIndex &	dofs,
		VectorDouble &	ent_field_data,
		VectorDofs &	ent_field_dofs
	)		const

Definition at line 683 of file ForcesAndSourcesCore.cpp.

                                                                     {
   MoFEMFunctionBeginHot;
   auto dit = dofs.get<FieldName_mi_tag>().lower_bound(field_name);
   auto hi_dit = dofs.get<FieldName_mi_tag>().upper_bound(field_name);
   int size = std::distance(dit, hi_dit);
   ent_field_data.resize(size, false);
   ent_field_dofs.resize(size, false);
   for (; dit != hi_dit; dit++) {
     int idx = (*dit)->getDofCoeffIdx();
     ent_field_data[idx] = (*dit)->getFieldData();
     ent_field_dofs[idx] = *dit;
   }
   MoFEMFunctionReturnHot(0);
 }

◆ getNoFieldFieldData() [2/2]

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldFieldData	(	DataForcesAndSourcesCore &	data,
		const boost::string_ref	field_name
	)		const

Definition at line 700 of file ForcesAndSourcesCore.cpp.

                                                                             {
   MoFEMFunctionBegin;
   if (data.dataOnEntities[MBENTITYSET].size() == 0) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
   }
   CHKERR getNoFieldFieldData(
       field_name,
       const_cast<FEDofEntity_multiIndex &>(
           numeredEntFiniteElementPtr->getDataDofs()),
       data.dataOnEntities[MBENTITYSET][0].getFieldData(),
       data.dataOnEntities[MBENTITYSET][0].getFieldDofs());
   MoFEMFunctionReturn(0);
 }

◆ getNoFieldIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldIndices	(	const std::string &	field_name,
		FENumeredDofEntity_multiIndex &	dofs,
		VectorInt &	nodes_indices
	)		const

get NoField indices

Definition at line 377 of file ForcesAndSourcesCore.cpp.

                                                                   {
   MoFEMFunctionBeginHot;
   auto dit = dofs.get<FieldName_mi_tag>().lower_bound(field_name);
   auto hi_dit = dofs.get<FieldName_mi_tag>().upper_bound(field_name);
   indices.resize(std::distance(dit, hi_dit));
   for (; dit != hi_dit; dit++) {
     int idx = (*dit)->getPetscGlobalDofIdx();
     indices[(*dit)->getDofCoeffIdx()] = idx;
   }
   MoFEMFunctionReturnHot(0);
 }

◆ getNoFieldRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNoFieldRowIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get col NoField indices

Definition at line 391 of file ForcesAndSourcesCore.cpp.

                                                                        {
   MoFEMFunctionBegin;
   if (data.dataOnEntities[MBENTITYSET].size() == 0) {
     SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "data inconsistency");
   }
   CHKERR getNoFieldIndices(field_name,
                            const_cast<FENumeredDofEntity_multiIndex &>(
                                numeredEntFiniteElementPtr->getRowsDofs()),
                            data.dataOnEntities[MBENTITYSET][0].getIndices());
   MoFEMFunctionReturn(0);
 }

◆ getNumberOfNodes()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getNumberOfNodes ( int & num_nodes ) const

Definition at line 111 of file ForcesAndSourcesCore.cpp.

                                                                           {
   MoFEMFunctionBeginHot;
 
   EntityHandle ent = numeredEntFiniteElementPtr->getEnt();
   switch (mField.get_moab().type_from_handle(ent)) {
   case MBVERTEX:
     num_nodes = 1;
     break;
   case MBEDGE:
     num_nodes = 2;
     break;
   case MBTRI:
     num_nodes = 3;
     break;
   case MBQUAD:
     num_nodes = 4;
     break;
   case MBTET:
     num_nodes = 4;
     break;
   case MBPRISM:
     num_nodes = 6;
     break;
   default:
     SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "not implemented");
   }
 
   MoFEMFunctionReturnHot(0);
 }

◆ getOpPtrVector()

boost::ptr_vector<UserDataOperator>& MoFEM::ForcesAndSourcesCore::getOpPtrVector ( )

Use to push back operator for row operator.

It can be used to calculate nodal forces or other quantities on the mesh.

Examples:: elasticity_mixed_formulation.cpp, forces_and_sources_testing_edge_element.cpp, and MagneticElement.hpp.

Definition at line 619 of file ForcesAndSourcesCore.hpp.

619 { return opPtrVector; }

MoFEM::ForcesAndSourcesCore::opPtrVector

boost::ptr_vector< UserDataOperator > opPtrVector

Vector of finite element users data operators.

Definition: ForcesAndSourcesCore.hpp:612

◆ getProblemNodesColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemNodesColIndices	(	const std::string &	field_name,
		VectorInt &	nodes_indices
	)		const

Definition at line 512 of file ForcesAndSourcesCore.cpp.

                                                                  {
   return getProblemNodesIndices(field_name, *(problemPtr->numeredDofsCols),
                                 nodes_indices);
 }

◆ getProblemNodesIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemNodesIndices	(	const std::string &	field_name,
		const NumeredDofEntity_multiIndex &	dofs,
		VectorInt &	nodes_indices
	)		const

get indices of nodal indices which are declared for problem but not this particular element

Definition at line 419 of file ForcesAndSourcesCore.cpp.

                                     {
   MoFEMFunctionBeginHot;
   nodes_indices.resize(0);
   auto &side_table = const_cast<SideNumber_multiIndex &>(
       numeredEntFiniteElementPtr->getSideNumberTable());
   auto siit = side_table.get<1>().lower_bound(boost::make_tuple(MBVERTEX, 0));
   auto hi_siit =
       side_table.get<1>().lower_bound(boost::make_tuple(MBVERTEX, 10000));
 
   int nn = 0;
   for (; siit != hi_siit; siit++, nn++) {
 
     if (siit->get()->side_number == -1)
       continue;
 
     const EntityHandle ent = siit->get()->ent;
     auto dit =
         dofs.get<Composite_Name_And_Ent_And_EntDofIdx_mi_tag>().lower_bound(
             boost::make_tuple(field_name, ent, 0));
     auto hi_dit =
         dofs.get<Composite_Name_And_Ent_And_EntDofIdx_mi_tag>().upper_bound(
             boost::make_tuple(field_name, ent, 10000)); /// very large number
     if (dit != hi_dit) {
       if (!nn) {
         nodes_indices.resize((*dit)->getNbOfCoeffs() *
                              std::distance(siit, hi_siit));
       }
       for (; dit != hi_dit; dit++) {
         nodes_indices[siit->get()->side_number * (*dit)->getNbOfCoeffs() +
                       (*dit)->getDofCoeffIdx()] =
             (*dit)->getPetscGlobalDofIdx();
       }
     }
   }
 
   MoFEMFunctionReturnHot(0);
 }

◆ getProblemNodesRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemNodesRowIndices	(	const std::string &	field_name,
		VectorInt &	nodes_indices
	)		const

Definition at line 498 of file ForcesAndSourcesCore.cpp.

                                                                  {
   return getProblemNodesIndices(field_name, *(problemPtr->numeredDofsRows),
                                 nodes_indices);
 }

◆ getProblemTypeColIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemTypeColIndices	(	const std::string &	field_name,
		EntityType	type,
		int	side_number,
		VectorInt &	indices
	)		const

Definition at line 519 of file ForcesAndSourcesCore.cpp.

                                                                          {
   return getProblemTypeIndices(field_name, *(problemPtr->numeredDofsCols), type,
                                side_number, indices);
 }

◆ getProblemTypeIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemTypeIndices	(	const std::string &	field_name,
		const NumeredDofEntity_multiIndex &	dofs,
		EntityType	type,
		int	side_number,
		VectorInt &	indices
	)		const

get indices by type (generic function) which are declared for problem but not this particular element

Definition at line 459 of file ForcesAndSourcesCore.cpp.

                                                                 {
   MoFEMFunctionBeginHot;
 
   indices.resize(0);
 
   auto &side_table = const_cast<SideNumber_multiIndex &>(
       numeredEntFiniteElementPtr->getSideNumberTable());
   auto siit =
       side_table.get<1>().lower_bound(boost::make_tuple(type, side_number));
   auto hi_siit =
       side_table.get<1>().upper_bound(boost::make_tuple(type, side_number));
 
   for (; siit != hi_siit; siit++) {
 
     if (siit->get()->side_number == -1)
       continue;
 
     const EntityHandle ent = siit->get()->ent;
     NumeredDofEntity_multiIndex::index<
         Composite_Name_And_Ent_And_EntDofIdx_mi_tag>::type::iterator dit,
         hi_dit;
     dit = dofs.get<Composite_Name_And_Ent_And_EntDofIdx_mi_tag>().lower_bound(
         boost::make_tuple(field_name, ent, 0));
     hi_dit =
         dofs.get<Composite_Name_And_Ent_And_EntDofIdx_mi_tag>().upper_bound(
             boost::make_tuple(field_name, ent, 10000)); /// very large number
 
     indices.resize(std::distance(dit, hi_dit));
     for (; dit != hi_dit; dit++) {
 
       indices[(*dit)->getEntDofIdx()] = (*dit)->getPetscGlobalDofIdx();
     }
   }
 
   MoFEMFunctionReturnHot(0);
 }

◆ getProblemTypeRowIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getProblemTypeRowIndices	(	const std::string &	field_name,
		EntityType	type,
		int	side_number,
		VectorInt &	indices
	)		const

Definition at line 505 of file ForcesAndSourcesCore.cpp.

                                                                          {
   return getProblemTypeIndices(field_name, *(problemPtr->numeredDofsRows), type,
                                side_number, indices);
 }

◆ getRowNodesIndices()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getRowNodesIndices	(	DataForcesAndSourcesCore &	data,
		const std::string &	field_name
	)		const

get row node indices from FENumeredDofEntity_multiIndex

Definition at line 303 of file ForcesAndSourcesCore.cpp.

                                                                             {
   return getNodesIndices(field_name,
                          const_cast<FENumeredDofEntity_multiIndex &>(
                              numeredEntFiniteElementPtr->getRowsDofs()),
                          data.dataOnEntities[MBVERTEX][0].getIndices(),
                          data.dataOnEntities[MBVERTEX][0].getLocalIndices());
 }

◆ getRule() [1/2]

virtual int MoFEM::ForcesAndSourcesCore::getRule ( int order )

virtual

set integration rule for finite element

This function is overloaded by the user. The integration rule is set such that specific operator implemented by the user is integrated accurately. For example if user implement bilinear operator

\[ b(u,v) = \int_\mathcal{T} \frac{\partial u_i}{\partial x_j}\frac{\partial v_i}{\partial x_j} \textrm{d}\mathcal{T} \]

then if \(u\) and \(v\) are polynomial of given order, then exact integral would be

int getRule(int order) { return 2*(order-1); };

The integration points and weights are set appropriately for given entity type and integration rule from quad.c

Method ForcesAndSourcesCore::getRule takes at argument takes maximal polynomial order set on the element on all fields defined on the element. If a user likes to have more control, another variant of this function can be called which distinguishing between field orders on rows, columns and data, the i.e. first argument of a bilinear form, the second argument of bilinear form and field coefficients on the element.

Note: If user set rule to -1 or any other negative integer, then method ForcesAndSourcesCore::setGaussPts is called. In that method user can implement own (specific) integration method.

Bug:: this function should be const

Reimplemented in EdgeSlidingConstrains::MyEdgeFE, PostProcFaceOnRefinedMesh, BoneRemodeling::DensityMapFe, MoFEM::VolumeElementForcesAndSourcesCoreOnSide, GelModule::Gel::GelFE, SurfaceSlidingConstrains::MyTriangleFE, SmallStrainPlasticity::MyVolumeFE, PostProcTemplateVolumeOnRefinedMesh< MoFEM::VolumeElementForcesAndSourcesCore >, KelvinVoigtDamper::DamperFE, NitscheMethod::MyVolumeFE, ThermalElement::MyTriFE, NeummanForcesSurface::MyTriangleFE, MixTransport::MixTransportElement::MyTriFE, NonlinearElasticElement::MyVolumeFE, MagneticElement::TriFE, MoFEM::NormElement::MyTriFE, NeummanForcesSurfaceComplexForLazy::MyTriangleSpatialFE, ThermalElement::MyVolumeFE, MixTransport::MixTransportElement::MyVolumeFE, MoFEM::NormElement::MyVolumeFE, OptimalMassTransport::TriFE, MagneticElement::VolumeFE, CellEngineering::FaceElement, BoneRemodeling::Remodeling::Fe, AnalyticalDirichletBC::ApproxField::MyTriFE, AnalyticalDirichletHelmholtzBC::ApproxField::MyTriFE, AnalyticalDirihletBC::ApproxField::MyTriFE, FluidPressure::MyTriangleFE, NitscheMethod::MyFace, OptimalMassTransport::VolumeFE, AnalyticalDirihletBC::ApproxField::MyTriFE, CohesiveElement::CohesiveInterfaceElement::MyPrism, EdgeForce::MyFE, BodyForceConstantField::MyVolumeFE, ThermalStressElement::MyVolumeFE, TestFE, and CellEngineering::DispMapFe.

Definition at line 275 of file ForcesAndSourcesCore.hpp.

275 { return 2 * order; }

◆ getRule() [2/2]

virtual int MoFEM::ForcesAndSourcesCore::getRule	(	int	order_row,
		int	order_col,
		int	order_data
	)

virtual

another variant of getRule

Parameters

order_row	order of base function on row
order_col	order of base function on columns
order_data	order of base function approximating data

Returns: integration rule

Bug:: this function should be const

Definition at line 298 of file ForcesAndSourcesCore.hpp.

                                                                     {
     return getRuleHook ? getRuleHook(order_row, order_col, order_data)
                        : getRule(order_data);
   }

◆ getSense()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getSense	(	EntityType	type,
		boost::ptr_vector< DataForcesAndSourcesCore::EntData > &	data
	)		const

get sense (orientation) of entity

Parameters

type	type of entity
data	entity data

Returns: error code

Definition at line 143 of file ForcesAndSourcesCore.cpp.

                                                                   {
   MoFEMFunctionBegin;
   auto &side_table = const_cast<SideNumber_multiIndex &>(
       numeredEntFiniteElementPtr->getSideNumberTable());
   if (PetscUnlikely(data.size() < side_table.get<2>().count(type))) {
     // prims has 9 edges, some of edges for "flat" prism are not active
     SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "wrong number of sides %u < %u", data.size(),
              side_table.get<2>().count(type));
   }
   const auto &st = side_table.get<2>();
   auto sit = st.lower_bound(type);
   auto hi_sit = st.upper_bound(type);
   for (; sit != hi_sit; sit++) {
     const auto &side = **sit;
     data[side.side_number].getSense() = sit->get()->sense;
     if (side.brother_side_number != -1) {
       if (PetscUnlikely(data.size() < (unsigned)side.brother_side_number)) {
         SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                 "data struture too small to keep data about brother node");
       }
       data[side.brother_side_number].getSense() = side.sense;
     }
   }
   MoFEMFunctionReturn(0);
 }

◆ getSpacesAndBaseOnEntities()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::getSpacesAndBaseOnEntities ( DataForcesAndSourcesCore & data ) const

Get field approximation space and base on entities.

Definition at line 811 of file ForcesAndSourcesCore.cpp.

                                           {
   MoFEMFunctionBeginHot;
   if (nInTheLoop == 0) {
     data.sPace.reset();
     data.bAse.reset();
     for (EntityType t = MBVERTEX; t != MBMAXTYPE; ++t) {
       data.spacesOnEntities[t].reset();
       data.basesOnEntities[t].reset();
     }
     for (int s = 0; s != LASTSPACE; ++s) {
       data.basesOnSpaces[s].reset();
     }
   }
   for (_IT_GET_FEDATA_DOFS_FOR_LOOP_(this, dof)) {
     if (dof->get()->getEntDofIdx() != 0)
       continue;
     const EntityType type = dof->get()->getEntType();
     const FieldSpace space = dof->get()->getSpace();
     const FieldApproximationBase approx = dof->get()->getApproxBase();
     data.sPace.set(space);
     data.bAse.set(approx);
     data.spacesOnEntities[type].set(space);
     data.basesOnEntities[type].set(approx);
     data.basesOnSpaces[space].set(approx);
   }
   MoFEMFunctionReturnHot(0);
 }

◆ getUserPolynomialBase()

auto& MoFEM::ForcesAndSourcesCore::getUserPolynomialBase ( )

Get the User Polynomial Base object.

Returns: boost::shared_ptr<BaseFunction>&

Definition at line 633 of file ForcesAndSourcesCore.hpp.

633 { return userPolynomialBasePtr; }

MoFEM::ForcesAndSourcesCore::userPolynomialBasePtr

boost::shared_ptr< BaseFunction > userPolynomialBasePtr

Pointer to user polynomail base.

Definition: ForcesAndSourcesCore.hpp:715

◆ loopOverOperators()

MoFEMErrorCode MoFEM::ForcesAndSourcesCore::loopOverOperators ( )

Iterate user data operators.

Returns: MoFEMErrorCode

Definition at line 954 of file ForcesAndSourcesCore.cpp.

                                                        {
   MoFEMFunctionBegin;
 
   const EntityType type = numeredEntFiniteElementPtr->getEntType();
   const UserDataOperator::OpType types[2] = {UserDataOperator::OPROW,
                                              UserDataOperator::OPCOL};
   std::vector<std::string> last_eval_field_name(2);
 
   boost::ptr_vector<UserDataOperator>::iterator oit, hi_oit;
   oit = opPtrVector.begin();
   hi_oit = opPtrVector.end();
 
   for (; oit != hi_oit; oit++) {
 
     oit->setPtrFE(this);
 
     if (oit->opType == UserDataOperator::OPLAST) {
 
       // Set field
       switch (oit->sPace) {
       case NOSPACE:
         SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "Unknown space");
       case NOFIELD:
       case H1:
       case HCURL:
       case HDIV:
       case L2:
         break;
       default:
         SETERRQ1(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
                  "Not implemented for this space", oit->sPace);
       }
 
       // Reseat all data which all field dependent
       dataOnElement[oit->sPace]->resetFieldDependentData();
       last_eval_field_name[0] = "";
 
       // Run operator
       try {
         CHKERR oit->opRhs(*dataOnElement[oit->sPace], oit->doVertices,
                           oit->doEdges, oit->doQuads, oit->doTris, oit->doTets,
                           false);
       }
       CATCH_OP_ERRORS(*oit);
 
     } else {
 
       boost::shared_ptr<DataForcesAndSourcesCore> op_data[2];
 
       for (int ss = 0; ss != 2; ss++) {
 
         const std::string field_name =
             !ss ? oit->rowFieldName : oit->colFieldName;
         const Field *field_struture = mField.get_field_structure(field_name);
         const BitFieldId data_id = field_struture->getId();
         const FieldSpace space = field_struture->getSpace();
         op_data[ss] = !ss ? dataOnElement[space] : derivedDataOnElement[space];
 
         if ((oit->getNumeredEntFiniteElementPtr()->getBitFieldIdData() &
              data_id)
                 .none()) {
           SETERRQ2(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                    "no data field < %s > on finite element < %s >",
                    field_name.c_str(), feName.c_str());
         }
 
         if (oit->getOpType() & types[ss] ||
             oit->getOpType() & UserDataOperator::OPROWCOL) {
 
           switch (space) {
           case NOSPACE:
             SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                     "unknown space");
             break;
           case NOFIELD:
           case H1:
           case HCURL:
           case HDIV:
           case L2:
             break;
           default:
             SETERRQ1(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
                      "Not implemented for this space", space);
           }
 
           if (last_eval_field_name[ss] != field_name) {
 
             CHKERR getEntityFieldData(*op_data[ss], field_name, MBEDGE);
             if (!ss)
               CHKERR getEntityRowIndices(*op_data[ss], field_name, MBEDGE);
             else
               CHKERR getEntityColIndices(*op_data[ss], field_name, MBEDGE);
 
             switch (space) {
             case NOSPACE:
               SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                       "unknown space");
               break;
             case H1:
               if (!ss)
                 CHKERR getRowNodesIndices(*op_data[ss], field_name);
               else
                 CHKERR getColNodesIndices(*op_data[ss], field_name);
               CHKERR getNodesFieldData(*op_data[ss], field_name);
               break;
             case HCURL:
             case HDIV:
               break;
             case L2:
               switch (type) {
               case MBVERTEX:
                 CHKERR getNodesFieldData(*op_data[ss], field_name);
                 break;
               default:
                 break;
               }
               break;
             case NOFIELD:
               if (!getNinTheLoop()) {
                 // NOFIELD data are the same for each element, can be
                 // retrieved only once
                 if (!ss) {
                   CHKERR getNoFieldRowIndices(*op_data[ss], field_name);
                 } else {
                   CHKERR getNoFieldColIndices(*op_data[ss], field_name);
                 }
                 CHKERR getNoFieldFieldData(*op_data[ss], field_name);
               }
               break;
             default:
               SETERRQ1(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED,
                        "Not implemented for this space", space);
             }
             last_eval_field_name[ss] = field_name;
 
           }
 
         }
       }
 
       if (oit->getOpType() & UserDataOperator::OPROW) {
         try {
           CHKERR oit->opRhs(*op_data[0], false);
         }
         CATCH_OP_ERRORS(*oit);
       }
 
       if (oit->getOpType() & UserDataOperator::OPCOL) {
         try {
           CHKERR oit->opRhs(*op_data[1], false);
         }
         CATCH_OP_ERRORS(*oit);
       }
 
       if (oit->getOpType() & UserDataOperator::OPROWCOL) {
         try {
           CHKERR oit->opLhs(*op_data[0], *op_data[1], oit->sYmm);
         }
         CATCH_OP_ERRORS(*oit);
       }
     }
   }
   MoFEMFunctionReturn(0);
 }

◆ operator()()

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::operator() ( )

virtual

function is run for every finite element

It is used to calculate element local matrices and assembly. It can be used for post-processing.

Reimplemented from MoFEM::BasicMethod.

Reimplemented in MoFEM::FaceElementForcesAndSourcesCore, MoFEM::VolumeElementForcesAndSourcesCore, MoFEM::FatPrismElementForcesAndSourcesCore, MoFEM::FlatPrismElementForcesAndSourcesCore, MoFEM::EdgeElementForcesAndSourcesCore, NeummanForcesSurfaceComplexForLazy::MyTriangleSpatialFE, and MoFEM::VertexElementForcesAndSourcesCore.

Definition at line 655 of file ForcesAndSourcesCore.hpp.

                                       {
     MoFEMFunctionBeginHot;
     if (operatorHook) {
       ierr = operatorHook();
       CHKERRG(ierr);
     }
     MoFEMFunctionReturnHot(0);
   }

◆ postProcess()

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::postProcess ( )

virtual

function is run at the end of loop

It is used to assembly matrices and vectors, calculating global variables, f.e. total internal energy, ect.

Iterating over dofs: Example1 iterating over dofs in row by name of the field for(IT_GET_FEROW_BY_NAME_DOFS_FOR_LOOP(this,"DISPLACEMENT",it)) { ... }

Reimplemented from MoFEM::BasicMethod.

Reimplemented in PostProcFaceOnRefinedMesh, PostProcFatPrismOnRefinedMesh, SolidShellModule::SolidShellPrismElement::PostProcFatPrismOnTriangleOnRefinedMesh, PostProcTemplateVolumeOnRefinedMesh< MoFEM::VolumeElementForcesAndSourcesCore >, GelModule::Gel::GelFE, MoFEM::FatPrismElementForcesAndSourcesCore, KelvinVoigtDamper::DamperFE, SolidShellModule::SolidShellPrismElement::SolidShellError, SolidShellModule::SolidShellPrismElement::SolidShell, Smoother::MyVolumeFE, and NonlinearElasticElement::MyVolumeFE.

Definition at line 663 of file ForcesAndSourcesCore.hpp.

                                        {
     MoFEMFunctionBeginHot;
     if (postProcessHook) {
       ierr = postProcessHook();
       CHKERRG(ierr);
     }
     MoFEMFunctionReturnHot(0);
   }

◆ preProcess()

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::preProcess ( )

virtual

function is run at the beginning of loop

It is used to zeroing matrices and vectors, calculation of shape functions on reference element, preprocessing boundary conditions, etc.

Reimplemented from MoFEM::BasicMethod.

Reimplemented in EdgeSlidingConstrains::MyEdgeFE, PostProcFaceOnRefinedMesh, PostProcFatPrismOnRefinedMesh, PostProcTemplateVolumeOnRefinedMesh< MoFEM::VolumeElementForcesAndSourcesCore >, GelModule::Gel::GelFE, SurfaceSlidingConstrains::MyTriangleFE, MoFEM::FatPrismElementForcesAndSourcesCore, KelvinVoigtDamper::DamperFE, SolidShellModule::SolidShellPrismElement::SolidShellError, SolidShellModule::SolidShellPrismElement::SolidShell, NeummanForcesSurfaceComplexForLazy::MyTriangleSpatialFE, Smoother::MyVolumeFE, NonlinearElasticElement::MyVolumeFE, and FluidPressure::MyTriangleFE.

Definition at line 647 of file ForcesAndSourcesCore.hpp.

                                       {
     MoFEMFunctionBeginHot;
     if (preProcessHook) {
       ierr = preProcessHook();
       CHKERRG(ierr);
     }
     MoFEMFunctionReturnHot(0);
   }

◆ setGaussPts() [1/2]

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::setGaussPts ( int order )

virtual

It will be removed in the future use other variant.

Reimplemented in PostProcFaceOnRefinedMesh, PostProcTemplateVolumeOnRefinedMesh< MoFEM::VolumeElementForcesAndSourcesCore >, MoFEM::VolumeElementForcesAndSourcesCoreOnSide, and NitscheMethod::MyVolumeFE.

Definition at line 305 of file ForcesAndSourcesCore.hpp.

                                                 {
     MoFEMFunctionBeginHot;
     SETERRQ(PETSC_COMM_SELF, MOFEM_NOT_IMPLEMENTED, "sorry, not implemented");
     MoFEMFunctionReturnHot(0);
   }

◆ setGaussPts() [2/2]

virtual MoFEMErrorCode MoFEM::ForcesAndSourcesCore::setGaussPts	(	int	order_row,
		int	order_col,
		int	order_data
	)

virtual

set user specific integration rule

This function allows for user defined integration rule. The key is to called matrix gaussPts, which is used by other MoFEM procedures. Matrix has number of rows equal to problem dimension plus one, where last index is used to store weight values. Number of columns is equal to number of integration points.

Note: This function is called if method ForcesAndSourcesCore::getRule is returning integer -1 or any other negative integer.

User sets

MatrixDouble gaussPts;

where

gaussPts.resize(dim+1,nb_gauss_pts);

number rows represents local coordinates of integration points in reference element, where last index in row is for integration weight.

Definition at line 334 of file ForcesAndSourcesCore.hpp.

                                                      {
 
     MoFEMFunctionBeginHot;
     ierr = setGaussPts(order_data);
     CHKERRG(ierr);
     MoFEMFunctionReturnHot(0);
   }

Member Data Documentation

◆ dataH1

DataForcesAndSourcesCore& MoFEM::ForcesAndSourcesCore::dataH1

Definition at line 57 of file ForcesAndSourcesCore.hpp.

◆ dataHcurl

DataForcesAndSourcesCore& MoFEM::ForcesAndSourcesCore::dataHcurl

Definition at line 58 of file ForcesAndSourcesCore.hpp.

◆ dataHdiv

DataForcesAndSourcesCore& MoFEM::ForcesAndSourcesCore::dataHdiv

Definition at line 59 of file ForcesAndSourcesCore.hpp.

◆ dataL2

DataForcesAndSourcesCore& MoFEM::ForcesAndSourcesCore::dataL2

Definition at line 60 of file ForcesAndSourcesCore.hpp.

◆ dataNoField

DataForcesAndSourcesCore& MoFEM::ForcesAndSourcesCore::dataNoField

Definition at line 56 of file ForcesAndSourcesCore.hpp.

◆ dataOnElement

const boost::shared_ptr<DataForcesAndSourcesCore> MoFEM::ForcesAndSourcesCore::dataOnElement[LASTSPACE]

Entity data on element entity rows fields.

FIXME: that should be moved to private class data and acessed only by member function

Definition at line 45 of file ForcesAndSourcesCore.hpp.

◆ derivedDataOnElement

const boost::shared_ptr<DataForcesAndSourcesCore> MoFEM::ForcesAndSourcesCore::derivedDataOnElement[LASTSPACE]

Entity data on element entity columns fields.

FIXME: that should be moved to private class data and acessed only by member function

Definition at line 54 of file ForcesAndSourcesCore.hpp.

◆ elementPolynomialBasePtr

boost::shared_ptr<BaseFunction> MoFEM::ForcesAndSourcesCore::elementPolynomialBasePtr

private

Pointer to entity polynomial base.

Definition at line 710 of file ForcesAndSourcesCore.hpp.

◆ gaussPts

MatrixDouble MoFEM::ForcesAndSourcesCore::gaussPts

Matrix of integration points.

Columns is equal to number of integration points, numver of rows depends on dimension of finite element entity, for example for tetrahedron rows are x,y,z,weight. Last row is integration weight.

FIXME: that should be moved to private class data and acessed only by member function

Definition at line 645 of file ForcesAndSourcesCore.hpp.

◆ getRuleHook

RuleHookFun MoFEM::ForcesAndSourcesCore::getRuleHook

Hook to get rule.

Todo:: check preferred format how works with gcc and clang, see http://www.boost.org/doc/libs/1_64_0/doc/html/function/tutorial.html#idp247873024

Definition at line 287 of file ForcesAndSourcesCore.hpp.

◆ lastEvaluatedElementEntityType

EntityType MoFEM::ForcesAndSourcesCore::lastEvaluatedElementEntityType

private

Last evaluated type of element entity.

Definition at line 704 of file ForcesAndSourcesCore.hpp.

◆ mField

Interface& MoFEM::ForcesAndSourcesCore::mField

Definition at line 36 of file ForcesAndSourcesCore.hpp.

◆ opPtrVector

boost::ptr_vector<UserDataOperator> MoFEM::ForcesAndSourcesCore::opPtrVector

Vector of finite element users data operators.

FIXME: that should be moved to private class data and acessed only by member function

Definition at line 612 of file ForcesAndSourcesCore.hpp.

◆ userPolynomialBasePtr

boost::shared_ptr<BaseFunction> MoFEM::ForcesAndSourcesCore::userPolynomialBasePtr

private

Pointer to user polynomail base.

Definition at line 715 of file ForcesAndSourcesCore.hpp.

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

src/finite_elements/ForcesAndSourcesCore.hpp
src/finite_elements/impl/ForcesAndSourcesCore.cpp

Classes

Public Types

Public Member Functions

Public Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ RuleHookFun

Constructor & Destructor Documentation

◆ ForcesAndSourcesCore()

◆ ~ForcesAndSourcesCore()

Member Function Documentation

◆ calculateBaseFunctionsOnElement() [1/2]

◆ calculateBaseFunctionsOnElement() [2/2]

◆ createDataOnElement()

◆ getColNodesIndices()

◆ getDataOrder()

◆ getElementPolynomialBase()

◆ getEntityColIndices()

◆ getEntityFieldData()

◆ getEntityFieldDataOrder()

◆ getEntityIndices()

◆ getEntityRowIndices()

◆ getEntitySense()

◆ getFaceTriNodes()

◆ getMaxColOrder()

◆ getMaxDataOrder()

◆ getMaxRowOrder()

◆ getNodesFieldData() [1/2]

◆ getNodesFieldData() [2/2]

◆ getNodesIndices()

◆ getNoFieldColIndices()

◆ getNoFieldFieldData() [1/2]

◆ getNoFieldFieldData() [2/2]

◆ getNoFieldIndices()

◆ getNoFieldRowIndices()

◆ getNumberOfNodes()

◆ getOpPtrVector()

◆ getProblemNodesColIndices()

◆ getProblemNodesIndices()

◆ getProblemNodesRowIndices()

◆ getProblemTypeColIndices()

◆ getProblemTypeIndices()

◆ getProblemTypeRowIndices()

◆ getRowNodesIndices()

◆ getRule() [1/2]

◆ getRule() [2/2]

◆ getSense()

◆ getSpacesAndBaseOnEntities()

◆ getUserPolynomialBase()

◆ loopOverOperators()

◆ operator()()

◆ postProcess()

◆ preProcess()

◆ setGaussPts() [1/2]

◆ setGaussPts() [2/2]

Member Data Documentation

◆ dataH1

◆ dataHcurl

◆ dataHdiv

◆ dataL2

◆ dataNoField

◆ dataOnElement

◆ derivedDataOnElement

◆ elementPolynomialBasePtr

◆ gaussPts

◆ getRuleHook

◆ lastEvaluatedElementEntityType

◆ mField

◆ opPtrVector

◆ userPolynomialBasePtr