![]() |
MFEM
v4.5.2
Finite element discretization library
|
Arbitrary order Nedelec elements in 3D on a cube. More...
#include <fe_nd.hpp>
Public Member Functions | |
ND_HexahedronElement (const int p, const int cb_type=BasisType::GaussLobatto, const int ob_type=BasisType::GaussLegendre) | |
Construct the ND_HexahedronElement of order p and closed and open BasisType cb_type and ob_type. More... | |
virtual void | CalcVShape (const IntegrationPoint &ip, DenseMatrix &shape) const |
Evaluate the values of all shape functions of a vector finite element in reference space at the given point ip. More... | |
virtual void | CalcVShape (ElementTransformation &Trans, DenseMatrix &shape) const |
Evaluate the values of all shape functions of a vector finite element in physical space at the point described by Trans. More... | |
virtual void | CalcCurlShape (const IntegrationPoint &ip, DenseMatrix &curl_shape) const |
Evaluate the curl of all shape functions of a vector finite element in reference space at the given point ip. More... | |
virtual void | GetLocalInterpolation (ElementTransformation &Trans, DenseMatrix &I) const |
Return the local interpolation matrix I (Dof x Dof) where the fine element is the image of the base geometry under the given transformation. More... | |
virtual void | GetLocalRestriction (ElementTransformation &Trans, DenseMatrix &R) const |
Return a local restriction matrix R (Dof x Dof) mapping fine dofs to coarse dofs. More... | |
virtual void | GetTransferMatrix (const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &I) const |
Return interpolation matrix, I, which maps dofs from a coarse element, fe, to the fine dofs on this finite element. More... | |
virtual void | Project (VectorCoefficient &vc, ElementTransformation &Trans, Vector &dofs) const |
Given a vector coefficient and a transformation, compute its projection (approximation) in the local finite dimensional space in terms of the degrees of freedom. (VectorFiniteElements) More... | |
virtual void | ProjectFromNodes (Vector &vc, ElementTransformation &Trans, Vector &dofs) const |
Given a vector of values at the finite element nodes and a transformation, compute its projection (approximation) in the local finite dimensional space in terms of the degrees of freedom. Valid for VectorFiniteElements. More... | |
virtual void | ProjectMatrixCoefficient (MatrixCoefficient &mc, ElementTransformation &T, Vector &dofs) const |
Given a matrix coefficient and a transformation, compute an approximation ("projection") in the local finite dimensional space in terms of the degrees of freedom. For VectorFiniteElements, the rows of the coefficient are projected in the vector space. More... | |
virtual void | Project (const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &I) const |
Compute the embedding/projection matrix from the given FiniteElement onto 'this' FiniteElement. The ElementTransformation is included to support cases when the projection depends on it. More... | |
virtual void | ProjectGrad (const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &grad) const |
Compute the discrete gradient matrix from the given FiniteElement onto 'this' FiniteElement. The ElementTransformation is included to support cases when the matrix depends on it. More... | |
virtual void | ProjectCurl (const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &curl) const |
Compute the discrete curl matrix from the given FiniteElement onto 'this' FiniteElement. The ElementTransformation is included to support cases when the matrix depends on it. More... | |
virtual void | Project (Coefficient &coeff, ElementTransformation &Trans, Vector &dofs) const |
Given a coefficient and a transformation, compute its projection (approximation) in the local finite dimensional space in terms of the degrees of freedom. More... | |
virtual void | Project (VectorCoefficient &vc, ElementTransformation &Trans, Vector &dofs) const |
Given a vector coefficient and a transformation, compute its projection (approximation) in the local finite dimensional space in terms of the degrees of freedom. (VectorFiniteElements) More... | |
virtual void | Project (const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &I) const |
Compute the embedding/projection matrix from the given FiniteElement onto 'this' FiniteElement. The ElementTransformation is included to support cases when the projection depends on it. More... | |
![]() | |
VectorTensorFiniteElement (const int dims, const int d, const int p, const int cbtype, const int obtype, const int M, const DofMapType dmtype) | |
VectorTensorFiniteElement (const int dims, const int d, const int p, const int obtype, const int M, const DofMapType dmtype) | |
const DofToQuad & | GetDofToQuad (const IntegrationRule &ir, DofToQuad::Mode mode) const override |
Return a DofToQuad structure corresponding to the given IntegrationRule using the given DofToQuad::Mode. More... | |
const DofToQuad & | GetDofToQuadOpen (const IntegrationRule &ir, DofToQuad::Mode mode) const |
virtual | ~VectorTensorFiniteElement () |
![]() | |
VectorFiniteElement (int D, Geometry::Type G, int Do, int O, int M, int F=FunctionSpace::Pk) | |
![]() | |
FiniteElement (int D, Geometry::Type G, int Do, int O, int F=FunctionSpace::Pk) | |
Construct FiniteElement with given. More... | |
int | GetDim () const |
Returns the reference space dimension for the finite element. More... | |
int | GetVDim () const |
Returns the vector dimension for vector-valued finite elements. More... | |
int | GetCurlDim () const |
Returns the dimension of the curl for vector-valued finite elements. More... | |
Geometry::Type | GetGeomType () const |
Returns the Geometry::Type of the reference element. More... | |
int | GetDof () const |
Returns the number of degrees of freedom in the finite element. More... | |
int | GetOrder () const |
Returns the order of the finite element. In the case of anisotropic orders, returns the maximum order. More... | |
bool | HasAnisotropicOrders () const |
Returns true if the FiniteElement basis may be using different orders/degrees in different spatial directions. More... | |
const int * | GetAnisotropicOrders () const |
Returns an array containing the anisotropic orders/degrees. More... | |
int | Space () const |
Returns the type of FunctionSpace on the element. More... | |
int | GetRangeType () const |
Returns the FiniteElement::RangeType of the element, one of {SCALAR, VECTOR}. More... | |
int | GetDerivRangeType () const |
Returns the FiniteElement::RangeType of the element derivative, either SCALAR or VECTOR. More... | |
int | GetMapType () const |
Returns the FiniteElement::MapType of the element describing how reference functions are mapped to physical space, one of {VALUE, INTEGRAL H_DIV, H_CURL}. More... | |
int | GetDerivType () const |
Returns the FiniteElement::DerivType of the element describing the spatial derivative method implemented, one of {NONE, GRAD, DIV, CURL}. More... | |
int | GetDerivMapType () const |
Returns the FiniteElement::DerivType of the element describing how reference function derivatives are mapped to physical space, one of {VALUE, INTEGRAL, H_DIV, H_CURL}. More... | |
void | CalcPhysShape (ElementTransformation &Trans, Vector &shape) const |
Evaluate the values of all shape functions of a scalar finite element in physical space at the point described by Trans. More... | |
void | CalcPhysDShape (ElementTransformation &Trans, DenseMatrix &dshape) const |
Evaluate the gradients of all shape functions of a scalar finite element in physical space at the point described by Trans. More... | |
const IntegrationRule & | GetNodes () const |
Get a const reference to the nodes of the element. More... | |
void | CalcPhysVShape (ElementTransformation &Trans, DenseMatrix &shape) const |
Equivalent to the CalcVShape() method with the same arguments. More... | |
virtual void | CalcDivShape (const IntegrationPoint &ip, Vector &divshape) const |
Evaluate the divergence of all shape functions of a vector finite element in reference space at the given point ip. More... | |
void | CalcPhysDivShape (ElementTransformation &Trans, Vector &divshape) const |
Evaluate the divergence of all shape functions of a vector finite element in physical space at the point described by Trans. More... | |
virtual void | CalcPhysCurlShape (ElementTransformation &Trans, DenseMatrix &curl_shape) const |
Evaluate the curl of all shape functions of a vector finite element in physical space at the point described by Trans. More... | |
virtual void | GetFaceDofs (int face, int **dofs, int *ndofs) const |
Get the dofs associated with the given face. *dofs is set to an internal array of the local dofc on the face, while *ndofs is set to the number of dofs on that face. More... | |
virtual void | CalcHessian (const IntegrationPoint &ip, DenseMatrix &Hessian) const |
Evaluate the Hessians of all shape functions of a scalar finite element in reference space at the given point ip. More... | |
virtual void | CalcPhysHessian (ElementTransformation &Trans, DenseMatrix &Hessian) const |
Evaluate the Hessian of all shape functions of a scalar finite element in reference space at the given point ip. More... | |
virtual void | CalcPhysLaplacian (ElementTransformation &Trans, Vector &Laplacian) const |
Evaluate the Laplacian of all shape functions of a scalar finite element in reference space at the given point ip. More... | |
virtual void | CalcPhysLinLaplacian (ElementTransformation &Trans, Vector &Laplacian) const |
virtual void | Project (Coefficient &coeff, ElementTransformation &Trans, Vector &dofs) const |
Given a coefficient and a transformation, compute its projection (approximation) in the local finite dimensional space in terms of the degrees of freedom. More... | |
virtual void | ProjectDelta (int vertex, Vector &dofs) const |
Project a delta function centered on the given vertex in the local finite dimensional space represented by the dofs. More... | |
virtual void | ProjectDiv (const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &div) const |
Compute the discrete divergence matrix from the given FiniteElement onto 'this' FiniteElement. The ElementTransformation is included to support cases when the matrix depends on it. More... | |
virtual | ~FiniteElement () |
Deconstruct the FiniteElement. More... | |
![]() | |
TensorBasisElement (const int dims, const int p, const int btype, const DofMapType dmtype) | |
int | GetBasisType () const |
const Poly_1D::Basis & | GetBasis1D () const |
const Array< int > & | GetDofMap () const |
Get an Array<int> that maps lexicographically ordered indices to the indices of the respective nodes/dofs/basis functions. If the dofs are ordered lexicographically, i.e. the mapping is identity, the returned Array will be empty. More... | |
Protected Member Functions | |
void | ProjectIntegrated (VectorCoefficient &vc, ElementTransformation &Trans, Vector &dofs) const |
![]() | |
void | SetDerivMembers () |
void | CalcVShape_RT (ElementTransformation &Trans, DenseMatrix &shape) const |
void | CalcVShape_ND (ElementTransformation &Trans, DenseMatrix &shape) const |
void | Project_RT (const double *nk, const Array< int > &d2n, VectorCoefficient &vc, ElementTransformation &Trans, Vector &dofs) const |
Project a vector coefficient onto the RT basis functions. More... | |
void | Project_RT (const double *nk, const Array< int > &d2n, Vector &vc, ElementTransformation &Trans, Vector &dofs) const |
Projects the vector of values given at FE nodes to RT space. More... | |
void | ProjectMatrixCoefficient_RT (const double *nk, const Array< int > &d2n, MatrixCoefficient &mc, ElementTransformation &T, Vector &dofs) const |
Project the rows of the matrix coefficient in an RT space. More... | |
void | Project_RT (const double *nk, const Array< int > &d2n, const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &I) const |
Project vector-valued basis functions onto the RT basis functions. More... | |
void | ProjectGrad_RT (const double *nk, const Array< int > &d2n, const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &grad) const |
void | ProjectCurl_ND (const double *tk, const Array< int > &d2t, const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &curl) const |
void | ProjectCurl_RT (const double *nk, const Array< int > &d2n, const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &curl) const |
void | Project_ND (const double *tk, const Array< int > &d2t, VectorCoefficient &vc, ElementTransformation &Trans, Vector &dofs) const |
Project a vector coefficient onto the ND basis functions. More... | |
void | Project_ND (const double *tk, const Array< int > &d2t, Vector &vc, ElementTransformation &Trans, Vector &dofs) const |
Projects the vector of values given at FE nodes to ND space. More... | |
void | ProjectMatrixCoefficient_ND (const double *tk, const Array< int > &d2t, MatrixCoefficient &mc, ElementTransformation &T, Vector &dofs) const |
Project the rows of the matrix coefficient in an ND space. More... | |
void | Project_ND (const double *tk, const Array< int > &d2t, const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &I) const |
Project vector-valued basis functions onto the ND basis functions. More... | |
void | ProjectGrad_ND (const double *tk, const Array< int > &d2t, const FiniteElement &fe, ElementTransformation &Trans, DenseMatrix &grad) const |
void | LocalL2Projection_RT (const VectorFiniteElement &cfe, ElementTransformation &Trans, DenseMatrix &I) const |
void | LocalInterpolation_RT (const VectorFiniteElement &cfe, const double *nk, const Array< int > &d2n, ElementTransformation &Trans, DenseMatrix &I) const |
void | LocalL2Projection_ND (const VectorFiniteElement &cfe, ElementTransformation &Trans, DenseMatrix &I) const |
void | LocalInterpolation_ND (const VectorFiniteElement &cfe, const double *tk, const Array< int > &d2t, ElementTransformation &Trans, DenseMatrix &I) const |
void | LocalRestriction_RT (const double *nk, const Array< int > &d2n, ElementTransformation &Trans, DenseMatrix &R) const |
void | LocalRestriction_ND (const double *tk, const Array< int > &d2t, ElementTransformation &Trans, DenseMatrix &R) const |
Additional Inherited Members | |
![]() | |
enum | RangeType { UNKNOWN_RANGE_TYPE = -1, SCALAR, VECTOR } |
Enumeration for range_type and deriv_range_type. More... | |
enum | MapType { UNKNOWN_MAP_TYPE = -1, VALUE, INTEGRAL, H_DIV, H_CURL } |
Enumeration for MapType: defines how reference functions are mapped to physical space. More... | |
enum | DerivType { NONE, GRAD, DIV, CURL } |
Enumeration for DerivType: defines which derivative method is implemented. More... | |
![]() | |
enum | DofMapType { L2_DOF_MAP = 0, H1_DOF_MAP = 1, Sr_DOF_MAP = 2 } |
![]() | |
static bool | IsClosedType (int b_type) |
Return true if the BasisType of b_type is closed (has Quadrature1D points on the boundary). More... | |
static bool | IsOpenType (int b_type) |
Return true if the BasisType of b_type is open (doesn't have Quadrature1D points on the boundary). More... | |
static int | VerifyClosed (int b_type) |
Ensure that the BasisType of b_type is closed (has Quadrature1D points on the boundary). More... | |
static int | VerifyOpen (int b_type) |
Ensure that the BasisType of b_type is open (doesn't have Quadrature1D points on the boundary). More... | |
static int | VerifyNodal (int b_type) |
Ensure that the BasisType of b_type nodal (satisfies the interpolation property). More... | |
![]() | |
static Geometry::Type | GetTensorProductGeometry (int dim) |
static int | Pow (int base, int dim) |
Return base raised to the power dim. More... | |
static const DofToQuad & | GetTensorDofToQuad (const FiniteElement &fe, const IntegrationRule &ir, DofToQuad::Mode mode, const Poly_1D::Basis &basis, bool closed, Array< DofToQuad *> &dof2quad_array) |
![]() | |
static const VectorFiniteElement & | CheckVectorFE (const FiniteElement &fe) |
![]() | |
Poly_1D::Basis & | obasis1d |
![]() | |
bool | is_nodal |
DenseMatrix | JtJ |
DenseMatrix | curlshape |
DenseMatrix | curlshape_J |
![]() | |
int | dim |
Dimension of reference space. More... | |
int | vdim |
Vector dimension of vector-valued basis functions. More... | |
int | cdim |
Dimension of curl for vector-valued basis functions. More... | |
Geometry::Type | geom_type |
Geometry::Type of the reference element. More... | |
int | func_space |
int | range_type |
int | map_type |
int | deriv_type |
int | deriv_range_type |
int | deriv_map_type |
int | dof |
Number of degrees of freedom. More... | |
int | order |
Order/degree of the shape functions. More... | |
int | orders [Geometry::MaxDim] |
Anisotropic orders. More... | |
IntegrationRule | Nodes |
DenseMatrix | vshape |
Array< DofToQuad * > | dof2quad_array |
Container for all DofToQuad objects created by the FiniteElement. More... | |
![]() | |
int | b_type |
Array< int > | dof_map |
Poly_1D::Basis & | basis1d |
Array< int > | inv_dof_map |
mfem::ND_HexahedronElement::ND_HexahedronElement | ( | const int | p, |
const int | cb_type = BasisType::GaussLobatto , |
||
const int | ob_type = BasisType::GaussLegendre |
||
) |
Construct the ND_HexahedronElement of order p and closed and open BasisType cb_type and ob_type.
|
virtual |
Evaluate the curl of all shape functions of a vector finite element in reference space at the given point ip.
Each row of the result DenseMatrix curl_shape contains the components of the curl of one vector shape function. The size (dof x CDim) of curl_shape must be set in advance, where CDim = 3 for dim = 3 and CDim = 1 for dim = 2.
Reimplemented from mfem::FiniteElement.
|
virtual |
Evaluate the values of all shape functions of a vector finite element in reference space at the given point ip.
Each row of the result DenseMatrix shape contains the components of one vector shape function. The size (dof x dim) of shape must be set in advance.
Reimplemented from mfem::FiniteElement.
|
inlinevirtual |
Evaluate the values of all shape functions of a vector finite element in physical space at the point described by Trans.
Each row of the result DenseMatrix shape contains the components of one vector shape function. The size (dof x SDim) of shape must be set in advance, where SDim >= dim is the physical space dimension as described by Trans.
Reimplemented from mfem::FiniteElement.
|
inlinevirtual |
Return the local interpolation matrix I (Dof x Dof) where the fine element is the image of the base geometry under the given transformation.
Reimplemented from mfem::FiniteElement.
|
inlinevirtual |
Return a local restriction matrix R (Dof x Dof) mapping fine dofs to coarse dofs.
The fine element is the image of the base geometry under the given transformation, Trans.
The assumption in this method is that a subset of the coarse dofs can be expressed only in terms of the dofs of the given fine element.
Rows in R corresponding to coarse dofs that cannot be expressed in terms of the fine dofs will be marked as invalid by setting the first entry (column 0) in the row to infinity().
This method assumes that the dimensions of R are set before it is called.
Reimplemented from mfem::FiniteElement.
|
inlinevirtual |
Return interpolation matrix, I, which maps dofs from a coarse element, fe, to the fine dofs on this finite element.
Trans represents the mapping from the reference element of this element into a subset of the reference space of the element fe, thus allowing the "coarse" FiniteElement to be different from the "fine" FiniteElement as when h-refinement is combined with p-refinement or p-derefinement. It is assumed that both finite elements use the same FiniteElement::MapType.
Reimplemented from mfem::FiniteElement.
void mfem::FiniteElement::Project |
Given a vector coefficient and a transformation, compute its projection (approximation) in the local finite dimensional space in terms of the degrees of freedom. (VectorFiniteElements)
The approximation used to project is usually local interpolation of degrees of freedom. The derived class could use other methods not implemented yet, e.g. local L2 projection.
Definition at line 131 of file fe_base.cpp.
void mfem::FiniteElement::Project |
Compute the embedding/projection matrix from the given FiniteElement onto 'this' FiniteElement. The ElementTransformation is included to support cases when the projection depends on it.
Definition at line 154 of file fe_base.cpp.
void mfem::FiniteElement::Project |
Given a coefficient and a transformation, compute its projection (approximation) in the local finite dimensional space in terms of the degrees of freedom.
The approximation used to project is usually local interpolation of degrees of freedom. The derived class could use other methods not implemented yet, e.g. local L2 projection.
Definition at line 125 of file fe_base.cpp.
|
inlinevirtual |
Given a vector coefficient and a transformation, compute its projection (approximation) in the local finite dimensional space in terms of the degrees of freedom. (VectorFiniteElements)
The approximation used to project is usually local interpolation of degrees of freedom. The derived class could use other methods not implemented yet, e.g. local L2 projection.
Reimplemented from mfem::FiniteElement.
|
inlinevirtual |
Compute the embedding/projection matrix from the given FiniteElement onto 'this' FiniteElement. The ElementTransformation is included to support cases when the projection depends on it.
Reimplemented from mfem::FiniteElement.
|
inlinevirtual |
Compute the discrete curl matrix from the given FiniteElement onto 'this' FiniteElement. The ElementTransformation is included to support cases when the matrix depends on it.
Reimplemented from mfem::FiniteElement.
|
inlinevirtual |
Given a vector of values at the finite element nodes and a transformation, compute its projection (approximation) in the local finite dimensional space in terms of the degrees of freedom. Valid for VectorFiniteElements.
Reimplemented from mfem::FiniteElement.
|
inlinevirtual |
Compute the discrete gradient matrix from the given FiniteElement onto 'this' FiniteElement. The ElementTransformation is included to support cases when the matrix depends on it.
Reimplemented from mfem::FiniteElement.
|
protected |
|
inlinevirtual |
Given a matrix coefficient and a transformation, compute an approximation ("projection") in the local finite dimensional space in terms of the degrees of freedom. For VectorFiniteElements, the rows of the coefficient are projected in the vector space.
Reimplemented from mfem::FiniteElement.