Classes and functions used to evaluate fields at integration pts, jacobians, etc.. More...

Collaboration diagram for Forms Integrators:

Classes
struct	MoFEM::FormsIntegrators< EleOp >::Assembly< A >::LinearForm< I >
	Linear form. More...

struct	MoFEM::FormsIntegrators< EleOp >::Assembly< A >::BiLinearForm< I >
	Bi linear form. More...

struct	MoFEM::FormsIntegrators< EleOp >::Assembly< A >::TriLinearForm< I >
	Tri linear form. More...

Typedefs
template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM>
using	MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradGrad = OpGradGradImpl< BASE_DIM, FIELD_DIM, SPACE_DIM, I, OpBase >
	Integrate $(v_{,i},\beta(\mathbf{x}) u_{,j}))_\Omega$. More...

template<int BASE_DIM, int FIELD_DIM>
using	MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpMass = OpMassImpl< BASE_DIM, FIELD_DIM, I, OpBase >
	Integrate $(v_i,\beta(\mathbf{x}) u_j)_\Omega$. More...

template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM, int S = 0>
using	MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradSymTensorGrad = OpGradSymTensorGradImpl< BASE_DIM, FIELD_DIM, SPACE_DIM, S, I, OpBase >
	Integrate $(v_k,D_{ijkl} u_{,l})_\Omega$. More...

template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM, int S = 0>
using	MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradGradSymTensorGradGrad = OpGradGradSymTensorGradGradImpl< BASE_DIM, FIELD_DIM, SPACE_DIM, S, I, OpBase >
	Integrate $(v_{,ij},D_{ijkl} u_{,kl})_\Omega$. More...

template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM, int S = 0>
using	MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradTensorGrad = OpGradTensorGradImpl< BASE_DIM, FIELD_DIM, SPACE_DIM, S, I, OpBase >
	Integrate $(v_{,j},D_{ijkl} u_{,l})_\Omega$. More...

using	MoFEM::ScalarFun = boost::function< double(const double, const double, const double)>
	Scalar function type. More...

template<int DIM>
using	MoFEM::VectorFun = boost::function< FTensor::Tensor1< double, DIM >(const double, const double, const double)>
	Vector function type. More...

template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM, int S = 1>
using	MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradTimesTensor = OpGradTimesTensorImpl< BASE_DIM, FIELD_DIM, SPACE_DIM, S, I, OpBase >
	Integrate $(v,f(\mathbf{x}))_\Omega$, f is a scalar. More...

Enumerations
enum	MoFEM::AssemblyType { MoFEM::PETSC, MoFEM::SCHUR, MoFEM::BLOCK_MAT, MoFEM::BLOCK_SCHUR, MoFEM::BLOCK_PRECONDITIONER_SCHUR, MoFEM::USER_ASSEMBLE, MoFEM::LAST_ASSEMBLE }
	[Storage and set boundary conditions] More...

enum	MoFEM::IntegrationType { MoFEM::GAUSS, MoFEM::USER_INTEGRATION, MoFEM::LAST_INTEGRATION }
	Form integrator integration types. More...

Functions
template<>
MoFEMErrorCode	MoFEM::VecSetValues< EssentialBcStorage > (Vec V, const EntitiesFieldData::EntData &data, const double *ptr, InsertMode iora)
	Set values to vector in operator. More...

Detailed Description

Classes and functions used to evaluate fields at integration pts, jacobians, etc..

Typedef Documentation

◆ OpGradGrad

template<typename EleOp >

template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM>

using MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradGrad = OpGradGradImpl<BASE_DIM, FIELD_DIM, SPACE_DIM, I, OpBase>

Integrate $(v_{,i},\beta(\mathbf{x}) u_{,j}))_\Omega$.

Template Parameters

SPACE_DIM

Definition at line 38 of file BiLinearFormsIntegrators.hpp.

◆ OpGradGradSymTensorGradGrad

template<typename EleOp >

template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM, int S = 0>

using MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradGradSymTensorGradGrad = OpGradGradSymTensorGradGradImpl<BASE_DIM, FIELD_DIM, SPACE_DIM, S, I, OpBase>

Integrate $(v_{,ij},D_{ijkl} u_{,kl})_\Omega$.

Note: $D_{ijkl}$ is * tensor with no symmetries

Template Parameters

SPACE_DIM
S

Definition at line 87 of file BiLinearFormsIntegrators.hpp.

◆ OpGradSymTensorGrad

template<typename EleOp >

template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM, int S = 0>

using MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradSymTensorGrad = OpGradSymTensorGradImpl<BASE_DIM, FIELD_DIM, SPACE_DIM, S, I, OpBase>

Integrate $(v_k,D_{ijkl} u_{,l})_\Omega$.

Note: $D_{ijkl}$ is * tensor with minor & major symmetry

Template Parameters

SPACE_DIM
S

Definition at line 72 of file BiLinearFormsIntegrators.hpp.

◆ OpGradTensorGrad

template<typename EleOp >

template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM, int S = 0>

using MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradTensorGrad = OpGradTensorGradImpl<BASE_DIM, FIELD_DIM, SPACE_DIM, S, I, OpBase>

Integrate $(v_{,j},D_{ijkl} u_{,l})_\Omega$.

Note: $D_{ijkl}$ is * tensor with no symmetries

Template Parameters

SPACE_DIM
S

Definition at line 101 of file BiLinearFormsIntegrators.hpp.

◆ OpGradTimesTensor

template<typename EleOp >

template<int BASE_DIM, int FIELD_DIM, int SPACE_DIM, int S = 1>

using MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpGradTimesTensor = OpGradTimesTensorImpl<BASE_DIM, FIELD_DIM, SPACE_DIM, S, I, OpBase>

Integrate $(v,f(\mathbf{x}))_\Omega$, f is a scalar.

Note: $f(\mathbf{x})\$ is scalar function of coordinates @tparam BASE_DIM @tparam FIELD_DIM */ template <int BASE_DIM, int FIELD_DIM, typename S = SourceFunctionSpecialization> using OpSource = OpSourceImpl<BASE_DIM, FIELD_DIM, I, typename S::template S<OpBase>>; /** @brief Vector field integrator $(v_i,f)_\Omega\f$, f is a vector

Template Parameters

BASE_DIM */ template <int BASE_DIM, int S = 1> using OpBaseTimesScalar = OpBaseTimesScalarImpl<BASE_DIM, S, I, OpBase>;

/**

Deprecated:: use instead OpBaseTimesScalar */ template <int BASE_DIM, int S = 1> using OpBaseTimesScalarField = OpBaseTimesScalar<BASE_DIM, S>;

/**

Vector field integrator $(v,f_i)_\Omega$, f is a vector

Template Parameters

BASE_DIM
FIELD_DIM
0	*/ template <int BASE_DIM, int FIELD_DIM, int S> using OpBaseTimesVector = OpBaseTimesVectorImpl<BASE_DIM, FIELD_DIM, S, I, OpBase>; //! [Source operator]

//! [Grad times tensor]

/**

Integrate $(v_{,i},f_{ij})_\Omega$, f is a vector

Note

Definition at line 82 of file LinearFormsIntegrators.hpp.

◆ OpMass

template<typename EleOp >

template<int BASE_DIM, int FIELD_DIM>

using MoFEM::FormsIntegrators< EleOp >::Assembly< A >::OpMass = OpMassImpl<BASE_DIM, FIELD_DIM, I, OpBase>

Integrate $(v_i,\beta(\mathbf{x}) u_j)_\Omega$.

Note: That FIELD_DIM = 4 or 9 is assumed that OpMass is for tensorial field 2x2 or 3x3

Todo:: It should be considered another template parameter RANK which will allow to distinguish between scalar, vectorial and tensorial fields

Template Parameters

BASE_DIM	dimension of base
FIELD_DIM	dimension of field

Definition at line 54 of file BiLinearFormsIntegrators.hpp.

◆ ScalarFun

using MoFEM::ScalarFun = typedef boost::function<double(const double, const double, const double)>

Scalar function type.

Examples: PoissonDiscontinousGalerkin.hpp.

Definition at line 144 of file FormsIntegrators.hpp.

◆ VectorFun

template<int DIM>

using MoFEM::VectorFun = typedef boost::function<FTensor::Tensor1<double, DIM>( const double, const double, const double)>

Vector function type.

Template Parameters

DIM	dimension of the return

Definition at line 176 of file FormsIntegrators.hpp.

Enumeration Type Documentation

◆ AssemblyType

enum MoFEM::AssemblyType

[Storage and set boundary conditions]

Form integrator assembly types

Enumerator
PETSC
SCHUR
BLOCK_MAT
BLOCK_SCHUR
BLOCK_PRECONDITIONER_SCHUR
USER_ASSEMBLE
LAST_ASSEMBLE

Definition at line 104 of file FormsIntegrators.hpp.

                   {
   PETSC,
   SCHUR,
   BLOCK_MAT,
   BLOCK_SCHUR,
   BLOCK_PRECONDITIONER_SCHUR,
   USER_ASSEMBLE,
   LAST_ASSEMBLE
 };

◆ IntegrationType

enum MoFEM::IntegrationType

Form integrator integration types.

Enumerator
GAUSS
USER_INTEGRATION
LAST_INTEGRATION

Definition at line 136 of file FormsIntegrators.hpp.

136 { GAUSS, USER_INTEGRATION, LAST_INTEGRATION };

Function Documentation

◆ VecSetValues< EssentialBcStorage >()

template<>

MoFEMErrorCode MoFEM::VecSetValues< EssentialBcStorage >	(	Vec	V,
		const EntitiesFieldData::EntData &	data,
		const double *	ptr,
		InsertMode	iora
	)

Set values to vector in operator.

MoFEM::FieldEntity provides MoFEM::FieldEntity::getWeakStoragePtr() storage function which allows to transfer data between FEs or operators processing the same entities.

When MoFEM::OpSetBc is pushed in weak storage indices taking in account indices which are skip to take boundary conditions are stored. Those entities are used by VecSetValues.

Parameters

V
data
ptr
iora

Returns: MoFEMErrorCode

Parameters

V
data
ptr
iora

Returns: MoFEMErrorCode

Examples: test_cache_on_entities.cpp.

Definition at line 76 of file FormsIntegrators.cpp.

                                                                      {
   MoFEMFunctionBegin;
  
   if(!V)
     SETERRQ(PETSC_COMM_SELF, MOFEM_INVALID_DATA,
             "Pointer to PETSc vector is null");
  
   CHKERR VecSetOption(V, VEC_IGNORE_NEGATIVE_INDICES, PETSC_TRUE);
   if (!data.getFieldEntities().empty()) {
     if (auto e_ptr = data.getFieldEntities()[0]) {
       if (auto stored_data_ptr =
               e_ptr->getSharedStoragePtr<EssentialBcStorage>()) {
         return ::VecSetValues(V, stored_data_ptr->entityIndices.size(),
                               &*stored_data_ptr->entityIndices.begin(), ptr,
                               iora);
       }
     }
   }
   return ::VecSetValues(V, data.getIndices().size(),
                         &*data.getIndices().begin(), ptr, iora);
   MoFEMFunctionReturn(0);
 }

Classes

Typedefs

Enumerations

Functions

Detailed Description

Typedef Documentation

◆ OpGradGrad

◆ OpGradGradSymTensorGradGrad

◆ OpGradSymTensorGrad

◆ OpGradTensorGrad

◆ OpGradTimesTensor

◆ OpMass

◆ ScalarFun

◆ VectorFun

Enumeration Type Documentation

◆ AssemblyType

◆ IntegrationType

Function Documentation

◆ VecSetValues< EssentialBcStorage >()