Intrepid2
Intrepid2_HCURL_WEDGE_I1_FEM.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_HCURL_WEDGE_I1_FEM_HPP__
50 #define __INTREPID2_HCURL_WEDGE_I1_FEM_HPP__
51 
52 #include "Intrepid2_Basis.hpp"
53 
54 namespace Intrepid2 {
55 
109  namespace Impl {
110 
115  public:
116  typedef struct Wedge<6> cell_topology_type;
117 
121  template<EOperator opType>
122  struct Serial {
123  template<typename outputViewType,
124  typename inputViewType>
125  KOKKOS_INLINE_FUNCTION
126  static void
127  getValues( outputViewType output,
128  const inputViewType input );
129 
130  };
131 
132  template<typename ExecSpaceType,
133  typename outputValueValueType, class ...outputValueProperties,
134  typename inputPointValueType, class ...inputPointProperties>
135  static void
136  getValues( Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValues,
137  const Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPoints,
138  const EOperator operatorType);
139 
143  template<typename outputValueViewType,
144  typename inputPointViewType,
145  EOperator opType>
146  struct Functor {
147  outputValueViewType _outputValues;
148  const inputPointViewType _inputPoints;
149 
150  KOKKOS_INLINE_FUNCTION
151  Functor( outputValueViewType outputValues_,
152  inputPointViewType inputPoints_ )
153  : _outputValues(outputValues_), _inputPoints(inputPoints_) {}
154 
155  KOKKOS_INLINE_FUNCTION
156  void operator()(const ordinal_type pt) const {
157  switch (opType) {
158  case OPERATOR_VALUE : {
159  auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), pt, Kokkos::ALL() );
160  const auto input = Kokkos::subview( _inputPoints, pt, Kokkos::ALL() );
161  Serial<opType>::getValues( output, input );
162  break;
163  }
164  case OPERATOR_CURL : {
165  auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), pt, Kokkos::ALL() );
166  const auto input = Kokkos::subview( _inputPoints, pt, Kokkos::ALL() );
167  Serial<opType>::getValues( output, input );
168  break;
169  }
170  default: {
171  INTREPID2_TEST_FOR_ABORT( opType != OPERATOR_VALUE &&
172  opType != OPERATOR_CURL,
173  ">>> ERROR: (Intrepid2::Basis_HCURL_WEDGE_I1_FEM::Serial::getValues) operator is not supported");
174  }
175  }
176  }
177  };
178  };
179  }
180 
181  template<typename ExecSpaceType = void,
182  typename outputValueType = double,
183  typename pointValueType = double>
184  class Basis_HCURL_WEDGE_I1_FEM : public Basis<ExecSpaceType, outputValueType, pointValueType> {
185  public:
189 
193 
197 
199 
200  virtual
201  void
202  getValues( outputViewType outputValues,
203  const pointViewType inputPoints,
204  const EOperator operatorType = OPERATOR_VALUE ) const {
205 #ifdef HAVE_INTREPID2_DEBUG
206  // Verify arguments
207  Intrepid2::getValues_HCURL_Args(outputValues,
208  inputPoints,
209  operatorType,
210  this->getBaseCellTopology(),
211  this->getCardinality() );
212 #endif
213  Impl::Basis_HCURL_WEDGE_I1_FEM::
214  getValues<ExecSpaceType>( outputValues,
215  inputPoints,
216  operatorType );
217  }
218 
219  virtual
220  void
221  getDofCoords( scalarViewType dofCoords ) const {
222 #ifdef HAVE_INTREPID2_DEBUG
223  // Verify rank of output array.
224  INTREPID2_TEST_FOR_EXCEPTION( dofCoords.rank() != 2, std::invalid_argument,
225  ">>> ERROR: (Intrepid2::Basis_HCURL_WEDGE_I1_FEM::getDofCoords) rank = 2 required for dofCoords array");
226  // Verify 0th dimension of output array.
227  INTREPID2_TEST_FOR_EXCEPTION( static_cast<ordinal_type>(dofCoords.extent(0)) != this->basisCardinality_, std::invalid_argument,
228  ">>> ERROR: (Intrepid2::Basis_HCURL_WEDGE_I1_FEM::getDofCoords) mismatch in number of dof and 0th dimension of dofCoords array");
229  // Verify 1st dimension of output array.
230  INTREPID2_TEST_FOR_EXCEPTION( dofCoords.extent(1) != this->basisCellTopology_.getDimension(), std::invalid_argument,
231  ">>> ERROR: (Intrepid2::Basis_HCURL_WEDGE_I1_FEM::getDofCoords) incorrect reference cell (1st) dimension in dofCoords array");
232 #endif
233  Kokkos::deep_copy(dofCoords, this->dofCoords_);
234  }
235 
236  virtual
237  const char*
238  getName() const {
239  return "Intrepid2_HCURL_WEDGE_I1_FEM";
240  }
241 
242  virtual
243  bool
245  return true;
246  }
247 
248  };
249 }// namespace Intrepid2
250 
252 
253 #endif
Implementation of the default H(curl)-compatible FEM basis of degree 1 on Wedge cell.
Kokkos::DynRankView< scalarType, Kokkos::LayoutStride, ExecSpaceType > scalarViewType
View type for scalars.
ordinal_type basisCardinality_
Cardinality of the basis, i.e., the number of basis functions/degrees-of-freedom. ...
See Intrepid2::Basis_HCURL_WEDGE_I1_FEM.
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...
Kokkos::DynRankView< scalarType, ExecSpaceType > dofCoords_
Coordinates of degrees-of-freedom for basis functions defined in physical space.
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.
Kokkos::View< ordinal_type *,typename ExecSpaceType::array_layout, Kokkos::HostSpace > ordinal_type_array_1d_host
View type for 1d host array.
virtual bool requireOrientation() const
True if orientation is required.
virtual void getDofCoords(scalarViewType dofCoords) const
Returns spatial locations (coordinates) of degrees of freedom on the reference cell.
Kokkos::DynRankView< pointValueType, Kokkos::LayoutStride, ExecSpaceType > pointViewType
View type for input points.
virtual void getValues(outputViewType outputValues, const pointViewType inputPoints, const EOperator operatorType=OPERATOR_VALUE) const
Evaluation of a FEM basis on a reference cell.
virtual const char * getName() const
Returns basis name.
Definition file for FEM basis functions of degree 1 for H(curl) functions on WEDGE cells...
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.
shards::CellTopology basisCellTopology_
Base topology of the cells for which the basis is defined. See the Shards package for definition of b...