50 #ifndef __INTREPID2_HGRAD_TET_Cn_FEM_ORTH_HPP__
51 #define __INTREPID2_HGRAD_TET_Cn_FEM_ORTH_HPP__
69 template<
typename OutputViewType,
70 typename inputViewType,
71 typename workViewType,
72 bool hasDeriv, ordinal_type n>
74 KOKKOS_INLINE_FUNCTION
76 generate( OutputViewType output,
77 const inputViewType input,
79 const ordinal_type p );
84 template<
typename OutputViewType,
85 typename inputViewType,
86 typename workViewType,
89 KOKKOS_INLINE_FUNCTION
91 generate( OutputViewType output,
92 const inputViewType input,
94 const ordinal_type p );
99 template<
typename OutputViewType,
100 typename inputViewType,
101 typename workViewType,
104 KOKKOS_INLINE_FUNCTION
106 generate( OutputViewType output,
107 const inputViewType input,
109 const ordinal_type p );
121 template<EOperator opType>
123 template<
typename OutputViewType,
124 typename inputViewType,
125 typename workViewType>
126 KOKKOS_INLINE_FUNCTION
128 getValues( OutputViewType output,
129 const inputViewType input,
131 const ordinal_type order);
134 template<
typename DeviceType, ordinal_type numPtsPerEval,
135 typename outputValueValueType,
class ...outputValueProperties,
136 typename inputPointValueType,
class ...inputPointProperties>
138 getValues(
const typename DeviceType::execution_space& space,
139 Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValues,
140 const Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPoints,
141 const ordinal_type order,
142 const EOperator operatorType );
147 template<
typename outputValueViewType,
148 typename inputPointViewType,
149 typename workViewType,
151 ordinal_type numPtsEval>
153 outputValueViewType _outputValues;
154 const inputPointViewType _inputPoints;
156 const ordinal_type _order;
158 KOKKOS_INLINE_FUNCTION
159 Functor( outputValueViewType outputValues_,
160 inputPointViewType inputPoints_,
162 const ordinal_type order_ )
163 : _outputValues(outputValues_), _inputPoints(inputPoints_), _work(work_),_order(order_){}
165 KOKKOS_INLINE_FUNCTION
166 void operator()(
const size_type iter)
const {
170 const auto ptRange = Kokkos::pair<ordinal_type,ordinal_type>(ptBegin, ptEnd);
171 const auto input = Kokkos::subview( _inputPoints, ptRange, Kokkos::ALL() );
174 case OPERATOR_VALUE : {
175 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange );
182 const auto work = Kokkos::subview( _work, Kokkos::ALL(), ptRange, Kokkos::ALL() );
183 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
197 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
202 INTREPID2_TEST_FOR_ABORT(
true,
203 ">>> ERROR: (Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH::Functor) operator is not supported");
214 template<
typename DeviceType = void,
215 typename outputValueType = double,
216 typename pointValueType =
double>
218 :
public Basis<DeviceType,outputValueType,pointValueType> {
220 typedef double value_type;
244 const EOperator operatorType = OPERATOR_VALUE )
const override {
245 #ifdef HAVE_INTREPID2_DEBUG
246 Intrepid2::getValues_HGRAD_Args(outputValues,
253 Impl::Basis_HGRAD_TET_Cn_FEM_ORTH::
254 getValues<DeviceType,numPtsPerEval>(space,
Kokkos::View< ordinal_type *, typename ExecutionSpace::array_layout, Kokkos::HostSpace > OrdinalTypeArray1DHost
View type for 1d host array.
See Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH.
ordinal_type getCardinality() const
Returns cardinality of the basis.
See Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH.
virtual void getValues(const ExecutionSpace &, OutputViewType, const PointViewType, const EOperator=OPERATOR_VALUE) const
Evaluation of a FEM basis on a reference cell.
virtual void getValues(const ExecutionSpace &space, OutputViewType outputValues, const PointViewType inputPoints, const EOperator operatorType=OPERATOR_VALUE) const override
Evaluation of a FEM basis on a reference cell.
Kokkos::DynRankView< OutputValueType, Kokkos::LayoutStride, DeviceType > OutputViewType
View type for basis value output.
An abstract base class that defines interface for concrete basis implementations for Finite Element (...
shards::CellTopology getBaseCellTopology() const
Returns the base cell topology for which the basis is defined. See Shards documentation https://trili...
Kokkos::DynRankView< scalarType, Kokkos::LayoutStride, DeviceType > ScalarViewType
View type for scalars.
typename DeviceType::execution_space ExecutionSpace
(Kokkos) Execution space for basis.
See Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH.
Kokkos::DynRankView< PointValueType, Kokkos::LayoutStride, DeviceType > PointViewType
View type for input points.
Implementation of the default H(grad)-compatible orthogonal basis of arbitrary degree on tetrahedron...
Kokkos::View< ordinal_type **, typename ExecutionSpace::array_layout, Kokkos::HostSpace > OrdinalTypeArray2DHost
View type for 2d host array.
static constexpr ordinal_type MaxNumPtsPerBasisEval
The maximum number of points to eval in serial mode.
Basis_HGRAD_TET_Cn_FEM_ORTH(const ordinal_type order)
Constructor.
Kokkos::DynRankView< OutputValueType, Kokkos::LayoutStride, DeviceType > OutputViewType
View type for basis value output.
Definition file for FEM orthogonal basis functions of degree n for H(grad) functions on TET cells...
Kokkos::View< ordinal_type ***, typename ExecutionSpace::array_layout, Kokkos::HostSpace > OrdinalTypeArray3DHost
View type for 3d host array.
ordinal_type getDegree() const
Returns the degree of the basis.
Kokkos::DynRankView< PointValueType, Kokkos::LayoutStride, DeviceType > PointViewType
View type for input points.
Header file for the abstract base class Intrepid2::Basis.
typename DeviceType::execution_space ExecutionSpace
(Kokkos) Execution space for basis.
See Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH.