Intrepid2
Intrepid2_DerivedBasis_HVOL_HEX.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),
38 // Mauro Perego (mperego@sandia.gov), or
39 // Nate Roberts (nvrober@sandia.gov)
40 //
41 // ************************************************************************
42 // @HEADER
43 
57 #ifndef Intrepid2_DerivedBasis_HVOL_HEX_h
58 #define Intrepid2_DerivedBasis_HVOL_HEX_h
59 
62 
63 namespace Intrepid2
64 {
69  template<class HVOL_LINE>
71  :
72  public Basis_TensorBasis<typename HVOL_LINE::BasisBase>
73  // TODO: make this a subclass of TensorBasis3 instead, following what we've done for H(curl) and H(div)
74  {
75  std::string name_;
76  ordinal_type polyOrder_x_, polyOrder_y_, polyOrder_z_;
77  EPointType pointType_;
78  public:
79  using ExecutionSpace = typename HVOL_LINE::ExecutionSpace;
80  using OutputValueType = typename HVOL_LINE::OutputValueType;
81  using PointValueType = typename HVOL_LINE::PointValueType;
82 
83  using OutputViewType = typename HVOL_LINE::OutputViewType;
84  using PointViewType = typename HVOL_LINE::PointViewType ;
85  using ScalarViewType = typename HVOL_LINE::ScalarViewType;
86 
87  using BasisBase = typename HVOL_LINE::BasisBase;
88 
89  using LineBasis = HVOL_LINE;
92 
99  Basis_Derived_HVOL_HEX(int polyOrder_x, int polyOrder_y, int polyOrder_z, const EPointType pointType=POINTTYPE_DEFAULT)
100  :
101  TensorBasis(Teuchos::rcp(new QuadBasis(polyOrder_x,polyOrder_y,pointType)),
102  Teuchos::rcp(new LineBasis(polyOrder_z,pointType))),
103  polyOrder_x_(polyOrder_x),
104  polyOrder_y_(polyOrder_y),
105  polyOrder_z_(polyOrder_z),
106  pointType_(pointType)
107  {
108  this->functionSpace_ = FUNCTION_SPACE_HVOL;
109 
110  std::ostringstream basisName;
111  basisName << "HVOL_HEX (" << this->TensorBasis::getName() << ")";
112  name_ = basisName.str();
113  }
114 
119  Basis_Derived_HVOL_HEX(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT) : Basis_Derived_HVOL_HEX(polyOrder, polyOrder, polyOrder,pointType) {}
120 
125  virtual
126  const char*
127  getName() const override {
128  return name_.c_str();
129  }
130 
133  virtual bool requireOrientation() const override {
134  return false;
135  }
136 
137  virtual OperatorTensorDecomposition getSimpleOperatorDecomposition(const EOperator operatorType) const override
138  {
139  const EOperator VALUE = Intrepid2::OPERATOR_VALUE;
140 
141  if (operatorType == VALUE)
142  {
143  return OperatorTensorDecomposition(VALUE, VALUE);
144  }
145  else
146  {
147  INTREPID2_TEST_FOR_EXCEPTION(true,std::invalid_argument,"operator not yet supported");
148  }
149  }
150 
152 
160  virtual void getValues(OutputViewType outputValues, const EOperator operatorType,
161  const PointViewType inputPoints1, const PointViewType inputPoints2,
162  bool tensorPoints) const override
163  {
164  Intrepid2::EOperator op1, op2;
165  if (operatorType == Intrepid2::OPERATOR_VALUE)
166  {
167  op1 = Intrepid2::OPERATOR_VALUE;
168  op2 = Intrepid2::OPERATOR_VALUE;
169 
170  this->TensorBasis::getValues(outputValues,
171  inputPoints1, op1,
172  inputPoints2, op2, tensorPoints);
173  }
174  else
175  {
176  INTREPID2_TEST_FOR_EXCEPTION(true,std::invalid_argument,"operator not yet supported");
177  }
178  }
179 
185  getHostBasis() const override {
187  return Teuchos::rcp( new HostBasisType(polyOrder_x_, polyOrder_y_, polyOrder_z_, pointType_) );
188  }
189  };
190 } // end namespace Intrepid2
191 
192 #endif /* Intrepid2_DerivedBasis_HVOL_HEX_h */
Teuchos::RCP< Basis< DeviceType, OutputType, PointType > > BasisPtr
Basis Pointer.
Implementation of bases that are tensor products of two or three component bases. ...
Implementation of H(vol) basis on the quadrilateral that is templated on H(vol) on the line...
virtual void getValues(OutputViewType outputValues, const EOperator operatorType, const PointViewType inputPoints1, const PointViewType inputPoints2, bool tensorPoints) const override
multi-component getValues() method (required/called by TensorBasis)
For a multi-component tensor basis, specifies the operators to be applied to the components to produc...
Basis_Derived_HVOL_HEX(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
virtual bool requireOrientation() const override
True if orientation is required.
virtual OperatorTensorDecomposition getSimpleOperatorDecomposition(const EOperator operatorType) const override
Returns a simple decomposition of the specified operator: what operator(s) should be applied to basis...
virtual const char * getName() const override
Returns basis name.
EOperator
Enumeration of primitive operators available in Intrepid. Primitive operators act on reconstructed fu...
Implementation of H(vol) basis on the quadrilateral that is templated on H(vol) on the line...
virtual BasisPtr< typename Kokkos::HostSpace::device_type, typename BasisBase::OutputValueType, typename BasisBase::PointValueType > getHostBasis() const override
Creates and returns a Basis object whose DeviceType template argument is Kokkos::HostSpace::device_ty...
virtual const char * getName() const override
Returns basis name.
EPointType
Enumeration of types of point distributions in Intrepid.
Implementation of H(vol) basis on the quadrilateral that is templated on H(vol) on the line...
virtual void getValues(BasisValues< OutputValueType, DeviceType > outputValues, const TensorPoints< PointValueType, DeviceType > inputPoints, const EOperator operatorType=OPERATOR_VALUE) const override
Evaluation of a FEM basis on a reference cell, using point and output value containers that allow pre...
Basis_Derived_HVOL_HEX(int polyOrder_x, int polyOrder_y, int polyOrder_z, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
Basis defined as the tensor product of two component bases.