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( Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValues,
139 const Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPoints,
140 const ordinal_type order,
146 template<
typename outputValueViewType,
147 typename inputPointViewType,
148 typename workViewType,
150 ordinal_type numPtsEval>
152 outputValueViewType _outputValues;
153 const inputPointViewType _inputPoints;
155 const ordinal_type _order;
157 KOKKOS_INLINE_FUNCTION
158 Functor( outputValueViewType outputValues_,
159 inputPointViewType inputPoints_,
161 const ordinal_type order_ )
162 : _outputValues(outputValues_), _inputPoints(inputPoints_), _work(work_),_order(order_){}
164 KOKKOS_INLINE_FUNCTION
165 void operator()(
const size_type iter)
const {
169 const auto ptRange = Kokkos::pair<ordinal_type,ordinal_type>(ptBegin, ptEnd);
170 const auto input = Kokkos::subview( _inputPoints, ptRange, Kokkos::ALL() );
173 case OPERATOR_VALUE : {
174 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange );
181 const auto work = Kokkos::subview( _work, Kokkos::ALL(), ptRange, Kokkos::ALL() );
182 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
196 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
201 INTREPID2_TEST_FOR_ABORT(
true,
202 ">>> ERROR: (Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH::Functor) operator is not supported");
213 template<
typename DeviceType = void,
214 typename outputValueType = double,
215 typename pointValueType =
double>
217 :
public Basis<DeviceType,outputValueType,pointValueType> {
219 typedef double value_type;
236 getValues( OutputViewType outputValues,
237 const PointViewType inputPoints,
238 const EOperator operatorType = OPERATOR_VALUE )
const override {
239 #ifdef HAVE_INTREPID2_DEBUG 247 Impl::Basis_HGRAD_TET_Cn_FEM_ORTH::
248 getValues<DeviceType,numPtsPerEval>( outputValues,
See Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH.
See Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH.
An abstract base class that defines interface for concrete basis implementations for Finite Element (...
See Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH.
ordinal_type getCardinality() const
Returns cardinality of the basis.
Implementation of the default H(grad)-compatible orthogonal basis of arbitrary degree on tetrahedron...
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.
EOperator
Enumeration of primitive operators available in Intrepid. Primitive operators act on reconstructed fu...
Definition file for FEM orthogonal basis functions of degree n for H(grad) functions on TET cells...
void getValues_HGRAD_Args(const outputValueViewType outputValues, const inputPointViewType inputPoints, const EOperator operatorType, const shards::CellTopology cellTopo, const ordinal_type basisCard)
Runtime check of the arguments for the getValues method in an HGRAD-conforming FEM basis...
ordinal_type getDegree() const
Returns the degree of the basis.
shards::CellTopology getBaseCellTopology() const
Returns the base cell topology for which the basis is defined. See Shards documentation https://trili...
Header file for the abstract base class Intrepid2::Basis.
See Intrepid2::Basis_HGRAD_TET_Cn_FEM_ORTH.