49 #ifndef __INTREPID2_HGRAD_HEX_C1_FEM_DEF_HPP__
50 #define __INTREPID2_HGRAD_HEX_C1_FEM_DEF_HPP__
57 template<EOperator opType>
58 template<
typename outputViewType,
59 typename inputViewType>
60 KOKKOS_INLINE_FUNCTION
62 Basis_HGRAD_HEX_C1_FEM::Serial<opType>::
63 getValues( outputViewType output,
64 const inputViewType input ) {
66 case OPERATOR_VALUE : {
67 const auto x = input(0);
68 const auto y = input(1);
69 const auto z = input(2);
72 output.access(0) = (1.0 - x)*(1.0 - y)*(1.0 - z)/8.0;
73 output.access(1) = (1.0 + x)*(1.0 - y)*(1.0 - z)/8.0;
74 output.access(2) = (1.0 + x)*(1.0 + y)*(1.0 - z)/8.0;
75 output.access(3) = (1.0 - x)*(1.0 + y)*(1.0 - z)/8.0;
77 output.access(4) = (1.0 - x)*(1.0 - y)*(1.0 + z)/8.0;
78 output.access(5) = (1.0 + x)*(1.0 - y)*(1.0 + z)/8.0;
79 output.access(6) = (1.0 + x)*(1.0 + y)*(1.0 + z)/8.0;
80 output.access(7) = (1.0 - x)*(1.0 + y)*(1.0 + z)/8.0;
83 case OPERATOR_GRAD : {
84 const auto x = input(0);
85 const auto y = input(1);
86 const auto z = input(2);
89 output.access(0, 0) = -(1.0 - y)*(1.0 - z)/8.0;
90 output.access(0, 1) = -(1.0 - x)*(1.0 - z)/8.0;
91 output.access(0, 2) = -(1.0 - x)*(1.0 - y)/8.0;
93 output.access(1, 0) = (1.0 - y)*(1.0 - z)/8.0;
94 output.access(1, 1) = -(1.0 + x)*(1.0 - z)/8.0;
95 output.access(1, 2) = -(1.0 + x)*(1.0 - y)/8.0;
97 output.access(2, 0) = (1.0 + y)*(1.0 - z)/8.0;
98 output.access(2, 1) = (1.0 + x)*(1.0 - z)/8.0;
99 output.access(2, 2) = -(1.0 + x)*(1.0 + y)/8.0;
101 output.access(3, 0) = -(1.0 + y)*(1.0 - z)/8.0;
102 output.access(3, 1) = (1.0 - x)*(1.0 - z)/8.0;
103 output.access(3, 2) = -(1.0 - x)*(1.0 + y)/8.0;
105 output.access(4, 0) = -(1.0 - y)*(1.0 + z)/8.0;
106 output.access(4, 1) = -(1.0 - x)*(1.0 + z)/8.0;
107 output.access(4, 2) = (1.0 - x)*(1.0 - y)/8.0;
109 output.access(5, 0) = (1.0 - y)*(1.0 + z)/8.0;
110 output.access(5, 1) = -(1.0 + x)*(1.0 + z)/8.0;
111 output.access(5, 2) = (1.0 + x)*(1.0 - y)/8.0;
113 output.access(6, 0) = (1.0 + y)*(1.0 + z)/8.0;
114 output.access(6, 1) = (1.0 + x)*(1.0 + z)/8.0;
115 output.access(6, 2) = (1.0 + x)*(1.0 + y)/8.0;
117 output.access(7, 0) = -(1.0 + y)*(1.0 + z)/8.0;
118 output.access(7, 1) = (1.0 - x)*(1.0 + z)/8.0;
119 output.access(7, 2) = (1.0 - x)*(1.0 + y)/8.0;
123 const auto x = input(0);
124 const auto y = input(1);
125 const auto z = input(2);
128 output.access(0, 0) = 0.0;
129 output.access(0, 1) = (1.0 - z)/8.0;
130 output.access(0, 2) = (1.0 - y)/8.0;
131 output.access(0, 3) = 0.0;
132 output.access(0, 4) = (1.0 - x)/8.0;
133 output.access(0, 5) = 0.0;
135 output.access(1, 0) = 0.0;
136 output.access(1, 1) = -(1.0 - z)/8.0;
137 output.access(1, 2) = -(1.0 - y)/8.0;
138 output.access(1, 3) = 0.0;
139 output.access(1, 4) = (1.0 + x)/8.0;
140 output.access(1, 5) = 0.0;
142 output.access(2, 0) = 0.0;
143 output.access(2, 1) = (1.0 - z)/8.0;
144 output.access(2, 2) = -(1.0 + y)/8.0;
145 output.access(2, 3) = 0.0;
146 output.access(2, 4) = -(1.0 + x)/8.0;
147 output.access(2, 5) = 0.0;
149 output.access(3, 0) = 0.0;
150 output.access(3, 1) = -(1.0 - z)/8.0;
151 output.access(3, 2) = (1.0 + y)/8.0;
152 output.access(3, 3) = 0.0;
153 output.access(3, 4) = -(1.0 - x)/8.0;
154 output.access(3, 5) = 0.0;
157 output.access(4, 0) = 0.0;
158 output.access(4, 1) = (1.0 + z)/8.0;
159 output.access(4, 2) = -(1.0 - y)/8.0;
160 output.access(4, 3) = 0.0;
161 output.access(4, 4) = -(1.0 - x)/8.0;
162 output.access(4, 5) = 0.0;
164 output.access(5, 0) = 0.0;
165 output.access(5, 1) = -(1.0 + z)/8.0;
166 output.access(5, 2) = (1.0 - y)/8.0;
167 output.access(5, 3) = 0.0;
168 output.access(5, 4) = -(1.0 + x)/8.0;
169 output.access(5, 5) = 0.0;
171 output.access(6, 0) = 0.0;
172 output.access(6, 1) = (1.0 + z)/8.0;
173 output.access(6, 2) = (1.0 + y)/8.0;
174 output.access(6, 3) = 0.0;
175 output.access(6, 4) = (1.0 + x)/8.0;
176 output.access(6, 5) = 0.0;
178 output.access(7, 0) = 0.0;
179 output.access(7, 1) = -(1.0 + z)/8.0;
180 output.access(7, 2) = -(1.0 + y)/8.0;
181 output.access(7, 3) = 0.0;
182 output.access(7, 4) = (1.0 - x)/8.0;
183 output.access(7, 5) = 0.0;
186 case OPERATOR_MAX : {
187 const ordinal_type jend = output.extent(1);
188 const ordinal_type iend = output.extent(0);
190 for (ordinal_type j=0;j<jend;++j)
191 for (ordinal_type i=0;i<iend;++i)
192 output.access(i, j) = 0.0;
196 INTREPID2_TEST_FOR_ABORT( opType != OPERATOR_VALUE &&
197 opType != OPERATOR_GRAD &&
198 opType != OPERATOR_CURL &&
199 opType != OPERATOR_D2 &&
200 opType != OPERATOR_MAX,
201 ">>> ERROR: (Intrepid2::Basis_HGRAD_HEX_C1_FEM::Serial::getValues) operator is not supported");
207 template<
typename SpT,
208 typename outputValueValueType,
class ...outputValueProperties,
209 typename inputPointValueType,
class ...inputPointProperties>
211 Basis_HGRAD_HEX_C1_FEM::
212 getValues( Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValues,
213 const Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPoints,
214 const EOperator operatorType ) {
215 typedef Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValueViewType;
216 typedef Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPointViewType;
217 typedef typename ExecSpace<typename inputPointViewType::execution_space,SpT>::ExecSpaceType ExecSpaceType;
220 const auto loopSize = inputPoints.extent(0);
221 Kokkos::RangePolicy<ExecSpaceType,Kokkos::Schedule<Kokkos::Static> > policy(0, loopSize);
223 switch (operatorType) {
225 case OPERATOR_VALUE: {
226 typedef Functor<outputValueViewType,inputPointViewType,OPERATOR_VALUE> FunctorType;
227 Kokkos::parallel_for( policy, FunctorType(outputValues, inputPoints) );
232 typedef Functor<outputValueViewType,inputPointViewType,OPERATOR_GRAD> FunctorType;
233 Kokkos::parallel_for( policy, FunctorType(outputValues, inputPoints) );
236 case OPERATOR_CURL: {
237 INTREPID2_TEST_FOR_EXCEPTION( operatorType == OPERATOR_CURL, std::invalid_argument,
238 ">>> ERROR (Basis_HGRAD_HEX_C1_FEM): CURL is invalid operator for rank-0 (scalar) functions in 3D");
243 INTREPID2_TEST_FOR_EXCEPTION( (operatorType == OPERATOR_DIV), std::invalid_argument,
244 ">>> ERROR (Basis_HGRAD_HEX_C1_FEM): DIV is invalid operator for rank-0 (scalar) functions in 3D");
249 typedef Functor<outputValueViewType,inputPointViewType,OPERATOR_D2> FunctorType;
250 Kokkos::parallel_for( policy, FunctorType(outputValues, inputPoints) );
261 typedef Functor<outputValueViewType,inputPointViewType,OPERATOR_MAX> FunctorType;
262 Kokkos::parallel_for( policy, FunctorType(outputValues, inputPoints) );
266 INTREPID2_TEST_FOR_EXCEPTION( !( Intrepid2::isValidOperator(operatorType) ), std::invalid_argument,
267 ">>> ERROR (Basis_HGRAD_HEX_C1_FEM): Invalid operator type");
275 template<
typename SpT,
typename OT,
typename PT>
278 this->basisCardinality_ = 8;
279 this->basisDegree_ = 1;
280 this->basisCellTopology_ = shards::CellTopology(shards::getCellTopologyData<shards::Hexahedron<8> >() );
281 this->basisType_ = BASIS_FEM_DEFAULT;
282 this->basisCoordinates_ = COORDINATES_CARTESIAN;
287 const ordinal_type tagSize = 4;
288 const ordinal_type posScDim = 0;
289 const ordinal_type posScOrd = 1;
290 const ordinal_type posDfOrd = 2;
293 ordinal_type tags[32] = { 0, 0, 0, 1,
308 this->setOrdinalTagData(this->tagToOrdinal_,
311 this->basisCardinality_,
325 Kokkos::DynRankView<typename scalarViewType::value_type,typename SpT::array_layout,Kokkos::HostSpace>
326 dofCoords(
"dofCoordsHost", this->basisCardinality_,this->basisCellTopology_.getDimension());
328 dofCoords(0,0) = -1.0; dofCoords(0,1) = -1.0; dofCoords(0,2) = -1.0;
329 dofCoords(1,0) = 1.0; dofCoords(1,1) = -1.0; dofCoords(1,2) = -1.0;
330 dofCoords(2,0) = 1.0; dofCoords(2,1) = 1.0; dofCoords(2,2) = -1.0;
331 dofCoords(3,0) = -1.0; dofCoords(3,1) = 1.0; dofCoords(3,2) = -1.0;
332 dofCoords(4,0) = -1.0; dofCoords(4,1) = -1.0; dofCoords(4,2) = 1.0;
333 dofCoords(5,0) = 1.0; dofCoords(5,1) = -1.0; dofCoords(5,2) = 1.0;
334 dofCoords(6,0) = 1.0; dofCoords(6,1) = 1.0; dofCoords(6,2) = 1.0;
335 dofCoords(7,0) = -1.0; dofCoords(7,1) = 1.0; dofCoords(7,2) = 1.0;
337 this->dofCoords_ = Kokkos::create_mirror_view(
typename SpT::memory_space(), dofCoords);
338 Kokkos::deep_copy(this->dofCoords_, dofCoords);
Kokkos::View< ordinal_type *,typename ExecSpaceType::array_layout, Kokkos::HostSpace > ordinal_type_array_1d_host
View type for 1d host array.
Basis_HGRAD_HEX_C1_FEM()
Constructor.