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
63namespace 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 this->setShardsTopologyAndTags();
115 }
116
121 Basis_Derived_HVOL_HEX(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT) : Basis_Derived_HVOL_HEX(polyOrder, polyOrder, polyOrder,pointType) {}
122
127 virtual
128 const char*
129 getName() const override {
130 return name_.c_str();
131 }
132
135 virtual bool requireOrientation() const override {
136 return false;
137 }
138
139 virtual OperatorTensorDecomposition getSimpleOperatorDecomposition(const EOperator &operatorType) const override
140 {
141 const EOperator VALUE = Intrepid2::OPERATOR_VALUE;
142
143 if (operatorType == VALUE)
144 {
145 return OperatorTensorDecomposition(VALUE, VALUE);
146 }
147 else
148 {
149 INTREPID2_TEST_FOR_EXCEPTION(true,std::invalid_argument,"operator not yet supported");
150 }
151 }
152
154
160 getHostBasis() const override {
162 return Teuchos::rcp( new HostBasisType(polyOrder_x_, polyOrder_y_, polyOrder_z_, pointType_) );
163 }
164 };
165} // end namespace Intrepid2
166
167#endif /* Intrepid2_DerivedBasis_HVOL_HEX_h */
Teuchos::RCP< Basis< DeviceType, OutputType, PointType > > BasisPtr
Basis Pointer.
Implementation of H(vol) basis on the quadrilateral that is templated on H(vol) on the line.
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_HEX(int polyOrder_x, int polyOrder_y, int polyOrder_z, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
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.
Basis_Derived_HVOL_HEX(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
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 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...