Nektar++
StdSegExp.h
Go to the documentation of this file.
1 ///////////////////////////////////////////////////////////////////////////////
2 //
3 // File StdSegExp.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 file for Standard Segment Expansions
32 //
33 ///////////////////////////////////////////////////////////////////////////////
34 
35 #ifndef NEKTAR_LIBS_STDREGIONS_STDSEGEXP_H
36 #define NEKTAR_LIBS_STDREGIONS_STDSEGEXP_H
37 
40 
41 namespace Nektar
42 {
43 namespace StdRegions
44 {
45 
46 class StdSegExp;
47 typedef std::shared_ptr<StdSegExp> StdSegExpSharedPtr;
48 
49 /// Class representing a segment element in reference space
50 
51 /// All interface of this class sits in StdExpansion class
52 
53 class StdSegExp : virtual public StdExpansion1D
54 {
55 public:
57 
59 
61 
63 
64 protected:
65  //----------------------------
66  // Integration Methods
67  //----------------------------
68 
70  const Array<OneD, const NekDouble> &inarray);
71 
72  //-----------------------------
73  // Differentiation Methods
74  //-----------------------------
75 
76  STD_REGIONS_EXPORT virtual void v_PhysDeriv(
77  const Array<OneD, const NekDouble> &inarray,
78  Array<OneD, NekDouble> &out_d0,
81 
82  STD_REGIONS_EXPORT virtual void v_PhysDeriv(
83  const int dir, const Array<OneD, const NekDouble> &inarray,
84  Array<OneD, NekDouble> &outarray);
85 
87  const Array<OneD, const NekDouble> &inarray,
88  Array<OneD, NekDouble> &out_d0,
91 
93  const int dir, const Array<OneD, const NekDouble> &inarray,
94  Array<OneD, NekDouble> &outarray);
95 
96  //-----------------------------
97  // Transforms
98  //-----------------------------
99 
100  STD_REGIONS_EXPORT virtual void v_BwdTrans(
101  const Array<OneD, const NekDouble> &inarray,
102  Array<OneD, NekDouble> &outarray);
103 
104  STD_REGIONS_EXPORT virtual void v_FwdTrans(
105  const Array<OneD, const NekDouble> &inarray,
106  Array<OneD, NekDouble> &outarray);
107 
109  const Array<OneD, const NekDouble> &inarray,
110  Array<OneD, NekDouble> &outarray);
111 
113  const Array<OneD, const NekDouble> &inarray,
114  Array<OneD, NekDouble> &outarray);
115 
116  //----------------------------
117  // Inner product functions
118  //----------------------------
119 
121  const Array<OneD, const NekDouble> &inarray,
122  Array<OneD, NekDouble> &outarray);
123 
125  const Array<OneD, const NekDouble> &base,
126  const Array<OneD, const NekDouble> &inarray,
127  Array<OneD, NekDouble> &outarray, int coll_check);
128 
130  const Array<OneD, const NekDouble> &inarray,
131  Array<OneD, NekDouble> &outarray, bool multiplybyweights = true);
132 
134  const int dir, const Array<OneD, const NekDouble> &inarray,
135  Array<OneD, NekDouble> &outarray);
136 
138  const int dir, const Array<OneD, const NekDouble> &inarray,
139  Array<OneD, NekDouble> &outarray);
140 
141  //----------------------------
142  // Evaluations Methods
143  //---------------------------
146 
149 
151  const Array<OneD, const NekDouble> &coords, int mode) final;
152 
154  const Array<OneD, const NekDouble> &inarray,
155  Array<OneD, NekDouble> &outarray, const StdMatrixKey &mkey);
156 
158  const Array<OneD, const NekDouble> &inarray,
159  Array<OneD, NekDouble> &outarray, const StdMatrixKey &mkey);
160 
162  Array<OneD, NekDouble> &array, const StdMatrixKey &mkey);
163 
165  Array<OneD, NekDouble> &array, const NekDouble alpha,
166  const NekDouble exponent, const NekDouble cutoff);
167 
169  const Array<OneD, const NekDouble> &inarray,
170  Array<OneD, NekDouble> &outarray);
171 
172  STD_REGIONS_EXPORT virtual void v_FillMode(
173  const int mode, Array<OneD, NekDouble> &outarray);
174 
175  STD_REGIONS_EXPORT virtual void v_GetCoords(
176  Array<OneD, NekDouble> &coords_0, Array<OneD, NekDouble> &coords_1,
177  Array<OneD, NekDouble> &coords_2);
178 
179  //----------------------------
180  // Public Mappings
181  //---------------------------
182 
184  Array<OneD, unsigned int> &outarray);
185 
187  Array<OneD, unsigned int> &outarray);
188 
189  STD_REGIONS_EXPORT virtual int v_GetVertexMap(int localVertexId,
190  bool useCoeffPacking = false);
191 
192  //----------------------------
193  // Helper functions
194  //---------------------------
195  STD_REGIONS_EXPORT virtual int v_GetNverts() const;
196  STD_REGIONS_EXPORT virtual int v_GetNtraces() const;
197  STD_REGIONS_EXPORT virtual int v_GetTraceNcoeffs(const int i) const;
198  STD_REGIONS_EXPORT virtual int v_GetTraceNumPoints(const int i) const;
199  STD_REGIONS_EXPORT virtual int v_NumBndryCoeffs() const;
200  STD_REGIONS_EXPORT virtual int v_NumDGBndryCoeffs() const;
203  const std::vector<unsigned int> &nummodes, int &modes_offset);
205 
206  //----------------------------
207  // Wrapper functions
208  //---------------------------
209 
211  const StdMatrixKey &mkey);
212 
214  const StdMatrixKey &mkey);
215 
216  //---------------------------------------
217  // Output interpolation functions
218  //---------------------------------------
220  Array<OneD, int> &conn, bool standard = true);
221 
222  // Operator evaluation functions
223  //---------------------------------------
224 
226  int numMin, const Array<OneD, const NekDouble> &inarray,
227  Array<OneD, NekDouble> &outarray);
228 
230  const unsigned int traceid, Array<OneD, unsigned int> &maparray);
231 
233  const unsigned int eid, Array<OneD, unsigned int> &maparray,
234  Array<OneD, int> &signarray, Orientation edgeOrient, int P, int Q);
235 
237  const int tid, Array<OneD, unsigned int> &maparray,
238  Array<OneD, int> &signarray, Orientation edgeOrient, int P, int Q);
239 
240 private:
241 };
242 
243 } // namespace StdRegions
244 } // namespace Nektar
245 
246 #endif // STDSEGEXP_H
#define STD_REGIONS_EXPORT
Describes the specification for a Basis.
Definition: Basis.h:50
Class representing a segment element in reference space.
Definition: StdSegExp.h:54
virtual int v_GetTraceNumPoints(const int i) const
Definition: StdSegExp.cpp:675
virtual void v_IProductWRTDerivBase_SumFac(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdSegExp.cpp:451
virtual void v_GetCoords(Array< OneD, NekDouble > &coords_0, Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2)
Definition: StdSegExp.cpp:647
virtual void v_ExponentialFilter(Array< OneD, NekDouble > &array, const NekDouble alpha, const NekDouble exponent, const NekDouble cutoff)
Definition: StdSegExp.cpp:596
virtual void v_StdPhysDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1=NullNekDouble1DArray, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
Definition: StdSegExp.cpp:162
virtual void v_FillMode(const int mode, Array< OneD, NekDouble > &outarray)
Definition: StdSegExp.cpp:498
virtual void v_GetElmtTraceToTraceMap(const unsigned int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation edgeOrient, int P, int Q)
Definition: StdSegExp.cpp:922
virtual int v_GetNverts() const
Definition: StdSegExp.cpp:659
virtual void v_GetBoundaryMap(Array< OneD, unsigned int > &outarray)
Definition: StdSegExp.cpp:776
virtual void v_LocCoordToLocCollapsed(const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta)
Definition: StdSegExp.cpp:486
virtual void v_HelmholtzMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdSegExp.cpp:533
NekDouble v_PhysEvaluateBasis(const Array< OneD, const NekDouble > &coords, int mode) final
Definition: StdSegExp.cpp:509
virtual DNekMatSharedPtr v_CreateStdMatrix(const StdMatrixKey &mkey)
Definition: StdSegExp.cpp:767
virtual void v_GetTraceToElementMap(const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation edgeOrient, int P, int Q)
Definition: StdSegExp.cpp:903
virtual void v_IProductWRTBase_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)
Definition: StdSegExp.cpp:460
virtual void v_PhysDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1=NullNekDouble1DArray, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
Evaluate the derivative at the physical quadrature points given by inarray and return in outarray.
Definition: StdSegExp.cpp:143
virtual void v_LaplacianMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdSegExp.cpp:515
virtual void v_MultiplyByStdQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdSegExp.cpp:636
virtual int v_NumDGBndryCoeffs() const
Definition: StdSegExp.cpp:686
virtual int v_GetVertexMap(int localVertexId, bool useCoeffPacking=false)
Definition: StdSegExp.cpp:838
virtual LibUtilities::ShapeType v_DetShapeType() const
Return Shape of region, using ShapeType enum list. i.e. Segment.
Definition: StdSegExp.cpp:79
virtual void v_ReduceOrderCoeffs(int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdSegExp.cpp:223
StdSegExp()
Default constructor.
Definition: StdSegExp.cpp:49
virtual void v_GetTraceCoeffMap(const unsigned int traceid, Array< OneD, unsigned int > &maparray)
Get the map of the coefficient location to teh local trace coefficients.
Definition: StdSegExp.cpp:874
virtual void v_LocCollapsedToLocCoord(const Array< OneD, const NekDouble > &eta, Array< OneD, NekDouble > &xi)
Definition: StdSegExp.cpp:492
virtual int v_GetNtraces() const
Definition: StdSegExp.cpp:664
virtual int v_NumBndryCoeffs() const
Definition: StdSegExp.cpp:681
virtual void v_BwdTrans_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdSegExp.cpp:368
virtual NekDouble v_Integral(const Array< OneD, const NekDouble > &inarray)
Integrate the physical point list inarray over region and return the value.
Definition: StdSegExp.cpp:114
virtual void v_BwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Backward transform from coefficient space given in inarray and evaluate at the physical quadrature po...
Definition: StdSegExp.cpp:205
virtual void v_FwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Forward transform from physical quadrature space stored in inarray and evaluate the expansion coeffic...
Definition: StdSegExp.cpp:272
virtual void v_SVVLaplacianFilter(Array< OneD, NekDouble > &array, const StdMatrixKey &mkey)
Definition: StdSegExp.cpp:555
virtual DNekMatSharedPtr v_GenMatrix(const StdMatrixKey &mkey)
Definition: StdSegExp.cpp:704
virtual void v_FwdTransBndConstrained(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdSegExp.cpp:294
virtual void v_IProductWRTDerivBase(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdSegExp.cpp:444
virtual bool v_IsBoundaryInteriorExpansion()
Definition: StdSegExp.cpp:84
virtual void v_GetInteriorMap(Array< OneD, unsigned int > &outarray)
Definition: StdSegExp.cpp:805
virtual void v_IProductWRTBase(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Inner product of inarray over region with respect to the expansion basis (this)->m_base[0] and return...
Definition: StdSegExp.cpp:438
virtual int v_GetTraceNcoeffs(const int i) const
Definition: StdSegExp.cpp:669
virtual void v_GetSimplexEquiSpacedConnectivity(Array< OneD, int > &conn, bool standard=true)
Definition: StdSegExp.cpp:855
virtual int v_CalcNumberOfCoefficients(const std::vector< unsigned int > &nummodes, int &modes_offset)
Definition: StdSegExp.cpp:691
std::shared_ptr< StdSegExp > StdSegExpSharedPtr
Definition: StdSegExp.h:46
The above copyright notice and this permission notice shall be included.
Definition: CoupledSolver.h:1
static Array< OneD, NekDouble > NullNekDouble1DArray
std::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:75
double NekDouble