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
58namespace 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 this->setShardsTopologyAndTags();
107 }
108
113 Basis_Derived_HVOL_QUAD(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT) : Basis_Derived_HVOL_QUAD(polyOrder,polyOrder,pointType) {}
114
119 virtual
120 const char*
121 getName() const override {
122 return name_.c_str();
123 }
124
127 virtual bool requireOrientation() const override {
128 return false;
129 }
130
131 virtual OperatorTensorDecomposition getSimpleOperatorDecomposition(const EOperator &operatorType) const override
132 {
133 const EOperator VALUE = Intrepid2::OPERATOR_VALUE;
134
135 if (operatorType == VALUE)
136 {
137 return OperatorTensorDecomposition(VALUE, VALUE);
138 }
139 else
140 {
141 INTREPID2_TEST_FOR_EXCEPTION(true,std::invalid_argument,"operator not yet supported");
142 }
143 }
144
146
152 getHostBasis() const override {
154 return Teuchos::rcp( new HostBasisType(polyOrder_x_, polyOrder_y_, pointType_) );
155 }
156 };
157} // end namespace Intrepid2
158
159#endif /* Intrepid2_DerivedBasis_HVOL_QUAD_h */
Teuchos::RCP< Basis< DeviceType, OutputType, PointType > > BasisPtr
Basis Pointer.
H(vol) basis on the line based on Legendre polynomials.
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 bool requireOrientation() const override
True if orientation is required.
Basis_Derived_HVOL_QUAD(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
virtual const char * getName() const override
Returns basis name.
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 OperatorTensorDecomposition getSimpleOperatorDecomposition(const EOperator &operatorType) const override
Returns a simple decomposition of the specified operator: what operator(s) should be applied to basis...
Basis_Derived_HVOL_QUAD(int polyOrder_x, int polyOrder_y, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
Basis defined as the tensor product of two component bases.
virtual const char * getName() const override
Returns basis name.
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...
For a multi-component tensor basis, specifies the operators to be applied to the components to produc...