Nektar++
StdTriExp.h
Go to the documentation of this file.
1///////////////////////////////////////////////////////////////////////////////
2//
3// File: StdTriExp.h
4//
5// For more information, please see: http://www.nektar.info
6//
7// The MIT License
8//
9// Copyright (c) 2006 Division of Applied Mathematics, Brown University (USA),
10// Department of Aeronautics, Imperial College London (UK), and Scientific
11// Computing and Imaging Institute, University of Utah (USA).
12//
13// Permission is hereby granted, free of charge, to any person obtaining a
14// copy of this software and associated documentation files (the "Software"),
15// to deal in the Software without restriction, including without limitation
16// the rights to use, copy, modify, merge, publish, distribute, sublicense,
17// and/or sell copies of the Software, and to permit persons to whom the
18// Software is furnished to do so, subject to the following conditions:
19//
20// The above copyright notice and this permission notice shall be included
21// in all copies or substantial portions of the Software.
22//
23// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
24// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
25// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
26// THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
27// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
28// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
29// DEALINGS IN THE SOFTWARE.
30//
31// Description: Header field for triangle routines built upon StdExpansion2D
32//
33///////////////////////////////////////////////////////////////////////////////
34
35#ifndef NEKTAR_LIB_STDREGIONS_STDTRIEXP_H
36#define NEKTAR_LIB_STDREGIONS_STDTRIEXP_H
37
39#include <StdRegions/StdSegExp.h> // for StdSegExp, etc
40
41namespace Nektar::StdRegions
42{
43class StdTriExp : virtual public StdExpansion2D
44{
45public:
47 const LibUtilities::BasisKey &Bb);
50 STD_REGIONS_EXPORT ~StdTriExp() override = default;
51
52protected:
53 //-------------------------------
54 // Integration Methods
55 //-------------------------------
57 v_Integral(const Array<OneD, const NekDouble> &inarray) override;
58
59 //----------------------------
60 // Differentiation Methods
61 //----------------------------
63 const Array<OneD, const NekDouble> &inarray,
67 const int dir, const Array<OneD, const NekDouble> &inarray,
68 Array<OneD, NekDouble> &outarray) override;
70 const Array<OneD, const NekDouble> &inarray,
74 const int dir, const Array<OneD, const NekDouble> &inarray,
75 Array<OneD, NekDouble> &outarray) override;
76
77 //---------------------------------------
78 // Transforms
79 //---------------------------------------
81 const Array<OneD, const NekDouble> &inarray,
82 Array<OneD, NekDouble> &outarray) override;
84 const Array<OneD, const NekDouble> &inarray,
85 Array<OneD, NekDouble> &outarray) override;
89 const Array<OneD, const NekDouble> &inarray,
91 bool doCheckCollDir0, bool doCheckCollDir1) override;
93 const Array<OneD, const NekDouble> &inarray,
94 Array<OneD, NekDouble> &outarray) override;
96 const Array<OneD, const NekDouble> &inarray,
97 Array<OneD, NekDouble> &outarray) override;
98
99 //---------------------------------------
100 // Inner product functions
101 //---------------------------------------
103 const Array<OneD, const NekDouble> &inarray,
104 Array<OneD, NekDouble> &outarray) override;
106 const Array<OneD, const NekDouble> &inarray,
107 Array<OneD, NekDouble> &outarray,
108 bool multiplybyweights = true) override;
110 const Array<OneD, const NekDouble> &base0,
111 const Array<OneD, const NekDouble> &base1,
112 const Array<OneD, const NekDouble> &inarray,
114 bool doCheckCollDir0, bool doCheckCollDir1) override;
116 const int dir, const Array<OneD, const NekDouble> &inarray,
117 Array<OneD, NekDouble> &outarray) override;
119 const int dir, const Array<OneD, const NekDouble> &inarray,
120 Array<OneD, NekDouble> &outarray) override;
121
122 //---------------------------------------
123 // Evaluation functions
124 //---------------------------------------
127 Array<OneD, NekDouble> &eta) override;
130 Array<OneD, NekDouble> &xi) override;
132 const int mode, Array<OneD, NekDouble> &outarray) override;
134 const Array<OneD, const NekDouble> &coords, int mode) final;
137 const Array<OneD, const NekDouble> &inarray,
138 std::array<NekDouble, 3> &firstOrderDerivs) override;
139
140 //---------------------------
141 // Helper functions
142 //---------------------------
143 STD_REGIONS_EXPORT int v_GetNverts() const final;
144 STD_REGIONS_EXPORT int v_GetNtraces() const final;
145 STD_REGIONS_EXPORT LibUtilities::ShapeType v_DetShapeType() const final;
146 STD_REGIONS_EXPORT int v_NumBndryCoeffs() const override;
147 STD_REGIONS_EXPORT int v_NumDGBndryCoeffs() const override;
148 STD_REGIONS_EXPORT int v_GetTraceNcoeffs(const int i) const override;
149 STD_REGIONS_EXPORT int v_GetTraceIntNcoeffs(const int i) const override;
150 STD_REGIONS_EXPORT int v_GetTraceNumPoints(const int i) const override;
152 const std::vector<unsigned int> &nummodes, int &modes_offset) override;
154 Array<OneD, NekDouble> &coords_x, Array<OneD, NekDouble> &coords_y,
155 Array<OneD, NekDouble> &coords_z) override;
158 const int i, const int j, bool UseGLL = false) const override;
159
160 //--------------------------
161 // Mappings
162 //--------------------------
164 int localVertexId, bool useCoeffPacking = false) override;
166 Array<OneD, unsigned int> &outarray) override;
168 Array<OneD, unsigned int> &outarray) override;
170 const unsigned int traceid,
171 Array<OneD, unsigned int> &maparray) override;
173 const int eid, Array<OneD, unsigned int> &maparray,
174 Array<OneD, int> &signarray,
175 const Orientation edgeOrient = eForwards) override;
176
177 //---------------------------------------
178 // Wrapper functions
179 //---------------------------------------
181 v_GenMatrix(const StdMatrixKey &mkey) override;
183 v_CreateStdMatrix(const StdMatrixKey &mkey) override;
184
185 //---------------------------------------
186 // Operator evaluation functions
187 //---------------------------------------
189 const Array<OneD, const NekDouble> &inarray,
190 Array<OneD, NekDouble> &outarray, const StdMatrixKey &mkey) override;
192 const Array<OneD, const NekDouble> &inarray,
193 Array<OneD, NekDouble> &outarray, const StdMatrixKey &mkey) override;
195 Array<OneD, NekDouble> &array, const StdMatrixKey &mkey) override;
197 int numMin, const Array<OneD, const NekDouble> &inarray,
198 Array<OneD, NekDouble> &outarray) override;
200 const int k1, const int k2, const Array<OneD, const NekDouble> &inarray,
201 Array<OneD, NekDouble> &outarray, const StdMatrixKey &mkey) override;
203 const int i, const Array<OneD, const NekDouble> &inarray,
204 Array<OneD, NekDouble> &outarray, const StdMatrixKey &mkey) override;
206 const Array<OneD, const NekDouble> &inarray,
207 Array<OneD, NekDouble> &outarray, const StdMatrixKey &mkey) override;
209 const Array<OneD, const NekDouble> &inarray,
210 Array<OneD, NekDouble> &outarray) override;
211
212 //---------------------------------------
213 // Output interpolation functions
214 //---------------------------------------
216 Array<OneD, int> &conn, bool standard = true) override;
217};
218typedef std::shared_ptr<StdTriExp> StdTriExpSharedPtr;
219} // namespace Nektar::StdRegions
220#endif // STDTRIEXP_H
#define STD_REGIONS_EXPORT
Describes the specification for a Basis.
Definition: Basis.h:45
void v_WeakDerivMatrixOp(const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
Definition: StdTriExp.cpp:1328
LibUtilities::ShapeType v_DetShapeType() const final
Definition: StdTriExp.cpp:728
int v_GetTraceNumPoints(const int i) const override
Definition: StdTriExp.cpp:781
void v_BwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Backward tranform for triangular elements.
Definition: StdTriExp.cpp:217
void v_IProductWRTDerivBase(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Definition: StdTriExp.cpp:491
void v_LocCollapsedToLocCoord(const Array< OneD, const NekDouble > &eta, Array< OneD, NekDouble > &xi) override
Definition: StdTriExp.cpp:596
void v_GetCoords(Array< OneD, NekDouble > &coords_x, Array< OneD, NekDouble > &coords_y, Array< OneD, NekDouble > &coords_z) override
Definition: StdTriExp.cpp:805
int v_CalcNumberOfCoefficients(const std::vector< unsigned int > &nummodes, int &modes_offset) override
Definition: StdTriExp.cpp:795
const LibUtilities::BasisKey v_GetTraceBasisKey(const int i, const int j, bool UseGLL=false) const override
Definition: StdTriExp.cpp:831
void v_GetSimplexEquiSpacedConnectivity(Array< OneD, int > &conn, bool standard=true) override
Definition: StdTriExp.cpp:1543
void v_MassMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
Definition: StdTriExp.cpp:1306
void v_GetBoundaryMap(Array< OneD, unsigned int > &outarray) override
Definition: StdTriExp.cpp:1064
int v_GetNtraces() const final
Definition: StdTriExp.cpp:723
void v_StdPhysDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray) override
Definition: StdTriExp.cpp:193
void v_GetTraceInteriorToElementMap(const int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, const Orientation edgeOrient=eForwards) override
Definition: StdTriExp.cpp:1151
int v_GetVertexMap(int localVertexId, bool useCoeffPacking=false) override
Definition: StdTriExp.cpp:971
DNekMatSharedPtr v_CreateStdMatrix(const StdMatrixKey &mkey) override
Definition: StdTriExp.cpp:1297
void v_BwdTrans_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Definition: StdTriExp.cpp:223
DNekMatSharedPtr v_GenMatrix(const StdMatrixKey &mkey) override
Definition: StdTriExp.cpp:1223
~StdTriExp() override=default
void v_IProductWRTBase(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Calculate the inner product of inarray with respect to the basis B=base0[p]*base1[pq] and put into ou...
Definition: StdTriExp.cpp:420
int v_GetTraceNcoeffs(const int i) const override
Definition: StdTriExp.cpp:753
void v_PhysDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray) override
Calculate the derivative of the physical points.
Definition: StdTriExp.cpp:106
void v_GetTraceCoeffMap(const unsigned int traceid, Array< OneD, unsigned int > &maparray) override
Definition: StdTriExp.cpp:1096
void v_BwdTrans_SumFacKernel(const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1) override
Definition: StdTriExp.cpp:233
bool v_IsBoundaryInteriorExpansion() const override
Definition: StdTriExp.cpp:825
void v_FillMode(const int mode, Array< OneD, NekDouble > &outarray) override
Definition: StdTriExp.cpp:603
void v_ReduceOrderCoeffs(int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Definition: StdTriExp.cpp:1445
void v_FwdTransBndConstrained(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Definition: StdTriExp.cpp:289
void v_IProductWRTDerivBase_SumFac(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Definition: StdTriExp.cpp:498
StdTriExp(const StdTriExp &T)=default
int v_GetTraceIntNcoeffs(const int i) const override
Definition: StdTriExp.cpp:767
int v_NumDGBndryCoeffs() const override
Definition: StdTriExp.cpp:743
void v_LocCoordToLocCollapsed(const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta) override
Definition: StdTriExp.cpp:577
NekDouble v_PhysEvaluateBasis(const Array< OneD, const NekDouble > &coords, int mode) final
Definition: StdTriExp.cpp:648
void v_HelmholtzMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
Definition: StdTriExp.cpp:1336
NekDouble v_Integral(const Array< OneD, const NekDouble > &inarray) override
Integrates the specified function over the domain.
Definition: StdTriExp.cpp:58
NekDouble v_PhysEvalFirstDeriv(const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) override
Definition: StdTriExp.cpp:672
void v_FwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Transform a given function from physical quadrature space to coefficient space.
Definition: StdTriExp.cpp:273
void v_GetInteriorMap(Array< OneD, unsigned int > &outarray) override
Definition: StdTriExp.cpp:1034
int v_GetNverts() const final
Definition: StdTriExp.cpp:718
void v_LaplacianMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
Definition: StdTriExp.cpp:1313
void v_SVVLaplacianFilter(Array< OneD, NekDouble > &array, const StdMatrixKey &mkey) override
Definition: StdTriExp.cpp:1343
void v_MultiplyByStdQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Definition: StdTriExp.cpp:1503
int v_NumBndryCoeffs() const override
Definition: StdTriExp.cpp:733
void v_IProductWRTBase_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true) override
Definition: StdTriExp.cpp:426
void v_IProductWRTBase_SumFacKernel(const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1) override
Definition: StdTriExp.cpp:453
std::shared_ptr< StdTriExp > StdTriExpSharedPtr
Definition: StdTriExp.h:218
static Array< OneD, NekDouble > NullNekDouble1DArray
std::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:75
double NekDouble
STL namespace.