Data on single entity (This is passed as argument to DataOperator::doWork) More...

#include "src/finite_elements/EntitiesFieldData.hpp"

Inheritance diagram for MoFEM::EntitiesFieldData::EntData:

Collaboration diagram for MoFEM::EntitiesFieldData::EntData:

Public Types
enum	BaseDerivatives { ZeroDerivative = 0 , FirstDerivative , SecondDerivative , ThirdDerivative , ForthDerivative , LastDerivative }

Bernstein-Bezier base only functions
static constexpr size_t	MaxBernsteinBezierOrder = BITFEID_SIZE

VectorInt &	getBBNodeOrder ()
	Get orders at the nodes.

MatrixInt &	getBBAlphaIndices ()
	Get file BB indices.

virtual boost::shared_ptr< MatrixInt > &	getBBAlphaIndicesSharedPtr (const std::string &field_name)

virtual boost::shared_ptr< MatrixDouble > &	getBBNSharedPtr (const std::string &field_name)

virtual const boost::shared_ptr< MatrixDouble > &	getBBNSharedPtr (const std::string &field_name) const

virtual boost::shared_ptr< MatrixDouble > &	getBBDiffNSharedPtr (const std::string &field_name)

virtual const boost::shared_ptr< MatrixDouble > &	getBBDiffNSharedPtr (const std::string &field_name) const

virtual std::map< std::string, boost::shared_ptr< MatrixInt > > &	getBBAlphaIndicesMap ()

virtual std::map< std::string, boost::shared_ptr< MatrixDouble > > &	getBBNMap ()
	get hash map of base function for BB base, key is a field name

virtual std::map< std::string, boost::shared_ptr< MatrixDouble > > &	getBBDiffNMap ()
	get hash map of derivatives base function for BB base, key is a field name

virtual boost::shared_ptr< MatrixInt > &	getBBAlphaIndicesByOrderSharedPtr (const size_t o)
	get ALpha indices for BB base by order

virtual boost::shared_ptr< MatrixDouble > &	getBBNByOrderSharedPtr (const size_t o)
	get BB base by order

virtual boost::shared_ptr< MatrixDouble > &	getBBDiffNByOrderSharedPtr (const size_t o)
	get BB base derivative by order

virtual std::array< boost::shared_ptr< MatrixInt >, MaxBernsteinBezierOrder > &	getBBAlphaIndicesByOrderArray ()

virtual std::array< boost::shared_ptr< MatrixDouble >, MaxBernsteinBezierOrder > &	getBBNByOrderArray ()

virtual std::array< boost::shared_ptr< MatrixDouble >, MaxBernsteinBezierOrder > &	getBBDiffNByOrderArray ()

virtual MoFEMErrorCode	baseSwap (const std::string &field_name, const FieldApproximationBase base)
	Swap bases functions.

Broken spaces functions
int	sEnse
	Entity sense (orientation)

ApproximationOrder	oRder
	Entity order.

FieldSpace	sPace
	Entity space.

FieldApproximationBase	bAse
	Field approximation base.

VectorInt	iNdices
	Global indices on entity.

VectorInt	localIndices
	Local indices on entity.

VectorDofs	dOfs
	DoFs on entity.

VectorFieldEntities	fieldEntities
	Field entities.

VectorDouble	fieldData
	Field data on entity.

std::vector< BitRefLevel >	entDataBitRefLevel
	Bit ref level in entity.

std::vector< int >	dofBrokenSideVec
	Map side to dofs number.

std::vector< EntityType >	dofBrokenTypeVec
	Map type of entity to dof number.

std::array< std::array< boost::shared_ptr< MatrixDouble >, LASTBASE >, LastDerivative >	baseFunctionsAndBaseDerivatives

std::array< boost::shared_ptr< MatrixDouble >, LASTBASE > &	N
	Base functions.

std::array< boost::shared_ptr< MatrixDouble >, LASTBASE > &	diffN
	Derivatives of base functions.

std::string	bbFieldName
	field name

VectorInt	bbNodeOrder
	order of nodes

std::map< std::string, boost::shared_ptr< MatrixDouble > >	bbN

std::map< std::string, boost::shared_ptr< MatrixDouble > >	bbDiffN

std::map< std::string, boost::shared_ptr< MatrixInt > >	bbAlphaIndices
	Indices for Bernstein-Bezier (BB) base.

std::array< boost::shared_ptr< MatrixDouble >, MaxBernsteinBezierOrder >	bbNByOrder
	BB base functions by order.

std::array< boost::shared_ptr< MatrixDouble >, MaxBernsteinBezierOrder >	bbDiffNByOrder
	BB base functions derivatives by order.

std::array< boost::shared_ptr< MatrixInt >, MaxBernsteinBezierOrder >	bbAlphaIndicesByOrder
	BB alpha indices by order.

boost::shared_ptr< MatrixDouble >	swapBaseNPtr
	Used by Bernstein base to keep temporary pointer.

boost::shared_ptr< MatrixDouble >	swapBaseDiffNPtr
	Used by Bernstein base to keep temporary pointer.

struct	OpAddParentEntData

template<typename OpBase >
struct	OpGetBrokenBaseSideData

Auxiliary functions
std::ostream &	operator<< (std::ostream &os, const EntitiesFieldData::EntData &e)

MoFEMErrorCode	resetFieldDependentData ()

Constructor and destructor
	EntData (const bool allocate_base_matrices=true)
	Construct EntData object.

virtual	~EntData ()=default

Sense, order and indices
virtual int	getSense () const
	Get entity sense for conforming approximation fields.

ApproximationOrder	getOrder () const
	Get approximation order.

const VectorInt &	getIndices () const
	Get global indices of degrees of freedom on entity.

const VectorIntAdaptor	getIndicesUpToOrder (int order)
	Get global indices of DOFs on entity up to given order.

const VectorInt &	getLocalIndices () const
	Get local indices of degrees of freedom on entity.

const VectorIntAdaptor	getLocalIndicesUpToOrder (int order)
	Get local indices of DOFs on entity up to given order.

int &	getSense ()
	Get reference to entity sense for modification.

ApproximationOrder &	getOrder ()
	Get reference to approximation order for modification.

VectorInt &	getIndices ()
	Get reference to global indices for modification.

VectorInt &	getLocalIndices ()
	Get reference to local indices for modification.

Data on entity
const VectorDouble &	getFieldData () const
	Get DOF values on entity.

const VectorAdaptor	getFieldDataUpToOrder (int order)
	Get DOF values on entity up to given order.

const VectorDofs &	getFieldDofs () const
	Get DOF data structures (const version)

VectorDouble &	getFieldData ()
	Get DOF values on entity for modification.

const VectorFieldEntities &	getFieldEntities () const
	Get field entities (const version)

VectorFieldEntities &	getFieldEntities ()
	Get field entities for modification.

virtual std::vector< BitRefLevel > &	getEntDataBitRefLevel ()
	Get entity bit reference level.

template<int Tensor_Dim>
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim >	getFTensor1FieldData ()
	Return FTensor of rank 1, i.e. vector from field data coefficients.

template<int Tensor_Dim0, int Tensor_Dim1>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 >	getFTensor2FieldData ()
	Return FTensor rank 2, i.e. matrix from field data coefficients.

template<int Tensor_Dim>
FTensor::Tensor2_symmetric< FTensor::PackPtr< double ,(Tensor_Dim (Tensor_Dim+1))/2 >, Tensor_Dim >	getFTensor2SymmetricFieldData ()
	Return symmetric FTensor rank 2, i.e. matrix from field data coefficients.

FTensor::Tensor0< FTensor::PackPtr< double *, 1 > >	getFTensor0FieldData ()
	Return scalar files as a FTensor of rank 0.

VectorDofs &	getFieldDofs ()

Base and space
FieldApproximationBase &	getBase ()
	Get approximation base.

FieldSpace &	getSpace ()
	Get field space.

virtual boost::shared_ptr< MatrixDouble > &	getNSharedPtr (const FieldApproximationBase base, const BaseDerivatives derivative)
	Get shared pointer to base functions with derivatives.

virtual boost::shared_ptr< MatrixDouble > &	getNSharedPtr (const FieldApproximationBase base)
	Get shared pointer to base functions (zero derivative)

virtual boost::shared_ptr< MatrixDouble > &	getDiffNSharedPtr (const FieldApproximationBase base)
	Get shared pointer to derivatives of base functions.

Get base functions for H1/L2
MatrixDouble &	getN (const FieldApproximationBase base)
	get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb. of columns is equal to number of base functions on this entity.

MatrixDouble &	getN (const std::string &field_name)
	Get base functions for specific field (Bernstein-Bezier)

MatrixDouble &	getN ()
	Get base functions using default base.

MatrixDouble &	getDiffN (const FieldApproximationBase base)
	get derivatives of base functions

MatrixDouble &	getN (const FieldApproximationBase base, const BaseDerivatives derivative)
	Get base function derivatives with specific derivative order.

MatrixDouble &	getDiffN (const std::string &field_name)
	Get derivatives of base functions for specific field.

MatrixDouble &	getDiffN ()
	Get derivatives of base functions using default base.

MatrixDouble &	getN (const BaseDerivatives derivative)
	Get base function derivatives using default base.

const VectorAdaptor	getN (const FieldApproximationBase base, const int gg)
	Get base functions at specific Gauss point.

const VectorAdaptor	getN (const int gg)
	Get base functions at specific Gauss point using default base.

const MatrixAdaptor	getDiffN (const FieldApproximationBase base, const int gg)
	Get derivative of base functions at specific Gauss point.

const MatrixAdaptor	getDiffN (const int gg)
	Get derivative of base functions at specific Gauss point using default base.

const VectorAdaptor	getN (const FieldApproximationBase base, const int gg, const int nb_base_functions)
	Get base functions at Gauss point with limited number of functions.

const VectorAdaptor	getN (const int gg, const int nb_base_functions)
	Get base functions at Gauss point with limited number (default base)

const MatrixAdaptor	getDiffN (const FieldApproximationBase base, const int gg, const int nb_base_functions)
	Get derivatives of base functions at Gauss point with limited number.

const MatrixAdaptor	getDiffN (const int gg, const int nb_base_functions)
	Get derivatives of base functions at Gauss point with limited number (default base)

Get base functions for vectorial approximation bases, i.e.
Hdiv/Hcurl
template<int DIM>
const MatrixAdaptor	getVectorN (const FieldApproximationBase base, const int gg)
	get Hdiv of base functions at Gauss pts

template<int DIM>
const MatrixAdaptor	getVectorN (const int gg)
	get Hdiv of base functions at Gauss pts

template<int DIM0, int DIM1>
const MatrixAdaptor	getVectorDiffN (FieldApproximationBase base, const int gg)
	get DiffHdiv of base functions at Gauss pts

template<int DIM0, int DIM1>
const MatrixAdaptor	getVectorDiffN (const int gg)
	get DiffHdiv of base functions at Gauss pts

template<int DIM0, int DIM1>
const MatrixAdaptor	getVectorDiffN (const FieldApproximationBase base, const int dof, const int gg)
	get DiffHdiv of base functions at Gauss pts

template<int DIM0, int DIM1>
const MatrixAdaptor	getVectorDiffN (const int dof, const int gg)
	get DiffHdiv of base functions at Gauss pts

Get base functions with FTensor
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > >	getFTensor0N (const FieldApproximationBase base)
	Get base function as Tensor0.

FTensor::Tensor0< FTensor::PackPtr< double *, 1 > >	getFTensor0N ()
	Get base function as Tensor0.

FTensor::Tensor0< FTensor::PackPtr< double *, 1 > >	getFTensor0N (const FieldApproximationBase base, const int gg, const int bb)
	Get base function as Tensor0 (Loop by integration points)

FTensor::Tensor0< FTensor::PackPtr< double *, 1 > >	getFTensor0N (const int gg, const int bb)
	Get base function as Tensor0 (Loop by integration points)

template<int Tensor_Dim>
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim >	getFTensor1DiffN (const FieldApproximationBase base)
	Get derivatives of base functions.

template<int Tensor_Dim>
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim >	getFTensor1DiffN ()
	Get derivatives of base functions.

template<int Tensor_Dim>
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim >	getFTensor1DiffN (const FieldApproximationBase base, const int gg, const int bb)
	Get derivatives of base functions (Loop by integration points)

template<int Tensor_Dim>
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim >	getFTensor1DiffN (const int gg, const int bb)
	Get derivatives of base functions (Loop by integration points)

template<int Tensor_Dim>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim Tensor_Dim >, Tensor_Dim, Tensor_Dim >	getFTensor2DiffN2 (const FieldApproximationBase base)
	Get second derivatives of scalar base functions.

template<int Tensor_Dim>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim Tensor_Dim >, Tensor_Dim, Tensor_Dim >	getFTensor2DiffN2 (const FieldApproximationBase base, const int gg, const int bb)
	Get second derivatives of scalar base functions at integration pts.

template<int Tensor_Dim>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim Tensor_Dim >, Tensor_Dim, Tensor_Dim >	getFTensor2DiffN2 ()
	Get second derivatives of scalar base functions.

template<int Tensor_Dim>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim Tensor_Dim >, Tensor_Dim, Tensor_Dim >	getFTensor2DiffN2 (const int gg, const int bb)
	Get second derivatives of scalar base functions at integration pts.

template<int Tensor_Dim>
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim >	getFTensor1N (FieldApproximationBase base)
	Get base functions for Hdiv/Hcurl spaces.

template<int Tensor_Dim>
auto	getFTensor1N ()
	Get base functions for Hdiv space.

template<int Tensor_Dim0, int Tensor_Dim1>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 >	getFTensor2DiffN (FieldApproximationBase base)
	Get derivatives of base functions for Hdiv space.

template<int Tensor_Dim0, int Tensor_Dim1>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 >	getFTensor2DiffN (FieldApproximationBase base, const int gg, const int bb)
	Get derivatives of base functions for Hdiv space at integration pts.

template<int Tensor_Dim0, int Tensor_Dim1>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 >	getFTensor2DiffN ()
	Get derivatives of base functions for Hdiv space.

template<int Tensor_Dim0, int Tensor_Dim1>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 >	getFTensor2DiffN (const int gg, const int bb)
	Get derivatives of base functions for Hdiv space at integration pts.

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>
FTensor::Tensor3< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 >	getFTensor3Diff2N (FieldApproximationBase base)
	Get second derivatives of base functions for Hvec space.

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>
FTensor::Tensor3< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 >	getFTensor3Diff2N ()
	Get second derivatives of base functions for Hvec space.

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>
FTensor::Tensor3< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 >	getFTensor3DiffN (FieldApproximationBase base)
	Get derivatives of base functions for tonsorial Hdiv space.

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>
FTensor::Tensor3< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 >	getFTensor3DiffN (FieldApproximationBase base, const int gg, const int bb)
	Get derivatives of base functions for tonsorial Hdiv space at integration pts.

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>
FTensor::Tensor3< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 >	getFTensor3DiffN ()
	Get derivatives of base functions for tonsorial Hdiv space.

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>
FTensor::Tensor3< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 >	getFTensor3DiffN (const int gg, const int bb)
	Get derivatives of base functions for tonsorial Hdiv space at integration pts.

template<int Tensor_Dim>
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim >	getFTensor1N (FieldApproximationBase base, const int gg, const int bb)
	Get Hdiv base functions at integration point.

template<int Tensor_Dim>
auto	getFTensor1N (const int gg, const int bb)
	Get Hdiv base functions at integration point.

template<int Tensor_Dim0, int Tensor_Dim1>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 >	getFTensor2N (FieldApproximationBase base)
	Get base functions for Hdiv/Hcurl spaces.

template<int Tensor_Dim0, int Tensor_Dim1>
auto	getFTensor2N ()
	Get base functions for Hdiv space.

template<int Tensor_Dim0, int Tensor_Dim1>
FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 >	getFTensor2N (FieldApproximationBase base, const int gg, const int bb)
	Get base functions for tensor Hdiv/Hcurl spaces.

template<int Tensor_Dim0, int Tensor_Dim1>
auto	getFTensor2N (const int gg, const int bb)
	Get base functions for Hdiv space.

Specializations for tensor base function
template<>
FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 >	getFTensor1N (FieldApproximationBase base)

template<>
FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 >	getFTensor1N (FieldApproximationBase base, const int gg, const int bb)

template<>
FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 >	getFTensor2N (FieldApproximationBase base)

Specializations for derivatives of base functions
template<>
FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 >	getFTensor1DiffN (const FieldApproximationBase base)
	Get spatial derivative of base function tensor for dimension 3d.

template<>
FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 >	getFTensor1DiffN ()
	Get spatial derivative of base function tensor for dimension 3d.

template<>
FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 >	getFTensor2DiffN (FieldApproximationBase base)

template<>
FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 >	getFTensor2DiffN (FieldApproximationBase base, const int gg, const int bb)

template<>
FTensor::Tensor3< FTensor::PackPtr< double *, 12 >, 3, 2, 2 >	getFTensor3Diff2N (FieldApproximationBase base)

template<>
FTensor::Tensor3< FTensor::PackPtr< double *, 27 >, 3, 3, 3 >	getFTensor3DiffN (FieldApproximationBase base)

template<>
FTensor::Tensor3< FTensor::PackPtr< double *, 27 >, 3, 3, 3 >	getFTensor3DiffN (FieldApproximationBase base, const int gg, const int bb)

template<>
FTensor::Tensor2< FTensor::PackPtr< double *, 4 >, 2, 2 >	getFTensor2DiffN2 (const FieldApproximationBase base)

template<>
FTensor::Tensor2< FTensor::PackPtr< double *, 4 >, 2, 2 >	getFTensor2DiffN2 (const FieldApproximationBase base, const int gg, const int bb)

template<>
FTensor::Tensor2< FTensor::PackPtr< double *, 4 >, 2, 2 >	getFTensor2DiffN2 ()

template<>
FTensor::Tensor2< FTensor::PackPtr< double *, 4 >, 2, 2 >	getFTensor2DiffN2 (const int gg, const int bb)

Specializations for field data
template<>
FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 >	getFTensor1FieldData ()

template<>
FTensor::Tensor2< FTensor::PackPtr< double *, 1 >, 1, 1 >	getFTensor2FieldData ()

template<>
FTensor::Tensor2_symmetric< FTensor::PackPtr< double *, 6 >, 3 >	getFTensor2SymmetricFieldData ()

Specializations for H1/L2
template<>
FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 >	getFTensor1DiffN (const FieldApproximationBase base, const int gg, const int bb)
	Get spatial derivative of base function tensor for dimension 3d.

template<>
FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 >	getFTensor1DiffN (const int gg, const int bb)
	Get spatial derivative of base function tensor for dimension 3d.

Specializations for HDiv/HCrul
template<>
FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 >	getFTensor2N (FieldApproximationBase base, const int gg, const int bb)

Detailed Description

Data on single entity (This is passed as argument to DataOperator::doWork)

Todo:: Hdiv and Hcurl functions should be accessed through common interface.

Examples: /home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianPlasticity.cpp, ContactOps.hpp, HenckyOps.hpp, PlasticOpsGeneric.hpp, ThermoElasticOps.hpp, mofem/atom_tests/boundary_marker.cpp, mofem/atom_tests/build_large_problem.cpp, mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/continuity_check_on_contact_prism_side_ele.cpp, mofem/atom_tests/continuity_check_on_skeleton_3d.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/dg_projection.cpp, mofem/atom_tests/field_evaluator.cpp, mofem/atom_tests/forces_and_sources_testing_edge_element.cpp, mofem/atom_tests/forces_and_sources_testing_flat_prism_element.cpp, mofem/atom_tests/forces_and_sources_testing_users_base.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_check_approx_in_2d.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hcurl_divergence_operator_2d.cpp, mofem/atom_tests/hdiv_check_approx_in_3d.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/higher_derivatives.cpp, mofem/atom_tests/operators_tests.cpp, mofem/atom_tests/prism_elements_from_surface.cpp, mofem/atom_tests/prism_polynomial_approximation.cpp, mofem/atom_tests/quad_polynomial_approximation.cpp, mofem/atom_tests/scalar_check_approximation.cpp, mofem/atom_tests/schur_test_diag_mat.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/simple_l2_only.cpp, mofem/atom_tests/tensor_divergence_operator.cpp, mofem/atom_tests/test_broken_space.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tutorials/adv-0_plasticity/plastic.cpp, mofem/tutorials/adv-0_plasticity/src/PlasticOpsGeneric.hpp, mofem/tutorials/adv-1_contact/src/ContactOps.hpp, mofem/tutorials/adv-3_level_set/level_set.cpp, mofem/tutorials/adv-4_dynamic_first_order_con_law/dynamic_first_order_con_law.cpp, mofem/tutorials/adv-5_poroelasticity/seepage.cpp, mofem/tutorials/adv-6_dg_projection/between_meshes_dg_projection.cpp, mofem/tutorials/cor-2to5_analytical_poisson/src/PoissonOperators.hpp, mofem/tutorials/cor-6_simple_elasticity/simple_elasticity.cpp, mofem/tutorials/cor-9_reaction_diffusion/reaction_diffusion.cpp, mofem/tutorials/fun-0_hello_world/hello_world.cpp, mofem/tutorials/fun-1_integration/integration.cpp, mofem/tutorials/fun-2_plot_base/plot_base.cpp, mofem/tutorials/max-1_lorentz_force//lorentz_force.cpp, mofem/tutorials/mix-0_mixed_poisson/mixed_poisson.cpp, mofem/tutorials/mix-1_light_intensity_equation/phase.cpp, mofem/tutorials/scl-0_least_squares/approximaton.cpp, mofem/tutorials/scl-1_poisson_2d_homogeneous/src/poisson_2d_homogeneous.hpp, mofem/tutorials/scl-6_heat_equation/heat_equation.cpp, mofem/tutorials/scl-7_wave_equation/wave_equation.cpp, mofem/tutorials/scl-8_heat_radiation/radiation.cpp, mofem/tutorials/scl-9_heat_method/heat_method.cpp, mofem/tutorials/vec-0_elasticity/elastic.cpp, mofem/tutorials/vec-10_elasticity_schur/schur_elastic.cpp, mofem/tutorials/vec-1_eigen_elasticity/eigen_elastic.cpp, mofem/tutorials/vec-2_nonlinear_elasticity/nonlinear_elastic.cpp, mofem/tutorials/vec-2_nonlinear_elasticity/src/HenckyOps.hpp, mofem/tutorials/vec-3_nonlinear_dynamics/nonlinear_dynamic_elastic.cpp, mofem/tutorials/vec-4_shallow_wave/approx_sphere.cpp, mofem/tutorials/vec-4_shallow_wave/shallow_wave.cpp, mofem/tutorials/vec-5_free_surface/free_surface.cpp, mofem/tutorials/vec-6_dg_kirchhoff_plate/plate.cpp, mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp, mofem/tutorials/vec-9_arc_length/arc_length.cpp, mofem/users_modules/adolc-plasticity/adolc_plasticity.cpp, mofem/users_modules/basic_finite_elements/src/HookeElement.hpp, mofem/users_modules/basic_finite_elements/src/NavierStokesElement.hpp, mofem/users_modules/basic_finite_elements/src/impl/HookeElement.cpp, mofem/users_modules/basic_finite_elements/src/impl/NavierStokesElement.cpp, mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp, nonlinear_elastic.cpp, plastic.cpp, thermo_elastic.cpp, and thermoplastic.cpp.

Definition at line 152 of file EntitiesFieldData.hpp.

Member Enumeration Documentation

◆ BaseDerivatives

enum MoFEM::EntitiesFieldData::EntData::BaseDerivatives

Enumerator
ZeroDerivative
FirstDerivative
SecondDerivative
ThirdDerivative
ForthDerivative
LastDerivative

Definition at line 154 of file EntitiesFieldData.hpp.

                       {
    ZeroDerivative = 0,
    FirstDerivative,
    SecondDerivative,
    ThirdDerivative,
    ForthDerivative,
    LastDerivative
  };

Constructor & Destructor Documentation

◆ EntData()

MoFEM::EntitiesFieldData::EntData::EntData ( const bool allocate_base_matrices = true )

Construct EntData object.

Parameters

allocate_base_matrices If true, allocate base function matrices (default: true)

Definition at line 10 of file EntitiesFieldData.cpp.

    : sEnse(0), oRder(0), bAse(NOBASE), entDataBitRefLevel(),
      N(baseFunctionsAndBaseDerivatives[ZeroDerivative]),
      diffN(baseFunctionsAndBaseDerivatives[FirstDerivative]) {
  if (allocate_base_matrices) {
 
    for (auto d = 0; d != LastDerivative; ++d) {
      for (int b = 0; b != LASTBASE; ++b) {
        baseFunctionsAndBaseDerivatives[d][b].reset(new MatrixDouble());
      }
    }
  } else {
 
    // set null to all base function pointers
    for (auto d = 0; d != LastDerivative; ++d) {
      std::fill(baseFunctionsAndBaseDerivatives[d].begin(),
                baseFunctionsAndBaseDerivatives[d].end(), nullptr);
    }
  }
}

◆ ~EntData()

virtual MoFEM::EntitiesFieldData::EntData::~EntData ( )

virtualdefault

Member Function Documentation

◆ baseSwap()

MoFEMErrorCode MoFEM::EntitiesFieldData::EntData::baseSwap	(	const std::string &	field_name,
		const FieldApproximationBase	base
	)

virtual

Swap bases functions.

Some base are not hierarchical and depend on approximation order. Such case demand special handling, that appropriate base order is set depending on field, such that is accessible in operator.

Note: Base is not swap on meshsets

Parameters

field_name
base

Returns: MoFEMErrorCode

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 143 of file EntitiesFieldData.cpp.

                                                                        {
  MoFEMFunctionBegin;
  auto make_swap = [](boost::shared_ptr<MatrixDouble> &ptr,
                      boost::shared_ptr<MatrixDouble> &ptrBB,
                      boost::shared_ptr<MatrixDouble> &swap_ptr) {
    if (swap_ptr) {
      ptr = swap_ptr;
      swap_ptr.reset();
    } else {
      swap_ptr = ptr;
      ptr = ptrBB;
    }
  };
  make_swap(getNSharedPtr(base), getBBNSharedPtr(field_name), swapBaseNPtr);
  make_swap(getDiffNSharedPtr(base), getBBDiffNSharedPtr(field_name),
            swapBaseDiffNPtr);
  MoFEMFunctionReturn(0);
}

◆ getBase()

FieldApproximationBase & MoFEM::EntitiesFieldData::EntData::getBase ( )

inline

Get approximation base.

Returns: Approximation base

Examples: mofem/atom_tests/hanging_node_approx.cpp.

Definition at line 1513 of file EntitiesFieldData.hpp.

1513{ return bAse; }

◆ getBBAlphaIndices()

MatrixInt & MoFEM::EntitiesFieldData::EntData::getBBAlphaIndices ( )

inline

Get file BB indices.

Returns: MatrixInt&

Definition at line 1771 of file EntitiesFieldData.hpp.

                                                       {
  return *getBBAlphaIndicesSharedPtr(bbFieldName);
}

◆ getBBAlphaIndicesByOrderArray()

std::array< boost::shared_ptr< MatrixInt >, EntitiesFieldData::EntData::MaxBernsteinBezierOrder > & MoFEM::EntitiesFieldData::EntData::getBBAlphaIndicesByOrderArray ( )

virtual

Definition at line 996 of file EntitiesFieldData.cpp.

                                                        {
  return bbAlphaIndicesByOrder;
}

◆ getBBAlphaIndicesByOrderSharedPtr()

boost::shared_ptr< MatrixInt > & MoFEM::EntitiesFieldData::EntData::getBBAlphaIndicesByOrderSharedPtr ( const size_t o )

virtual

get ALpha indices for BB base by order

Parameters

o	approximation order

Returns: boost::shared_ptr<MatrixInt>&

Definition at line 980 of file EntitiesFieldData.cpp.

                                                                          {
  return bbAlphaIndicesByOrder[o];
}

◆ getBBAlphaIndicesMap()

std::map< std::string, boost::shared_ptr< MatrixInt > > & MoFEM::EntitiesFieldData::EntData::getBBAlphaIndicesMap ( )

virtual

Definition at line 965 of file EntitiesFieldData.cpp.

                                               {
  return bbAlphaIndices;
}

◆ getBBAlphaIndicesSharedPtr()

boost::shared_ptr< MatrixInt > & MoFEM::EntitiesFieldData::EntData::getBBAlphaIndicesSharedPtr ( const std::string & field_name )

virtual

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 928 of file EntitiesFieldData.cpp.

                                 {
  return bbAlphaIndices[field_name];
}

◆ getBBDiffNByOrderArray()

std::array< boost::shared_ptr< MatrixDouble >, EntitiesFieldData::EntData::MaxBernsteinBezierOrder > & MoFEM::EntitiesFieldData::EntData::getBBDiffNByOrderArray ( )

virtual

Definition at line 1008 of file EntitiesFieldData.cpp.

                                                 {
  return bbDiffNByOrder;
}

◆ getBBDiffNByOrderSharedPtr()

boost::shared_ptr< MatrixDouble > & MoFEM::EntitiesFieldData::EntData::getBBDiffNByOrderSharedPtr ( const size_t o )

virtual

get BB base derivative by order

Parameters

o

Returns: boost::shared_ptr<MatrixDouble>&

Definition at line 990 of file EntitiesFieldData.cpp.

                                                                   {
  return bbDiffNByOrder[o];
}

◆ getBBDiffNMap()

std::map< std::string, boost::shared_ptr< MatrixDouble > > & MoFEM::EntitiesFieldData::EntData::getBBDiffNMap ( )

virtual

get hash map of derivatives base function for BB base, key is a field name

Returns: std::map<std::string, boost::shared_ptr<MatrixDouble>>&

Definition at line 975 of file EntitiesFieldData.cpp.

                                        {
  return bbDiffN;
}

◆ getBBDiffNSharedPtr() [1/2]

boost::shared_ptr< MatrixDouble > & MoFEM::EntitiesFieldData::EntData::getBBDiffNSharedPtr ( const std::string & field_name )

virtual

Get shared pointer to BB derivatives of base base functions

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 951 of file EntitiesFieldData.cpp.

                                                                         {
  return bbDiffN[field_name];
}

◆ getBBDiffNSharedPtr() [2/2]

const boost::shared_ptr< MatrixDouble > & MoFEM::EntitiesFieldData::EntData::getBBDiffNSharedPtr ( const std::string & field_name ) const

virtual

Get shared pointer to derivatives of BB base base functions

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 959 of file EntitiesFieldData.cpp.

                                       {
  return bbDiffN.at(field_name);
}

◆ getBBNByOrderArray()

std::array< boost::shared_ptr< MatrixDouble >, EntitiesFieldData::EntData::MaxBernsteinBezierOrder > & MoFEM::EntitiesFieldData::EntData::getBBNByOrderArray ( )

virtual

Definition at line 1002 of file EntitiesFieldData.cpp.

                                             {
  return bbNByOrder;
}

◆ getBBNByOrderSharedPtr()

boost::shared_ptr< MatrixDouble > & MoFEM::EntitiesFieldData::EntData::getBBNByOrderSharedPtr ( const size_t o )

virtual

get BB base by order

Parameters

o

Returns: boost::shared_ptr<MatrixDouble>&

Definition at line 985 of file EntitiesFieldData.cpp.

                                                               {
  return bbNByOrder[o];
}

◆ getBBNMap()

std::map< std::string, boost::shared_ptr< MatrixDouble > > & MoFEM::EntitiesFieldData::EntData::getBBNMap ( )

virtual

get hash map of base function for BB base, key is a field name

Returns: std::map<std::string, boost::shared_ptr<MatrixDouble>>&

Definition at line 970 of file EntitiesFieldData.cpp.

                                    {
  return bbN;
}

◆ getBBNodeOrder()

VectorInt & MoFEM::EntitiesFieldData::EntData::getBBNodeOrder ( )

inline

Get orders at the nodes.

Returns: VectorInt&

Definition at line 1769 of file EntitiesFieldData.hpp.

1769{ return bbNodeOrder; }

MoFEM::EntitiesFieldData::EntData::bbNodeOrder

VectorInt bbNodeOrder

order of nodes

Definition EntitiesFieldData.hpp:1293

◆ getBBNSharedPtr() [1/2]

boost::shared_ptr< MatrixDouble > & MoFEM::EntitiesFieldData::EntData::getBBNSharedPtr ( const std::string & field_name )

virtual

Get shared pointer to BB base base functions

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 934 of file EntitiesFieldData.cpp.

                                                                     {
  return bbN[field_name];
}

◆ getBBNSharedPtr() [2/2]

const boost::shared_ptr< MatrixDouble > & MoFEM::EntitiesFieldData::EntData::getBBNSharedPtr ( const std::string & field_name ) const

virtual

Get shared pointer to BB base base functions

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 942 of file EntitiesFieldData.cpp.

                                       {
  return bbN.at(field_name);
}

◆ getDiffN() [1/7]

MatrixDouble & MoFEM::EntitiesFieldData::EntData::getDiffN ( )

inline

Get derivatives of base functions using default base.

Returns: MatrixDouble& Reference to derivative matrix

Definition at line 1574 of file EntitiesFieldData.hpp.

1574{ return getDiffN(bAse); }

MoFEM::EntitiesFieldData::EntData::getDiffN

MatrixDouble & getDiffN()

Get derivatives of base functions using default base.

Definition EntitiesFieldData.hpp:1574

◆ getDiffN() [2/7]

MatrixDouble & MoFEM::EntitiesFieldData::EntData::getDiffN ( const FieldApproximationBase base )

inline

get derivatives of base functions

Matrix at rows has nb. of Gauss pts, at columns it has derivative of base functions. Columns are structured as follows, [ dN1/dx, dN1/dy, dN1/dz, dN2/dx, dN2/dy, dN2/dz, ... ]

Scalar base functions: Note that base functions are calculated in file H1.c Above description not apply for derivatives of nodal functions, since derivative of nodal functions in case of simplexes, EDGES, TRIANGLES and TETS are constant. So that matrix rows represents nb. of base functions, columns are derivatives. Nb. of columns depend on element dimension, for EDGES is one, for TRIS is 2 and TETS is 3.

Note: Note that for node element this function make no sense.

Tensorial base functions:

Note: Note: In rows are integration pts, columns are formatted that that components of vectors and then derivatives, for example row for given integration points is formatted in array
\[ t_{0,0}, t_{1,0}, t_{1,0}, t_{0,1}, t_{1,1}, t_{1,1}, t_{0,2}, t_{1,2}, t_{1,2} \]
where comma express derivative, i.e. \(t_{2,1} = \frac{\partial t_2}{\partial \xi_1}\)

Examples: mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, and mofem/atom_tests/scalar_check_approximation.cpp.

Definition at line 1545 of file EntitiesFieldData.hpp.

                                                                    {
  return *(getDiffNSharedPtr(base));
}

◆ getDiffN() [3/7]

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getDiffN	(	const FieldApproximationBase	base,
		const int	gg
	)

inline

Get derivative of base functions at specific Gauss point.

Returned matrix has base functions in rows and derivatives in columns.

Parameters

base	Approximation base type
gg	Gauss point number

Returns: const MatrixAdaptor Adaptor to derivative matrix at Gauss point

Definition at line 1594 of file EntitiesFieldData.hpp.

                                                   {
  // FIXME: That is bug, it will not work if number of integration pts is
  // equal to number of nodes on entity.  User who not implementing low
  // level DataOperator will not experience this.
  if (getN(base).size1() == getDiffN(base).size1()) {
    int size = getN(base).size2();
    int dim = getDiffN(base).size2() / size;
    double *data = &getDiffN(base)(gg, 0);
    return MatrixAdaptor(
        getN(base).size2(), dim,
        ublas::shallow_array_adaptor<double>(getDiffN(base).size2(), data));
  } else {
    // in some cases, f.e. for derivatives of nodal base functions at only
    // one gauss point is needed
    return MatrixAdaptor(
        getN(base).size1(), getN(base).size2(),
        ublas::shallow_array_adaptor<double>(getDiffN(base).data().size(),
                                             &getDiffN(base).data()[0]));
  }
}

◆ getDiffN() [4/7]

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getDiffN	(	const FieldApproximationBase	base,
		const int	gg,
		const int	nb_base_functions
	)

inline

Get derivatives of base functions at Gauss point with limited number.

Note that for multi-field elements, different fields can have different approximation orders. Since hierarchical approximation basis is used, base functions are calculated once for element using maximal order.

Parameters

base	Approximation base type
gg	Gauss point number
nb_base_functions	Number of base functions

Returns: const MatrixAdaptor Adaptor to limited derivative matrix

Definition at line 1636 of file EntitiesFieldData.hpp.

                                                                  {
  // FIXME: That is bug, it will not work if number of integration pts is
  // equal to number of nodes on entity.  User who not implementing low
  // level DataOperator will not experience this.
  if (getN(base).size1() == getDiffN(base).size1()) {
    (void)getN(base)(gg,
                     nb_base_functions -
                         1); // throw error if nb_base_functions is to big
    int dim = getDiffN(base).size2() / getN(base).size2();
    double *data = &getDiffN(base)(gg, 0);
    return MatrixAdaptor(
        nb_base_functions, dim,
        ublas::shallow_array_adaptor<double>(dim * nb_base_functions, data));
  } else {
    // in some cases, f.e. for derivatives of nodal base functions only one
    // gauss point is needed
    return MatrixAdaptor(
        getN(base).size1(), getN(base).size2(),
        ublas::shallow_array_adaptor<double>(getDiffN(base).data().size(),
                                             &getDiffN(base).data()[0]));
  }
}

◆ getDiffN() [5/7]

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getDiffN ( const int gg )

inline

Get derivative of base functions at specific Gauss point using default base.

Parameters

gg	Gauss point number

Returns: const MatrixAdaptor Adaptor to derivative matrix at Gauss point

Definition at line 1616 of file EntitiesFieldData.hpp.

                                                                   {
  return getDiffN(bAse, gg);
}

◆ getDiffN() [6/7]

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getDiffN	(	const int	gg,
		const int	nb_base_functions
	)

inline

Get derivatives of base functions at Gauss point with limited number (default base)

Parameters

gg	Gauss point number
nb_base_functions	Number of base functions

Returns: const MatrixAdaptor Adaptor to limited derivative matrix

Definition at line 1662 of file EntitiesFieldData.hpp.

                                                                  {
  return getDiffN(bAse, gg, nb_base_functions);
}

◆ getDiffN() [7/7]

MatrixDouble & MoFEM::EntitiesFieldData::EntData::getDiffN ( const std::string & field_name )

inline

Get derivatives of base functions for specific field.

Parameters

field_name Name of the field

Returns: MatrixDouble& Reference to derivative matrix

Definition at line 1564 of file EntitiesFieldData.hpp.

                                                              {
  #ifndef NDEBUG
  if (!getBBDiffNSharedPtr(field_name)) {
    CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                      "Null pointer to base functions derivative");
  }
  #endif // NDEBUG
  return *(getBBDiffNSharedPtr(field_name));
}

◆ getDiffNSharedPtr()

boost::shared_ptr< MatrixDouble > & MoFEM::EntitiesFieldData::EntData::getDiffNSharedPtr ( const FieldApproximationBase base )

virtual

Get shared pointer to derivatives of base functions.

Parameters

base	Approximation base type

Returns: boost::shared_ptr<MatrixDouble>& Shared pointer to derivative matrix

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 44 of file EntitiesFieldData.cpp.

                                       {
  return diffN[base];
}

◆ getEntDataBitRefLevel()

std::vector< BitRefLevel > & MoFEM::EntitiesFieldData::EntData::getEntDataBitRefLevel ( )

virtual

Get entity bit reference level.

Returns: std::vector<BitRefLevel>& Reference to bit reference level vector

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 1050 of file EntitiesFieldData.cpp.

                                                                        {
  return entDataBitRefLevel;
}

◆ getFieldData() [1/2]

VectorDouble & MoFEM::EntitiesFieldData::EntData::getFieldData ( )

inline

Get DOF values on entity for modification.

Returns: VectorDouble& Reference to vector of DOF values

Definition at line 1475 of file EntitiesFieldData.hpp.

1475{ return fieldData; }

MoFEM::EntitiesFieldData::EntData::fieldData

VectorDouble fieldData

Field data on entity.

Definition EntitiesFieldData.hpp:1277

◆ getFieldData() [2/2]

const VectorDouble & MoFEM::EntitiesFieldData::EntData::getFieldData ( ) const

inline

Get DOF values on entity.

Returns: const VectorDouble& Reference to vector of DOF values

Examples: mofem/atom_tests/continuity_check_on_skeleton_3d.cpp, mofem/atom_tests/forces_and_sources_testing_flat_prism_element.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/prism_elements_from_surface.cpp, mofem/tutorials/cor-2to5_analytical_poisson/src/PoissonOperators.hpp, mofem/tutorials/max-1_lorentz_force//lorentz_force.cpp, and mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 1454 of file EntitiesFieldData.hpp.

                                                                 {
  return fieldData;
}

◆ getFieldDataUpToOrder()

const VectorAdaptor MoFEM::EntitiesFieldData::EntData::getFieldDataUpToOrder ( int order )

inline

Get DOF values on entity up to given order.

Parameters

order Maximum approximation order to include

Returns: const VectorAdaptor Adaptor to subset of DOF values

Definition at line 1459 of file EntitiesFieldData.hpp.

                                                         {
  unsigned int size = 0;
  if (auto dof = dOfs[0]) {
    size = dof->getOrderNbDofs(order) * dof->getNbOfCoeffs();
    size = size < fieldData.size() ? size : fieldData.size();
  }
  double *data = &*fieldData.data().begin();
  return getVectorAdaptor(data, size);
}

◆ getFieldDofs() [1/2]

VectorDofs & MoFEM::EntitiesFieldData::EntData::getFieldDofs ( )

inline

Definition at line 1473 of file EntitiesFieldData.hpp.

1473{ return dOfs; }

◆ getFieldDofs() [2/2]

const VectorDofs & MoFEM::EntitiesFieldData::EntData::getFieldDofs ( ) const

inline

Get DOF data structures (const version)

Returns: const VectorDofs& Reference to vector of FEDofEntity pointers

Examples: mofem/atom_tests/prism_elements_from_surface.cpp, mofem/tutorials/cor-2to5_analytical_poisson/src/PoissonOperators.hpp, and mofem/users_modules/basic_finite_elements/src/impl/HookeElement.cpp.

Definition at line 1469 of file EntitiesFieldData.hpp.

                                                               {
  return dOfs;
}

◆ getFieldEntities() [1/2]

VectorFieldEntities & MoFEM::EntitiesFieldData::EntData::getFieldEntities ( )

inline

Get field entities for modification.

Returns: VectorFieldEntities& Reference to vector of FieldEntity pointers

Definition at line 1477 of file EntitiesFieldData.hpp.

                                                                {
  return fieldEntities;
}

◆ getFieldEntities() [2/2]

const VectorFieldEntities & MoFEM::EntitiesFieldData::EntData::getFieldEntities ( ) const

inline

Get field entities (const version)

Returns: const VectorFieldEntities& Reference to vector of FieldEntity pointers

Examples: mofem/atom_tests/test_cache_on_entities.cpp, mofem/tutorials/adv-6_dg_projection/between_meshes_dg_projection.cpp, and mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 1482 of file EntitiesFieldData.hpp.

                                                 {
  return fieldEntities;
}

◆ getFTensor0FieldData()

FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > MoFEM::EntitiesFieldData::EntData::getFTensor0FieldData ( )

Return scalar files as a FTensor of rank 0.

Returns: FTensor::Tensor0<FTensor::PackPtr<double *,1> >

Examples: mofem/atom_tests/continuity_check_on_skeleton_3d.cpp, and mofem/atom_tests/scalar_check_approximation.cpp.

Definition at line 513 of file EntitiesFieldData.cpp.

                                               {
#ifndef NDEBUG
  for (auto &d : dOfs) {
    if (d) {
      if (d->getNbOfCoeffs() != 1) {
        std::stringstream s;
        s << "Wrong number of coefficients is " << d->getNbOfCoeffs();
        s << " but expected scalar field, tensor of rank 0";
        CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY, s.str());
      }
      break;
    }
  }
#endif
  return FTensor::Tensor0<FTensor::PackPtr<double *, 1>>(
      &*fieldData.data().begin());
}

◆ getFTensor0N() [1/4]

FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > MoFEM::EntitiesFieldData::EntData::getFTensor0N ( )

inline

Get base function as Tensor0.

Return base functions for field base

Returns: Tensor0

Definition at line 1730 of file EntitiesFieldData.hpp.

                                       {
  return getFTensor0N(bAse);
};

◆ getFTensor0N() [2/4]

FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > MoFEM::EntitiesFieldData::EntData::getFTensor0N ( const FieldApproximationBase base )

inline

Get base function as Tensor0.

Parameters

base

Returns: Tensor0

Examples: mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/prism_polynomial_approximation.cpp, mofem/atom_tests/quad_polynomial_approximation.cpp, mofem/atom_tests/scalar_check_approximation.cpp, mofem/tutorials/adv-3_level_set/level_set.cpp, mofem/tutorials/cor-2to5_analytical_poisson/src/PoissonOperators.hpp, mofem/tutorials/cor-6_simple_elasticity/simple_elasticity.cpp, mofem/tutorials/cor-9_reaction_diffusion/reaction_diffusion.cpp, mofem/tutorials/scl-0_least_squares/approximaton.cpp, mofem/tutorials/scl-1_poisson_2d_homogeneous/src/poisson_2d_homogeneous.hpp, mofem/tutorials/scl-8_heat_radiation/radiation.cpp, mofem/tutorials/vec-4_shallow_wave/approx_sphere.cpp, mofem/tutorials/vec-4_shallow_wave/shallow_wave.cpp, mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp, mofem/users_modules/basic_finite_elements/src/impl/HookeElement.cpp, mofem/users_modules/basic_finite_elements/src/impl/NavierStokesElement.cpp, mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp, thermo_elastic.cpp, and thermoplastic.cpp.

Definition at line 1724 of file EntitiesFieldData.hpp.

                                                                        {
  double *ptr = &*getN(base).data().begin();
  return FTensor::Tensor0<FTensor::PackPtr<double *, 1>>(ptr);
};

◆ getFTensor0N() [3/4]

FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > MoFEM::EntitiesFieldData::EntData::getFTensor0N	(	const FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

inline

Get base function as Tensor0 (Loop by integration points)

Parameters

base
gg	integration points
bb	base function

Returns: Tensor0

Note that:

t0 = data.getFTensor0N(base,bb);

++t0

Increment in above code will move pointer to base function in next integration point.

Definition at line 1735 of file EntitiesFieldData.hpp.

                                                                     {
  double *ptr = &getN(base)(gg, bb);
  return FTensor::Tensor0<FTensor::PackPtr<double *, 1>>(ptr);
};

◆ getFTensor0N() [4/4]

FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > MoFEM::EntitiesFieldData::EntData::getFTensor0N	(	const int	gg,
		const int	bb
	)

inline

Get base function as Tensor0 (Loop by integration points)

Return base functions for field base

Parameters

bb	base function

Returns: Tensor0

Definition at line 1742 of file EntitiesFieldData.hpp.

                                                                 {
  return getFTensor0N(bAse, gg, bb);
};

◆ getFTensor1DiffN() [1/8]

template<int Tensor_Dim>

FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor1DiffN ( )

Get derivatives of base functions.

For volume element like tetrahedral or prism,

Tensor1<double*,3> diff_base = data.getFTensor1DiffN<3>();

Tensor1

Definition single.cpp:41

For face element like triangle or quad

Tensor1<double*,2> diff_base = data.getFTensor1DiffN<2>();

Returns: Tensor rank 1 (vector)

Definition at line 543 of file EntitiesFieldData.cpp.

                                           {
  return getFTensor1DiffN<Tensor_Dim>(bAse);
}

◆ getFTensor1DiffN() [2/8]

FTensor::Tensor1< FTensor::PackPtr< double *, 2 >, 2 > MoFEM::EntitiesFieldData::EntData::getFTensor1DiffN< 2 > ( )

Get spatial derivative of base function tensor for dimension 3d.

Get spatial derivative of base function tensor for dimension 2d.

Definition at line 543 of file EntitiesFieldData.cpp.

                                              {
  return getFTensor1DiffN<3>(bAse);
}

◆ getFTensor1DiffN() [3/8]

template<int Tensor_Dim>

FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor1DiffN ( const FieldApproximationBase base )

Get derivatives of base functions.

For volume element like tetrahedral or prism,

Tensor1<double*,3> diff_base = data.getFTensor1DiffN<3>();

For face element like triangle or quad

Tensor1<double*,2> diff_base = data.getFTensor1DiffN<2>();

Parameters

base functions

Returns: Tensor rank 1 (vector)

Examples: mofem/atom_tests/scalar_check_approximation.cpp, mofem/tutorials/adv-3_level_set/level_set.cpp, mofem/tutorials/cor-2to5_analytical_poisson/src/PoissonOperators.hpp, mofem/tutorials/cor-6_simple_elasticity/simple_elasticity.cpp, mofem/tutorials/cor-9_reaction_diffusion/reaction_diffusion.cpp, mofem/tutorials/scl-1_poisson_2d_homogeneous/src/poisson_2d_homogeneous.hpp, mofem/tutorials/scl-9_heat_method/heat_method.cpp, mofem/tutorials/vec-4_shallow_wave/shallow_wave.cpp, mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp, mofem/users_modules/basic_finite_elements/src/HookeElement.hpp, mofem/users_modules/basic_finite_elements/src/impl/HookeElement.cpp, mofem/users_modules/basic_finite_elements/src/impl/NavierStokesElement.cpp, and mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 533 of file EntitiesFieldData.cpp.

                                       {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor1<double *, Tensor_Dim>();
}

◆ getFTensor1DiffN() [4/8]

FTensor::Tensor1< FTensor::PackPtr< double *, 2 >, 2 > MoFEM::EntitiesFieldData::EntData::getFTensor1DiffN< 2 > ( const FieldApproximationBase base )

Get spatial derivative of base function tensor for dimension 3d.

Get spatial derivative of base function tensor for dimension 2d.

Definition at line 543 of file EntitiesFieldData.cpp.

                                       {
  double *ptr = &*getDiffN(base).data().begin();
  return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>(ptr, &ptr[1],
                                                            &ptr[2]);
}

◆ getFTensor1DiffN() [5/8]

template<int Tensor_Dim>

FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor1DiffN	(	const FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Get derivatives of base functions (Loop by integration points)

For volume element like tetrahedral or prism,

Tensor1<double*,3> diff_base = data.getFTensor1DiffN<3>(base,gg,bb);

where bb is base function and gg is integration pt. Operator ++diff_base will move tensor pointer to next integration point.

For face element like triangle or quad

Tensor1<double*,2> diff_base = data.getFTensor1DiffN<2>(base,gg,bb);

Returns: Tensor rank 1 (vector)

Definition at line 590 of file EntitiesFieldData.cpp.

                                                                         {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor1<FTensor::PackPtr<double *, Tensor_Dim>, Tensor_Dim>();
}

◆ getFTensor1DiffN() [6/8]

FTensor::Tensor1< FTensor::PackPtr< double *, 2 >, 2 > MoFEM::EntitiesFieldData::EntData::getFTensor1DiffN< 2 >	(	const FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Get spatial derivative of base function tensor for dimension 3d.

Get spatial derivative of base function tensor for dimension 2d.

Definition at line 590 of file EntitiesFieldData.cpp.

                                                                   {
  double *ptr = &getDiffN(base)(gg, 3 * bb);
  return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>(ptr, &ptr[1],
                                                            &ptr[2]);
}

◆ getFTensor1DiffN() [7/8]

template<int Tensor_Dim>

FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor1DiffN	(	const int	gg,
		const int	bb
	)

Get derivatives of base functions (Loop by integration points)

For volume element like tetrahedral or prism,

Tensor1<double*,3> diff_base = data.getFTensor1DiffN<3>(gg,bb);

where bb is base function and gg is integration pt. Operator ++diff_base will move tensor pointer to next base function.

For face element like triangle or quad

Tensor1<double*,2> diff_base = data.getFTensor1DiffN<2>(gg,bb);

Returns: Tensor rank 1 (vector)

◆ getFTensor1DiffN() [8/8]

FTensor::Tensor1< FTensor::PackPtr< double *, 2 >, 2 > MoFEM::EntitiesFieldData::EntData::getFTensor1DiffN< 2 >	(	const int	gg,
		const int	bb
	)

Get spatial derivative of base function tensor for dimension 3d.

Get spatial derivative of base function tensor for dimension 2d.

Definition at line 590 of file EntitiesFieldData.cpp.

                                                                        {
  return getFTensor1DiffN<3>(bAse, gg, bb);
}

◆ getFTensor1FieldData() [1/2]

template<int Tensor_Dim>

FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor1FieldData ( )

Return FTensor of rank 1, i.e. vector from field data coefficients.

auto t_vec = data.getFTensor1FieldData<3>();

Template Parameters

Tensor_Dim size of vector

Returns: FTensor::Tensor1<FTensor::PackPtr<double *, Tensor_Dim>, Tensor_Dim>

Definition at line 1488 of file EntitiesFieldData.hpp.

                                               {
  std::stringstream s;
  s << "Not implemented for this dimension dim = " << Tensor_Dim;
  THROW_MESSAGE(s.str());
}

◆ getFTensor1FieldData() [2/2]

FTensor::Tensor1< FTensor::PackPtr< double *, 1 >, 1 > MoFEM::EntitiesFieldData::EntData::getFTensor1FieldData< 1 > ( )

Definition at line 277 of file EntitiesFieldData.cpp.

                                                  {
#ifndef NDEBUG
  for (auto &d : dOfs) {
    if (d) {
      if (d->getNbOfCoeffs() != 3) {
        std::stringstream s;
        s << "Wrong number of coefficients is " << d->getNbOfCoeffs();
        s << " but you ask for tensor rank 1 dimension 3";
        CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY, s.str());
      }
      break;
    }
  }
#endif
  double *ptr = &*fieldData.data().begin();
  return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>(ptr, &ptr[1],
                                                            &ptr[2]);
}

◆ getFTensor1N() [1/6]

template<int Tensor_Dim>

auto MoFEM::EntitiesFieldData::EntData::getFTensor1N ( )

Get base functions for Hdiv space.

Example:

FTensor::Index<'i',3> i;
int nb_dofs = data.getFieldData().size();
auto t_n_hdiv = data.getFTensor1N<3>();
for(int gg = 0;gg!=nb_gauss_pts;gg++) {
  int ll = 0;
  for(;ll!=nb_dofs;ll++) {
    double dot_product = t_n_hdiv(i)*t_n_hdiv(i);
    ++t_n_hdiv;
  }
  for(;ll!=data.getVectorN().size2()/3;ll++) {
    ++t_n_hdiv;
  }
}

Definition at line 1746 of file EntitiesFieldData.hpp.

                                                                      {
  return getFTensor1N<Tensor_Dim>(bAse);
}

◆ getFTensor1N() [2/6]

template<int Tensor_Dim>

auto MoFEM::EntitiesFieldData::EntData::getFTensor1N	(	const int	gg,
		const int	bb
	)

inline

Get Hdiv base functions at integration point.

FTensor::Index<'i',3> i;
for(int gg = 0;gg!=nb_gauss_pts;gg++) {
 auto t_base = data.getFTensor1N(gg,0);
 for(int bb = 0;bb!=nb_base_functions;bb++) {
   double dot = t_base(i)*t_base(i);
 }
}

Definition at line 1751 of file EntitiesFieldData.hpp.

                                                                      {
  return getFTensor1N<Tensor_Dim>(bAse, gg, bb);
}

◆ getFTensor1N() [3/6]

template<int Tensor_Dim>

FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor1N ( FieldApproximationBase base )

Get base functions for Hdiv/Hcurl spaces.

Note: You probably like to use getFTensor1N(), in typical use base is set automatically based on base set to field.

Parameters

base	Approximation base

Example:

FTensor::Index<'i',3> i;
int nb_dofs = data.getFieldData().size();
auto t_n_hdiv = data.getFTensor1N<3>();
for(int gg = 0;gg!=nb_gauss_pts;gg++) {
  int ll = 0;
  for(;ll!=nb_dofs;ll++) {
    double dot_product = t_n_hdiv(i)*t_n_hdiv(i);
    ++t_n_hdiv;
  }
  for(;ll!=data.getVectorN().size2()/3;ll++) {
    ++t_n_hdiv;
  }
}

Examples: mofem/atom_tests/continuity_check_on_skeleton_3d.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hcurl_divergence_operator_2d.cpp, mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp, thermo_elastic.cpp, and thermoplastic.cpp.

Definition at line 647 of file EntitiesFieldData.cpp.

                                                                  {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor1<FTensor::PackPtr<double *, Tensor_Dim>, Tensor_Dim>();
}

◆ getFTensor1N() [4/6]

FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor1N< 3 > ( FieldApproximationBase base )

Definition at line 678 of file EntitiesFieldData.cpp.

                                                                     {
  double *t_n_ptr = &*getN(base).data().begin();
  return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>(t_n_ptr, // HVEC0
                                                            &t_n_ptr[HVEC1],
                                                            &t_n_ptr[HVEC2]);
}

◆ getFTensor1N() [5/6]

template<int Tensor_Dim>

FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor1N	(	FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Get Hdiv base functions at integration point.

FTensor::Index<'i',3> i;
for(int gg = 0;gg!=nb_gauss_pts;gg++) {
 auto t_base = data.getFTensor1N(base,gg,bb);
 for(int bb = 0;bb!=nb_base_functions;bb++) {
   auto dot = t_base(i)*t_base(i);
 }
}

Definition at line 656 of file EntitiesFieldData.cpp.

                                                                     {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor1<FTensor::PackPtr<double *, Tensor_Dim>, Tensor_Dim>();
}

◆ getFTensor1N() [6/6]

FTensor::Tensor1< FTensor::PackPtr< double *, 3 >, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor1N< 3 >	(	FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Definition at line 678 of file EntitiesFieldData.cpp.

                                                                        {
  double *t_n_ptr = &getN(base)(gg, 3 * bb);
  return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>(t_n_ptr, // HVEC0
                                                            &t_n_ptr[HVEC1],
                                                            &t_n_ptr[HVEC2]);
}

◆ getFTensor2DiffN() [1/6]

template<int Tensor_Dim0, int Tensor_Dim1>

FTensor::Tensor2< FTensor::PackPtr< double *, Tensor_Dim0 *Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN ( )

inline

Get derivatives of base functions for Hdiv space.

Definition at line 913 of file EntitiesFieldData.hpp.

                     {
    return getFTensor2DiffN<Tensor_Dim0, Tensor_Dim1>(bAse);
  }

◆ getFTensor2DiffN() [2/6]

template<int Tensor_Dim0, int Tensor_Dim1>

FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN	(	const int	gg,
		const int	bb
	)

inline

Get derivatives of base functions for Hdiv space at integration pts.

Definition at line 923 of file EntitiesFieldData.hpp.

                                               {
    return getFTensor2DiffN<Tensor_Dim0, Tensor_Dim1>(bAse, gg, bb);
  }

◆ getFTensor2DiffN() [3/6]

template<int Tensor_Dim0, int Tensor_Dim1>

FTensor::Tensor2< FTensor::PackPtr< double *, Tensor_Dim0 *Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN ( FieldApproximationBase base )

Get derivatives of base functions for Hdiv space.

Examples: mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hcurl_divergence_operator_2d.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, and mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 667 of file EntitiesFieldData.cpp.

                                                                      {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim0 << "x" << Tensor_Dim1
    << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor2<double *, Tensor_Dim0, Tensor_Dim1>();
}

◆ getFTensor2DiffN() [4/6]

FTensor::Tensor2< FTensor::PackPtr< double *, 6 >, 3, 2 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN< 3, 2 > ( FieldApproximationBase base )

Definition at line 678 of file EntitiesFieldData.cpp.

                                 {
  double *t_diff_n_ptr = &*getDiffN(base).data().begin();
  return FTensor::Tensor2<FTensor::PackPtr<double *, 9>, 3, 3>(
      t_diff_n_ptr, &t_diff_n_ptr[HVEC0_1], &t_diff_n_ptr[HVEC0_2],
      &t_diff_n_ptr[HVEC1_0], &t_diff_n_ptr[HVEC1_1], &t_diff_n_ptr[HVEC1_2],
      &t_diff_n_ptr[HVEC2_0], &t_diff_n_ptr[HVEC2_1], &t_diff_n_ptr[HVEC2_2]);
}

◆ getFTensor2DiffN() [5/6]

template<int Tensor_Dim0, int Tensor_Dim1>

FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN	(	FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Get derivatives of base functions for Hdiv space at integration pts.

Definition at line 678 of file EntitiesFieldData.cpp.

                                                                         {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim0 << "x" << Tensor_Dim1
    << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor2<double *, Tensor_Dim0, Tensor_Dim1>();
}

◆ getFTensor2DiffN() [6/6]

FTensor::Tensor2< FTensor::PackPtr< double *, 6 >, 3, 2 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN< 3, 2 >	(	FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Definition at line 678 of file EntitiesFieldData.cpp.

                                                                               {
  double *t_diff_n_ptr = &getDiffN(base)(gg, 9 * bb);
  return FTensor::Tensor2<FTensor::PackPtr<double *, 9>, 3, 3>(
      t_diff_n_ptr, &t_diff_n_ptr[HVEC0_1], &t_diff_n_ptr[HVEC0_2],
      &t_diff_n_ptr[HVEC1_0], &t_diff_n_ptr[HVEC1_1], &t_diff_n_ptr[HVEC1_2],
      &t_diff_n_ptr[HVEC2_0], &t_diff_n_ptr[HVEC2_1], &t_diff_n_ptr[HVEC2_2]);
}

◆ getFTensor2DiffN2() [1/8]

template<int Tensor_Dim>

FTensor::Tensor2< FTensor::PackPtr< double *, Tensor_Dim *Tensor_Dim >, Tensor_Dim, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN2 ( )

Get second derivatives of scalar base functions.

◆ getFTensor2DiffN2() [2/8]

FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN2< 3 > ( )

Definition at line 750 of file EntitiesFieldData.cpp.

                                               {
  return getFTensor2DiffN2<2>(bAse);
}

◆ getFTensor2DiffN2() [3/8]

template<int Tensor_Dim>

FTensor::Tensor2< FTensor::PackPtr< double *, Tensor_Dim *Tensor_Dim >, Tensor_Dim, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN2 ( const FieldApproximationBase base )

Get second derivatives of scalar base functions.

◆ getFTensor2DiffN2() [4/8]

FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN2< 3 > ( const FieldApproximationBase base )

Definition at line 750 of file EntitiesFieldData.cpp.

                                                                          {
  double *ptr = &(getN(base, BaseDerivatives::SecondDerivative))(0, 0);
  return FTensor::Tensor2<FTensor::PackPtr<double *, 4>, 2, 2>(
      ptr, &ptr[1], &ptr[2], &ptr[3]);
}

◆ getFTensor2DiffN2() [5/8]

template<int Tensor_Dim>

FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim Tensor_Dim >, Tensor_Dim, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN2	(	const FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Get second derivatives of scalar base functions at integration pts.

◆ getFTensor2DiffN2() [6/8]

FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN2< 3 >	(	const FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Definition at line 750 of file EntitiesFieldData.cpp.

                                                                             {
  double *ptr = &getN(base, BaseDerivatives::SecondDerivative)(gg, 4 * bb);
  return FTensor::Tensor2<FTensor::PackPtr<double *, 4>, 2, 2>(
      ptr, &ptr[1], &ptr[2], &ptr[3]);
}

◆ getFTensor2DiffN2() [7/8]

template<int Tensor_Dim>

FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim Tensor_Dim >, Tensor_Dim, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN2	(	const int	gg,
		const int	bb
	)

Get second derivatives of scalar base functions at integration pts.

◆ getFTensor2DiffN2() [8/8]

FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor2DiffN2< 3 >	(	const int	gg,
		const int	bb
	)

Definition at line 750 of file EntitiesFieldData.cpp.

                                                                         {
  return getFTensor2DiffN2<2>(bAse, gg, bb);
}

◆ getFTensor2FieldData() [1/2]

template<int Tensor_Dim0, int Tensor_Dim1>

FTensor::Tensor2< FTensor::PackPtr< double *, Tensor_Dim0 *Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 > MoFEM::EntitiesFieldData::EntData::getFTensor2FieldData ( )

Return FTensor rank 2, i.e. matrix from field data coefficients.

auto t_mat = data.getFTensor2FieldData<3,3>();

Template Parameters

Tensor_Dim0	First dimension of tensor
Tensor_Dim1	Second dimension of tensor

Returns: FTensor::Tensor2<FTensor::PackPtr<double *, Tensor_Dim0 * Tensor_Dim1>, Tensor_Dim0, Tensor_Dim1> Rank 2 tensor

Definition at line 1497 of file EntitiesFieldData.hpp.

                                               {
  std::stringstream s;
  s << "Not implemented for this dimension dim0 = " << Tensor_Dim0;
  s << " and dim1 " << Tensor_Dim1;
  THROW_MESSAGE(s.str());
}

◆ getFTensor2FieldData() [2/2]

FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor2FieldData< 3, 3 > ( )

Definition at line 277 of file EntitiesFieldData.cpp.

                                                     {
#ifndef NDEBUG
  for (auto &d : dOfs) {
    if (d) {
      if (d->getNbOfCoeffs() != 1) {
        std::stringstream s;
        s << "Wrong number of coefficients is " << d->getNbOfCoeffs();
        s << " but you ask for tensor rank 2 dimensions 1 by 1 so 1 "
             "coefficients "
             "is expected";
        CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY, s.str());
      }
      break;
    }
  }
#endif
  double *ptr = &*fieldData.data().begin();
  return FTensor::Tensor2<FTensor::PackPtr<double *, 1>, 1, 1>(ptr);
}

◆ getFTensor2N() [1/6]

template<int Tensor_Dim0, int Tensor_Dim1>

auto MoFEM::EntitiesFieldData::EntData::getFTensor2N ( )

Get base functions for Hdiv space.

Example:

FTensor::Index<'i',3> i;
FTensor::Index<'j',3> j;
 
int nb_dofs = data.getFieldData().size();
auto t_n_hdiv = data.getFTensor2N<3,3>();
for(int gg = 0;gg!=nb_gauss_pts;gg++) {
  int ll = 0;
  for(;ll!=nb_dofs;ll++) {
    double dot_product = t_n_hdiv(i,j)*t_n_hdiv(i,j);
    ++t_n_hdiv;
  }
  for(;ll!=data.getVectorN().size2()/3;ll++) {
    ++t_n_hdiv;
  }
}

Definition at line 1756 of file EntitiesFieldData.hpp.

                                            {
  return getFTensor2N<Tensor_Dim0, Tensor_Dim1>(bAse);
}

◆ getFTensor2N() [2/6]

template<int Tensor_Dim0, int Tensor_Dim1>

auto MoFEM::EntitiesFieldData::EntData::getFTensor2N	(	const int	gg,
		const int	bb
	)

Get base functions for Hdiv space.

Example:

FTensor::Index<'i',3> i;
FTensor::Index<'j',3> j;
int nb_dofs = data.getFieldData().size();
for(int gg = 0;gg!=nb_gauss_pts;++gg) {
  int ll = 0;
  auto t_n_hdiv = data.getFTensor2N<3,3>(gg,0);
  for(;ll!=nb_dofs;ll++) {
    double dot_product = t_n_hdiv(i)*t_n_hdiv(i);
    ++t_n_hdiv;
  }
  for(;ll!=data.getVectorN().size2()/3;ll++) {
    ++t_n_hdiv;
  }
}

Definition at line 1761 of file EntitiesFieldData.hpp.

                                                                      {
  return getFTensor2N<Tensor_Dim0, Tensor_Dim1>(bAse, gg, bb);
}

◆ getFTensor2N() [3/6]

template<int Tensor_Dim0, int Tensor_Dim1>

FTensor::Tensor2< FTensor::PackPtr< double *, Tensor_Dim0 *Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 > MoFEM::EntitiesFieldData::EntData::getFTensor2N ( FieldApproximationBase base )

Get base functions for Hdiv/Hcurl spaces.

Note: You probably like to use getFTensor1N(), in typical use base is set automatically based on base set to field.

Parameters

base	Approximation base

Example:

FTensor::Index<'i',3> i;
FTensor::Index<'i',3> j;
int nb_dofs = data.getFieldData().size();
auto t_n_hdiv = data.getFTensor2N<3,3>();
for(int gg = 0;gg!=nb_gauss_pts;gg++) {
  int ll = 0;
  for(;ll!=nb_dofs;ll++) {
    double dot_product = t_n_hdiv(i,j)*t_n_hdiv(i,j);
    ++t_n_hdiv;
  }
  for(;ll!=data.getVectorN().size2()/3;ll++) {
    ++t_n_hdiv;
  }
}

Examples: mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 829 of file EntitiesFieldData.cpp.

                                                                  {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim0 << ", " << Tensor_Dim1
    << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor2<FTensor::PackPtr<double *, Tensor_Dim0 * Tensor_Dim1>,
                          Tensor_Dim0, Tensor_Dim1>();
}

◆ getFTensor2N() [4/6]

FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor2N< 3, 3 > ( FieldApproximationBase base )

Definition at line 829 of file EntitiesFieldData.cpp.

                                                                        {
  double *t_n_ptr = &*(getN(base).data().begin());
  return FTensor::Tensor2<FTensor::PackPtr<double *, 9>, 3, 3>(
 
      &t_n_ptr[HVEC0], &t_n_ptr[HVEC1], &t_n_ptr[HVEC2],
 
      &t_n_ptr[3 + HVEC0], &t_n_ptr[3 + HVEC1], &t_n_ptr[3 + HVEC2],
 
      &t_n_ptr[6 + HVEC0], &t_n_ptr[6 + HVEC1], &t_n_ptr[6 + HVEC2]
 
  );
}

◆ getFTensor2N() [5/6]

template<int Tensor_Dim0, int Tensor_Dim1>

FTensor::Tensor2< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 >, Tensor_Dim0, Tensor_Dim1 > MoFEM::EntitiesFieldData::EntData::getFTensor2N	(	FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Get base functions for tensor Hdiv/Hcurl spaces.

Note: You probably like to use getFTensor2N(), in typical use base is set automatically based on base set to field.

Parameters

base	Approximation base

Example:

FTensor::Index<'i',3> i;
FTensor::Index<'j',3> i;
int nb_dofs = data.getFieldData().size();
for(int gg = 0;gg!=nb_gauss_pts;gg++) {
  auto t_n_hdiv = data.getFTensor2N<3>(base,gg,bb);
  int ll = 0;
  for(;ll!=nb_dofs;ll++) {
    double dot_product = t_n_hdiv(i,j)*t_n_hdiv(i,j);
    ++t_n_hdiv;
  }
  for(;ll!=data.getVectorN().size2()/3;ll++) {
    ++t_n_hdiv;
  }
}

Definition at line 856 of file EntitiesFieldData.cpp.

                                                                     {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim0 << ", " << Tensor_Dim1
    << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor2<FTensor::PackPtr<double *, Tensor_Dim0 * Tensor_Dim1>,
                          Tensor_Dim0, Tensor_Dim1>();
}

◆ getFTensor2N() [6/6]

template<>

FTensor::Tensor2< FTensor::PackPtr< double *, 9 >, 3, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor2N	(	FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Definition at line 856 of file EntitiesFieldData.cpp.

                                                                           {
  double *t_n_ptr = &getN(base)(gg, 9 * bb);
  return FTensor::Tensor2<FTensor::PackPtr<double *, 9>, 3, 3>(
 
      &t_n_ptr[HVEC0], &t_n_ptr[HVEC1], &t_n_ptr[HVEC2],
 
      &t_n_ptr[3 + HVEC0], &t_n_ptr[3 + HVEC1], &t_n_ptr[3 + HVEC2],
 
      &t_n_ptr[6 + HVEC0], &t_n_ptr[6 + HVEC1], &t_n_ptr[6 + HVEC2]
 
  );
}

◆ getFTensor2SymmetricFieldData() [1/2]

template<int Tensor_Dim>

FTensor::Tensor2_symmetric< FTensor::PackPtr< double *,(Tensor_Dim *(Tensor_Dim+1))/2 >, Tensor_Dim > MoFEM::EntitiesFieldData::EntData::getFTensor2SymmetricFieldData ( )

Return symmetric FTensor rank 2, i.e. matrix from field data coefficients.

auto t_mat = data.getFTensor2SymmetricFieldData<3>();

Template Parameters

Tensor_Dim Dimension of the symmetric tensor

Returns

FTensor::Tensor2_symmetric<FTensor::PackPtr<double *, (Tensor_Dim

(Tensor_Dim + 1)) / 2>, Tensor_Dim> Symmetric rank 2 tensor

Definition at line 1507 of file EntitiesFieldData.hpp.

                                                        {
  std::stringstream s;
  s << "Not implemented for this dimension dim = " << Tensor_Dim;
  THROW_MESSAGE(s.str());
}

◆ getFTensor2SymmetricFieldData() [2/2]

FTensor::Tensor2_symmetric< FTensor::PackPtr< double *, 3 >, 2 > MoFEM::EntitiesFieldData::EntData::getFTensor2SymmetricFieldData< 2 > ( )

Definition at line 277 of file EntitiesFieldData.cpp.

                                                           {
#ifndef NDEBUG
  for (auto &d : dOfs) {
    if (d) {
      if (d->getNbOfCoeffs() != 6) {
        std::stringstream s;
        s << "Wrong number of coefficients is " << d->getNbOfCoeffs();
        s << " but you ask for symmetric tensor rank 2 dimensions 3 by 3 so 6 "
             "coefficients "
             "is expected";
        CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY, s.str());
      }
      break;
    }
  }
#endif
  double *ptr = &*fieldData.data().begin();
  return FTensor::Tensor2_symmetric<FTensor::PackPtr<double *, 6>, 3>(
      ptr, &ptr[1], &ptr[2], &ptr[3], &ptr[4], &ptr[5]);
}

◆ getFTensor3Diff2N() [1/3]

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>

FTensor::Tensor3< FTensor::PackPtr< double *, Tensor_Dim0 *Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 > MoFEM::EntitiesFieldData::EntData::getFTensor3Diff2N ( )

inline

Get second derivatives of base functions for Hvec space.

Definition at line 941 of file EntitiesFieldData.hpp.

                      {
    return getFTensor3Diff2N<Tensor_Dim0, Tensor_Dim1, Tensor_Dim2>(bAse);
  }

◆ getFTensor3Diff2N() [2/3]

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>

FTensor::Tensor3< FTensor::PackPtr< double *, Tensor_Dim0 *Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 > MoFEM::EntitiesFieldData::EntData::getFTensor3Diff2N ( FieldApproximationBase base )

Get second derivatives of base functions for Hvec space.

◆ getFTensor3Diff2N() [3/3]

FTensor::Tensor3< FTensor::PackPtr< double *, 12 >, 3, 2, 2 > MoFEM::EntitiesFieldData::EntData::getFTensor3Diff2N ( FieldApproximationBase base )

Definition at line 750 of file EntitiesFieldData.cpp.

                                                                       {
  double *ptr = &(getN(base, BaseDerivatives::SecondDerivative))(0, 0);
  return FTensor::Tensor3<FTensor::PackPtr<double *, 12>, 3, 2, 2>{
 
      &ptr[2 * HVEC0_0 + 0], &ptr[2 * HVEC0_0 + 1], &ptr[2 * HVEC0_1 + 0],
      &ptr[2 * HVEC0_1 + 1],
 
      &ptr[2 * HVEC1_0 + 0], &ptr[2 * HVEC1_0 + 1], &ptr[2 * HVEC1_1 + 0],
      &ptr[2 * HVEC1_1 + 1],
 
      &ptr[2 * HVEC2_0 + 0], &ptr[2 * HVEC2_0 + 1], &ptr[2 * HVEC2_1 + 0],
      &ptr[2 * HVEC2_1 + 1]
 
  };
}

◆ getFTensor3DiffN() [1/6]

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>

FTensor::Tensor3< FTensor::PackPtr< double *, Tensor_Dim0 *Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 > MoFEM::EntitiesFieldData::EntData::getFTensor3DiffN ( )

inline

Get derivatives of base functions for tonsorial Hdiv space.

Definition at line 968 of file EntitiesFieldData.hpp.

                     {
    return getFTensor3DiffN<Tensor_Dim0, Tensor_Dim1, Tensor_Dim2>(bAse);
  }

◆ getFTensor3DiffN() [2/6]

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>

FTensor::Tensor3< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 > MoFEM::EntitiesFieldData::EntData::getFTensor3DiffN	(	const int	gg,
		const int	bb
	)

inline

Get derivatives of base functions for tonsorial Hdiv space at integration pts.

Definition at line 979 of file EntitiesFieldData.hpp.

                                               {
    return getFTensor3DiffN<Tensor_Dim0, Tensor_Dim1, Tensor_Dim2>(bAse, gg,
                                                                   bb);
  }

◆ getFTensor3DiffN() [3/6]

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>

FTensor::Tensor3< FTensor::PackPtr< double *, Tensor_Dim0 *Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 > MoFEM::EntitiesFieldData::EntData::getFTensor3DiffN ( FieldApproximationBase base )

Get derivatives of base functions for tonsorial Hdiv space.

Definition at line 886 of file EntitiesFieldData.cpp.

                                                                      {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim0 << "x" << Tensor_Dim1
    << "x" << Tensor_Dim2 << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor3<double *, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2>();
}

◆ getFTensor3DiffN() [4/6]

FTensor::Tensor3< FTensor::PackPtr< double *, 27 >, 3, 3, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor3DiffN ( FieldApproximationBase base )

Definition at line 908 of file EntitiesFieldData.cpp.

                                                                      {
  double *ptr = &*getDiffN(base).data().begin();
  return getFTensor3FromPtr<3, 3, 3>(ptr);
};

◆ getFTensor3DiffN() [5/6]

template<int Tensor_Dim0, int Tensor_Dim1, int Tensor_Dim2>

FTensor::Tensor3< FTensor::PackPtr< double , Tensor_Dim0 Tensor_Dim1 *Tensor_Dim2 >, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2 > MoFEM::EntitiesFieldData::EntData::getFTensor3DiffN	(	FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Get derivatives of base functions for tonsorial Hdiv space at integration pts.

Definition at line 897 of file EntitiesFieldData.cpp.

                                                                         {
  std::stringstream s;
  s << "Template for tensor dimension " << Tensor_Dim0 << "x" << Tensor_Dim1
    << " not implemented";
  THROW_MESSAGE(s.str());
  return FTensor::Tensor3<double *, Tensor_Dim0, Tensor_Dim1, Tensor_Dim2>();
}

◆ getFTensor3DiffN() [6/6]

FTensor::Tensor3< FTensor::PackPtr< double *, 27 >, 3, 3, 3 > MoFEM::EntitiesFieldData::EntData::getFTensor3DiffN	(	FieldApproximationBase	base,
		const int	gg,
		const int	bb
	)

Definition at line 915 of file EntitiesFieldData.cpp.

                                                                         {
  double *ptr = &getDiffN(base)(gg, 27 * bb);
  return getFTensor3FromPtr<3, 3, 3>(ptr);
};

◆ getIndices() [1/2]

VectorInt & MoFEM::EntitiesFieldData::EntData::getIndices ( )

inline

Get reference to global indices for modification.

Returns: VectorInt& Reference to global indices vector

Definition at line 1448 of file EntitiesFieldData.hpp.

1448{ return iNdices; }

MoFEM::EntitiesFieldData::EntData::iNdices

VectorInt iNdices

Global indices on entity.

Definition EntitiesFieldData.hpp:1273

◆ getIndices() [2/2]

const VectorInt & MoFEM::EntitiesFieldData::EntData::getIndices ( ) const

inline

Get global indices of degrees of freedom on entity.

Returns: const VectorInt& Reference to vector of global DOF indices

Examples: mofem/atom_tests/boundary_marker.cpp, mofem/atom_tests/child_and_parent.cpp, mofem/atom_tests/forces_and_sources_testing_users_base.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_divergence_operator_2d.cpp, mofem/atom_tests/prism_elements_from_surface.cpp, mofem/atom_tests/quad_polynomial_approximation.cpp, mofem/atom_tests/scalar_check_approximation.cpp, mofem/atom_tests/simple_interface.cpp, mofem/atom_tests/test_cache_on_entities.cpp, mofem/tutorials/adv-3_level_set/level_set.cpp, mofem/tutorials/adv-6_dg_projection/between_meshes_dg_projection.cpp, mofem/tutorials/cor-2to5_analytical_poisson/src/PoissonOperators.hpp, mofem/tutorials/cor-6_simple_elasticity/simple_elasticity.cpp, mofem/tutorials/cor-9_reaction_diffusion/reaction_diffusion.cpp, mofem/tutorials/fun-0_hello_world/hello_world.cpp, mofem/tutorials/scl-0_least_squares/approximaton.cpp, mofem/tutorials/scl-1_poisson_2d_homogeneous/src/poisson_2d_homogeneous.hpp, mofem/tutorials/scl-9_heat_method/heat_method.cpp, mofem/tutorials/vec-6_dg_kirchhoff_plate/plate.cpp, mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp, mofem/users_modules/basic_finite_elements/src/impl/HookeElement.cpp, mofem/users_modules/basic_finite_elements/src/impl/NavierStokesElement.cpp, and mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 1414 of file EntitiesFieldData.hpp.

                                                            {
  return iNdices;
}

◆ getIndicesUpToOrder()

const VectorIntAdaptor MoFEM::EntitiesFieldData::EntData::getIndicesUpToOrder ( int order )

inline

Get global indices of DOFs on entity up to given order.

Parameters

order Maximum approximation order to include

Returns: const VectorIntAdaptor Adaptor to subset of global indices

Definition at line 1419 of file EntitiesFieldData.hpp.

                                                       {
  unsigned int size = 0;
  if (auto dof = dOfs[0]) {
    size = dof->getOrderNbDofs(order) * dof->getNbOfCoeffs();
    size = size < iNdices.size() ? size : iNdices.size();
  }
  int *data = &*iNdices.data().begin();
  return VectorIntAdaptor(size, ublas::shallow_array_adaptor<int>(size, data));
}

◆ getLocalIndices() [1/2]

VectorInt & MoFEM::EntitiesFieldData::EntData::getLocalIndices ( )

inline

Get reference to local indices for modification.

Returns: VectorInt& Reference to local indices vector

Definition at line 1450 of file EntitiesFieldData.hpp.

                                                     {
  return localIndices;
}

◆ getLocalIndices() [2/2]

const VectorInt & MoFEM::EntitiesFieldData::EntData::getLocalIndices ( ) const

inline

Get local indices of degrees of freedom on entity.

Returns: const VectorInt& Reference to vector of local DOF indices

Examples: mofem/atom_tests/boundary_marker.cpp, and mofem/atom_tests/prism_elements_from_surface.cpp.

Definition at line 1429 of file EntitiesFieldData.hpp.

                                                                 {
  return localIndices;
}

◆ getLocalIndicesUpToOrder()

const VectorIntAdaptor MoFEM::EntitiesFieldData::EntData::getLocalIndicesUpToOrder ( int order )

inline

Get local indices of DOFs on entity up to given order.

Parameters

order Maximum approximation order to include

Returns: const VectorIntAdaptor Adaptor to subset of local indices

Definition at line 1434 of file EntitiesFieldData.hpp.

                                                            {
  unsigned int size = 0;
  if (auto dof = dOfs[0]) {
    size = dof->getOrderNbDofs(order) * dof->getNbOfCoeffs();
    size = size < localIndices.size() ? size : localIndices.size();
  }
  int *data = &*localIndices.data().begin();
  return VectorIntAdaptor(size, ublas::shallow_array_adaptor<int>(size, data));
}

◆ getN() [1/9]

MatrixDouble & MoFEM::EntitiesFieldData::EntData::getN ( )

inline

Get base functions using default base.

Returns: MatrixDouble& Reference to base function matrix

Definition at line 1542 of file EntitiesFieldData.hpp.

1542{ return getN(bAse); }

◆ getN() [2/9]

MatrixDouble & MoFEM::EntitiesFieldData::EntData::getN ( const BaseDerivatives derivative )

inline

Get base function derivatives using default base.

Parameters

derivative Derivative order

Returns: MatrixDouble& Reference to derivative matrix

Definition at line 1577 of file EntitiesFieldData.hpp.

                                                               {
  return getN(bAse, derivative);
}

◆ getN() [3/9]

MatrixDouble & MoFEM::EntitiesFieldData::EntData::getN ( const FieldApproximationBase base )

inline

get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb. of columns is equal to number of base functions on this entity.

Note: Note that for vectorial base, like Hdiv or Hcurl, in columns are vectorial base functions. For tonsorial would be base functions. Interpretation depends on type of base, scalar, vectorial or tonsorial and dimension of problem.

Examples: mofem/atom_tests/continuity_check_on_contact_prism_side_ele.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_h1.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hcurl.cpp, mofem/atom_tests/continuity_check_on_skeleton_with_simple_2d_for_hdiv.cpp, mofem/atom_tests/forces_and_sources_testing_edge_element.cpp, mofem/atom_tests/forces_and_sources_testing_users_base.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_check_approx_in_2d.cpp, mofem/atom_tests/hcurl_curl_operator.cpp, mofem/atom_tests/hcurl_divergence_operator_2d.cpp, mofem/atom_tests/hdiv_divergence_operator.cpp, mofem/atom_tests/prism_elements_from_surface.cpp, mofem/atom_tests/prism_polynomial_approximation.cpp, mofem/atom_tests/quad_polynomial_approximation.cpp, mofem/atom_tests/scalar_check_approximation.cpp, mofem/tutorials/adv-3_level_set/level_set.cpp, mofem/tutorials/cor-6_simple_elasticity/simple_elasticity.cpp, mofem/tutorials/scl-9_heat_method/heat_method.cpp, mofem/tutorials/vec-6_dg_kirchhoff_plate/plate.cpp, mofem/tutorials/vec-7_shape_optimisation/adjoint.cpp, mofem/users_modules/basic_finite_elements/src/HookeElement.hpp, mofem/users_modules/basic_finite_elements/src/impl/HookeElement.cpp, mofem/users_modules/basic_finite_elements/src/impl/NavierStokesElement.cpp, and mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 1518 of file EntitiesFieldData.hpp.

                                                                {
#ifndef NDEBUG
  if (!getNSharedPtr(base)) {
    MOFEM_LOG_C("SELF", Sev::error, "Ptr to base %s functions is null",
                ApproximationBaseNames[base]);
    CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                      "Null pointer to base functions");
  }
#endif // NDEBUG  
  return *(getNSharedPtr(base));
}

◆ getN() [4/9]

MatrixDouble & MoFEM::EntitiesFieldData::EntData::getN	(	const FieldApproximationBase	base,
		const BaseDerivatives	derivative
	)

inline

Get base function derivatives with specific derivative order.

Parameters

base	Approximation base type
derivative	Derivative order

Returns: MatrixDouble& Reference to derivative matrix

Definition at line 1550 of file EntitiesFieldData.hpp.

                                                                   {
#ifndef NDEBUG
  if (!getNSharedPtr(base, derivative)) {
    MOFEM_LOG_C("SELF", Sev::error,
                "Ptr to base %s functions derivative %d is null",
                ApproximationBaseNames[base], derivative);
    THROW_MESSAGE("Null pointer");
  }
#endif
  return *(getNSharedPtr(base, derivative));
}

◆ getN() [5/9]

const VectorAdaptor MoFEM::EntitiesFieldData::EntData::getN	(	const FieldApproximationBase	base,
		const int	gg
	)

inline

Get base functions at specific Gauss point.

Parameters

base	Approximation base type
gg	Gauss point number

Returns: const VectorAdaptor Adaptor to base functions at Gauss point

Definition at line 1582 of file EntitiesFieldData.hpp.

                                               {
  int size = getN(base).size2();
  double *data = &getN(base)(gg, 0);
  return VectorAdaptor(size, ublas::shallow_array_adaptor<double>(size, data));
}

◆ getN() [6/9]

const VectorAdaptor MoFEM::EntitiesFieldData::EntData::getN	(	const FieldApproximationBase	base,
		const int	gg,
		const int	nb_base_functions
	)

inline

Get base functions at Gauss point with limited number of functions.

Note that for multi-field elements, different fields can have different approximation orders. Since hierarchical approximation basis is used, base functions are calculated once for element using maximal order.

Parameters

base	Approximation base type
gg	Gauss point number
nb_base_functions	Number of base functions to return

Returns: const VectorAdaptor Adaptor to limited base functions

Definition at line 1621 of file EntitiesFieldData.hpp.

                                                                            {
  (void)getN()(gg, nb_base_functions -
                       1); // throw error if nb_base_functions is to big
  double *data = &getN(base)(gg, 0);
  return VectorAdaptor(nb_base_functions, ublas::shallow_array_adaptor<double>(
                                              nb_base_functions, data));
}

◆ getN() [7/9]

const VectorAdaptor MoFEM::EntitiesFieldData::EntData::getN ( const int gg )

inline

Get base functions at specific Gauss point using default base.

Parameters

gg	Gauss point number

Returns: const VectorAdaptor Adaptor to base functions at Gauss point

Definition at line 1589 of file EntitiesFieldData.hpp.

                                                               {
  return getN(bAse, gg);
}

◆ getN() [8/9]

const VectorAdaptor MoFEM::EntitiesFieldData::EntData::getN	(	const int	gg,
		const int	nb_base_functions
	)

inline

Get base functions at Gauss point with limited number (default base)

Parameters

gg	Gauss point number
nb_base_functions	Number of base functions to return

Returns: const VectorAdaptor Adaptor to limited base functions

Definition at line 1631 of file EntitiesFieldData.hpp.

                                                                        {
  return getN(bAse, gg, nb_base_functions);
}

◆ getN() [9/9]

MatrixDouble & MoFEM::EntitiesFieldData::EntData::getN ( const std::string & field_name )

inline

Get base functions for specific field (Bernstein-Bezier)

Parameters

field_name Name of the field

Returns: MatrixDouble& Reference to base function matrix

Definition at line 1530 of file EntitiesFieldData.hpp.

                                                                        {
#ifndef NDEBUG
  if (!getBBNSharedPtr(field_name)) {
    MOFEM_LOG_C("SELF", Sev::error, "Ptr to field %s base functions is null ",
                field_name.c_str());
    CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                      "Null pointer to base functions");
  }
#endif // NDEBUG
  return *(getBBNSharedPtr(field_name));
}

◆ getNSharedPtr() [1/2]

boost::shared_ptr< MatrixDouble > & MoFEM::EntitiesFieldData::EntData::getNSharedPtr ( const FieldApproximationBase base )

virtual

Get shared pointer to base functions (zero derivative)

Parameters

base	Approximation base type

Returns: boost::shared_ptr<MatrixDouble>& Shared pointer to base function matrix

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 40 of file EntitiesFieldData.cpp.

                                                                         {
  return N[base];
}

◆ getNSharedPtr() [2/2]

boost::shared_ptr< MatrixDouble > & MoFEM::EntitiesFieldData::EntData::getNSharedPtr	(	const FieldApproximationBase	base,
		const BaseDerivatives	derivative
	)

virtual

Get shared pointer to base functions with derivatives.

Parameters

base	Approximation base type
derivative	Derivative order (ZeroDerivative, FirstDerivative, etc.)

Returns: boost::shared_ptr<MatrixDouble>& Shared pointer to base function matrix

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 34 of file EntitiesFieldData.cpp.

                                                                            {
  return baseFunctionsAndBaseDerivatives[derivative][base];
}

◆ getOrder() [1/2]

ApproximationOrder & MoFEM::EntitiesFieldData::EntData::getOrder ( )

inline

Get reference to approximation order for modification.

Returns: ApproximationOrder& Reference to approximation order

Definition at line 1446 of file EntitiesFieldData.hpp.

1446{ return oRder; }

◆ getOrder() [2/2]

ApproximationOrder MoFEM::EntitiesFieldData::EntData::getOrder ( ) const

inline

Get approximation order.

Returns: ApproximationOrder Approximation order for this entity

Definition at line 1410 of file EntitiesFieldData.hpp.

                                                            {
  return oRder;
}

◆ getSense() [1/2]

int & MoFEM::EntitiesFieldData::EntData::getSense ( )

inline

Get reference to entity sense for modification.

Returns: int& Reference to entity sense

Definition at line 1444 of file EntitiesFieldData.hpp.

1444{ return sEnse; }

◆ getSense() [2/2]

int MoFEM::EntitiesFieldData::EntData::getSense ( ) const

virtual

Get entity sense for conforming approximation fields.

Entity sense is needed to calculate base functions with conforming approximation fields.

Returns: int Entity sense (orientation)

Reimplemented in MoFEM::DerivedEntitiesFieldData::DerivedEntData.

Definition at line 31 of file EntitiesFieldData.cpp.

31{ return sEnse; }

◆ getSpace()

FieldSpace & MoFEM::EntitiesFieldData::EntData::getSpace ( )

inline

Get field space.

Returns: Field space

Examples: mofem/atom_tests/continuity_check_on_contact_prism_side_ele.cpp, mofem/atom_tests/continuity_check_on_skeleton_3d.cpp, and mofem/atom_tests/hanging_node_approx.cpp.

Definition at line 1515 of file EntitiesFieldData.hpp.

1515{ return sPace; }

MoFEM::EntitiesFieldData::EntData::sPace

FieldSpace sPace

Entity space.

Definition EntitiesFieldData.hpp:1271

◆ getVectorDiffN() [1/4]

template<int DIM0, int DIM1>

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getVectorDiffN	(	const FieldApproximationBase	base,
		const int	dof,
		const int	gg
	)

inline

get DiffHdiv of base functions at Gauss pts

Parameters

base	Approximation base
gg	nb. of Gauss point
number	of of base functions

Definition at line 1709 of file EntitiesFieldData.hpp.

                                                                        {
  double *data =
      &EntitiesFieldData::EntData::getDiffN(base)(gg, DIM0 * DIM1 * dof);
  return MatrixAdaptor(DIM0, DIM1,
                       ublas::shallow_array_adaptor<double>(DIM0 * DIM1, data));
}

◆ getVectorDiffN() [2/4]

template<int DIM0, int DIM1>

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getVectorDiffN	(	const int	dof,
		const int	gg
	)

inline

get DiffHdiv of base functions at Gauss pts

Parameters

gg	nb. of Gauss point
number	of of base functions

Definition at line 1718 of file EntitiesFieldData.hpp.

                                                                             {
  return getVectorDiffN<DIM0, DIM1>(bAse, dof, gg);
}

◆ getVectorDiffN() [3/4]

template<int DIM0, int DIM1>

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getVectorDiffN ( const int gg )

inline

get DiffHdiv of base functions at Gauss pts

Parameters

gg	nb. of Gauss point
number	of of base functions

Definition at line 1703 of file EntitiesFieldData.hpp.

                                                                         {
  return getVectorDiffN<DIM0, DIM1>(bAse, gg);
}

◆ getVectorDiffN() [4/4]

template<int DIM0, int DIM1>

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getVectorDiffN	(	FieldApproximationBase	base,
		const int	gg
	)

inline

get DiffHdiv of base functions at Gauss pts

Parameters

base	Approximation base
gg	nb. of Gauss point
number	of of base functions

Definition at line 1689 of file EntitiesFieldData.hpp.

                                                         {
  if (PetscUnlikely(getDiffN(base).size2() % (DIM0 * DIM1))) {
    THROW_MESSAGE("Wrong dimension");
  }
 
  const int nb_base_functions = getN(base).size2() / (DIM0 * DIM1);
  double *data = &getN(base)(gg, 0);
  return MatrixAdaptor(nb_base_functions, DIM0 * DIM1,
                       ublas::shallow_array_adaptor<double>(
                           DIM0 * DIM1 * nb_base_functions, data));
}

◆ getVectorN() [1/2]

template<int DIM>

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getVectorN	(	const FieldApproximationBase	base,
		const int	gg
	)

inline

get Hdiv of base functions at Gauss pts

Parameters

base	Approximation base
gg	nb. of Gauss point

Examples: mofem/atom_tests/hdiv_divergence_operator.cpp.

Definition at line 1669 of file EntitiesFieldData.hpp.

                                                     {
  if (PetscUnlikely(getN(base).size2() % DIM)) {
    THROW_MESSAGE("Wrong dimension");
  }
 
  const int nb_base_functions = getN(base).size2() / DIM;
  double *data = &getN(base)(gg, 0);
  return MatrixAdaptor(
      nb_base_functions, DIM,
      ublas::shallow_array_adaptor<double>(DIM * nb_base_functions, data));
}

◆ getVectorN() [2/2]

template<int DIM>

const MatrixAdaptor MoFEM::EntitiesFieldData::EntData::getVectorN ( const int gg )

inline

get Hdiv of base functions at Gauss pts

Parameters

gg	nb. of Gauss point
number	of of base functions

Definition at line 1683 of file EntitiesFieldData.hpp.

                                                                     {
  return getVectorN<DIM>(bAse, gg);
}

◆ resetFieldDependentData()

MoFEMErrorCode MoFEM::EntitiesFieldData::EntData::resetFieldDependentData ( )

inline

Reset data associated with particular field name

Returns: error code

Definition at line 122 of file EntitiesFieldData.cpp.

                                                                 {
  MoFEMFunctionBeginHot;
  sPace = NOSPACE;
  bAse = NOBASE;
  fieldEntities.resize(0, false);
  iNdices.resize(0, false);
  localIndices.resize(0, false);
  dOfs.resize(0, false);
  fieldData.resize(0, false);
  MoFEMFunctionReturnHot(0);
}

Friends And Related Symbol Documentation

◆ OpAddParentEntData

friend struct OpAddParentEntData

friend

Definition at line 1321 of file EntitiesFieldData.hpp.

◆ operator<<

std::ostream & operator<<	(	std::ostream &	os,
		const EntitiesFieldData::EntData &	e
	)

friend

Definition at line 247 of file EntitiesFieldData.cpp.

                                                            {
  os << "sEnse: " << e.getSense() << std::endl
     << "oRder: " << e.getOrder() << std::endl
     << "global indices: " << e.getIndices() << std::endl
     << "local indices: " << e.getLocalIndices() << std::endl;
  // FIXME: precision should not be set here
  os << "fieldData: " << std::fixed << std::setprecision(2) << e.getFieldData()
     << std::endl;
  MatrixDouble base = const_cast<EntitiesFieldData::EntData &>(e).getN();
  MatrixDouble diff_base =
      const_cast<EntitiesFieldData::EntData &>(e).getDiffN();
  const double eps = 1e-6;
  for (unsigned int ii = 0; ii != base.size1(); ii++) {
    for (unsigned int jj = 0; jj != base.size2(); jj++) {
      if (fabs(base(ii, jj)) < eps)
        base(ii, jj) = 0;
    }
  }
  for (unsigned int ii = 0; ii != diff_base.size1(); ii++) {
    for (unsigned int jj = 0; jj != diff_base.size2(); jj++) {
      if (fabs(diff_base(ii, jj)) < eps)
        diff_base(ii, jj) = 0;
    }
  }
  os << "N: " << std::fixed << base << std::endl
     << "diffN: " << std::fixed << diff_base;
  return os;
}

◆ OpGetBrokenBaseSideData

template<typename OpBase >

friend struct OpGetBrokenBaseSideData

friend

Definition at line 1323 of file EntitiesFieldData.hpp.

Member Data Documentation

◆ bAse

FieldApproximationBase MoFEM::EntitiesFieldData::EntData::bAse

protected

Field approximation base.

Definition at line 1272 of file EntitiesFieldData.hpp.

◆ baseFunctionsAndBaseDerivatives

std::array<std::array<boost::shared_ptr<MatrixDouble>, LASTBASE>, LastDerivative> MoFEM::EntitiesFieldData::EntData::baseFunctionsAndBaseDerivatives

protected

Definition at line 1286 of file EntitiesFieldData.hpp.

◆ bbAlphaIndices

std::map<std::string, boost::shared_ptr<MatrixInt> > MoFEM::EntitiesFieldData::EntData::bbAlphaIndices

protected

Indices for Bernstein-Bezier (BB) base.

Definition at line 1297 of file EntitiesFieldData.hpp.

◆ bbAlphaIndicesByOrder

std::array<boost::shared_ptr<MatrixInt>, MaxBernsteinBezierOrder> MoFEM::EntitiesFieldData::EntData::bbAlphaIndicesByOrder

protected

BB alpha indices by order.

Definition at line 1304 of file EntitiesFieldData.hpp.

◆ bbDiffN

std::map<std::string, boost::shared_ptr<MatrixDouble> > MoFEM::EntitiesFieldData::EntData::bbDiffN

protected

Definition at line 1295 of file EntitiesFieldData.hpp.

◆ bbDiffNByOrder

std::array<boost::shared_ptr<MatrixDouble>, MaxBernsteinBezierOrder> MoFEM::EntitiesFieldData::EntData::bbDiffNByOrder

protected

BB base functions derivatives by order.

Definition at line 1302 of file EntitiesFieldData.hpp.

◆ bbFieldName

std::string MoFEM::EntitiesFieldData::EntData::bbFieldName

protected

field name

Definition at line 1292 of file EntitiesFieldData.hpp.

◆ bbN

std::map<std::string, boost::shared_ptr<MatrixDouble> > MoFEM::EntitiesFieldData::EntData::bbN

protected

Definition at line 1294 of file EntitiesFieldData.hpp.

◆ bbNByOrder

std::array<boost::shared_ptr<MatrixDouble>, MaxBernsteinBezierOrder> MoFEM::EntitiesFieldData::EntData::bbNByOrder

protected

BB base functions by order.

Definition at line 1300 of file EntitiesFieldData.hpp.

◆ bbNodeOrder

VectorInt MoFEM::EntitiesFieldData::EntData::bbNodeOrder

protected

order of nodes

Definition at line 1293 of file EntitiesFieldData.hpp.

◆ diffN

std::array<boost::shared_ptr<MatrixDouble>, LASTBASE>& MoFEM::EntitiesFieldData::EntData::diffN

protected

Derivatives of base functions.

Definition at line 1290 of file EntitiesFieldData.hpp.

◆ dofBrokenSideVec

std::vector<int> MoFEM::EntitiesFieldData::EntData::dofBrokenSideVec

protected

Map side to dofs number.

Definition at line 1280 of file EntitiesFieldData.hpp.

◆ dofBrokenTypeVec

std::vector<EntityType> MoFEM::EntitiesFieldData::EntData::dofBrokenTypeVec

protected

Map type of entity to dof number.

Definition at line 1282 of file EntitiesFieldData.hpp.

◆ dOfs

VectorDofs MoFEM::EntitiesFieldData::EntData::dOfs

protected

DoFs on entity.

Definition at line 1275 of file EntitiesFieldData.hpp.

◆ entDataBitRefLevel

std::vector<BitRefLevel> MoFEM::EntitiesFieldData::EntData::entDataBitRefLevel

protected

Bit ref level in entity.

Definition at line 1278 of file EntitiesFieldData.hpp.

◆ fieldData

VectorDouble MoFEM::EntitiesFieldData::EntData::fieldData

protected

Field data on entity.

Definition at line 1277 of file EntitiesFieldData.hpp.

◆ fieldEntities

VectorFieldEntities MoFEM::EntitiesFieldData::EntData::fieldEntities

protected

Field entities.

Definition at line 1276 of file EntitiesFieldData.hpp.

◆ iNdices

VectorInt MoFEM::EntitiesFieldData::EntData::iNdices

protected

Global indices on entity.

Definition at line 1273 of file EntitiesFieldData.hpp.

◆ localIndices

VectorInt MoFEM::EntitiesFieldData::EntData::localIndices

protected

Local indices on entity.

Definition at line 1274 of file EntitiesFieldData.hpp.

◆ MaxBernsteinBezierOrder

constexpr size_t MoFEM::EntitiesFieldData::EntData::MaxBernsteinBezierOrder = BITFEID_SIZE

staticconstexpr

Definition at line 1237 of file EntitiesFieldData.hpp.

◆ N

std::array<boost::shared_ptr<MatrixDouble>, LASTBASE>& MoFEM::EntitiesFieldData::EntData::N

protected

Base functions.

Definition at line 1288 of file EntitiesFieldData.hpp.

◆ oRder

ApproximationOrder MoFEM::EntitiesFieldData::EntData::oRder

protected

Entity order.

Definition at line 1270 of file EntitiesFieldData.hpp.

◆ sEnse

int MoFEM::EntitiesFieldData::EntData::sEnse

protected

Entity sense (orientation)

Definition at line 1269 of file EntitiesFieldData.hpp.

◆ sPace

FieldSpace MoFEM::EntitiesFieldData::EntData::sPace

protected

Entity space.

Definition at line 1271 of file EntitiesFieldData.hpp.

◆ swapBaseDiffNPtr

boost::shared_ptr<MatrixDouble> MoFEM::EntitiesFieldData::EntData::swapBaseDiffNPtr

protected

Used by Bernstein base to keep temporary pointer.

Swap approximation base.

Bernstein-Bezier (BB) base is not hierarchical, and is calculated for particular field, since it all shape functions change with the order. BB base is precalculated for every field, and when user push operator with particular field using BB base, pointers to shape functions and BaseDerivatives of shape functions are set to particular location, once operator is executed, pointers are switch back to its original position.

getNSharedPtr(base) <=== getBBNSharedPtr(field_name); // DO OPERATOR WORK getNSharedPtr(base) ==> getBBNSharedPtr(field_name);

Parameters

field_name	Name of the field for which to swap base
base	Approximation base type to swap

Returns: MoFEMErrorCode Error code (0 on success)

Definition at line 1319 of file EntitiesFieldData.hpp.

◆ swapBaseNPtr

boost::shared_ptr<MatrixDouble> MoFEM::EntitiesFieldData::EntData::swapBaseNPtr

protected

Used by Bernstein base to keep temporary pointer.

Swap approximation base.

Bernstein-Bezier (BB) base is not hierarchical, and is calculated for particular field, since it all shape functions change with the order. BB base is precalculated for every field, and when user push operator with particular field using BB base, pointers to shape functions and BaseDerivatives of shape functions are set to particular location, once operator is executed, pointers are switch back to its original position.

getNSharedPtr(base) <=== getBBNSharedPtr(field_name); // DO OPERATOR WORK getNSharedPtr(base) ==> getBBNSharedPtr(field_name);

Parameters

field_name	Name of the field for which to swap base
base	Approximation base type to swap

Returns: MoFEMErrorCode Error code (0 on success)

Definition at line 1312 of file EntitiesFieldData.hpp.

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

src/finite_elements/EntitiesFieldData.hpp
src/finite_elements/impl/EntitiesFieldData.cpp

Public Types

Bernstein-Bezier base only functions

Broken spaces functions

Auxiliary functions

Constructor and destructor

Sense, order and indices

Data on entity

Base and space

Get base functions for H1/L2

Get base functions for vectorial approximation bases, i.e.

Get base functions with FTensor

Specializations for tensor base function

Specializations for derivatives of base functions

Specializations for field data

Specializations for H1/L2

Specializations for HDiv/HCrul

Detailed Description

Member Enumeration Documentation

◆ BaseDerivatives

Constructor & Destructor Documentation

◆ EntData()

◆ ~EntData()

Member Function Documentation

◆ baseSwap()

◆ getBase()

◆ getBBAlphaIndices()

◆ getBBAlphaIndicesByOrderArray()

◆ getBBAlphaIndicesByOrderSharedPtr()

◆ getBBAlphaIndicesMap()

◆ getBBAlphaIndicesSharedPtr()

◆ getBBDiffNByOrderArray()

◆ getBBDiffNByOrderSharedPtr()

◆ getBBDiffNMap()

◆ getBBDiffNSharedPtr() [1/2]

◆ getBBDiffNSharedPtr() [2/2]

◆ getBBNByOrderArray()

◆ getBBNByOrderSharedPtr()

◆ getBBNMap()

◆ getBBNodeOrder()

◆ getBBNSharedPtr() [1/2]

◆ getBBNSharedPtr() [2/2]

◆ getDiffN() [1/7]

◆ getDiffN() [2/7]

◆ getDiffN() [3/7]

◆ getDiffN() [4/7]

◆ getDiffN() [5/7]

◆ getDiffN() [6/7]

◆ getDiffN() [7/7]

◆ getDiffNSharedPtr()

◆ getEntDataBitRefLevel()

◆ getFieldData() [1/2]

◆ getFieldData() [2/2]

◆ getFieldDataUpToOrder()

◆ getFieldDofs() [1/2]

◆ getFieldDofs() [2/2]

◆ getFieldEntities() [1/2]

◆ getFieldEntities() [2/2]

◆ getFTensor0FieldData()

◆ getFTensor0N() [1/4]

◆ getFTensor0N() [2/4]

◆ getFTensor0N() [3/4]

◆ getFTensor0N() [4/4]

◆ getFTensor1DiffN() [1/8]

◆ getFTensor1DiffN() [2/8]

◆ getFTensor1DiffN() [3/8]

◆ getFTensor1DiffN() [4/8]

◆ getFTensor1DiffN() [5/8]

◆ getFTensor1DiffN() [6/8]

◆ getFTensor1DiffN() [7/8]

◆ getFTensor1DiffN() [8/8]

◆ getFTensor1FieldData() [1/2]

◆ getFTensor1FieldData() [2/2]

◆ getFTensor1N() [1/6]

◆ getFTensor1N() [2/6]

◆ getFTensor1N() [3/6]

◆ getFTensor1N() [4/6]

◆ getFTensor1N() [5/6]

◆ getFTensor1N() [6/6]

◆ getFTensor2DiffN() [1/6]

◆ getFTensor2DiffN() [2/6]