Intrepid2
|
Public Types | |
using | BasisPtr = Teuchos::RCP<BasisBase> |
![]() | |
using | BasisBase = BasisBaseClass |
using | BasisPtr = Teuchos::RCP<BasisBase> |
using | DeviceType = typename BasisBase::DeviceType |
using | ExecutionSpace = typename BasisBase::ExecutionSpace |
using | OutputValueType = typename BasisBase::OutputValueType |
using | PointValueType = typename BasisBase::PointValueType |
using | OrdinalTypeArray1DHost = typename BasisBase::OrdinalTypeArray1DHost |
using | OrdinalTypeArray2DHost = typename BasisBase::OrdinalTypeArray2DHost |
using | OutputViewType = typename BasisBase::OutputViewType |
using | PointViewType = typename BasisBase::PointViewType |
using | TensorBasis = Basis_TensorBasis<BasisBaseClass> |
Public Member Functions | |
Basis_TensorBasis3 (BasisPtr basis1, BasisPtr basis2, BasisPtr basis3, const bool useShardsCellTopologyAndTags=false) | |
virtual OperatorTensorDecomposition | getOperatorDecomposition (const EOperator operatorType) const override |
Returns a full decomposition of the specified operator. (Full meaning that all TensorBasis components are expanded into their non-TensorBasis components.) | |
virtual void | getValues (OutputViewType outputValues, const EOperator operatorType, const PointViewType inputPoints12, const PointViewType inputPoints3, bool tensorPoints) const override |
Evaluation of a tensor FEM basis on a reference cell. | |
virtual void | getValues (OutputViewType outputValues, const EOperator operatorType, const PointViewType inputPoints1, const PointViewType inputPoints2, const PointViewType inputPoints3, bool tensorPoints) const =0 |
Evaluation of a tensor FEM basis on a reference cell; subclasses should override this. | |
void | getValues (OutputViewType outputValues, const PointViewType inputPoints1, const EOperator operatorType1, const PointViewType inputPoints2, const EOperator operatorType2, const PointViewType inputPoints3, const EOperator operatorType3, bool tensorPoints, double weight=1.0) const |
Evaluation of a tensor FEM basis on a reference cell; subclasses should override this. | |
virtual HostBasisPtr< OutputValueType, PointValueType > | getHostBasis () const override |
Creates and returns a Basis object whose DeviceType template argument is Kokkos::HostSpace::device_type, but is otherwise identical to this. | |
virtual void | getValues (BasisValues< OutputValueType, DeviceType > outputValues, const TensorPoints< PointValueType, DeviceType > inputPoints, const EOperator operatorType=OPERATOR_VALUE) const override |
Evaluation of a FEM basis on a reference cell, using point and output value containers that allow preservation of tensor-product structure. | |
void | getValues (OutputViewType outputValues, const PointViewType inputPoints, const EOperator operatorType=OPERATOR_VALUE) const override |
Evaluation of a FEM basis on a reference cell. | |
virtual void | getValues (OutputViewType outputValues, const EOperator operatorType, const PointViewType inputPoints1, const PointViewType inputPoints2, bool tensorPoints) const |
Evaluation of a tensor FEM basis on a reference cell; subclasses should override this. | |
void | getValues (OutputViewType outputValues, const PointViewType inputPoints1, const EOperator operatorType1, const PointViewType inputPoints2, const EOperator operatorType2, bool tensorPoints, double weight=1.0) const |
Evaluation of a tensor FEM basis on a reference cell. | |
![]() | |
Basis_TensorBasis (BasisPtr basis1, BasisPtr basis2, EFunctionSpace functionSpace=FUNCTION_SPACE_MAX, const bool useShardsCellTopologyAndTags=false) | |
Constructor. | |
void | setShardsTopologyAndTags () |
virtual int | getNumTensorialExtrusions () const override |
ordinal_type | getTensorDkEnumeration (ordinal_type dkEnum1, ordinal_type operatorOrder1, ordinal_type dkEnum2, ordinal_type operatorOrder2) const |
Given "Dk" enumeration indices for the component bases, returns a Dk enumeration index for the composite basis. | |
virtual OperatorTensorDecomposition | getSimpleOperatorDecomposition (const EOperator &operatorType) const |
Returns a simple decomposition of the specified operator: what operator(s) should be applied to basis1, and what operator(s) to basis2. A one-element OperatorTensorDecomposition corresponds to a single TensorData entry; a multiple-element OperatorTensorDecomposition corresponds to a VectorData object with axialComponents = false. | |
virtual BasisValues< OutputValueType, DeviceType > | allocateBasisValues (TensorPoints< PointValueType, DeviceType > points, const EOperator operatorType=OPERATOR_VALUE) const override |
Allocate BasisValues container suitable for passing to the getValues() variant that takes a TensorPoints container as argument. | |
void | getComponentPoints (const PointViewType inputPoints, const bool attemptTensorDecomposition, PointViewType &inputPoints1, PointViewType &inputPoints2, bool &tensorDecompositionSucceeded) const |
Method to extract component points from composite points. | |
virtual void | getDofCoords (typename BasisBase::ScalarViewType dofCoords) const override |
Fills in spatial locations (coordinates) of degrees of freedom (nodes) on the reference cell. | |
virtual void | getDofCoeffs (typename BasisBase::ScalarViewType dofCoeffs) const override |
Fills in coefficients of degrees of freedom on the reference cell. | |
virtual const char * | getName () const override |
Returns basis name. | |
std::vector< BasisPtr > | getTensorBasisComponents () const |
void | getValues (OutputViewType outputValues, const PointViewType inputPoints, const EOperator operatorType=OPERATOR_VALUE) const override |
Evaluation of a FEM basis on a reference cell. | |
void | getValues (OutputViewType outputValues, const PointViewType inputPoints1, const EOperator operatorType1, const PointViewType inputPoints2, const EOperator operatorType2, bool tensorPoints, double weight=1.0) const |
Evaluation of a tensor FEM basis on a reference cell. | |
Protected Attributes | |
BasisPtr | basis1_ |
BasisPtr | basis2_ |
BasisPtr | basis3_ |
![]() | |
BasisPtr | basis1_ |
BasisPtr | basis2_ |
std::vector< BasisPtr > | tensorComponents_ |
std::string | name_ |
int | numTensorialExtrusions_ |
Private Types | |
using | BasisBase = BasisBaseClass |
using | TensorBasis = Basis_TensorBasis<BasisBase> |
Definition at line 2125 of file Intrepid2_TensorBasis.hpp.
|
private |
Definition at line 2128 of file Intrepid2_TensorBasis.hpp.
using Intrepid2::Basis_TensorBasis3< BasisBaseClass >::BasisPtr = Teuchos::RCP<BasisBase> |
Definition at line 2140 of file Intrepid2_TensorBasis.hpp.
|
private |
Definition at line 2129 of file Intrepid2_TensorBasis.hpp.
|
inline |
Definition at line 2146 of file Intrepid2_TensorBasis.hpp.
|
inlineoverridevirtual |
Creates and returns a Basis object whose DeviceType template argument is Kokkos::HostSpace::device_type, but is otherwise identical to this.
Reimplemented from Intrepid2::Basis_TensorBasis< BasisBaseClass >.
Definition at line 2384 of file Intrepid2_TensorBasis.hpp.
|
inlineoverridevirtual |
Returns a full decomposition of the specified operator. (Full meaning that all TensorBasis components are expanded into their non-TensorBasis components.)
The outer container has a length that corresponds to the number of vector components in output; the inner container has one entry per basis component.
A one-element outer container corresponds to a single TensorData entry; a multiple-element outer container corresponds to a VectorData object with axialComponents = false.
Reimplemented from Intrepid2::Basis_TensorBasis< BasisBaseClass >.
Definition at line 2162 of file Intrepid2_TensorBasis.hpp.
References Intrepid2::OperatorTensorDecomposition::expandedDecomposition(), and Intrepid2::Basis_TensorBasis< BasisBaseClass >::getSimpleOperatorDecomposition().
|
inlineoverridevirtual |
Evaluation of a FEM basis on a reference cell, using point and output value containers that allow preservation of tensor-product structure.
Returns values of operatorType acting on FEM basis functions for a set of points in the reference cell for which the basis is defined.
outputValues | [out] - variable rank array with the basis values. Should be allocated using Basis::allocateBasisValues(). |
inputPoints | [in] - rank-2 array (P,D) with the evaluation points. This should be allocated using Cubature::allocateCubaturePoints() and filled using Cubature::getCubature(). |
operatorType | [in] - the operator acting on the basis function |
This is the preferred getValues() method for TensorBasis and DirectSumBasis and their subclasses. It allows a reduced memory footprint and optimized integration, etc.
Reimplemented from Intrepid2::Basis_TensorBasis< BasisBaseClass >.
Definition at line 1459 of file Intrepid2_TensorBasis.hpp.
|
inlinevirtual |
Evaluation of a tensor FEM basis on a reference cell; subclasses should override this.
Returns values of operatorType acting on FEM basis functions for a set of points in the reference cell for which the basis is defined.
outputValues | [out] - variable rank array with the basis values |
operatorType | [in] - the operator acting on the basis functions |
inputPoints1 | [in] - rank-2 array (P1,D1) with the evaluation points for basis1 |
inputPoints2 | [in] - rank-2 array (P2,D2) with the evaluation points for basis2 |
tensorPoints | [in] - whether the points should be interpreted as tensor components of the evaluation points, or in a one-to-one correspondence |
Subclasses should override this method; this gives them an opportunity to specify how operatorType should be decomposed into operators on the component bases.
If tensorPoints is true, then the points dimension of outputValues should be (P1*P2). If tensorPoints is false, then P1 should equal P2, and these should match the points dimension of outputValues.
There are three variants of getValues:
The intent is that subclasses implement this five-argument version; in that implementation, they need to do little else than call the seven-argument version below.
Note that the three-argument implementation handles the OPERATOR_Dn operators directly; that is, subclasses can omit any consideration of OPERATOR_Dn operators in their implementation of the five-argument version.
Reimplemented from Intrepid2::Basis_TensorBasis< BasisBaseClass >.
Definition at line 1756 of file Intrepid2_TensorBasis.hpp.
|
pure virtual |
Evaluation of a tensor FEM basis on a reference cell; subclasses should override this.
Returns values of operatorType acting on FEM basis functions for a set of points in the reference cell for which the basis is defined.
outputValues | [out] - variable rank array with the basis values |
operatorType | [in] - the operator acting on the basis functions |
inputPoints1 | [in] - rank-2 array (P1,D1) with the evaluation points for basis1 |
inputPoints1 | [in] - rank-2 array (P2,D2) with the evaluation points for basis2 |
inputPoints3 | [in] - rank-2 array (P3,D3) with the evaluation points for basis3 |
tensorPoints | [in] - whether the points should be interpreted as tensor components of the evaluation points, or in a one-to-one correspondence |
Subclasses should override this method; this gives them an opportunity to specify how operatorType should be decomposed into operators on the component bases.
If tensorPoints is true, then the points dimension of outputValues should be (P1*P2*P3). If tensorPoints is false, then P1 should equal P2 and P2 should equal P3, and these should match the points dimension of outputValues.
There are four variants of getValues:
The intent is that subclasses implement this six-argument version; in that implementation, they need to do little else than call the nine-argument version below.
Note that the three-argument implementation handles the OPERATOR_Dn operators directly; that is, subclasses can omit any consideration of OPERATOR_Dn operators in their implementation of the five-argument version.
|
inlineoverridevirtual |
Evaluation of a tensor FEM basis on a reference cell.
Returns values of operatorType acting on FEM basis functions for a set of points in the reference cell for which the basis is defined.
outputValues | [out] - variable rank array with the basis values |
operatorType | [in] - the operator acting on the basis functions |
inputPoints12 | [in] - rank-2 array (P12,D12) with the evaluation points for basis12 |
inputPoints3 | [in] - rank-2 array (P3,D3) with the evaluation points for basis3 |
tensorPoints | [in] - whether the points should be interpreted as tensor components of the evaluation points, or in a one-to-one correspondence |
If tensorPoints is true, then the points dimension of outputValues should be (P12*P3). If tensorPoints is false, then P12 should equal P3, and these should match the points dimension of outputValues.
There are four variants of getValues:
The intent is that subclasses implement the six-argument version; in that implementation, they need to do little else than call the nine-argument version below.
Note that the three-argument implementation handles the OPERATOR_Dn operators directly; that is, subclasses can omit any consideration of OPERATOR_Dn operators in their implementation of the six-argument version.
Reimplemented from Intrepid2::Basis_TensorBasis< BasisBaseClass >.
Definition at line 2195 of file Intrepid2_TensorBasis.hpp.
References Intrepid2::Basis_TensorBasis3< BasisBaseClass >::getValues().
Referenced by Intrepid2::Basis_Derived_HCURL_Family1_HEX< HGRAD_LINE, HVOL_LINE >::getValues(), Intrepid2::Basis_Derived_HCURL_Family2_HEX< HGRAD_LINE, HVOL_LINE >::getValues(), Intrepid2::Basis_Derived_HCURL_Family3_HEX< HGRAD_LINE, HVOL_LINE >::getValues(), Intrepid2::Basis_Derived_HDIV_Family1_HEX< HGRAD_LINE, HVOL_LINE >::getValues(), Intrepid2::Basis_Derived_HDIV_Family2_HEX< HGRAD_LINE, HVOL_LINE >::getValues(), Intrepid2::Basis_Derived_HDIV_Family3_HEX< HGRAD_LINE, HVOL_LINE >::getValues(), and Intrepid2::Basis_TensorBasis3< BasisBaseClass >::getValues().
|
inlineoverride |
Evaluation of a FEM basis on a reference cell.
Returns values of operatorType acting on FEM basis functions for a set of points in the reference cell for which the basis is defined.
outputValues | [out] - variable rank array with the basis values |
inputPoints | [in] - rank-2 array (P,D) with the evaluation points |
operatorType | [in] - the operator acting on the basis functions |
Definition at line 1614 of file Intrepid2_TensorBasis.hpp.
|
inline |
Evaluation of a tensor FEM basis on a reference cell.
Returns values of operatorType acting on FEM basis functions for a set of points in the reference cell for which the basis is defined.
outputValues | [out] - variable rank array with the basis values |
inputPoints1 | [in] - rank-2 array (P1,D1) with the evaluation points for basis1 |
operatorType1 | [in] - the operator acting on basis1 |
inputPoints2 | [in] - rank-2 array (P2,D2) with the evaluation points for basis2 |
operatorType2 | [in] - the operator acting on basis2 |
tensorPoints | [in] - whether the points should be interpreted as tensor components of the evaluation points, or in a one-to-one correspondence |
weight | [in] - optional weight (typically 1.0 or -1.0) |
If tensorPoints is true, then the points dimension of outputValues should be (P1*P2). If tensorPoints is false, then P1 should equal P2, and these should match the points dimension of outputValues.
There are three variants of getValues:
Subclasses should override the five-argument version above; in their implementation, they need to do little else than call this seven-argument version.
Definition at line 1786 of file Intrepid2_TensorBasis.hpp.
|
inline |
Evaluation of a tensor FEM basis on a reference cell; subclasses should override this.
Returns values of operatorType acting on FEM basis functions for a set of points in the reference cell for which the basis is defined.
outputValues | [out] - variable rank array with the basis values |
inputPoints1 | [in] - rank-2 array (P1,D1) with the evaluation points for basis1 |
operatorType1 | [in] - the operator acting on basis1 |
inputPoints2 | [in] - rank-2 array (P2,D2) with the evaluation points for basis2 |
operatorType2 | [in] - the operator acting on basis2 |
inputPoints3 | [in] - rank-2 array (P3,D3) with the evaluation points for basis3 |
operatorType3 | [in] - the operator acting on basis3 |
tensorPoints | [in] - whether the points should be interpreted as tensor components of the evaluation points, or in a one-to-one correspondence |
If tensorPoints is true, then the points dimension of outputValues should be (P1*P2*P3). If tensorPoints is false, then P1 should equal P2 and P2 should equal P3, and these should match the points dimension of outputValues.
There are four variants of getValues:
The intent is that subclasses implement this six-argument version; in that implementation, they need to do little else than call the nine-argument version below.
Note that the three-argument implementation handles the OPERATOR_Dn operators directly; that is, subclasses can omit any consideration of OPERATOR_Dn operators in their implementation of the five-argument version.
Definition at line 2282 of file Intrepid2_TensorBasis.hpp.
References Intrepid2::getMatchingViewWithLabel().
|
protected |
Definition at line 2142 of file Intrepid2_TensorBasis.hpp.
|
protected |
Definition at line 2143 of file Intrepid2_TensorBasis.hpp.
|
protected |
Definition at line 2144 of file Intrepid2_TensorBasis.hpp.