Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
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;
59 
60 
61  class Expansion1D: virtual public Expansion,
62  virtual public StdRegions::StdExpansion1D
63  {
64  public:
66  Geometry1DSharedPtr pGeom)
67  : Expansion(pGeom),
69  {
70  m_elementEdgeLeft = -1;
71  m_elementEdgeRight = -1;
72  }
73 
75 
77  const int dir,
79  Array<OneD,NekDouble> &outarray);
80 
81  inline Expansion2DSharedPtr GetLeftAdjacentElementExp() const;
82 
83  inline Expansion2DSharedPtr GetRightAdjacentElementExp() const;
84 
85  inline int GetLeftAdjacentElementEdge() const;
86 
87  inline int GetRightAdjacentElementEdge() const;
88 
89  inline void SetAdjacentElementExp(
90  int edge,
91  Expansion2DSharedPtr &e);
92 
94  const NekDouble tau,
95  const Array<OneD, const NekDouble> &inarray,
96  Array<OneD, NekDouble> &outarray);
97 
99 
100  protected:
102  const StdRegions::StdMatrixKey &mkey);
103 
104  virtual void v_AddRobinMassMatrix(
105  const int vert,
106  const Array<OneD, const NekDouble > &primCoeffs,
107  DNekMatSharedPtr &inoutmat);
108 
109  virtual void v_AddRobinEdgeContribution(
110  const int vert,
111  const Array<OneD, const NekDouble > &primCoeffs,
112  Array<OneD, NekDouble> &coeffs);
113 
114  private:
115  Expansion2DWeakPtr m_elementLeft;
116  Expansion2DWeakPtr m_elementRight;
119 
120  };
121 
122  inline Expansion2DSharedPtr Expansion1D::
124  {
125  ASSERTL1(m_elementLeft.lock().get(),
126  "Left adjacent element not set.");
127  return m_elementLeft.lock();
128  }
129 
130  inline Expansion2DSharedPtr Expansion1D::
132  {
133  ASSERTL1(m_elementLeft.lock().get(),
134  "Right adjacent element not set.");
135 
136  return m_elementRight.lock();
137  }
138 
140  {
141  return m_elementEdgeLeft;
142  }
143 
145  {
146  return m_elementEdgeRight;
147  }
148 
150  int edge,
151  Expansion2DSharedPtr &e)
152  {
153  if (m_elementLeft.lock().get())
154  {
155  ASSERTL1(!m_elementRight.lock().get(),
156  "Both adjacent elements already set.");
157 
158  m_elementRight = e;
159  m_elementEdgeRight = edge;
160  }
161  else
162  {
163  m_elementLeft = e;
164  m_elementEdgeLeft = edge;
165  }
166  }
167 
170  {
171  return boost::dynamic_pointer_cast<SpatialDomains
172  ::Geometry1D>(m_geom);
173  }
174  } //end of namespace
175 } //end of namespace
176 
177 #endif
178 
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:45
Expansion2DWeakPtr m_elementRight
Definition: Expansion1D.h:116
std::vector< Expansion1DSharedPtr > Expansion1DVector
Definition: Expansion1D.h:56
Expansion2DWeakPtr m_elementLeft
Definition: Expansion1D.h:115
boost::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:70
Expansion2DSharedPtr GetLeftAdjacentElementExp() const
Definition: Expansion1D.h:123
boost::weak_ptr< Expansion1D > Expansion1DWeakPtr
Definition: Expansion1D.h:55
Expansion2DSharedPtr GetRightAdjacentElementExp() const
Definition: Expansion1D.h:131
void SetAdjacentElementExp(int edge, Expansion2DSharedPtr &e)
Definition: Expansion1D.h:149
double NekDouble
Expansion1D(SpatialDomains::Geometry1DSharedPtr pGeom)
Definition: Expansion1D.h:65
std::vector< Expansion1DSharedPtr >::iterator Expansion1DVectorIter
Definition: Expansion1D.h:58
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:169
#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:218
boost::shared_ptr< Expansion1D > Expansion1DSharedPtr
Definition: Expansion1D.h:53
boost::shared_ptr< Expansion2D > Expansion2DSharedPtr
Definition: Expansion1D.h:49