Intrepid2
Intrepid2_HGRAD_TRI_Cn_FEM_ORTH.hpp
Go to the documentation of this file.
1 // @HEADER
2 // ************************************************************************
3 //
4 // Intrepid2 Package
5 // Copyright (2007) Sandia Corporation
6 //
7 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
8 // license for use of this work by or on behalf of the U.S. Government.
9 //
10 // Redistribution and use in source and binary forms, with or without
11 // modification, are permitted provided that the following conditions are
12 // met:
13 //
14 // 1. Redistributions of source code must retain the above copyright
15 // notice, this list of conditions and the following disclaimer.
16 //
17 // 2. Redistributions in binary form must reproduce the above copyright
18 // notice, this list of conditions and the following disclaimer in the
19 // documentation and/or other materials provided with the distribution.
20 //
21 // 3. Neither the name of the Corporation nor the names of the
22 // contributors may be used to endorse or promote products derived from
23 // this software without specific prior written permission.
24 //
25 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36 //
37 // Questions? Contact Kyungjoo Kim (kyukim@sandia.gov), or
38 // Mauro Perego (mperego@sandia.gov)
39 //
40 // ************************************************************************
41 // @HEADER
42 
49 #ifndef __INTREPID2_HGRAD_TRI_CN_FEM_ORTH_HPP__
50 #define __INTREPID2_HGRAD_TRI_CN_FEM_ORTH_HPP__
51 
52 #include "Intrepid2_Basis.hpp"
53 
54 namespace Intrepid2 {
55 
65 namespace Impl {
66 
70 template<typename outputViewType,
71 typename inputViewType,
72 typename workViewType,
73 bool hasDeriv, ordinal_type n>
75  KOKKOS_INLINE_FUNCTION
76  static void
77  generate( outputViewType output,
78  const inputViewType input,
79  workViewType work,
80  const ordinal_type p );
81 };
82 
86 template<typename outputViewType,
87 typename inputViewType,
88 typename workViewType,
89 bool hasDeriv>
90 struct OrthPolynomialTri<outputViewType,inputViewType,workViewType,hasDeriv,0> {
91  KOKKOS_INLINE_FUNCTION
92  static void
93  generate( outputViewType output,
94  const inputViewType input,
95  workViewType work,
96  const ordinal_type p );
97 };
98 
102 template<typename outputViewType,
103 typename inputViewType,
104 typename workViewType,
105 bool hasDeriv>
106 struct OrthPolynomialTri<outputViewType,inputViewType,workViewType,hasDeriv,1> {
107  KOKKOS_INLINE_FUNCTION
108  static void
109  generate( outputViewType output,
110  const inputViewType input,
111  workViewType work,
112  const ordinal_type p );
113 };
114 
119 public:
120 
124  template<EOperator opType>
125  struct Serial {
126  template<typename outputViewType,
127  typename inputViewType,
128  typename workViewType>
129  KOKKOS_INLINE_FUNCTION
130  static void
131  getValues( outputViewType output,
132  const inputViewType input,
133  workViewType work,
134  const ordinal_type order);
135  };
136 
137  template<typename ExecSpaceType, ordinal_type numPtsPerEval,
138  typename outputValueValueType, class ...outputValueProperties,
139  typename inputPointValueType, class ...inputPointProperties>
140  static void
141  getValues( Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValues,
142  const Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPoints,
143  const ordinal_type order,
144  const EOperator operatorType );
145 
149  template<typename outputValueViewType,
150  typename inputPointViewType,
151  typename workViewType,
152  EOperator opType,
153  ordinal_type numPtsEval>
154  struct Functor {
155  outputValueViewType _outputValues;
156  const inputPointViewType _inputPoints;
157  workViewType _work;
158  const ordinal_type _order;
159 
160  KOKKOS_INLINE_FUNCTION
161  Functor( outputValueViewType outputValues_,
162  inputPointViewType inputPoints_,
163  workViewType work_,
164  const ordinal_type order_ )
165  : _outputValues(outputValues_), _inputPoints(inputPoints_), _work(work_),_order(order_){}
166 
167  KOKKOS_INLINE_FUNCTION
168  void operator()(const size_type iter) const {
169  const auto ptBegin = Util<ordinal_type>::min(iter*numPtsEval, _inputPoints.extent(0));
170  const auto ptEnd = Util<ordinal_type>::min(ptBegin+numPtsEval, _inputPoints.extent(0));
171 
172  const auto ptRange = Kokkos::pair<ordinal_type,ordinal_type>(ptBegin, ptEnd);
173  const auto input = Kokkos::subview( _inputPoints, ptRange, Kokkos::ALL() );
174 
175  switch (opType) {
176  case OPERATOR_VALUE : {
177  auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange );
178  Serial<opType>::getValues( output, input, _work, _order ); //here work is not used
179  break;
180  }
181  case OPERATOR_GRAD :
182  case OPERATOR_D1 :
183  {
184  const auto work = Kokkos::subview( _work, Kokkos::ALL(), ptRange, Kokkos::ALL() );
185  auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
186  Serial<opType>::getValues( output, input, work, _order);
187  break;
188  }
189  case OPERATOR_D2 :
190  case OPERATOR_D3 :
191  case OPERATOR_D4 :
192  case OPERATOR_D5 :
193  case OPERATOR_D6 :
194  case OPERATOR_D7 :
195  case OPERATOR_D8 :
196  case OPERATOR_D9 :
197  case OPERATOR_D10 :
198  {
199  auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), ptRange, Kokkos::ALL() );
200  Serial<opType>::getValues( output, input, _work, _order); //here work is not used
201  break;
202  }
203  default: {
204  INTREPID2_TEST_FOR_ABORT( true,
205  ">>> ERROR: (Intrepid2::Basis_HGRAD_TRI_Cn_FEM_ORTH::Functor) operator is not supported");
206 
207  }
208  }
209  }
210  };
211 
212 };
213 
214 }
215 
216 template<typename ExecSpaceType = void,
217  typename outputValueType = double,
218  typename pointValueType = double>
220  : public Basis<ExecSpaceType,outputValueType,pointValueType> {
221  public:
222  typedef double value_type;
226 
229  Basis_HGRAD_TRI_Cn_FEM_ORTH( const ordinal_type order );
230 
233 
235 
236  virtual
237  void
238  getValues( outputViewType outputValues,
239  const pointViewType inputPoints,
240  const EOperator operatorType = OPERATOR_VALUE ) const {
241  #ifdef HAVE_INTREPID2_DEBUG
242  Intrepid2::getValues_HGRAD_Args(outputValues,
243  inputPoints,
244  operatorType,
245  this->getBaseCellTopology(),
246  this->getCardinality() );
247  #endif
248  constexpr ordinal_type numPtsPerEval = Parameters::MaxNumPtsPerBasisEval;
249  Impl::Basis_HGRAD_TRI_Cn_FEM_ORTH::
250  getValues<ExecSpaceType,numPtsPerEval>( outputValues,
251  inputPoints,
252  this->getDegree(),
253  operatorType );
254  }
255 };
256 
257 }// namespace Intrepid2
258 
260 
261 #endif
262 
small utility functions
An abstract base class that defines interface for concrete basis implementations for Finite Element (...
Kokkos::DynRankView< outputValueType, Kokkos::LayoutStride, ExecSpaceType > outputViewType
View type for basis value output.
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...
ordinal_type getCardinality() const
Returns cardinality of the basis.
Kokkos::View< ordinal_type ***, typename ExecSpaceType::array_layout, Kokkos::HostSpace > ordinal_type_array_3d_host
View type for 3d host array.
See Intrepid2::Basis_HGRAD_TRI_Cn_FEM_ORTH.
See Intrepid2::Basis_HGRAD_TRI_Cn_FEM_ORTH.
static constexpr ordinal_type MaxNumPtsPerBasisEval
The maximum number of points to eval in serial mode.
Kokkos::View< ordinal_type *,typename ExecSpaceType::array_layout, Kokkos::HostSpace > ordinal_type_array_1d_host
View type for 1d host array.
Definition file for FEM orthogonal basis functions of arbitrary degree for H(grad) functions on TRI...
Kokkos::DynRankView< pointValueType, Kokkos::LayoutStride, ExecSpaceType > pointViewType
View type for input points.
ordinal_type getDegree() const
Returns the degree of the basis.
virtual void getValues(outputViewType outputValues, const pointViewType inputPoints, const EOperator operatorType=OPERATOR_VALUE) const
Evaluation of a FEM basis on a reference cell.
Header file for the abstract base class Intrepid2::Basis.
Kokkos::View< ordinal_type **,typename ExecSpaceType::array_layout, Kokkos::HostSpace > ordinal_type_array_2d_host
View type for 2d host array.
Basis_HGRAD_TRI_Cn_FEM_ORTH(const ordinal_type order)
Constructor.