Intrepid2
Intrepid2_DerivedBasis_HVOL_QUAD.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 
52 #ifndef Intrepid2_DerivedBasis_HVOL_QUAD_h
53 #define Intrepid2_DerivedBasis_HVOL_QUAD_h
54 
57 
58 namespace Intrepid2
59 {
64  template<class HVOL_LINE>
66  :
67  public Basis_TensorBasis<typename HVOL_LINE::BasisBase>
68  {
69  protected:
70  std::string name_;
71  using LineBasis = HVOL_LINE;
73 
74  ordinal_type polyOrder_x_, polyOrder_y_;
75  EPointType pointType_;
76  public:
77  using ExecutionSpace = typename HVOL_LINE::ExecutionSpace;
78  using OutputValueType = typename HVOL_LINE::OutputValueType;
79  using PointValueType = typename HVOL_LINE::PointValueType;
80 
81  using OutputViewType = typename HVOL_LINE::OutputViewType;
82  using PointViewType = typename HVOL_LINE::PointViewType ;
83  using ScalarViewType = typename HVOL_LINE::ScalarViewType;
84 
85  using BasisBase = typename HVOL_LINE::BasisBase;
86 
92  Basis_Derived_HVOL_QUAD(int polyOrder_x, int polyOrder_y, const EPointType pointType=POINTTYPE_DEFAULT)
93  :
94  TensorBasis(Teuchos::rcp( new LineBasis(polyOrder_x, pointType)),
95  Teuchos::rcp( new LineBasis(polyOrder_y, pointType))),
96  polyOrder_x_(polyOrder_x),
97  polyOrder_y_(polyOrder_y),
98  pointType_(pointType)
99  {
100  this->functionSpace_ = FUNCTION_SPACE_HVOL;
101 
102  std::ostringstream basisName;
103  basisName << "HVOL_QUAD (" << this->TensorBasis::getName() << ")";
104  name_ = basisName.str();
105  }
106 
111  Basis_Derived_HVOL_QUAD(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT) : Basis_Derived_HVOL_QUAD(polyOrder,polyOrder,pointType) {}
112 
117  virtual
118  const char*
119  getName() const override {
120  return name_.c_str();
121  }
122 
125  virtual bool requireOrientation() const override {
126  return false;
127  }
128 
129  virtual OperatorTensorDecomposition getSimpleOperatorDecomposition(const EOperator operatorType) const override
130  {
131  const EOperator VALUE = Intrepid2::OPERATOR_VALUE;
132 
133  if (operatorType == VALUE)
134  {
135  return OperatorTensorDecomposition(VALUE, VALUE);
136  }
137  else
138  {
139  INTREPID2_TEST_FOR_EXCEPTION(true,std::invalid_argument,"operator not yet supported");
140  }
141  }
142 
144 
152  virtual void getValues(OutputViewType outputValues, const EOperator operatorType,
153  const PointViewType inputPoints1, const PointViewType inputPoints2,
154  bool tensorPoints) const override
155  {
156  Intrepid2::EOperator op1, op2;
157  if (operatorType == Intrepid2::OPERATOR_VALUE)
158  {
159  op1 = Intrepid2::OPERATOR_VALUE;
160  op2 = Intrepid2::OPERATOR_VALUE;
161 
162  this->TensorBasis::getValues(outputValues,
163  inputPoints1, op1,
164  inputPoints2, op2, tensorPoints);
165  }
166  else
167  {
168  INTREPID2_TEST_FOR_EXCEPTION(true,std::invalid_argument,"operator not yet supported");
169  }
170  }
171 
177  getHostBasis() const override {
179  return Teuchos::rcp( new HostBasisType(polyOrder_x_, polyOrder_y_, pointType_) );
180  }
181  };
182 } // end namespace Intrepid2
183 
184 #endif /* Intrepid2_DerivedBasis_HVOL_QUAD_h */
Teuchos::RCP< Basis< DeviceType, OutputType, PointType > > BasisPtr
Basis Pointer.
Basis_Derived_HVOL_QUAD(int polyOrder_x, int polyOrder_y, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
Implementation of bases that are tensor products of two or three component bases. ...
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...
Implementation of H(vol) basis on the quadrilateral that is templated on H(vol) on the line...
Basis_Derived_HVOL_QUAD(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
For a multi-component tensor basis, specifies the operators to be applied to the components to produc...
virtual const char * getName() const override
Returns basis name.
virtual OperatorTensorDecomposition getSimpleOperatorDecomposition(const EOperator operatorType) const override
Returns a simple decomposition of the specified operator: what operator(s) should be applied to basis...
EOperator
Enumeration of primitive operators available in Intrepid. Primitive operators act on reconstructed fu...
H(vol) basis on the line based on Legendre polynomials.
virtual const char * getName() const override
Returns basis name.
EPointType
Enumeration of types of point distributions in Intrepid.
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)
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 defined as the tensor product of two component bases.
virtual bool requireOrientation() const override
True if orientation is required.