Intrepid2
Intrepid2_DerivedBasis_HGRAD_WEDGE.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
51#ifndef Intrepid2_DerivedBasis_HGRAD_WEDGE_h
52#define Intrepid2_DerivedBasis_HGRAD_WEDGE_h
53
55
56namespace Intrepid2
57{
58 template<class HGRAD_TRI, class HGRAD_LINE>
60 : public Basis_TensorBasis<typename HGRAD_LINE::BasisBase>
61 {
62 protected:
63 std::string name_;
64 ordinal_type order_xy_, order_z_;
65 EPointType pointType_;
66 public:
67 using ExecutionSpace = typename HGRAD_LINE::ExecutionSpace;
68 using OutputValueType = typename HGRAD_LINE::OutputValueType;
69 using PointValueType = typename HGRAD_LINE::PointValueType;
70
71 using OutputViewType = typename HGRAD_LINE::OutputViewType;
72 using PointViewType = typename HGRAD_LINE::PointViewType ;
73 using ScalarViewType = typename HGRAD_LINE::ScalarViewType;
74
75 using BasisBase = typename HGRAD_LINE::BasisBase;
76 using TriBasis = HGRAD_TRI;
77 using LineBasis = HGRAD_LINE;
79
85 Basis_Derived_HGRAD_WEDGE(int polyOrder_xy, int polyOrder_z, const EPointType pointType=POINTTYPE_DEFAULT)
86 :
87 TensorBasis(Teuchos::rcp( new TriBasis(polyOrder_xy, pointType)),
88 Teuchos::rcp( new LineBasis(polyOrder_z, pointType)))
89 {
90 this->functionSpace_ = FUNCTION_SPACE_HGRAD;
91
92 std::ostringstream basisName;
93 basisName << "HGRAD_WEDGE (" << this->TensorBasis::getName() << ")";
94 name_ = basisName.str();
95
96 order_xy_= polyOrder_xy;
97 order_z_ = polyOrder_z;
98 pointType_ = pointType;
99
100 this->setShardsTopologyAndTags();
101 }
102
107 Basis_Derived_HGRAD_WEDGE(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT) : Basis_Derived_HGRAD_WEDGE(polyOrder,polyOrder,pointType) {}
108
109
112 virtual bool requireOrientation() const override {
113 return (this->getDofCount(1,0) > 1); //if it has more than 1 DOF per edge, than it needs orientations
114 }
115
116 using BasisBase::getValues;
117
122 virtual
123 const char*
124 getName() const override {
125 return name_.c_str();
126 }
127
137 Teuchos::RCP<BasisBase>
138 getSubCellRefBasis(const ordinal_type subCellDim, const ordinal_type subCellOrd) const override
139 {
140 INTREPID2_TEST_FOR_EXCEPTION(true,std::invalid_argument,"Method not yet implemented");
141 if(subCellDim == 2) {
142 switch(subCellOrd) {
143 // TODO: check this subcell numbering -- which are the tris, and which are the quads?
144 // (one way: use the shards CellTopology to look this up…)
145 case 0:
146 case 1:
147 return Teuchos::rcp( new TriBasis(order_xy_, pointType_) );
148 case 2:
149 case 3:
150 case 4:
151 // TODO: check what order the poly order arguments should go in...
152 // (one way: use the shards CellTopology to look this up…)
153 // TODO: define QuadBasis somehow, somewhere. (Probably use DerivedBasis_HGRAD_QUAD.)
154// return Teuchos::rcp( new QuadBasis(order_xy_, order_z_, pointType_) );
155 return Teuchos::null;
156 }
157 return Teuchos::null;
158 }
159
160 INTREPID2_TEST_FOR_EXCEPTION(true,std::invalid_argument,"Input parameters out of bounds");
161 }
162
168 getHostBasis() const override {
170
171 auto hostBasis = Teuchos::rcp(new HostBasis(order_xy_, order_z_, pointType_));
172
173 return hostBasis;
174 }
175 };
176} // end namespace Intrepid2
177
178#endif /* Intrepid2_DerivedBasis_HGRAD_WEDGE_h */
BasisPtr< typename Kokkos::HostSpace::device_type, OutputType, PointType > HostBasisPtr
Pointer to a Basis whose device type is on the host (Kokkos::HostSpace::device_type),...
Implementation of bases that are tensor products of two or three component bases.
Teuchos::RCP< BasisBase > getSubCellRefBasis(const ordinal_type subCellDim, const ordinal_type subCellOrd) const override
returns the basis associated to a subCell.
Basis_Derived_HGRAD_WEDGE(int polyOrder, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
virtual HostBasisPtr< OutputValueType, PointValueType > getHostBasis() const override
Creates and returns a Basis object whose DeviceType template argument is Kokkos::HostSpace::device_ty...
Basis_Derived_HGRAD_WEDGE(int polyOrder_xy, int polyOrder_z, const EPointType pointType=POINTTYPE_DEFAULT)
Constructor.
virtual const char * getName() const override
Returns basis name.
virtual bool requireOrientation() const override
True if orientation is required.
Basis defined as the tensor product of two component bases.
virtual const char * getName() const override
Returns basis name.