49 #ifndef __INTREPID2_HGRAD_LINE_CN_FEM_JACOBI_HPP__
50 #define __INTREPID2_HGRAD_LINE_CN_FEM_JACOBI_HPP__
142 template<EOperator opType>
144 template<
typename OutputViewType,
145 typename inputViewType>
146 KOKKOS_INLINE_FUNCTION
148 getValues( OutputViewType output,
149 const inputViewType input,
150 const ordinal_type order,
153 const ordinal_type opDn = 0 );
156 template<
typename DeviceType, ordinal_type numPtsPerEval,
157 typename outputValueValueType,
class ...outputValueProperties,
158 typename inputPointValueType,
class ...inputPointProperties>
160 getValues(
const typename DeviceType::execution_space& space,
161 Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValues,
162 const Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPoints,
163 const ordinal_type order,
166 const EOperator operatorType );
171 template<
typename outputValueViewType,
172 typename inputPointViewType,
174 ordinal_type numPtsEval>
176 outputValueViewType _outputValues;
177 const inputPointViewType _inputPoints;
178 const ordinal_type _order;
179 const double _alpha, _beta;
180 const ordinal_type _opDn;
182 KOKKOS_INLINE_FUNCTION
183 Functor( outputValueViewType outputValues_,
184 inputPointViewType inputPoints_,
185 const ordinal_type order_,
188 const ordinal_type opDn_ = 0 )
189 : _outputValues(outputValues_), _inputPoints(inputPoints_),
190 _order(order_), _alpha(alpha_), _beta(beta_), _opDn(opDn_) {}
192 KOKKOS_INLINE_FUNCTION
193 void operator()(
const size_type iter)
const {
197 const auto ptRange = Kokkos::pair<ordinal_type,ordinal_type>(ptBegin, ptEnd);
198 const auto input = Kokkos::subview( _inputPoints, ptRange, Kokkos::ALL() );
200 if (input.extent(0)) {
202 case OPERATOR_VALUE : {
203 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange );
207 case OPERATOR_GRAD : {
208 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
213 auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
218 INTREPID2_TEST_FOR_ABORT(
true,
219 ">>> ERROR: (Intrepid2::Basis_HGRAD_LINE_Cn_FEM_JACOBI::Functor) operator is not supported");
229 template<
typename DeviceType = void,
230 typename outputValueType = double,
231 typename pointValueType =
double>
233 :
public Basis<DeviceType,outputValueType,pointValueType> {
235 typedef double value_type;
251 const double alpha = 0,
252 const double beta = 0 );
261 const EOperator operatorType = OPERATOR_VALUE )
const override {
262 #ifdef HAVE_INTREPID2_DEBUG
263 Intrepid2::getValues_HGRAD_Args(outputValues,
269 constexpr ordinal_type numPtsPerEval = 1;
270 Impl::Basis_HGRAD_LINE_Cn_FEM_JACOBI::
271 getValues<DeviceType,numPtsPerEval>(space,
281 double alpha_, beta_;
Kokkos::DynRankView< OutputValueType, Kokkos::LayoutStride, DeviceType > OutputViewType
View type for basis value output.
Kokkos::View< ordinal_type *, typename ExecutionSpace::array_layout, Kokkos::HostSpace > OrdinalTypeArray1DHost
View type for 1d host array.
Kokkos::DynRankView< PointValueType, Kokkos::LayoutStride, DeviceType > PointViewType
View type for input points.
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.
Basis_HGRAD_LINE_Cn_FEM_JACOBI(const ordinal_type order, const double alpha=0, const double beta=0)
Constructor.
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.
Header file for Intrepid2::Polylib class providing orthogonal polynomial calculus and interpolation...
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.
See Intrepid2::Basis_HGRAD_LINE_Cn_FEM_JACOBI.
typename DeviceType::execution_space ExecutionSpace
(Kokkos) Execution space for basis.
See Intrepid2::Basis_HGRAD_LINE_Cn_FEM_JACOBI.
Kokkos::View< ordinal_type **, typename ExecutionSpace::array_layout, Kokkos::HostSpace > OrdinalTypeArray2DHost
View type for 2d host array.
Kokkos::DynRankView< OutputValueType, Kokkos::LayoutStride, DeviceType > OutputViewType
View type for basis value output.
Implementation of the locally H(grad)-compatible FEM basis of variable order on the [-1...
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.
See Intrepid2::Basis_HGRAD_LINE_Cn_FEM_JACOBI.
Definition file for FEM orthogonal basis functions of degree n for H(grad) functions on LINE...
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.