Nektar++
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Expansion3D.h
Go to the documentation of this file.
1 ///////////////////////////////////////////////////////////////////////////////
2 //
3 // File: Expansion3D.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 Expansion3D routines
32 //
33 ///////////////////////////////////////////////////////////////////////////////
34 
35 #ifndef EXPANSION3D_H
36 #define EXPANSION3D_H
37 
38 #include <LocalRegions/Expansion.h>
42 
43 namespace Nektar
44 {
45  namespace LocalRegions
46  {
47 
48  class Expansion2D;
49  typedef std::shared_ptr<Expansion2D> Expansion2DSharedPtr;
50  typedef std::weak_ptr<Expansion2D> Expansion2DWeakPtr;
51 
52  class Expansion3D;
53  typedef std::shared_ptr<Expansion3D> Expansion3DSharedPtr;
54  typedef std::weak_ptr<Expansion3D> Expansion3DWeakPtr;
55  typedef std::vector< Expansion3DSharedPtr > Expansion3DVector;
56 
57  class Expansion3D: virtual public Expansion,
58  virtual public StdRegions::StdExpansion3D
59  {
60  public:
63 
65  (Array<OneD, NekDouble> &inout);
67  (const int face, Array<OneD, NekDouble> &inout);
68  inline void AddHDGHelmholtzFaceTerms(
69  const NekDouble tau,
70  const int edge,
71  Array<OneD, NekDouble> &facePhys,
72  const StdRegions::VarCoeffMap &dirForcing,
73  Array<OneD, NekDouble> &outarray);
74 
75  inline void AddNormTraceInt(
76  const int dir,
78  Array<OneD, Array<OneD, NekDouble> > &faceCoeffs,
79  Array<OneD,NekDouble> &outarray);
80 
81  inline void AddNormTraceInt(
82  const int dir,
85  Array<OneD,NekDouble> &outarray,
86  const StdRegions::VarCoeffMap &varcoeffs);
87 
88  inline void AddFaceBoundaryInt(
89  const int face,
90  ExpansionSharedPtr &FaceExp,
91  Array<OneD, NekDouble> &facePhys,
92  Array<OneD, NekDouble> &outarray,
93  const StdRegions::VarCoeffMap &varcoeffs
95 
97 
99  const StdRegions::Orientation orient,
100  Array<OneD, int> &idmap,
101  const int nq0,
102  const int nq1);
103 
105  const Array<OneD, const Array<OneD, NekDouble> > &Fvec,
106  Array<OneD, NekDouble> &outarray);
107 
109  GetEdgeInverseBoundaryMap(int eid);
110 
114  faceOrient = StdRegions::eNoOrientation,
115  int P1=-1, int P2=-1);
116 
121 
122  protected:
123  std::map<int, NormalVector> m_faceNormals;
124 
125  virtual void v_DGDeriv(
126  const int dir,
127  const Array<OneD, const NekDouble> &incoeffs,
129  Array<OneD, Array<OneD, NekDouble> > &faceCoeffs,
130  Array<OneD, NekDouble> &out_d);
132  const StdRegions::StdMatrixKey &mkey);
133  virtual void v_AddFaceNormBoundaryInt(
134  const int face,
135  const ExpansionSharedPtr &FaceExp,
137  Array<OneD, NekDouble> &outarray);
138 
139  virtual void v_AddRobinMassMatrix(
140  const int face,
141  const Array<OneD, const NekDouble> &primCoeffs,
142  DNekMatSharedPtr &inoutmat);
143 
144  virtual StdRegions::Orientation v_GetTraceOrient(int face);
145 
146  virtual void v_GetTracePhysVals(
147  const int face,
148  const StdRegions::StdExpansionSharedPtr &FaceExp,
149  const Array<OneD, const NekDouble> &inarray,
150  Array<OneD, NekDouble> &outarray,
151  StdRegions::Orientation orient);
152 
153 
155  const int face,
156  ExpansionSharedPtr &FaceExp,
157  const Array<OneD, const NekDouble> &varcoeff,
158  Array<OneD,NekDouble> &outarray);
159 
161  const int dir,
162  const int face,
163  ExpansionSharedPtr &FaceExp_f,
164  const Array<OneD, const Array<OneD, NekDouble> > &normals,
165  const StdRegions::VarCoeffMap &varcoeffs);
166 
167  //-----------------------------
168  // Low Energy Basis functions
169  //-----------------------------
170 
172  const DNekScalMatSharedPtr &r_bnd,
173  const StdRegions::MatrixType matrixType);
174 
176  const DNekScalMatSharedPtr & transformationmatrix);
177 
179  const DNekScalMatSharedPtr &r_bnd);
180 
182  (const int face) const;
183 
184  private:
185  // Do not add members here since it may lead to conflicts.
186  // Only use this class for member functions
187 
188  std::vector<bool> m_requireNeg;
189  };
190 
192  {
193  return std::dynamic_pointer_cast<SpatialDomains::Geometry3D>(m_geom);
194  }
195  } //end of namespace
196 } //end of namespace
197 
198 #endif
#define LOCAL_REGIONS_EXPORT
Expansion3D(SpatialDomains::Geometry3DSharedPtr pGeom)
Definition: Expansion3D.h:61
void SetTraceToGeomOrientation(Array< OneD, NekDouble > &inout)
Align trace orientation with the geometry orientation.
void AddNormTraceInt(const int dir, Array< OneD, ExpansionSharedPtr > &FaceExp, Array< OneD, Array< OneD, NekDouble > > &faceCoeffs, Array< OneD, NekDouble > &outarray)
Array< OneD, unsigned int > GetTraceInverseBoundaryMap(int fid, StdRegions::Orientation faceOrient=StdRegions::eNoOrientation, int P1=-1, int P2=-1)
virtual const NormalVector & v_GetTraceNormal(const int face) const
virtual DNekMatSharedPtr v_BuildVertexMatrix(const DNekScalMatSharedPtr &r_bnd)
virtual void v_DGDeriv(const int dir, const Array< OneD, const NekDouble > &incoeffs, Array< OneD, ExpansionSharedPtr > &FaceExp, Array< OneD, Array< OneD, NekDouble > > &faceCoeffs, Array< OneD, NekDouble > &out_d)
Evaluate coefficients of weak deriviative in the direction dir given the input coefficicents incoeffs...
std::vector< bool > m_requireNeg
Definition: Expansion3D.h:188
void GetInverseBoundaryMaps(Array< OneD, unsigned int > &vmap, Array< OneD, Array< OneD, unsigned int > > &emap, Array< OneD, Array< OneD, unsigned int > > &fmap)
void v_ReOrientTracePhysMap(const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1)
void AddFaceBoundaryInt(const int face, ExpansionSharedPtr &FaceExp, Array< OneD, NekDouble > &facePhys, Array< OneD, NekDouble > &outarray, const StdRegions::VarCoeffMap &varcoeffs=StdRegions::NullVarCoeffMap)
void GetPhysFaceVarCoeffsFromElement(const int face, ExpansionSharedPtr &FaceExp, const Array< OneD, const NekDouble > &varcoeff, Array< OneD, NekDouble > &outarray)
virtual void v_AddFaceNormBoundaryInt(const int face, const ExpansionSharedPtr &FaceExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
Array< OneD, unsigned int > GetEdgeInverseBoundaryMap(int eid)
std::map< int, NormalVector > m_faceNormals
Definition: Expansion3D.h:123
virtual DNekMatSharedPtr v_GenMatrix(const StdRegions::StdMatrixKey &mkey)
SpatialDomains::Geometry3DSharedPtr GetGeom3D() const
Definition: Expansion3D.h:191
virtual StdRegions::Orientation v_GetTraceOrient(int face)
virtual DNekMatSharedPtr v_BuildInverseTransformationMatrix(const DNekScalMatSharedPtr &transformationmatrix)
Build inverse and inverse transposed transformation matrix: and .
void AddHDGHelmholtzFaceTerms(const NekDouble tau, const int edge, Array< OneD, NekDouble > &facePhys, const StdRegions::VarCoeffMap &dirForcing, Array< OneD, NekDouble > &outarray)
Definition: Expansion3D.cpp:52
virtual DNekMatSharedPtr v_BuildTransformationMatrix(const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType)
void v_NormVectorIProductWRTBase(const Array< OneD, const Array< OneD, NekDouble > > &Fvec, Array< OneD, NekDouble > &outarray)
void SetFaceToGeomOrientation(const int face, Array< OneD, NekDouble > &inout)
Align face orientation with the geometry orientation.
virtual Array< OneD, NekDouble > v_GetnFacecdotMF(const int dir, const int face, ExpansionSharedPtr &FaceExp_f, const Array< OneD, const Array< OneD, NekDouble > > &normals, const StdRegions::VarCoeffMap &varcoeffs)
virtual void v_AddRobinMassMatrix(const int face, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat)
virtual void v_GetTracePhysVals(const int face, const StdRegions::StdExpansionSharedPtr &FaceExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient)
Extract the physical values along face face from inarray into outarray following the local face orien...
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:272
std::weak_ptr< Expansion2D > Expansion2DWeakPtr
Definition: Expansion1D.h:49
std::vector< Expansion3DSharedPtr > Expansion3DVector
Definition: Expansion3D.h:55
std::shared_ptr< Expansion > ExpansionSharedPtr
Definition: Expansion.h:68
std::shared_ptr< Expansion2D > Expansion2DSharedPtr
Definition: Expansion1D.h:47
std::weak_ptr< Expansion3D > Expansion3DWeakPtr
Definition: Expansion2D.h:51
std::shared_ptr< Expansion3D > Expansion3DSharedPtr
Definition: Expansion2D.h:49
std::shared_ptr< Geometry3D > Geometry3DSharedPtr
Definition: Geometry3D.h:52
std::shared_ptr< StdExpansion > StdExpansionSharedPtr
std::map< StdRegions::VarCoeffType, Array< OneD, NekDouble > > VarCoeffMap
Definition: StdRegions.hpp:272
static VarCoeffMap NullVarCoeffMap
Definition: StdRegions.hpp:273
The above copyright notice and this permission notice shall be included.
Definition: CoupledSolver.h:1
std::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
std::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:69
double NekDouble