Nektar++
Public Member Functions | Protected Attributes | Private Member Functions | List of all members
Nektar::SpatialDomains::BoundaryConditions Class Reference

#include <Conditions.h>

Collaboration diagram for Nektar::SpatialDomains::BoundaryConditions:
Collaboration graph
[legend]

Public Member Functions

 BoundaryConditions (const LibUtilities::SessionReaderSharedPtr &pSession, const MeshGraphSharedPtr &meshGraph)
 
 BoundaryConditions (void)
 
 ~BoundaryConditions (void)
 
const BoundaryRegionCollection & GetBoundaryRegions (void) const
 
void AddBoundaryRegions (const int regionID, BoundaryRegionShPtr &bRegion)
 
const BoundaryConditionCollection & GetBoundaryConditions (void) const
 
void AddBoundaryConditions (const int regionID, BoundaryConditionMapShPtr &bCond)
 
const std::string GetVariable (unsigned int indx)
 

Protected Attributes

MeshGraphSharedPtr m_meshGraph
 The mesh graph to use for referencing geometry info. More...
 
LibUtilities::SessionReaderSharedPtr m_session
 
BoundaryRegionCollection m_boundaryRegions
 
BoundaryConditionCollection m_boundaryConditions
 

Private Member Functions

void Read (TiXmlElement *conditions)
 Read segments (and general MeshGraph) given TiXmlDocument. More...
 
void ReadBoundaryRegions (TiXmlElement *regions)
 
void ReadBoundaryConditions (TiXmlElement *conditions)
 

Detailed Description

Definition at line 217 of file Conditions.h.

Constructor & Destructor Documentation

Nektar::SpatialDomains::BoundaryConditions::BoundaryConditions ( const LibUtilities::SessionReaderSharedPtr &  pSession,
const MeshGraphSharedPtr &  meshGraph 
)

Definition at line 45 of file Conditions.cpp.

References m_session, and Read().

46  : m_meshGraph(meshGraph),
47  m_session (pSession)
48 
49  {
50  Read(m_session->GetElement("Nektar/Conditions"));
51  }
Nektar::SpatialDomains::BoundaryConditions::Read
void Read(TiXmlElement *conditions)
Read segments (and general MeshGraph) given TiXmlDocument.
Definition: Conditions.cpp:66
Nektar::SpatialDomains::BoundaryConditions::m_session
LibUtilities::SessionReaderSharedPtr m_session
Definition: Conditions.h:254
Nektar::SpatialDomains::BoundaryConditions::m_meshGraph
MeshGraphSharedPtr m_meshGraph
The mesh graph to use for referencing geometry info.
Definition: Conditions.h:253
Nektar::SpatialDomains::BoundaryConditions::BoundaryConditions ( void  )

Definition at line 54 of file Conditions.cpp.

55  {
56  }
Nektar::SpatialDomains::BoundaryConditions::~BoundaryConditions ( void  )

Definition at line 58 of file Conditions.cpp.

59  {
60  }

Member Function Documentation

void Nektar::SpatialDomains::BoundaryConditions::AddBoundaryConditions ( const int  regionID,
BoundaryConditionMapShPtr &  bCond 
)
inline

Definition at line 241 of file Conditions.h.

References m_boundaryConditions.

242  {
243  m_boundaryConditions[regionID] = bCond;
244  }
Nektar::SpatialDomains::BoundaryConditions::m_boundaryConditions
BoundaryConditionCollection m_boundaryConditions
Definition: Conditions.h:257
void Nektar::SpatialDomains::BoundaryConditions::AddBoundaryRegions ( const int  regionID,
BoundaryRegionShPtr &  bRegion 
)
inline

Definition at line 230 of file Conditions.h.

References m_boundaryRegions.

231  {
232  m_boundaryRegions[regionID] = bRegion;
233  }
Nektar::SpatialDomains::BoundaryConditions::m_boundaryRegions
BoundaryRegionCollection m_boundaryRegions
Definition: Conditions.h:256
const BoundaryConditionCollection& Nektar::SpatialDomains::BoundaryConditions::GetBoundaryConditions ( void  ) const
inline

Definition at line 235 of file Conditions.h.

References m_boundaryConditions.

Referenced by Nektar::MultiRegions::DisContField2D::FindPeriodicEdges(), Nektar::MultiRegions::DisContField3D::FindPeriodicFaces(), Nektar::MultiRegions::DisContField1D::FindPeriodicVertices(), Nektar::MultiRegions::DisContField3D::GenerateBoundaryConditionExpansion(), Nektar::MultiRegions::DisContField1D::GenerateBoundaryConditionExpansion(), Nektar::MultiRegions::DisContField2D::GenerateBoundaryConditionExpansion(), Nektar::MultiRegions::DisContField1D::GetDomainBCs(), Nektar::MultiRegions::DisContField1D::SetBoundaryConditionExpansion(), and Nektar::MultiRegions::DisContField3DHomogeneous1D::SetupBoundaryConditions().

236  {
237  return m_boundaryConditions;
238  }
Nektar::SpatialDomains::BoundaryConditions::m_boundaryConditions
BoundaryConditionCollection m_boundaryConditions
Definition: Conditions.h:257
const BoundaryRegionCollection& Nektar::SpatialDomains::BoundaryConditions::GetBoundaryRegions ( void  ) const
inline

Definition at line 225 of file Conditions.h.

References m_boundaryRegions.

Referenced by Extractlayerdata(), Nektar::MultiRegions::DisContField2D::FindPeriodicEdges(), Nektar::MultiRegions::DisContField3D::FindPeriodicFaces(), Nektar::MultiRegions::DisContField1D::FindPeriodicVertices(), Nektar::MultiRegions::DisContField3D::GenerateBoundaryConditionExpansion(), Nektar::MultiRegions::DisContField1D::GenerateBoundaryConditionExpansion(), Nektar::MultiRegions::DisContField2D::GenerateBoundaryConditionExpansion(), Nektar::MultiRegions::DisContField1D::GetDomainBCs(), main(), Nektar::Utilities::OutputFld::Process(), Nektar::MultiRegions::DisContField1D::SetBoundaryConditionExpansion(), Nektar::MultiRegions::DisContField3DHomogeneous2D::SetupBoundaryConditions(), Nektar::MultiRegions::DisContField3DHomogeneous1D::SetupBoundaryConditions(), Nektar::SolverUtils::FilterAeroForces::v_Initialise(), and Nektar::FilterMovingBody::v_Initialise().

226  {
227  return m_boundaryRegions;
228  }
Nektar::SpatialDomains::BoundaryConditions::m_boundaryRegions
BoundaryRegionCollection m_boundaryRegions
Definition: Conditions.h:256
const std::string Nektar::SpatialDomains::BoundaryConditions::GetVariable ( unsigned int  indx)
inline

Definition at line 246 of file Conditions.h.

References m_session.

247  {
248  return m_session->GetVariable(indx);
249  }
Nektar::SpatialDomains::BoundaryConditions::m_session
LibUtilities::SessionReaderSharedPtr m_session
Definition: Conditions.h:254
void Nektar::SpatialDomains::BoundaryConditions::Read ( TiXmlElement *  conditions)
private

Read segments (and general MeshGraph) given TiXmlDocument.

Definition at line 66 of file Conditions.cpp.

References ASSERTL0, ReadBoundaryConditions(), and ReadBoundaryRegions().

Referenced by BoundaryConditions().

67  {
68  ASSERTL0(conditions, "Unable to find CONDITIONS tag in file.");
69 
70  TiXmlElement *boundaryRegions = conditions->FirstChildElement("BOUNDARYREGIONS");
71 
72  if(boundaryRegions)
73  {
74  ReadBoundaryRegions(conditions);
75 
76  ReadBoundaryConditions(conditions);
77  }
78  }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::SpatialDomains::BoundaryConditions::ReadBoundaryRegions
void ReadBoundaryRegions(TiXmlElement *regions)
Definition: Conditions.cpp:84
Nektar::SpatialDomains::BoundaryConditions::ReadBoundaryConditions
void ReadBoundaryConditions(TiXmlElement *conditions)
Definition: Conditions.cpp:147
void Nektar::SpatialDomains::BoundaryConditions::ReadBoundaryConditions ( TiXmlElement *  conditions)
private

Definition at line 147 of file Conditions.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, ErrorUtil::ewarning, Nektar::StdRegions::find(), Nektar::ParseUtils::GenerateSeqVector(), Nektar::iterator, m_boundaryConditions, m_boundaryRegions, m_session, and NEKERROR.

Referenced by Read().

148  {
149  // Protect against multiple reads.
150  if(m_boundaryConditions.size() != 0)
151  {
152  return;
153  }
154 
155  // Read REGION tags
156  TiXmlElement *boundaryConditionsElement = conditions->FirstChildElement("BOUNDARYCONDITIONS");
157  ASSERTL0(boundaryConditionsElement, "Boundary conditions must be specified.");
158 
159  TiXmlElement *regionElement = boundaryConditionsElement->FirstChildElement("REGION");
160 
161  // Read R (Robin), D (Dirichlet), N (Neumann), P (Periodic) C(Cauchy) tags
162  while (regionElement)
163  {
164  BoundaryConditionMapShPtr boundaryConditions = MemoryManager<BoundaryConditionMap>::AllocateSharedPtr();
165 
166  int boundaryRegionID;
167  int err = regionElement->QueryIntAttribute("REF", &boundaryRegionID);
168  ASSERTL0(err == TIXML_SUCCESS, "Error reading boundary region reference.");
169 
170  ASSERTL0(m_boundaryConditions.count(boundaryRegionID) == 0,
171  "Boundary region '" + boost::lexical_cast<std::string>(boundaryRegionID)
172  + "' appears multiple times.");
173 
174  // Find the boundary region corresponding to this ID.
175  std::string boundaryRegionIDStr;
176  std::ostringstream boundaryRegionIDStrm(boundaryRegionIDStr);
177  boundaryRegionIDStrm << boundaryRegionID;
178 
179  ASSERTL0(m_boundaryRegions.count(boundaryRegionID) == 1,
180  "Boundary region " + boost::lexical_cast<
181  string>(boundaryRegionID)+ " not found");
182 
183  TiXmlElement *conditionElement = regionElement->FirstChildElement();
184  std::vector<std::string> vars = m_session->GetVariables();
185 
186  while (conditionElement)
187  {
188  // Check type.
189  std::string conditionType = conditionElement->Value();
190  std::string attrData;
191  bool isTimeDependent = false;
192 
193  // All have var specified, or else all variables are zero.
194  TiXmlAttribute *attr = conditionElement->FirstAttribute();
195 
196  std::vector<std::string>::iterator iter;
197  std::string attrName;
198 
199  attrData = conditionElement->Attribute("VAR");
200 
201  if (!attrData.empty())
202  {
203  iter = std::find(vars.begin(), vars.end(), attrData);
204  ASSERTL0(iter != vars.end(),
205  (std::string("Cannot find variable: ")
206  + attrData).c_str());
207  }
208 
209  if (conditionType == "N")
210  {
211  if (attrData.empty())
212  {
213  // All variables are Neumann and are set to zero.
214  for (std::vector<std::string>::iterator varIter = vars.begin();
215  varIter != vars.end(); ++varIter)
216  {
217  BoundaryConditionShPtr neumannCondition(MemoryManager<NeumannBoundaryCondition>::AllocateSharedPtr(m_session,"0.0"));
218  (*boundaryConditions)[*varIter] = neumannCondition;
219  }
220  }
221  else
222  {
223  // Use the iterator from above, which must point to the variable.
224  attr = attr->Next();
225 
226  if (attr)
227  {
228  std::string equation, userDefined, filename;
229 
230  while(attr)
231  {
232 
233  attrName = attr->Name();
234 
235  if (attrName=="USERDEFINEDTYPE")
236  {
237  // Do stuff for the user defined attribute
238  attrData = attr->Value();
239  ASSERTL0(!attrData.empty(),
240  "USERDEFINEDTYPE attribute must have associated value.");
241 
242  m_session->SubstituteExpressions(attrData);
243 
244  userDefined = attrData;
245  isTimeDependent = boost::iequals(attrData,"TimeDependent");
246  }
247  else if(attrName=="VALUE")
248  {
249  attrData = attr->Value();
250  ASSERTL0(!attrData.empty(),
251  "VALUE attribute must be specified.");
252 
253  m_session->SubstituteExpressions(attrData);
254 
255  equation = attrData;
256  }
257  else if(attrName=="FILE")
258  {
259  attrData = attr->Value();
260  ASSERTL0(!attrData.empty(), "FILE attribute must be specified.");
261 
262  m_session->SubstituteExpressions(attrData);
263 
264  filename = attrData;
265  }
266  else
267  {
268  ASSERTL0(false,
269  (std::string("Unknown boundary condition attribute: ") + attrName).c_str());
270  }
271  attr = attr->Next();
272  }
273 
274  BoundaryConditionShPtr neumannCondition(MemoryManager<NeumannBoundaryCondition>::AllocateSharedPtr(m_session, equation, userDefined, filename));
275  neumannCondition->SetIsTimeDependent(isTimeDependent);
276  (*boundaryConditions)[*iter] = neumannCondition;
277  }
278  else
279  {
280  // This variable's condition is zero.
281  BoundaryConditionShPtr neumannCondition(MemoryManager<NeumannBoundaryCondition>::AllocateSharedPtr(m_session, "0.0"));
282  (*boundaryConditions)[*iter] = neumannCondition;
283  }
284  }
285  }
286  else if (conditionType == "D")
287  {
288  if (attrData.empty())
289  {
290  // All variables are Dirichlet and are set to zero.
291  for (std::vector<std::string>::iterator varIter = vars.begin();
292  varIter != vars.end(); ++varIter)
293  {
294  BoundaryConditionShPtr dirichletCondition(MemoryManager<DirichletBoundaryCondition>::AllocateSharedPtr(m_session, "0"));
295  (*boundaryConditions)[*varIter] = dirichletCondition;
296  }
297  }
298  else
299  {
300  // Use the iterator from above, which must point to the variable.
301  attr = attr->Next();
302 
303  if (attr)
304  {
305  std::string equation, userDefined, filename;
306 
307  while(attr)
308  {
309  attrName = attr->Name();
310 
311  if (attrName=="USERDEFINEDTYPE") {
312 
313  // Do stuff for the user defined attribute
314  attrData = attr->Value();
315  ASSERTL0(!attrData.empty(), "USERDEFINEDTYPE attribute must have associated value.");
316 
317  m_session->SubstituteExpressions(attrData);
318 
319  userDefined = attrData;
320  isTimeDependent = boost::iequals(attrData,"TimeDependent");
321  }
322  else if(attrName=="VALUE")
323  {
324  attrData = attr->Value();
325  ASSERTL0(!attrData.empty(), "VALUE attribute must have associated value.");
326 
327  m_session->SubstituteExpressions(attrData);
328 
329  equation = attrData;
330  }
331  else if(attrName=="FILE")
332  {
333  attrData = attr->Value();
334  ASSERTL0(!attrData.empty(), "FILE attribute must be specified.");
335 
336  m_session->SubstituteExpressions(attrData);
337 
338  filename = attrData;
339  }
340  else
341  {
342  ASSERTL0(false,
343  (std::string("Unknown boundary condition attribute: ") + attrName).c_str());
344  }
345  attr = attr->Next();
346  }
347 
348  BoundaryConditionShPtr dirichletCondition(MemoryManager<DirichletBoundaryCondition>::AllocateSharedPtr(m_session, equation, userDefined, filename));
349  dirichletCondition->SetIsTimeDependent(isTimeDependent);
350  (*boundaryConditions)[*iter] = dirichletCondition;
351  }
352  else
353  {
354  // This variable's condition is zero.
355  BoundaryConditionShPtr dirichletCondition(MemoryManager<DirichletBoundaryCondition>::AllocateSharedPtr(m_session, "0"));
356  (*boundaryConditions)[*iter] = dirichletCondition;
357  }
358  }
359  }
360  else if (conditionType == "R") // Read du/dn + PRIMCOEFF u = VALUE
361  {
362  if (attrData.empty())
363  {
364  // All variables are Robin and are set to zero.
365  for (std::vector<std::string>::iterator varIter = vars.begin();
366  varIter != vars.end(); ++varIter)
367  {
368  BoundaryConditionShPtr robinCondition(MemoryManager<RobinBoundaryCondition>::AllocateSharedPtr(m_session, "0", "0"));
369  (*boundaryConditions)[*varIter] = robinCondition;
370  }
371  }
372  else
373  {
374  // Use the iterator from above, which must
375  // point to the variable. Read the A and
376  // B attributes.
377  attr = attr->Next();
378 
379  if (attr)
380  {
381  std::string attrName1;
382  std::string attrData1;
383  std::string equation1, equation2, userDefined;
384  std::string filename;
385 
386  while(attr){
387 
388  attrName1 = attr->Name();
389 
390  if (attrName1=="USERDEFINEDTYPE") {
391 
392  // Do stuff for the user defined attribute
393  attrData1 = attr->Value();
394  ASSERTL0(!attrData1.empty(), "USERDEFINEDTYPE attribute must have associated value.");
395 
396  m_session->SubstituteExpressions(attrData1);
397  userDefined = attrData1;
398  isTimeDependent = boost::iequals(attrData,"TimeDependent");
399  }
400  else if(attrName1 == "VALUE"){
401 
402  attrData1 = attr->Value();
403  ASSERTL0(!attrData1.empty(), "VALUE attributes must have associated values.");
404 
405  m_session->SubstituteExpressions(attrData1);
406 
407  equation1 = attrData1;
408 
409  attr = attr->Next();
410  ASSERTL0(attr, "Unable to read PRIMCOEFF attribute.");
411 
412  attrName1= attr->Name();
413  ASSERTL0(attrName1 == "PRIMCOEFF", (std::string("Unknown attribute: ") + attrName1).c_str());
414 
415  attrData1 = attr->Value();
416  ASSERTL0(!attrData1.empty(), "PRIMCOEFF attributes must have associated values.");
417 
418  m_session->SubstituteExpressions(attrData1);
419 
420  equation2 = attrData1;
421 
422  }
423  else if(attrName1=="FILE")
424  {
425  attrData1 = attr->Value();
426  ASSERTL0(!attrData1.empty(), "FILE attribute must be specified.");
427 
428  m_session->SubstituteExpressions(attrData1);
429 
430  filename = attrData1;
431  }
432  else
433  {
434  ASSERTL0(false, (std::string("Unknown boundary condition attribute: ") + attrName1).c_str());
435 
436  }
437  attr = attr->Next();
438  }
439 
440  BoundaryConditionShPtr robinCondition(MemoryManager<RobinBoundaryCondition>::AllocateSharedPtr(m_session, equation1, equation2, userDefined, filename));
441  (*boundaryConditions)[*iter] = robinCondition;
442  }
443  else
444  {
445  // This variable's condition is zero.
446  BoundaryConditionShPtr robinCondition(MemoryManager<RobinBoundaryCondition>::AllocateSharedPtr(m_session, "0", "0"));
447  robinCondition->SetIsTimeDependent(isTimeDependent);
448  (*boundaryConditions)[*iter] = robinCondition;
449  }
450  }
451  }
452  else if (conditionType == "P")
453  {
454  if (attrData.empty())
455  {
456  attr = attr->Next();
457 
458  if (attr)
459  {
460  attrName = attr->Name();
461 
462  ASSERTL0(attrName == "VALUE", (std::string("Unknown attribute: ") + attrName).c_str());
463 
464  attrData = attr->Value();
465  ASSERTL0(!attrData.empty(), "VALUE attribute must have associated value.");
466 
467  int beg = attrData.find_first_of("[");
468  int end = attrData.find_first_of("]");
469  std::string periodicBndRegionIndexStr = attrData.substr(beg+1,end-beg-1);
470  ASSERTL0(beg < end, (std::string("Error reading periodic boundary region definition for boundary region: ")
471  + boundaryRegionIDStrm.str()).c_str());
472 
473  vector<unsigned int> periodicBndRegionIndex;
474  bool parseGood = ParseUtils::GenerateSeqVector(periodicBndRegionIndexStr.c_str(), periodicBndRegionIndex);
475 
476  ASSERTL0(parseGood && (periodicBndRegionIndex.size()==1), (std::string("Unable to read periodic boundary condition for boundary region: ")
477  + boundaryRegionIDStrm.str()).c_str());
478 
479  BoundaryConditionShPtr periodicCondition(MemoryManager<PeriodicBoundaryCondition>::AllocateSharedPtr(periodicBndRegionIndex[0]));
480 
481  for (std::vector<std::string>::iterator varIter = vars.begin();
482  varIter != vars.end(); ++varIter)
483  {
484  (*boundaryConditions)[*varIter] = periodicCondition;
485  }
486  }
487  else
488  {
489  ASSERTL0(false, "Periodic boundary conditions should be explicitely defined");
490  }
491  }
492  else
493  {
494  // Use the iterator from above, which must point to the variable.
495  // Read the VALUE attribute. It is the next and only other attribute.
496  attr = attr->Next();
497 
498  if (attr)
499  {
500  attrName = attr->Name();
501 
502  ASSERTL0(attrName == "VALUE", (std::string("Unknown attribute: ") + attrName).c_str());
503 
504  attrData = attr->Value();
505  ASSERTL0(!attrData.empty(), "VALUE attribute must have associated value.");
506 
507  int beg = attrData.find_first_of("[");
508  int end = attrData.find_first_of("]");
509  std::string periodicBndRegionIndexStr = attrData.substr(beg+1,end-beg-1);
510  ASSERTL0(beg < end, (std::string("Error reading periodic boundary region definition for boundary region: ") + boundaryRegionIDStrm.str()).c_str());
511 
512  vector<unsigned int> periodicBndRegionIndex;
513  bool parseGood = ParseUtils::GenerateSeqVector(periodicBndRegionIndexStr.c_str(), periodicBndRegionIndex);
514 
515  ASSERTL0(parseGood && (periodicBndRegionIndex.size()==1), (std::string("Unable to read periodic boundary condition for boundary region: ") + boundaryRegionIDStrm.str()).c_str());
516 
517  BoundaryConditionShPtr periodicCondition(MemoryManager<PeriodicBoundaryCondition>::AllocateSharedPtr(periodicBndRegionIndex[0]));
518  (*boundaryConditions)[*iter] = periodicCondition;
519  }
520  else
521  {
522  ASSERTL0(false, "Periodic boundary conditions should be explicitely defined");
523  }
524  }
525  }
526  else if (conditionType == "C")
527  {
528  NEKERROR(ErrorUtil::ewarning, "Cauchy type boundary conditions not implemented.");
529  }
530 
531  conditionElement = conditionElement->NextSiblingElement();
532  }
533 
534  m_boundaryConditions[boundaryRegionID] = boundaryConditions;
535  regionElement = regionElement->NextSiblingElement("REGION");
536  }
537  }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
NEKERROR
#define NEKERROR(type, msg)
Assert Level 0 – Fundamental assert which is used whether in FULLDEBUG, DEBUG or OPT compilation mod...
Definition: ErrorUtil.hpp:158
Nektar::SpatialDomains::BoundaryConditions::m_boundaryRegions
BoundaryRegionCollection m_boundaryRegions
Definition: Conditions.h:256
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::ParseUtils::GenerateSeqVector
static bool GenerateSeqVector(const char *const str, std::vector< unsigned int > &vec)
Definition: ParseUtils.hpp:78
Nektar::iterator
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
Definition: MatrixOperations.cpp:96
Nektar::SpatialDomains::BoundaryConditionMapShPtr
boost::shared_ptr< BoundaryConditionMap > BoundaryConditionMapShPtr
Definition: Conditions.h:212
Nektar::SpatialDomains::BoundaryConditions::m_session
LibUtilities::SessionReaderSharedPtr m_session
Definition: Conditions.h:254
Nektar::SpatialDomains::BoundaryConditions::m_boundaryConditions
BoundaryConditionCollection m_boundaryConditions
Definition: Conditions.h:257
Nektar::StdRegions::find
InputIterator find(InputIterator first, InputIterator last, InputIterator startingpoint, const EqualityComparable &value)
Definition: StdRegions.hpp:312
Nektar::SpatialDomains::BoundaryConditionShPtr
boost::shared_ptr< BoundaryConditionBase > BoundaryConditionShPtr
Definition: Conditions.h:206
void Nektar::SpatialDomains::BoundaryConditions::ReadBoundaryRegions ( TiXmlElement *  regions)
private

All elements are of the form: "<B ID="#"> ... </B>", with ? being the element type.

Definition at line 84 of file Conditions.cpp.

References ASSERTL0, m_boundaryRegions, and m_meshGraph.

Referenced by Read().

85  {
86  // ensure boundary regions only read once per class definition
87  if(m_boundaryRegions.size() != 0)
88  {
89  return;
90  }
91 
92  TiXmlElement *boundaryRegions = conditions->FirstChildElement("BOUNDARYREGIONS");
93  ASSERTL0(boundaryRegions, "Unable to find BOUNDARYREGIONS block.");
94 
95  // See if we have boundary regions defined.
96  TiXmlElement *boundaryRegionsElement = boundaryRegions->FirstChildElement("B");
97 
98  while (boundaryRegionsElement)
99  {
100  /// All elements are of the form: "<B ID="#"> ... </B>", with
101  /// ? being the element type.
102  int indx;
103  int err = boundaryRegionsElement->QueryIntAttribute("ID", &indx);
104  ASSERTL0(err == TIXML_SUCCESS, "Unable to read attribute ID.");
105 
106  TiXmlNode* boundaryRegionChild = boundaryRegionsElement->FirstChild();
107  // This is primarily to skip comments that may be present.
108  // Comments appear as nodes just like elements.
109  // We are specifically looking for text in the body
110  // of the definition.
111  while(boundaryRegionChild && boundaryRegionChild->Type() != TiXmlNode::TINYXML_TEXT)
112  {
113  boundaryRegionChild = boundaryRegionChild->NextSibling();
114  }
115 
116  ASSERTL0(boundaryRegionChild, "Unable to read variable definition body.");
117  std::string boundaryRegionStr = boundaryRegionChild->ToText()->ValueStr();
118 
119  std::string::size_type indxBeg = boundaryRegionStr.find_first_of('[') + 1;
120  std::string::size_type indxEnd = boundaryRegionStr.find_last_of(']') - 1;
121 
122  ASSERTL0(indxBeg <= indxEnd, (std::string("Error reading boundary region definition:") + boundaryRegionStr).c_str());
123 
124  std::string indxStr = boundaryRegionStr.substr(indxBeg, indxEnd - indxBeg + 1);
125 
126  if (!indxStr.empty())
127  {
128  // Extract the composites from the string and return them in a list.
129  BoundaryRegionShPtr boundaryRegion(MemoryManager<BoundaryRegion>::AllocateSharedPtr());
130 
131  ASSERTL0(m_boundaryRegions.count(indx) == 0,
132  "Boundary region "+indxStr+ " defined more than "
133  "once!");
134 
135  m_meshGraph->GetCompositeList(indxStr, *boundaryRegion);
136  m_boundaryRegions[indx] = boundaryRegion;
137  }
138 
139  boundaryRegionsElement = boundaryRegionsElement->NextSiblingElement("B");
140  }
141  }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::SpatialDomains::BoundaryConditions::m_boundaryRegions
BoundaryRegionCollection m_boundaryRegions
Definition: Conditions.h:256
Nektar::SpatialDomains::BoundaryConditions::m_meshGraph
MeshGraphSharedPtr m_meshGraph
The mesh graph to use for referencing geometry info.
Definition: Conditions.h:253
Nektar::SpatialDomains::BoundaryRegionShPtr
boost::shared_ptr< BoundaryRegion > BoundaryRegionShPtr
Definition: Conditions.h:202

Member Data Documentation

BoundaryConditionCollection Nektar::SpatialDomains::BoundaryConditions::m_boundaryConditions
protected

Definition at line 257 of file Conditions.h.

Referenced by AddBoundaryConditions(), GetBoundaryConditions(), and ReadBoundaryConditions().

BoundaryRegionCollection Nektar::SpatialDomains::BoundaryConditions::m_boundaryRegions
protected

Definition at line 256 of file Conditions.h.

Referenced by AddBoundaryRegions(), GetBoundaryRegions(), ReadBoundaryConditions(), and ReadBoundaryRegions().

MeshGraphSharedPtr Nektar::SpatialDomains::BoundaryConditions::m_meshGraph
protected

The mesh graph to use for referencing geometry info.

Definition at line 253 of file Conditions.h.

Referenced by ReadBoundaryRegions().

LibUtilities::SessionReaderSharedPtr Nektar::SpatialDomains::BoundaryConditions::m_session
protected

Definition at line 254 of file Conditions.h.

Referenced by BoundaryConditions(), GetVariable(), and ReadBoundaryConditions().

Generated on Tue Dec 8 2015 23:45:40 for Nektar++ by   doxygen 1.8.9.1