Nektar++
Expansion1D.h
Go to the documentation of this file.
1 ///////////////////////////////////////////////////////////////////////////////
2 //
3 // File Expansion1D.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 // License for the specific language governing rights and limitations under
14 // Permission is hereby granted, free of charge, to any person obtaining a
15 // copy of this software and associated documentation files (the "Software"),
16 // to deal in the Software without restriction, including without limitation
17 // the rights to use, copy, modify, merge, publish, distribute, sublicense,
18 // and/or sell copies of the Software, and to permit persons to whom the
19 // Software is furnished to do so, subject to the following conditions:
20 //
21 // The above copyright notice and this permission notice shall be included
22 // in all copies or substantial portions of the Software.
23 //
24 // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
25 // OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
26 // FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
27 // THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
28 // LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
29 // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
30 // DEALINGS IN THE SOFTWARE.
31 //
32 // Description: Header file for Expansion1D routines
33 //
34 ///////////////////////////////////////////////////////////////////////////////
35 
36 #ifndef EXPANSION1D_H
37 #define EXPANSION1D_H
38 
39 #include <LocalRegions/Expansion.h>
42 #include <boost/weak_ptr.hpp>
44 
45 namespace Nektar
46 {
47  namespace LocalRegions
48  {
49  class Expansion2D;
50  typedef boost::shared_ptr<Expansion2D> Expansion2DSharedPtr;
51  typedef boost::weak_ptr<Expansion2D> Expansion2DWeakPtr;
52 
53  class Expansion1D;
54  typedef boost::shared_ptr<Expansion1D> Expansion1DSharedPtr;
55  typedef boost::weak_ptr<Expansion1D> Expansion1DWeakPtr;
56  typedef std::vector< Expansion1DSharedPtr > Expansion1DVector;
58 
59 
60  class Expansion1D: virtual public Expansion, virtual public StdRegions::StdExpansion1D
61  {
62  public:
64  {
65  m_elementEdgeLeft = -1;
66  m_elementEdgeRight = -1;
67  }
68 
70 
72  const int dir,
74  Array<OneD,NekDouble> &outarray);
75 
76  inline Expansion2DSharedPtr GetLeftAdjacentElementExp() const;
77 
78  inline Expansion2DSharedPtr GetRightAdjacentElementExp() const;
79 
80  inline int GetLeftAdjacentElementEdge() const;
81 
82  inline int GetRightAdjacentElementEdge() const;
83 
84  inline void SetAdjacentElementExp(
85  int edge,
86  Expansion2DSharedPtr &e);
87 
89  const NekDouble tau,
90  const Array<OneD,const NekDouble> &inarray,
91  Array<OneD,NekDouble> &outarray);
92 
94 
95  protected:
97 
98  virtual void v_AddRobinMassMatrix(
99  const int vert,
100  const Array<OneD, const NekDouble > &primCoeffs,
101  DNekMatSharedPtr &inoutmat);
102 
103  virtual void v_AddRobinEdgeContribution(
104  const int vert,
105  const Array<OneD, const NekDouble > &primCoeffs,
106  Array<OneD, NekDouble> &coeffs);
107 
108  private:
109  Expansion2DWeakPtr m_elementLeft;
110  Expansion2DWeakPtr m_elementRight;
113 
114  };
115 
116  inline Expansion2DSharedPtr Expansion1D::GetLeftAdjacentElementExp() const
117  {
118  ASSERTL1(m_elementLeft.lock().get(), "Left adjacent element not set.");
119  return m_elementLeft.lock();
120  }
121 
122  inline Expansion2DSharedPtr Expansion1D::GetRightAdjacentElementExp() const
123  {
124  ASSERTL1(m_elementLeft.lock().get(), "Right adjacent element not set.");
125  return m_elementRight.lock();
126  }
127 
129  {
130  return m_elementEdgeLeft;
131  }
132 
134  {
135  return m_elementEdgeRight;
136  }
137 
138  inline void Expansion1D::SetAdjacentElementExp(int edge, Expansion2DSharedPtr &e)
139  {
140  if (m_elementLeft.lock().get())
141  {
142  ASSERTL1(!m_elementRight.lock().get(),
143  "Both adjacent elements already set.");
144  m_elementRight = e;
145  m_elementEdgeRight = edge;
146  }
147  else
148  {
149  m_elementLeft = e;
150  m_elementEdgeLeft = edge;
151  }
152  }
153 
155  {
156  return boost::dynamic_pointer_cast<SpatialDomains::Geometry1D>(m_geom);
157  }
158  } //end of namespace
159 } //end of namespace
160 
161 #endif
162 
virtual void v_AddRobinMassMatrix(const int vert, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat)
virtual DNekMatSharedPtr v_GenMatrix(const StdRegions::StdMatrixKey &mkey)
Definition: Expansion1D.cpp:43
Expansion2DWeakPtr m_elementRight
Definition: Expansion1D.h:110
std::vector< Expansion1DSharedPtr > Expansion1DVector
Definition: Expansion1D.h:56
Expansion2DWeakPtr m_elementLeft
Definition: Expansion1D.h:109
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:125
boost::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:70
Expansion2DSharedPtr GetLeftAdjacentElementExp() const
Definition: Expansion1D.h:116
boost::weak_ptr< Expansion1D > Expansion1DWeakPtr
Definition: Expansion1D.h:55
Expansion2DSharedPtr GetRightAdjacentElementExp() const
Definition: Expansion1D.h:122
void SetAdjacentElementExp(int edge, Expansion2DSharedPtr &e)
Definition: Expansion1D.h:138
double NekDouble
Expansion1D(SpatialDomains::Geometry1DSharedPtr pGeom)
Definition: Expansion1D.h:63
std::vector< Expansion1DSharedPtr >::iterator Expansion1DVectorIter
Definition: Expansion1D.h:57
void AddNormTraceInt(const int dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
1D geometry information
Definition: Geometry1D.h:56
boost::shared_ptr< Geometry1D > Geometry1DSharedPtr
Definition: Geometry1D.h:48
boost::weak_ptr< Expansion2D > Expansion2DWeakPtr
Definition: Expansion1D.h:51
virtual void v_AddRobinEdgeContribution(const int vert, const Array< OneD, const NekDouble > &primCoeffs, Array< OneD, NekDouble > &coeffs)
SpatialDomains::Geometry1DSharedPtr GetGeom1D() const
Definition: Expansion1D.h:154
#define LOCAL_REGIONS_EXPORT
void AddHDGHelmholtzTraceTerms(const NekDouble tau, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
boost::shared_ptr< Expansion1D > Expansion1DSharedPtr
Definition: Expansion1D.h:53
boost::shared_ptr< Expansion2D > Expansion2DSharedPtr
Definition: Expansion1D.h:49