49 #ifndef __INTREPID2_HGRAD_TRI_CN_FEM_ORTH_HPP__
50 #define __INTREPID2_HGRAD_TRI_CN_FEM_ORTH_HPP__
70 template<
typename OutputViewType,
71 typename inputViewType,
72 typename workViewType,
73 bool hasDeriv, ordinal_type n>
75 KOKKOS_INLINE_FUNCTION
77 generate( OutputViewType output,
78 const inputViewType input,
80 const ordinal_type p );
86 template<
typename OutputViewType,
87 typename inputViewType,
88 typename workViewType,
91 KOKKOS_INLINE_FUNCTION
93 generate( OutputViewType output,
94 const inputViewType input,
96 const ordinal_type p );
102 template<
typename OutputViewType,
103 typename inputViewType,
104 typename workViewType,
107 KOKKOS_INLINE_FUNCTION
109 generate( OutputViewType output,
110 const inputViewType input,
112 const ordinal_type p );
124 template<EOperator opType>
126 template<
typename OutputViewType,
127 typename inputViewType,
128 typename workViewType>
129 KOKKOS_INLINE_FUNCTION
131 getValues( OutputViewType output,
132 const inputViewType input,
134 const ordinal_type order);
137 template<
typename DeviceType, ordinal_type numPtsPerEval,
138 typename outputValueValueType,
class ...outputValueProperties,
139 typename inputPointValueType,
class ...inputPointProperties>
141 getValues(
const typename DeviceType::execution_space& space,
142 Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValues,
143 const Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPoints,
144 const ordinal_type order,
145 const EOperator operatorType );
150 template<
typename outputValueViewType,
151 typename inputPointViewType,
152 typename workViewType,
154 ordinal_type numPtsEval>
156 outputValueViewType _outputValues;
157 const inputPointViewType _inputPoints;
159 const ordinal_type _order;
161 KOKKOS_INLINE_FUNCTION
162 Functor( outputValueViewType outputValues_,
163 inputPointViewType inputPoints_,
165 const ordinal_type order_ )
166 : _outputValues(outputValues_), _inputPoints(inputPoints_), _work(work_),_order(order_){}
168 KOKKOS_INLINE_FUNCTION
169 void operator()(
const size_type iter)
const {
173 const auto ptRange = Kokkos::pair<ordinal_type,ordinal_type>(ptBegin, ptEnd);
174 const auto input = Kokkos::subview( _inputPoints, ptRange, Kokkos::ALL() );
177 case OPERATOR_VALUE : {
178 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange );
185 const auto work = Kokkos::subview( _work, Kokkos::ALL(), ptRange, Kokkos::ALL() );
186 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
200 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
205 INTREPID2_TEST_FOR_ABORT(
true,
206 ">>> ERROR: (Intrepid2::Basis_HGRAD_TRI_Cn_FEM_ORTH::Functor) operator is not supported");
217 template<
typename DeviceType = void,
218 typename outputValueType = double,
219 typename pointValueType =
double>
221 :
public Basis<DeviceType,outputValueType,pointValueType> {
223 typedef double value_type;
247 const EOperator operatorType = OPERATOR_VALUE )
const override {
248 #ifdef HAVE_INTREPID2_DEBUG
249 Intrepid2::getValues_HGRAD_Args(outputValues,
256 Impl::Basis_HGRAD_TRI_Cn_FEM_ORTH::
257 getValues<DeviceType,numPtsPerEval>(space,
Kokkos::View< ordinal_type *, typename ExecutionSpace::array_layout, Kokkos::HostSpace > OrdinalTypeArray1DHost
View type for 1d host array.
ordinal_type getCardinality() const
Returns cardinality of the basis.
virtual void getValues(const ExecutionSpace &, OutputViewType, const PointViewType, const EOperator=OPERATOR_VALUE) const
Evaluation of a FEM basis on a reference cell.
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...
Implementation of the default H(grad)-compatible orthogonal basis (Dubiner) of arbitrary degree on tr...
Kokkos::DynRankView< scalarType, Kokkos::LayoutStride, DeviceType > ScalarViewType
View type for scalars.
Kokkos::DynRankView< PointValueType, Kokkos::LayoutStride, DeviceType > PointViewType
View type for input points.
See Intrepid2::Basis_HGRAD_TRI_Cn_FEM_ORTH.
See Intrepid2::Basis_HGRAD_TRI_Cn_FEM_ORTH.
Kokkos::View< ordinal_type **, typename ExecutionSpace::array_layout, Kokkos::HostSpace > OrdinalTypeArray2DHost
View type for 2d host array.
typename DeviceType::execution_space ExecutionSpace
(Kokkos) Execution space for basis.
static constexpr ordinal_type MaxNumPtsPerBasisEval
The maximum number of points to eval in serial mode.
Kokkos::DynRankView< OutputValueType, Kokkos::LayoutStride, DeviceType > OutputViewType
View type for basis value output.
See Intrepid2::Basis_HGRAD_TRI_Cn_FEM_ORTH.
Definition file for FEM orthogonal basis functions of arbitrary degree for H(grad) functions on TRI...
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.
Basis_HGRAD_TRI_Cn_FEM_ORTH(const ordinal_type order)
Constructor.
Kokkos::View< ordinal_type ***, typename ExecutionSpace::array_layout, Kokkos::HostSpace > OrdinalTypeArray3DHost
View type for 3d host array.
Kokkos::DynRankView< OutputValueType, Kokkos::LayoutStride, DeviceType > OutputViewType
View type for basis value output.
ordinal_type getDegree() const
Returns the degree of the basis.
Kokkos::DynRankView< PointValueType, Kokkos::LayoutStride, DeviceType > PointViewType
View type for input points.
See Intrepid2::Basis_HGRAD_TRI_Cn_FEM_ORTH.
Header file for the abstract base class Intrepid2::Basis.
typename DeviceType::execution_space ExecutionSpace
(Kokkos) Execution space for basis.