Intrepid2
Intrepid2_HDIV_TRI_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_HDIV_TRI_I1_FEM_HPP__
50 #define __INTREPID2_HDIV_TRI_I1_FEM_HPP__
51 
52 #include "Intrepid2_Basis.hpp"
53 
54 namespace Intrepid2 {
55 
101  namespace Impl {
102 
107  public:
108  typedef struct Triangle<3> cell_topology_type;
112  template<EOperator opType>
113  struct Serial {
114  template<typename outputViewType,
115  typename inputViewType>
116  KOKKOS_INLINE_FUNCTION
117  static void
118  getValues( outputViewType output,
119  const inputViewType input );
120 
121  };
122 
123  template<typename ExecSpaceType,
124  typename outputValueValueType, class ...outputValueProperties,
125  typename inputPointValueType, class ...inputPointProperties>
126  static void
127  getValues( Kokkos::DynRankView<outputValueValueType,outputValueProperties...> outputValues,
128  const Kokkos::DynRankView<inputPointValueType, inputPointProperties...> inputPoints,
129  const EOperator operatorType);
130 
134  template<typename outputValueViewType,
135  typename inputPointViewType,
136  EOperator opType>
137  struct Functor {
138  outputValueViewType _outputValues;
139  const inputPointViewType _inputPoints;
140 
141  KOKKOS_INLINE_FUNCTION
142  Functor( outputValueViewType outputValues_,
143  inputPointViewType inputPoints_ )
144  : _outputValues(outputValues_), _inputPoints(inputPoints_) {}
145 
146  KOKKOS_INLINE_FUNCTION
147  void operator()(const ordinal_type pt) const {
148  switch (opType) {
149  case OPERATOR_VALUE : {
150  auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), pt, Kokkos::ALL() );
151  const auto input = Kokkos::subview( _inputPoints, pt, Kokkos::ALL() );
152  Serial<opType>::getValues( output, input );
153  break;
154  }
155  case OPERATOR_DIV : {
156  auto output = Kokkos::subview( _outputValues, Kokkos::ALL(), pt );
157  const auto input = Kokkos::subview( _inputPoints, pt, Kokkos::ALL() );
158  Serial<opType>::getValues( output, input );
159  break;
160  }
161  default: {
162  INTREPID2_TEST_FOR_ABORT( opType != OPERATOR_VALUE &&
163  opType != OPERATOR_DIV,
164  ">>> ERROR: (Intrepid2::Basis_HDIV_TRI_I1_FEM::Serial::getValues) operator is not supported");
165  }
166  }
167  }
168  };
169  };
170  }
171 
172  template<typename ExecSpaceType = void,
173  typename outputValueType = double,
174  typename pointValueType = double>
175  class Basis_HDIV_TRI_I1_FEM: public Basis<ExecSpaceType,outputValueType,pointValueType> {
176  public:
180 
184 
188 
190 
191  virtual
192  void
193  getValues( outputViewType outputValues,
194  const pointViewType inputPoints,
195  const EOperator operatorType = OPERATOR_VALUE ) const {
196 #ifdef HAVE_INTREPID2_DEBUG
197  // Verify arguments
198  Intrepid2::getValues_HDIV_Args(outputValues,
199  inputPoints,
200  operatorType,
201  this->getBaseCellTopology(),
202  this->getCardinality() );
203 #endif
204  Impl::Basis_HDIV_TRI_I1_FEM::
205  getValues<ExecSpaceType>( outputValues,
206  inputPoints,
207  operatorType );
208  }
209 
210  virtual
211  void
212  getDofCoords( scalarViewType dofCoords ) const {
213 #ifdef HAVE_INTREPID2_DEBUG
214  // Verify rank of output array.
215  INTREPID2_TEST_FOR_EXCEPTION( dofCoords.rank() != 2, std::invalid_argument,
216  ">>> ERROR: (Intrepid2::Basis_HDIV_TRI_I1_FEM::getDofCoords) rank = 2 required for dofCoords array");
217  // Verify 0th dimension of output array.
218  INTREPID2_TEST_FOR_EXCEPTION( static_cast<ordinal_type>(dofCoords.extent(0)) != this->basisCardinality_, std::invalid_argument,
219  ">>> ERROR: (Intrepid2::Basis_HDIV_TRI_I1_FEM::getDofCoords) mismatch in number of dof and 0th dimension of dofCoords array");
220  // Verify 1st dimension of output array.
221  INTREPID2_TEST_FOR_EXCEPTION( dofCoords.extent(1) != this->basisCellTopology_.getDimension(), std::invalid_argument,
222  ">>> ERROR: (Intrepid2::Basis_HDIV_TRI_I1_FEM::getDofCoords) incorrect reference cell (1st) dimension in dofCoords array");
223 #endif
224  Kokkos::deep_copy(dofCoords, this->dofCoords_);
225  }
226 
227  virtual
228  void
229  getDofCoeffs( scalarViewType dofCoeffs ) const {
230 #ifdef HAVE_INTREPID2_DEBUG
231  // Verify rank of output array.
232  INTREPID2_TEST_FOR_EXCEPTION( dofCoeffs.rank() != 2, std::invalid_argument,
233  ">>> ERROR: (Intrepid2::Basis_HDIV_TRI_I1_FEM::getDofCoeffs) rank = 2 required for dofCoeffs array");
234  // Verify 0th dimension of output array.
235  INTREPID2_TEST_FOR_EXCEPTION( static_cast<ordinal_type>(dofCoeffs.extent(0)) != this->getCardinality(), std::invalid_argument,
236  ">>> ERROR: (Intrepid2::Basis_HDIV_TRI_I1_FEM::getDofCoeffs) mismatch in number of dof and 0th dimension of dofCoeffs array");
237  // Verify 1st dimension of output array.
238  INTREPID2_TEST_FOR_EXCEPTION( dofCoeffs.extent(1) != this->getBaseCellTopology().getDimension(), std::invalid_argument,
239  ">>> ERROR: (Intrepid2::Basis_HDIV_TRI_I1_FEM::getDofCoeffs) incorrect reference cell (1st) dimension in dofCoeffs array");
240 #endif
241  Kokkos::deep_copy(dofCoeffs, this->dofCoeffs_);
242  }
243 
244  virtual
245  const char*
246  getName() const {
247  return "Intrepid2_HDIV_TRI_I1_FEM";
248  }
249 
250  virtual
251  bool
253  return true;
254  }
255 
256  };
257 
258 }// namespace Intrepid2
259 
261 
262 #endif
virtual void getDofCoeffs(scalarViewType dofCoeffs) const
Coefficients for computing degrees of freedom for Lagrangian basis If P is an element of the space sp...
Kokkos::DynRankView< scalarType, Kokkos::LayoutStride, ExecSpaceType > scalarViewType
View type for scalars.
void getValues_HDIV_Args(const outputValueViewType outputValues, const inputPointViewType inputPoints, const EOperator operatorType, const shards::CellTopology cellTopo, const ordinal_type basisCard)
Runtime check of the arguments for the getValues method in an HDIV-conforming FEM basis...
ordinal_type basisCardinality_
Cardinality of the basis, i.e., the number of basis functions/degrees-of-freedom. ...
virtual void getValues(outputViewType outputValues, const pointViewType inputPoints, const EOperator operatorType=OPERATOR_VALUE) const
Evaluation of a FEM basis on a reference cell.
See Intrepid2::Basis_HDIV_TRI_I1_FEM.
ordinal_type getCardinality() const
Returns cardinality of the basis.
virtual void getDofCoords(scalarViewType dofCoords) const
Returns spatial locations (coordinates) of degrees of freedom on the reference cell.
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.
Kokkos::DynRankView< scalarType, ExecSpaceType > dofCoords_
Coordinates of degrees-of-freedom for basis functions defined in physical space.
shards::CellTopology getBaseCellTopology() const
Returns the base cell topology for which the basis is defined. See Shards documentation https://trili...
Kokkos::View< ordinal_type ***, typename ExecSpaceType::array_layout, Kokkos::HostSpace > ordinal_type_array_3d_host
View type for 3d host array.
Implementation of the default H(div)-compatible FEM basis of degree 1 on a Triangle cell...
EOperator
Enumeration of primitive operators available in Intrepid. Primitive operators act on reconstructed fu...
Kokkos::View< ordinal_type *,typename ExecSpaceType::array_layout, Kokkos::HostSpace > ordinal_type_array_1d_host
View type for 1d host array.
Kokkos::DynRankView< pointValueType, Kokkos::LayoutStride, ExecSpaceType > pointViewType
View type for input points.
virtual bool requireOrientation() const
True if orientation is required.
Kokkos::DynRankView< scalarType, ExecSpaceType > dofCoeffs_
Coefficients for computing degrees of freedom for Lagrangian basis If P is an element of the space sp...
Definition file for FEM basis functions of degree 1 for H(div) functions on TRI cells.
virtual const char * getName() const
Returns basis name.
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...