Intrepid2
Intrepid2_CubatureTensorPyr.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_CUBATURE_TENSOR_PYR_HPP__
50#define __INTREPID2_CUBATURE_TENSOR_PYR_HPP__
51
53
54namespace Intrepid2 {
55
59 template<typename DeviceType = void,
60 typename pointValueType = double,
61 typename weightValueType = double>
63 : public CubatureTensor<DeviceType,pointValueType,weightValueType> {
64 public:
65
66 template<typename cubPointViewType,
67 typename cubWeightViewType>
68 struct Functor {
69 cubPointViewType _cubPoints;
70 cubWeightViewType _cubWeights;
71
72 KOKKOS_INLINE_FUNCTION
73 Functor( cubPointViewType cubPoints_,
74 cubWeightViewType cubWeights_)
75 : _cubPoints(cubPoints_), _cubWeights(cubWeights_) {}
76
77 KOKKOS_INLINE_FUNCTION
78 void operator()(const ordinal_type i) const {
79 const auto tmp = 0.5*(1.0 - _cubPoints(i,2));
80 _cubPoints(i,0) *= tmp;
81 _cubPoints(i,1) *= tmp;
82 _cubPoints(i,2) = 1.0 - tmp;
83
84 _cubWeights(i) /= 8.0; // when line jacobi20 is used (exact)
85 //_cubWeights(i) *= (tmp*tmp*.5); // when gauss integration rule is used (need over-integration)
86 }
87 };
88
89 template<typename cubPointValueType, class ...cubPointProperties,
90 typename cubWeightValueType, class ...cubWeightProperties>
91 void
92 getCubatureImpl( Kokkos::DynRankView<cubPointValueType, cubPointProperties...> cubPoints,
93 Kokkos::DynRankView<cubWeightValueType,cubWeightProperties...> cubWeights ) const;
94
95 using typename Cubature<DeviceType,pointValueType,weightValueType>::PointViewType;
96 using typename Cubature<DeviceType,pointValueType,weightValueType>::weightViewType;
97
98 using CubatureTensor<DeviceType,pointValueType,weightValueType>::getCubature;
99 using typename CubatureTensor<DeviceType,pointValueType,weightValueType>::ExecSpaceType;
100
101 virtual
102 void
103 getCubature( PointViewType cubPoints,
104 weightViewType cubWeights ) const {
105 getCubatureImpl( cubPoints,
106 cubWeights );
107 }
108
110 : CubatureTensor<DeviceType,pointValueType,weightValueType>() {}
111
112 CubatureTensorPyr(const CubatureTensorPyr &b)
113 : CubatureTensor<DeviceType,pointValueType,weightValueType>(b) {}
114
115 template<typename CubatureLineType>
116 CubatureTensorPyr( const CubatureLineType line )
117 : CubatureTensor<DeviceType,pointValueType,weightValueType>(line, line, line) {}
118
119 template<typename CubatureLineType0,
120 typename CubatureLineType1,
121 typename CubatureLineType2>
122 CubatureTensorPyr( const CubatureLineType0 line0,
123 const CubatureLineType1 line1,
124 const CubatureLineType2 line2 )
125 : CubatureTensor<DeviceType,pointValueType,weightValueType>(line0, line1, line2) {}
126
127 };
128}
129
131
132#endif
Definition file for the Intrepid2::CubatureTensorPyr class.
Header file for the Intrepid2::CubatureTensor class.
Defines tensor-product cubature (integration) rules in Intrepid.
virtual void getCubature(PointViewType cubPoints, weightViewType cubWeights) const
Returns cubature points and weights (return arrays must be pre-sized/pre-allocated).
Defines tensor-product cubature (integration) rules in Intrepid.