Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Geometry.cpp
Go to the documentation of this file.
1 ////////////////////////////////////////////////////////////////////////////////
2 //
3 // File: Geometry.cpp
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: This file contains the base class implementation for the
33 // Geometry class.
34 //
35 //
36 ////////////////////////////////////////////////////////////////////////////////
37 #include <SpatialDomains/Geometry.h>
38 
39 namespace Nektar
40 {
41  namespace SpatialDomains
42  {
43 
44  // static class property
45  GeomFactorsVector Geometry::m_regGeomFactorsManager;
46 
47  Geometry::Geometry():
48  m_coordim(0),
49  m_geomFactorsState(eNotFilled),
50  m_state(eNotFilled),
51  m_shapeType(LibUtilities::eNoShapeType),
52  m_globalID(-1)
53  {
54  }
55 
56  Geometry::Geometry(const int coordim):
57  m_coordim(coordim),
58  m_geomFactorsState(eNotFilled),
59  m_state(eNotFilled),
60  m_shapeType(LibUtilities::eNoShapeType),
61  m_globalID(-1)
62  {
63  }
64 
65  Geometry::~Geometry()
66  {
67  }
68 
69  GeomFactorsSharedPtr Geometry::ValidateRegGeomFactor(
70  GeomFactorsSharedPtr geomFactor)
71  {
72  GeomFactorsSharedPtr returnval = geomFactor;
73 
74 /// \todo should this '#if 0' statement be removed?
75 #if 0
76  bool found = false;
77  if (geomFactor->GetGtype() == eRegular)
78  {
79  for (GeomFactorsVectorIter iter = m_regGeomFactorsManager.begin();
80  iter != m_regGeomFactorsManager.end();
81  ++iter)
82  {
83  if (**iter == *geomFactor)
84  {
85  returnval = *iter;
86  found = true;
87  break;
88  }
89  }
90 
91  if (!found)
92  {
93  m_regGeomFactorsManager.push_back(geomFactor);
94  returnval = geomFactor;
95  }
96  }
97 #endif
98  return returnval;
99  }
100 
101  bool SortByGlobalId(const boost::shared_ptr<Geometry>& lhs,
102  const boost::shared_ptr<Geometry>& rhs)
103  {
104  return lhs->GetGlobalID() < rhs->GetGlobalID();
105  }
106 
107  bool GlobalIdEquality(const boost::shared_ptr<Geometry>& lhs,
108  const boost::shared_ptr<Geometry>& rhs)
109  {
110  return lhs->GetGlobalID() == rhs->GetGlobalID();
111  }
112 
113  void Geometry::v_AddElmtConnected(int gvo_id, int locid)
114  {
115  NEKERROR(ErrorUtil::efatal,
116  "This function is only valid for shape type geometries");
117  }
118 
119  int Geometry::v_NumElmtConnected() const
120  {
121  NEKERROR(ErrorUtil::efatal,
122  "This function is only valid for shape type geometries");
123  return 0;
124  }
125 
126  bool Geometry::v_IsElmtConnected(int gvo_id, int locid) const
127  {
128  NEKERROR(ErrorUtil::efatal,
129  "This function is only valid for shape type geometries");
130  return false;
131  }
132 
133  int Geometry::v_GetVid(int i) const
134  {
135  NEKERROR(ErrorUtil::efatal,
136  "This function is only valid for shape type geometries");
137  return 0;
138  }
139 
140  int Geometry::v_GetEid(int i) const
141  {
142  NEKERROR(ErrorUtil::efatal,
143  "This function is only valid for shape type geometries");
144  return 0;
145  }
146 
147  int Geometry::v_GetFid(int i) const
148  {
149  NEKERROR(ErrorUtil::efatal,
150  "This function is only valid for expansion type geometries");
151  return 0;
152  }
153 
154  int Geometry::v_GetNumVerts() const
155  {
156  NEKERROR(ErrorUtil::efatal,
157  "This function is only valid for shape type geometries");
158  return 0;
159  }
160 
161  StdRegions::Orientation Geometry::v_GetEorient(const int i) const
162  {
163  NEKERROR(ErrorUtil::efatal,
164  "This function is not valid for this geometry.");
165  return StdRegions::eForwards;
166  }
167 
168  StdRegions::Orientation Geometry::v_GetPorient(const int i) const
169  {
170  NEKERROR(ErrorUtil::efatal,
171  "This function is not valid for this geometry.");
172  return StdRegions::eFwd;
173  }
174 
175  int Geometry::v_GetNumEdges() const
176  {
177  NEKERROR(ErrorUtil::efatal,
178  "This function is only valid for shape type geometries");
179  return 0;
180  }
181 
182  int Geometry::v_GetNumFaces() const
183  {
184  NEKERROR(ErrorUtil::efatal,
185  "This function is only valid for shape type geometries");
186  return 0;
187  }
188 
189  int Geometry::v_GetShapeDim() const
190  {
191  NEKERROR(ErrorUtil::efatal,
192  "This function is only valid for shape type geometries");
193  return 0;
194  }
195 
196  StdRegions::StdExpansionSharedPtr Geometry::v_GetXmap() const
197  {
198  return m_xmap;
199  }
200 
201  bool Geometry::v_ContainsPoint(
202  const Array<OneD, const NekDouble>& gloCoord,
203  NekDouble tol)
204  {
205  NEKERROR(ErrorUtil::efatal,
206  "This function has not been defined for this geometry");
207  return false;
208  }
209 
210  bool Geometry::v_ContainsPoint(
211  const Array<OneD, const NekDouble>& gloCoord,
212  Array<OneD, NekDouble>& locCoord,
213  NekDouble tol)
214  {
215  NEKERROR(ErrorUtil::efatal,
216  "This function has not been defined for this geometry");
217  return false;
218  }
219 
220  bool Geometry::v_ContainsPoint(
221  const Array<OneD, const NekDouble>& gloCoord,
222  Array<OneD, NekDouble>& locCoord,
223  NekDouble tol,
224  NekDouble &resid)
225  {
226  NEKERROR(ErrorUtil::efatal,
227  "This function has not been defined for this geometry");
228  return false;
229  }
230 
231  int Geometry::v_GetVertexEdgeMap(const int i, const int j) const
232  {
233  NEKERROR(ErrorUtil::efatal,
234  "This function has not been defined for this geometry");
235  return 0;
236  }
237 
238  int Geometry::v_GetVertexFaceMap(const int i, const int j) const
239  {
240  NEKERROR(ErrorUtil::efatal,
241  "This function has not been defined for this geometry");
242  return 0;
243  }
244 
245  int Geometry::v_GetEdgeFaceMap(const int i, const int j) const
246  {
247  NEKERROR(ErrorUtil::efatal,
248  "This function has not been defined for this geometry");
249  return 0;
250  }
251 
252  NekDouble Geometry::v_GetCoord(
253  const int i,
254  const Array<OneD,const NekDouble>& Lcoord)
255  {
256  NEKERROR(ErrorUtil::efatal,
257  "This function is only valid for expansion type geometries");
258  return 0.0;
259  }
260 
261  NekDouble Geometry::v_GetLocCoords(
262  const Array<OneD,const NekDouble> &coords,
263  Array<OneD,NekDouble> &Lcoords)
264  {
265  NEKERROR(ErrorUtil::efatal,
266  "This function is only valid for expansion type geometries");
267  return 0.0;
268  }
269 
270  void Geometry::v_FillGeom()
271  {
272  NEKERROR(ErrorUtil::efatal,
273  "This function is only valid for expansion type geometries");
274  }
275 
276  void Geometry::v_SetOwnData()
277  {
278  NEKERROR(ErrorUtil::efatal,
279  "This function is only valid for expansion type geometries");
280  }
281 
282  const LibUtilities::BasisSharedPtr Geometry::v_GetBasis(const int i)
283  {
284  NEKERROR(ErrorUtil::efatal,
285  "This function is only valid for shape type geometries");
286  LibUtilities::BasisSharedPtr returnval;
287  return returnval;
288  }
289 
290  int Geometry::v_GetCoordim() const
291  {
292  return m_coordim;
293  }
294 
295  }; //end of namespace
296 }; //end of namespace
297 
Generated on Sun Mar 15 2015 00:11:31 for Nektar++ by   doxygen 1.8.1.2