49 #ifndef __INTREPID2_HCURL_TET_I1_FEM_DEF_HPP__
50 #define __INTREPID2_HCURL_TET_I1_FEM_DEF_HPP__
57 template<EOperator opType>
58 template<
typename OutputViewType,
59 typename inputViewType>
60 KOKKOS_INLINE_FUNCTION
62 Basis_HCURL_TET_I1_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, 0) = 2.0*(1.0 - y - z);
73 output.access(0, 1) = 2.0*x;
74 output.access(0, 2) = 2.0*x;
76 output.access(1, 0) =-2.0*y;
77 output.access(1, 1) = 2.0*x;
78 output.access(1, 2) = 0.0;
80 output.access(2, 0) = -2.0*y;
81 output.access(2, 1) = 2.0*(-1.0 + x + z);
82 output.access(2, 2) = -2.0*y;
84 output.access(3, 0) = 2.0*z;
85 output.access(3, 1) = 2.0*z;
86 output.access(3, 2) = 2.0*(1.0 - x - y);
88 output.access(4, 0) =-2.0*z;
89 output.access(4, 1) = 0.0;
90 output.access(4, 2) = 2.0*x;
92 output.access(5, 0) = 0.0;
93 output.access(5, 1) =-2.0*z;
94 output.access(5, 2) = 2.0*y;
99 output.access(0, 0) = 0.0;
100 output.access(0, 1) =-4.0;
101 output.access(0, 2) = 4.0;
103 output.access(1, 0) = 0.0;
104 output.access(1, 1) = 0.0;
105 output.access(1, 2) = 4.0;
107 output.access(2, 0) =-4.0;
108 output.access(2, 1) = 0.0;
109 output.access(2, 2) = 4.0;
111 output.access(3, 0) =-4.0;
112 output.access(3, 1) = 4.0;
113 output.access(3, 2) = 0.0;
115 output.access(4, 0) = 0.0;
116 output.access(4, 1) =-4.0;
117 output.access(4, 2) = 0.0;
119 output.access(5, 0) = 4.0;
120 output.access(5, 1) = 0.0;
121 output.access(5, 2) = 0.0;
125 INTREPID2_TEST_FOR_ABORT( opType != OPERATOR_VALUE &&
126 opType != OPERATOR_CURL,
127 ">>> ERROR: (Intrepid2::Basis_HCURL_TET_I1_FEM::Serial::getValues) operator is not supported");
132 template<
typename DT,
133 typename outputValueValueType,
class ...outputValueProperties,
134 typename inputPointValueType,
class ...inputPointProperties>
136 Basis_HCURL_TET_I1_FEM::
137 getValues(
const typename DT::execution_space& space,
138 Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValues,
139 const Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPoints,
140 const EOperator operatorType ) {
141 typedef Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValueViewType;
142 typedef Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPointViewType;
143 typedef typename ExecSpace<typename inputPointViewType::execution_space,typename DT::execution_space>::ExecSpaceType ExecSpaceType;
146 const auto loopSize = inputPoints.extent(0);
147 Kokkos::RangePolicy<ExecSpaceType,Kokkos::Schedule<Kokkos::Static> > policy(space, 0, loopSize);
149 switch (operatorType) {
151 case OPERATOR_VALUE: {
152 typedef Functor<outputValueViewType, inputPointViewType, OPERATOR_VALUE> FunctorType;
153 Kokkos::parallel_for( policy, FunctorType(outputValues, inputPoints) );
156 case OPERATOR_CURL: {
157 typedef Functor<outputValueViewType, inputPointViewType, OPERATOR_CURL> FunctorType;
158 Kokkos::parallel_for( policy, FunctorType(outputValues, inputPoints) );
162 INTREPID2_TEST_FOR_EXCEPTION( (operatorType == OPERATOR_DIV), std::invalid_argument,
163 ">>> ERROR (Basis_HCURL_TET_I1_FEM): DIV is invalid operator for HCURL Basis Functions");
166 case OPERATOR_GRAD: {
167 INTREPID2_TEST_FOR_EXCEPTION( (operatorType == OPERATOR_GRAD), std::invalid_argument,
168 ">>> ERROR (Basis_HCURL_TET_I1_FEM): GRAD is invalid operator for HCURL Basis Functions");
181 INTREPID2_TEST_FOR_EXCEPTION( ( (operatorType == OPERATOR_D1) ||
182 (operatorType == OPERATOR_D2) ||
183 (operatorType == OPERATOR_D3) ||
184 (operatorType == OPERATOR_D4) ||
185 (operatorType == OPERATOR_D5) ||
186 (operatorType == OPERATOR_D6) ||
187 (operatorType == OPERATOR_D7) ||
188 (operatorType == OPERATOR_D8) ||
189 (operatorType == OPERATOR_D9) ||
190 (operatorType == OPERATOR_D10) ),
191 std::invalid_argument,
192 ">>> ERROR (Basis_HCURL_TET_I1_FEM): Invalid operator type");
196 INTREPID2_TEST_FOR_EXCEPTION( ( (operatorType != OPERATOR_VALUE) &&
197 (operatorType != OPERATOR_GRAD) &&
198 (operatorType != OPERATOR_CURL) &&
199 (operatorType != OPERATOR_DIV) &&
200 (operatorType != OPERATOR_D1) &&
201 (operatorType != OPERATOR_D2) &&
202 (operatorType != OPERATOR_D3) &&
203 (operatorType != OPERATOR_D4) &&
204 (operatorType != OPERATOR_D5) &&
205 (operatorType != OPERATOR_D6) &&
206 (operatorType != OPERATOR_D7) &&
207 (operatorType != OPERATOR_D8) &&
208 (operatorType != OPERATOR_D9) &&
209 (operatorType != OPERATOR_D10) ),
210 std::invalid_argument,
211 ">>> ERROR (Basis_HCURL_TET_I1_FEM): Invalid operator type");
220 template<
typename DT,
typename OT,
typename PT >
223 this->basisCardinality_ = 6;
224 this->basisDegree_ = 1;
225 this->basisCellTopology_ = shards::CellTopology(shards::getCellTopologyData<shards::Tetrahedron<4> >() );
226 this->basisType_ = BASIS_FEM_DEFAULT;
227 this->basisCoordinates_ = COORDINATES_CARTESIAN;
228 this->functionSpace_ = FUNCTION_SPACE_HCURL;
233 const ordinal_type tagSize = 4;
234 const ordinal_type posScDim = 0;
235 const ordinal_type posScOrd = 1;
236 const ordinal_type posDfOrd = 2;
239 ordinal_type tags[24] = {
251 this->setOrdinalTagData(this->tagToOrdinal_,
254 this->basisCardinality_,
261 Kokkos::DynRankView<typename ScalarViewType::value_type,typename DT::execution_space::array_layout,Kokkos::HostSpace>
262 dofCoords(
"dofCoordsHost", this->basisCardinality_,this->basisCellTopology_.getDimension());
264 dofCoords(0,0) = 0.5; dofCoords(0,1) = 0.0; dofCoords(0,2) = 0.0;
265 dofCoords(1,0) = 0.5; dofCoords(1,1) = 0.5; dofCoords(1,2) = 0.0;
266 dofCoords(2,0) = 0.0; dofCoords(2,1) = 0.5; dofCoords(2,2) = 0.0;
267 dofCoords(3,0) = 0.0; dofCoords(3,1) = 0.0; dofCoords(3,2) = 0.5;
268 dofCoords(4,0) = 0.5; dofCoords(4,1) = 0.0; dofCoords(4,2) = 0.5;
269 dofCoords(5,0) = 0.0; dofCoords(5,1) = 0.5; dofCoords(5,2) = 0.5;
271 this->dofCoords_ = Kokkos::create_mirror_view(
typename DT::memory_space(), dofCoords);
272 Kokkos::deep_copy(this->dofCoords_, dofCoords);
275 Kokkos::DynRankView<typename ScalarViewType::value_type,typename DT::execution_space::array_layout,Kokkos::HostSpace>
276 dofCoeffs(
"dofCoeffsHost", this->basisCardinality_,this->basisCellTopology_.getDimension());
278 dofCoeffs(0,0) = 0.5; dofCoeffs(0,1) = 0.0; dofCoeffs(0,2) = 0.0;
279 dofCoeffs(1,0) = -0.5; dofCoeffs(1,1) = 0.5; dofCoeffs(1,2) = 0.0;
280 dofCoeffs(2,0) = 0.0; dofCoeffs(2,1) = -0.5; dofCoeffs(2,2) = 0.0;
281 dofCoeffs(3,0) = 0.0; dofCoeffs(3,1) = 0.0; dofCoeffs(3,2) = 0.5;
282 dofCoeffs(4,0) = -0.5; dofCoeffs(4,1) = 0.0; dofCoeffs(4,2) = 0.5;
283 dofCoeffs(5,0) = 0.0; dofCoeffs(5,1) = -0.5; dofCoeffs(5,2) = 0.5;
285 this->dofCoeffs_ = Kokkos::create_mirror_view(
typename DT::memory_space(), dofCoeffs);
286 Kokkos::deep_copy(this->dofCoeffs_, dofCoeffs);
Basis_HCURL_TET_I1_FEM()
Constructor.
Kokkos::View< ordinal_type *, typename ExecutionSpace::array_layout, Kokkos::HostSpace > OrdinalTypeArray1DHost
View type for 1d host array.