Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Public Member Functions | Public Attributes | Private Attributes | List of all members
Nektar::FieldUtils::Field Struct Reference

#include <Field.hpp>

Collaboration diagram for Nektar::FieldUtils::Field:
Collaboration graph
[legend]

Public Member Functions

FIELD_UTILS_EXPORT Field ()
 
FIELD_UTILS_EXPORT ~Field ()
 
FIELD_UTILS_EXPORT
MultiRegions::ExpListSharedPtr 		SetUpFirstExpList (int NumHomogeneousDir, bool fldfilegiven=false)
 
FIELD_UTILS_EXPORT
LibUtilities::FieldIOSharedPtr 		FieldIOForFile (string filename)
 Construct a FieldIO object for the file filename. More...
 
FIELD_UTILS_EXPORT
MultiRegions::ExpListSharedPtr 		AppendExpList (int NumHomogeneousDir, string var="DefaultVar", bool NewField=false)
 

Public Attributes

bool m_verbose
 
vector
< LibUtilities::FieldDefinitionsSharedPtr > 		m_fielddef
 
vector< vector< double > > m_data
 
vector
< MultiRegions::ExpListSharedPtr > 		m_exp
 
bool m_declareExpansionAsContField
 
bool m_declareExpansionAsDisContField
 
bool m_requireBoundaryExpansion
 
bool m_useFFT
 
LibUtilities::CommSharedPtr m_comm
 
LibUtilities::SessionReaderSharedPtr m_session
 
SpatialDomains::MeshGraphSharedPtr m_graph
 
LibUtilities::PtsIOSharedPtr m_ptsIO
 
map< string, vector< string > > m_inputfiles
 
bool m_writeBndFld
 
vector< unsigned int > m_bndRegionsToWrite
 
bool m_fldToBnd
 
bool m_addNormals
 
bool m_setUpEquiSpacedFields
 
LibUtilities::PtsFieldSharedPtr m_fieldPts
 
MultiRegions::AssemblyMapCGSharedPtr m_locToGlobalMap
 
LibUtilities::FieldMetaDataMap m_fieldMetaDataMap
 

Private Attributes

map< string,
LibUtilities::FieldIOSharedPtr > 		m_fld
 Map to store FieldIO instances. Key is the reader type, value is the FieldIO object. More...
 

Detailed Description

Definition at line 71 of file Field.hpp.

Constructor & Destructor Documentation

FIELD_UTILS_EXPORT Nektar::FieldUtils::Field::Field ( )
inline

Definition at line 73 of file Field.hpp.

74  : m_verbose(false), m_declareExpansionAsContField(false),
75  m_declareExpansionAsDisContField(false),
76  m_requireBoundaryExpansion(false), m_writeBndFld(false),
77  m_fldToBnd(false), m_addNormals(false),
78  m_setUpEquiSpacedFields(false), m_fieldPts(LibUtilities::NullPtsField)
79  {
80  }
Nektar::FieldUtils::Field::m_declareExpansionAsDisContField
bool m_declareExpansionAsDisContField
Definition: Field.hpp:95
Nektar::FieldUtils::Field::m_declareExpansionAsContField
bool m_declareExpansionAsContField
Definition: Field.hpp:94
Nektar::LibUtilities::NullPtsField
static PtsFieldSharedPtr NullPtsField
Definition: PtsField.h:179
Nektar::FieldUtils::Field::m_fieldPts
LibUtilities::PtsFieldSharedPtr m_fieldPts
Definition: Field.hpp:114
FIELD_UTILS_EXPORT Nektar::FieldUtils::Field::~Field ( )
inline

Definition at line 82 of file Field.hpp.

83  {
84  if (m_comm)
85  {
86  m_comm->Finalise();
87  }
88  }
Nektar::FieldUtils::Field::m_comm
LibUtilities::CommSharedPtr m_comm
Definition: Field.hpp:101

Member Function Documentation

FIELD_UTILS_EXPORT MultiRegions::ExpListSharedPtr Nektar::FieldUtils::Field::AppendExpList ( int  NumHomogeneousDir,
string  var = "DefaultVar",
bool  NewField = false 
)
inline

Definition at line 442 of file Field.hpp.

References ASSERTL0.

444  {
445  if (var.compare("DefaultVar") == 0 && m_requireBoundaryExpansion)
446  {
447  if (m_session->GetVariables().size())
448  {
449  var = m_session->GetVariables()[0];
450  }
451  }
452  MultiRegions::ExpListSharedPtr tmp;
453  switch (m_graph->GetMeshDimension())
454  {
455  case 1:
456  {
457  if (NumHomogeneousDir == 1)
458  {
459  ASSERTL0(!(m_declareExpansionAsContField ||
460  m_declareExpansionAsDisContField),
461  "ContField2DHomogeneous1D or "
462  "DisContField2DHomogenenous1D has not been "
463  "implemented");
464 
465  MultiRegions::ExpList2DHomogeneous1DSharedPtr tmp2 =
466  boost::dynamic_pointer_cast<
467  MultiRegions::ExpList2DHomogeneous1D>(m_exp[0]);
468 
469  tmp = MemoryManager<MultiRegions::ExpList2DHomogeneous1D>::
470  AllocateSharedPtr(*tmp2);
471  }
472  else if (NumHomogeneousDir == 2)
473  {
474  if (m_declareExpansionAsContField)
475  {
476  MultiRegions::ContField3DHomogeneous2DSharedPtr tmp2 =
477  boost::dynamic_pointer_cast<
478  MultiRegions::ContField3DHomogeneous2D>(
479  m_exp[0]);
480 
481  tmp = MemoryManager<
482  MultiRegions::ContField3DHomogeneous2D>::
483  AllocateSharedPtr(*tmp2);
484  }
485  else if (m_declareExpansionAsDisContField)
486  {
487  MultiRegions::DisContField3DHomogeneous2DSharedPtr
488  tmp2 = boost::dynamic_pointer_cast<
489  MultiRegions::DisContField3DHomogeneous2D>(
490  m_exp[0]);
491 
492  tmp = MemoryManager<
493  MultiRegions::DisContField3DHomogeneous2D>::
494  AllocateSharedPtr(*tmp2);
495  }
496  else
497  {
498  MultiRegions::ExpList3DHomogeneous2DSharedPtr tmp2 =
499  boost::dynamic_pointer_cast<
500  MultiRegions::ExpList3DHomogeneous2D>(m_exp[0]);
501 
502  tmp = MemoryManager<
503  MultiRegions::ExpList3DHomogeneous2D>::
504  AllocateSharedPtr(*tmp2);
505  }
506  }
507  else
508  {
509  if (m_declareExpansionAsContField)
510  {
511  MultiRegions::ContField1DSharedPtr tmp2 =
512  boost::dynamic_pointer_cast<
513  MultiRegions::ContField1D>(m_exp[0]);
514 
515  tmp = MemoryManager<MultiRegions::ContField1D>::
516  AllocateSharedPtr(m_session, m_graph, var);
517  }
518  else if (m_declareExpansionAsDisContField)
519  {
520  MultiRegions::DisContField1DSharedPtr tmp2 =
521  boost::dynamic_pointer_cast<
522  MultiRegions::DisContField1D>(m_exp[0]);
523 
524  tmp = MemoryManager<MultiRegions::DisContField1D>::
525  AllocateSharedPtr(m_session, m_graph, var);
526  }
527  else
528  {
529  MultiRegions::ExpList1DSharedPtr tmp2 =
530  boost::dynamic_pointer_cast<
531  MultiRegions::ExpList1D>(m_exp[0]);
532 
533  tmp = MemoryManager<
534  MultiRegions::ExpList1D>::AllocateSharedPtr(*tmp2);
535  }
536  }
537  }
538  break;
539  case 2:
540  {
541  if (NumHomogeneousDir == 1)
542  {
543  if (m_declareExpansionAsContField)
544  {
545  if (NewField)
546  {
547  bool dealiasing = false;
548 
549  tmp = MemoryManager<
550  MultiRegions::ContField3DHomogeneous1D>::
551  AllocateSharedPtr(
552  m_session, m_exp[0]
553  ->GetHomogeneousBasis()
554  ->GetBasisKey(),
555  m_exp[0]->GetHomoLen(), m_useFFT,
556  dealiasing, m_graph, var);
557  }
558  else
559  {
560  MultiRegions::ContField3DHomogeneous1DSharedPtr
561  tmp2 = boost::dynamic_pointer_cast<
562  MultiRegions::ContField3DHomogeneous1D>(
563  m_exp[0]);
564 
565  ASSERTL0(tmp2, "Failed to type cast m_exp[0]");
566  tmp = MemoryManager<
567  MultiRegions::ContField3DHomogeneous1D>::
568  AllocateSharedPtr(*tmp2, m_graph, var);
569  }
570  }
571  else if (m_declareExpansionAsDisContField)
572  {
573  if (NewField)
574  {
575  bool dealiasing = false;
576 
577  tmp = MemoryManager<
578  MultiRegions::DisContField3DHomogeneous1D>::
579  AllocateSharedPtr(
580  m_session, m_exp[0]
581  ->GetHomogeneousBasis()
582  ->GetBasisKey(),
583  m_exp[0]->GetHomoLen(), m_useFFT,
584  dealiasing, m_graph, var);
585  }
586  else
587  {
588  MultiRegions::DisContField3DHomogeneous1DSharedPtr
589  tmp2 = boost::dynamic_pointer_cast<
590  MultiRegions::DisContField3DHomogeneous1D>(
591  m_exp[0]);
592  ASSERTL0(tmp2, "Failed to type cast m_exp[0]");
593 
594  tmp = MemoryManager<
595  MultiRegions::DisContField3DHomogeneous1D>::
596  AllocateSharedPtr(*tmp2);
597  }
598  }
599  else
600  {
601  if (NewField)
602  {
603  bool dealiasing = false;
604 
605  tmp = MemoryManager<
606  MultiRegions::ExpList3DHomogeneous1D>::
607  AllocateSharedPtr(
608  m_session, m_exp[0]
609  ->GetHomogeneousBasis()
610  ->GetBasisKey(),
611  m_exp[0]->GetHomoLen(), m_useFFT,
612  dealiasing, m_graph);
613  }
614  else
615  {
616  MultiRegions::ExpList3DHomogeneous1DSharedPtr tmp2 =
617  boost::dynamic_pointer_cast<
618  MultiRegions::ExpList3DHomogeneous1D>(
619  m_exp[0]);
620  ASSERTL0(tmp2, "Failed to type cast m_exp[0]");
621 
622  tmp = MemoryManager<
623  MultiRegions::ExpList3DHomogeneous1D>::
624  AllocateSharedPtr(*tmp2);
625  }
626  }
627  }
628  else
629  {
630  if (m_declareExpansionAsContField)
631  {
632  if (NewField)
633  {
634  tmp = MemoryManager<MultiRegions::ContField2D>::
635  AllocateSharedPtr(m_session, m_graph, var);
636  }
637  else // call copy constructor
638  {
639 
640  MultiRegions::ContField2DSharedPtr tmp2 =
641  boost::dynamic_pointer_cast<
642  MultiRegions::ContField2D>(m_exp[0]);
643 
644  tmp = MemoryManager<MultiRegions::ContField2D>::
645  AllocateSharedPtr(*tmp2, m_graph, var);
646  }
647  }
648  else if (m_declareExpansionAsDisContField)
649  {
650  if (NewField)
651  {
652  tmp = MemoryManager<MultiRegions::DisContField2D>::
653  AllocateSharedPtr(m_session, m_graph, var);
654  }
655  else // call copy constructor
656  {
657  MultiRegions::DisContField2DSharedPtr tmp2 =
658  boost::dynamic_pointer_cast<
659  MultiRegions::DisContField2D>(m_exp[0]);
660 
661  tmp = MemoryManager<MultiRegions::DisContField2D>::
662  AllocateSharedPtr(*tmp2, m_graph, var);
663  }
664  }
665  else
666  {
667  MultiRegions::ExpList2DSharedPtr tmp2 =
668  boost::dynamic_pointer_cast<
669  MultiRegions::ExpList2D>(m_exp[0]);
670 
671  tmp = MemoryManager<
672  MultiRegions::ExpList2D>::AllocateSharedPtr(*tmp2);
673  }
674  }
675  }
676  break;
677  case 3:
678  {
679  if (m_declareExpansionAsContField)
680  {
681  if (NewField)
682  {
683  tmp = MemoryManager<MultiRegions::ContField3D>::
684  AllocateSharedPtr(m_session, m_graph, var);
685  }
686  else
687  {
688  MultiRegions::ContField3DSharedPtr tmp2 =
689  boost::dynamic_pointer_cast<
690  MultiRegions::ContField3D>(m_exp[0]);
691 
692  tmp = MemoryManager<MultiRegions::ContField3D>::
693  AllocateSharedPtr(*tmp2, m_graph, var);
694 
695  m_locToGlobalMap = tmp2->GetLocalToGlobalMap();
696  }
697  }
698  else if (m_declareExpansionAsDisContField)
699  {
700  if (NewField)
701  {
702  tmp = MemoryManager<MultiRegions::DisContField3D>::
703  AllocateSharedPtr(m_session, m_graph, var);
704  }
705  else
706  {
707  MultiRegions::DisContField3DSharedPtr tmp2 =
708  boost::dynamic_pointer_cast<
709  MultiRegions::DisContField3D>(m_exp[0]);
710 
711  tmp = MemoryManager<MultiRegions::DisContField3D>::
712  AllocateSharedPtr(*tmp2, m_graph, var);
713  }
714  }
715  else
716  {
717  MultiRegions::ExpList3DSharedPtr tmp2 =
718  boost::dynamic_pointer_cast<MultiRegions::ExpList3D>(
719  m_exp[0]);
720 
721  tmp = MemoryManager<
722  MultiRegions::ExpList3D>::AllocateSharedPtr(*tmp2);
723  }
724  }
725  break;
726  default:
727  ASSERTL0(false, "Expansion dimension not recognised");
728  break;
729  }
730 
731  return tmp;
732  }
Nektar::FieldUtils::Field::m_exp
vector< MultiRegions::ExpListSharedPtr > m_exp
Definition: Field.hpp:92
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
Nektar::MultiRegions::ContField1DSharedPtr
boost::shared_ptr< ContField1D > ContField1DSharedPtr
Definition: ContField1D.h:236
Nektar::MultiRegions::ContField3DHomogeneous1DSharedPtr
boost::shared_ptr< ContField3DHomogeneous1D > ContField3DHomogeneous1DSharedPtr
Definition: ContField3DHomogeneous1D.h:106
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::FieldUtils::Field::m_declareExpansionAsDisContField
bool m_declareExpansionAsDisContField
Definition: Field.hpp:95
Nektar::MultiRegions::DisContField3DHomogeneous1DSharedPtr
boost::shared_ptr< DisContField3DHomogeneous1D > DisContField3DHomogeneous1DSharedPtr
Definition: DisContField3DHomogeneous1D.h:299
Nektar::FieldUtils::Field::m_session
LibUtilities::SessionReaderSharedPtr m_session
Definition: Field.hpp:102
Nektar::MultiRegions::ContField2DSharedPtr
boost::shared_ptr< ContField2D > ContField2DSharedPtr
Definition: ContField2D.h:287
Nektar::FieldUtils::Field::m_graph
SpatialDomains::MeshGraphSharedPtr m_graph
Definition: Field.hpp:103
Nektar::MultiRegions::ExpList3DHomogeneous2DSharedPtr
boost::shared_ptr< ExpList3DHomogeneous2D > ExpList3DHomogeneous2DSharedPtr
Shared pointer to an ExpList3DHomogeneous2D object.
Definition: ExpList3DHomogeneous2D.h:50
Nektar::FieldUtils::Field::m_declareExpansionAsContField
bool m_declareExpansionAsContField
Definition: Field.hpp:94
Nektar::MultiRegions::DisContField3DHomogeneous2DSharedPtr
boost::shared_ptr< DisContField3DHomogeneous2D > DisContField3DHomogeneous2DSharedPtr
Definition: DisContField3DHomogeneous2D.h:156
Nektar::MultiRegions::ContField3DHomogeneous2DSharedPtr
boost::shared_ptr< ContField3DHomogeneous2D > ContField3DHomogeneous2DSharedPtr
Definition: ContField3DHomogeneous2D.h:99
Nektar::MultiRegions::ExpList1DSharedPtr
boost::shared_ptr< ExpList1D > ExpList1DSharedPtr
Shared pointer to an ExpList1D object.
Definition: ExpList1D.h:50
Nektar::MultiRegions::ExpListSharedPtr
boost::shared_ptr< ExpList > ExpListSharedPtr
Shared pointer to an ExpList object.
Definition: LocTraceToTraceMap.h:56
Nektar::MultiRegions::DisContField2DSharedPtr
boost::shared_ptr< DisContField2D > DisContField2DSharedPtr
Definition: DisContField2D.h:301
Nektar::MultiRegions::DisContField1DSharedPtr
boost::shared_ptr< DisContField1D > DisContField1DSharedPtr
Definition: DisContField1D.h:277
Nektar::MultiRegions::ExpList2DSharedPtr
boost::shared_ptr< ExpList2D > ExpList2DSharedPtr
Shared pointer to an ExpList2D object.
Definition: ExpList2D.h:49
Nektar::MultiRegions::DisContField3DSharedPtr
boost::shared_ptr< DisContField3D > DisContField3DSharedPtr
Definition: DisContField3D.h:271
Nektar::MultiRegions::ExpList2DHomogeneous1DSharedPtr
boost::shared_ptr< ExpList2DHomogeneous1D > ExpList2DHomogeneous1DSharedPtr
Shared pointer to an ExpList2DHomogeneous1D object.
Definition: ExpList2DHomogeneous1D.h:50
Nektar::MultiRegions::ExpList3DSharedPtr
boost::shared_ptr< ExpList3D > ExpList3DSharedPtr
Shared pointer to an ExpList3D object.
Definition: ExpList3D.h:115
Nektar::FieldUtils::Field::m_locToGlobalMap
MultiRegions::AssemblyMapCGSharedPtr m_locToGlobalMap
Definition: Field.hpp:116
Nektar::MultiRegions::ExpList3DHomogeneous1DSharedPtr
boost::shared_ptr< ExpList3DHomogeneous1D > ExpList3DHomogeneous1DSharedPtr
Shared pointer to an ExpList3DHomogeneous1D object.
Definition: ExpList3DHomogeneous1D.h:50
Nektar::MultiRegions::ContField3DSharedPtr
boost::shared_ptr< ContField3D > ContField3DSharedPtr
Definition: ContField3D.h:206
FIELD_UTILS_EXPORT LibUtilities::FieldIOSharedPtr Nektar::FieldUtils::Field::FieldIOForFile ( string  filename)
inline

Construct a FieldIO object for the file filename.

This routine constructs an appropriate FieldIO object for a filename through the LibUtilities::FieldIO::GetFileType function to detect the file format. The result is then cached in Field::m_fld to avoid needing to repeatedly construct the object.

Parameters
filenameFilename to open.
Returns
Reader for filename.

Definition at line 420 of file Field.hpp.

References Nektar::LibUtilities::NekFactory< tKey, tBase, >::CreateInstance(), Nektar::LibUtilities::GetCommFactory(), Nektar::LibUtilities::GetFieldIOFactory(), and Nektar::iterator.

422  {
423  LibUtilities::CommSharedPtr c = m_session ? m_session->GetComm() :
424  LibUtilities::GetCommFactory().CreateInstance("Serial", 0, 0);
425  string fmt = LibUtilities::FieldIO::GetFileType(filename, c);
426  map<string, LibUtilities::FieldIOSharedPtr>::iterator it =
427  m_fld.find(fmt);
428 
429  if (it == m_fld.end())
430  {
431  LibUtilities::FieldIOSharedPtr fld =
432  LibUtilities::GetFieldIOFactory().CreateInstance(fmt, c, true);
433  m_fld[fmt] = fld;
434  return fld;
435  }
436  else
437  {
438  return it->second;
439  }
440  }
Nektar::LibUtilities::NekFactory::CreateInstance
tBaseSharedPtr CreateInstance(tKey idKey BOOST_PP_COMMA_IF(MAX_PARAM) BOOST_PP_ENUM_BINARY_PARAMS(MAX_PARAM, tParam, x))
Create an instance of the class referred to by idKey.
Definition: NekFactory.hpp:162
Nektar::FieldUtils::Field::m_session
LibUtilities::SessionReaderSharedPtr m_session
Definition: Field.hpp:102
Nektar::LibUtilities::GetFieldIOFactory
FieldIOFactory & GetFieldIOFactory()
Returns the FieldIO factory.
Definition: FieldIO.cpp:74
Nektar::LibUtilities::GetCommFactory
CommFactory & GetCommFactory()
Definition: Comm.cpp:61
Nektar::LibUtilities::CommSharedPtr
boost::shared_ptr< Comm > CommSharedPtr
Pointer to a Communicator object.
Definition: Comm.h:55
Nektar::LibUtilities::FieldIO::GetFileType
static const std::string GetFileType(const std::string &filename, CommSharedPtr comm)
Determine file type of given input file.
Definition: FieldIO.cpp:101
Nektar::LibUtilities::FieldIOSharedPtr
boost::shared_ptr< FieldIO > FieldIOSharedPtr
Definition: FieldIO.h:309
Nektar::FieldUtils::Field::m_fld
map< string, LibUtilities::FieldIOSharedPtr > m_fld
Map to store FieldIO instances. Key is the reader type, value is the FieldIO object.
Definition: Field.hpp:737
Nektar::iterator
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
Definition: MatrixOperations.cpp:96
FIELD_UTILS_EXPORT MultiRegions::ExpListSharedPtr Nektar::FieldUtils::Field::SetUpFirstExpList ( int  NumHomogeneousDir,
bool  fldfilegiven = false 
)
inline

Definition at line 120 of file Field.hpp.

References ASSERTL0, Nektar::LibUtilities::eFourier, Nektar::LibUtilities::eFourierEvenlySpaced, Nektar::LibUtilities::eFourierHalfModeRe, Nektar::LibUtilities::eFourierSingleMode, and Nektar::LibUtilities::ePolyEvenlySpaced.

122  {
123 
124  MultiRegions::ExpListSharedPtr exp;
125 
126  // Set up expansion list
127  int expdim = m_graph->GetMeshDimension();
128  bool dealiasing = false;
129 
130  m_session->MatchSolverInfo("USEFFT", "FFTW", m_useFFT, false);
131 
132  switch (expdim)
133  {
134  case 1:
135  {
136  ASSERTL0(NumHomogeneousDir <= 2,
137  "Quasi-3D approach is only set up for 1 or 2 "
138  "homogeneous directions");
139 
140  if (NumHomogeneousDir == 1)
141  {
142  MultiRegions::ExpList2DHomogeneous1DSharedPtr Exp2DH1;
143 
144  // Define Homogeneous expansion
145  int nplanes;
146  NekDouble ly;
147  LibUtilities::BasisType btype;
148 
149  if (fldfilegiven)
150  {
151  nplanes = m_fielddef[0]->m_numModes[1];
152  ly = m_fielddef[0]->m_homogeneousLengths[0];
153  btype = m_fielddef[0]->m_basis[1];
154  }
155  else
156  {
157  m_session->LoadParameter("HomModesZ", nplanes);
158  m_session->LoadParameter("LY", ly);
159  btype = LibUtilities::eFourier;
160  }
161 
162  // Choose points to be at evenly spaced points at
163  // nplanes points
164  const LibUtilities::PointsKey Pkey(
165  nplanes, LibUtilities::eFourierEvenlySpaced);
166 
167  const LibUtilities::BasisKey Bkey(btype, nplanes, Pkey);
168 
169  if (m_declareExpansionAsContField ||
170  m_declareExpansionAsDisContField)
171  {
172  ASSERTL0(false, "ContField2DHomogeneous1D or "
173  "DisContField2DHomogenenous1D has "
174  "not been implemented");
175  }
176 
177  Exp2DH1 =
178  MemoryManager<MultiRegions::ExpList2DHomogeneous1D>::
179  AllocateSharedPtr(m_session, Bkey, ly, m_useFFT,
180  dealiasing, m_graph);
181  exp = Exp2DH1;
182  }
183  else if (NumHomogeneousDir == 2)
184  {
185  MultiRegions::ExpList3DHomogeneous2DSharedPtr Exp3DH2;
186 
187  int nylines, nzlines;
188  NekDouble ly, lz;
189  LibUtilities::BasisType btype1, btype2;
190 
191  if (fldfilegiven)
192  {
193  nylines = m_fielddef[0]->m_numModes[1];
194  nzlines = m_fielddef[0]->m_numModes[2];
195  ly = m_fielddef[0]->m_homogeneousLengths[0];
196  lz = m_fielddef[0]->m_homogeneousLengths[1];
197  btype1 = m_fielddef[0]->m_basis[1];
198  btype2 = m_fielddef[0]->m_basis[2];
199  }
200  else
201  {
202  m_session->LoadParameter("HomModesY", nylines);
203  m_session->LoadParameter("HomModesZ", nzlines);
204  m_session->LoadParameter("LY", ly);
205  m_session->LoadParameter("LZ", lz);
206  btype1 = LibUtilities::eFourier;
207  btype2 = LibUtilities::eFourier;
208  }
209 
210  // Choose points to be at evenly spaced points at
211  // nplanes points
212  const LibUtilities::PointsKey PkeyY(
213  nylines, LibUtilities::eFourierEvenlySpaced);
214  const LibUtilities::BasisKey BkeyY(btype1, nylines, PkeyY);
215 
216  const LibUtilities::PointsKey PkeyZ(
217  nzlines, LibUtilities::eFourierEvenlySpaced);
218  const LibUtilities::BasisKey BkeyZ(btype2, nzlines, PkeyZ);
219 
220  if (m_declareExpansionAsContField)
221  {
222  Exp3DH2 = MemoryManager<
223  MultiRegions::ContField3DHomogeneous2D>::
224  AllocateSharedPtr(m_session, BkeyY, BkeyZ, ly, lz,
225  m_useFFT, dealiasing, m_graph,
226  m_session->GetVariable(0));
227  }
228  else if (m_declareExpansionAsDisContField)
229  {
230  Exp3DH2 = MemoryManager<
231  MultiRegions::DisContField3DHomogeneous2D>::
232  AllocateSharedPtr(m_session, BkeyY, BkeyZ, ly, lz,
233  m_useFFT, dealiasing, m_graph,
234  m_session->GetVariable(0));
235  }
236  else
237  {
238  Exp3DH2 = MemoryManager<
239  MultiRegions::ExpList3DHomogeneous2D>::
240  AllocateSharedPtr(m_session, BkeyY, BkeyZ, ly, lz,
241  m_useFFT, dealiasing, m_graph);
242  }
243 
244  exp = Exp3DH2;
245  }
246  else
247  {
248  MultiRegions::ExpList1DSharedPtr Exp1D;
249 
250  if (m_declareExpansionAsContField)
251  {
252  Exp1D = MemoryManager<MultiRegions::ContField1D>::
253  AllocateSharedPtr(m_session, m_graph,
254  m_session->GetVariable(0));
255  }
256  else if (m_declareExpansionAsDisContField)
257  {
258  Exp1D = MemoryManager<MultiRegions::DisContField1D>::
259  AllocateSharedPtr(m_session, m_graph,
260  m_session->GetVariable(0));
261  }
262  else
263  {
264  Exp1D = MemoryManager<MultiRegions::ExpList1D>::
265  AllocateSharedPtr(m_session, m_graph);
266  }
267 
268  exp = Exp1D;
269  }
270  }
271  break;
272  case 2:
273  {
274  ASSERTL0(NumHomogeneousDir <= 1,
275  "NumHomogeneousDir is only set up for 1");
276 
277  if (NumHomogeneousDir == 1)
278  {
279  MultiRegions::ExpList3DHomogeneous1DSharedPtr Exp3DH1;
280 
281  // Define Homogeneous expansion
282  int nplanes;
283  NekDouble lz;
284  LibUtilities::BasisType btype;
285  LibUtilities::PointsType ptype =
286  LibUtilities::eFourierEvenlySpaced;
287 
288  if (fldfilegiven)
289  {
290  nplanes = m_fielddef[0]->m_numModes[2];
291  lz = m_fielddef[0]->m_homogeneousLengths[0];
292  btype = m_fielddef[0]->m_basis[2];
293 
294  if (btype == LibUtilities::eFourierSingleMode)
295  {
296  btype = LibUtilities::eFourier;
297  m_fielddef[0]->m_basis[2] =
298  LibUtilities::eFourierSingleMode;
299  if (nplanes <= 2)
300  {
301  nplanes = 4;
302  }
303  }
304  else if (btype == LibUtilities::eFourierHalfModeRe &&
305  nplanes == 1)
306  {
307  ptype = LibUtilities::ePolyEvenlySpaced;
308  }
309  }
310  else
311  {
312  m_session->LoadParameter("HomModesZ", nplanes);
313  m_session->LoadParameter("LZ", lz);
314  btype = LibUtilities::eFourier;
315  }
316 
317  // Choose points to be at evenly spaced points at
318  // nplanes points
319  const LibUtilities::PointsKey Pkey(
320  nplanes, ptype);
321 
322  const LibUtilities::BasisKey Bkey(btype, nplanes, Pkey);
323 
324  if (m_declareExpansionAsContField)
325  {
326  Exp3DH1 = MemoryManager<
327  MultiRegions::ContField3DHomogeneous1D>::
328  AllocateSharedPtr(m_session, Bkey, lz, m_useFFT,
329  dealiasing, m_graph,
330  m_session->GetVariable(0));
331  }
332  else if (m_declareExpansionAsDisContField)
333  {
334  Exp3DH1 = MemoryManager<
335  MultiRegions::DisContField3DHomogeneous1D>::
336  AllocateSharedPtr(m_session, Bkey, lz, m_useFFT,
337  dealiasing, m_graph,
338  m_session->GetVariable(0));
339  }
340  else
341  {
342  Exp3DH1 = MemoryManager<
343  MultiRegions::ExpList3DHomogeneous1D>::
344  AllocateSharedPtr(m_session, Bkey, lz, m_useFFT,
345  dealiasing, m_graph);
346  }
347  exp = Exp3DH1;
348  }
349  else
350  {
351  MultiRegions::ExpList2DSharedPtr Exp2D;
352 
353  if (m_declareExpansionAsContField)
354  {
355  Exp2D = MemoryManager<MultiRegions::ContField2D>::
356  AllocateSharedPtr(m_session, m_graph,
357  m_session->GetVariable(0));
358  }
359  else if (m_declareExpansionAsDisContField)
360  {
361  Exp2D = MemoryManager<MultiRegions::DisContField2D>::
362  AllocateSharedPtr(m_session, m_graph,
363  m_session->GetVariable(0));
364  }
365  else
366  {
367  Exp2D = MemoryManager<MultiRegions::ExpList2D>::
368  AllocateSharedPtr(m_session, m_graph);
369  }
370 
371  exp = Exp2D;
372  }
373  }
374  break;
375  case 3:
376  {
377  MultiRegions::ExpList3DSharedPtr Exp3D;
378 
379  if (m_declareExpansionAsContField)
380  {
381  Exp3D = MemoryManager<MultiRegions::ContField3D>::
382  AllocateSharedPtr(m_session, m_graph,
383  m_session->GetVariable(0));
384  }
385  else if (m_declareExpansionAsDisContField)
386  {
387  Exp3D = MemoryManager<MultiRegions::DisContField3D>::
388  AllocateSharedPtr(m_session, m_graph,
389  m_session->GetVariable(0));
390  }
391  else
392  {
393  Exp3D = MemoryManager<
394  MultiRegions::ExpList3D>::AllocateSharedPtr(m_session,
395  m_graph);
396  }
397 
398  exp = Exp3D;
399  }
400  break;
401  default:
402  ASSERTL0(false, "Expansion dimension not recognised");
403  break;
404  }
405 
406  return exp;
407  };
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::FieldUtils::Field::m_declareExpansionAsDisContField
bool m_declareExpansionAsDisContField
Definition: Field.hpp:95
Nektar::FieldUtils::Field::m_session
LibUtilities::SessionReaderSharedPtr m_session
Definition: Field.hpp:102
Nektar::LibUtilities::eFourier
Fourier Expansion .
Definition: BasisType.h:52
Nektar::FieldUtils::Field::m_graph
SpatialDomains::MeshGraphSharedPtr m_graph
Definition: Field.hpp:103
Nektar::MultiRegions::ExpList3DHomogeneous2DSharedPtr
boost::shared_ptr< ExpList3DHomogeneous2D > ExpList3DHomogeneous2DSharedPtr
Shared pointer to an ExpList3DHomogeneous2D object.
Definition: ExpList3DHomogeneous2D.h:50
Nektar::FieldUtils::Field::m_declareExpansionAsContField
bool m_declareExpansionAsContField
Definition: Field.hpp:94
Nektar::LibUtilities::ePolyEvenlySpaced
1D Evenly-spaced points using Lagrange polynomial
Definition: PointsType.h:65
Nektar::LibUtilities::eFourierEvenlySpaced
1D Evenly-spaced points using Fourier Fit
Definition: PointsType.h:66
Nektar::LibUtilities::eFourierHalfModeRe
Fourier Modified expansions with just the real part of the first mode .
Definition: BasisType.h:59
Nektar::MultiRegions::ExpList1DSharedPtr
boost::shared_ptr< ExpList1D > ExpList1DSharedPtr
Shared pointer to an ExpList1D object.
Definition: ExpList1D.h:50
Nektar::MultiRegions::ExpListSharedPtr
boost::shared_ptr< ExpList > ExpListSharedPtr
Shared pointer to an ExpList object.
Definition: LocTraceToTraceMap.h:56
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:44
Nektar::MultiRegions::ExpList2DSharedPtr
boost::shared_ptr< ExpList2D > ExpList2DSharedPtr
Shared pointer to an ExpList2D object.
Definition: ExpList2D.h:49
Nektar::MultiRegions::ExpList2DHomogeneous1DSharedPtr
boost::shared_ptr< ExpList2DHomogeneous1D > ExpList2DHomogeneous1DSharedPtr
Shared pointer to an ExpList2DHomogeneous1D object.
Definition: ExpList2DHomogeneous1D.h:50
Nektar::MultiRegions::ExpList3DSharedPtr
boost::shared_ptr< ExpList3D > ExpList3DSharedPtr
Shared pointer to an ExpList3D object.
Definition: ExpList3D.h:115
Nektar::LibUtilities::eFourierSingleMode
Fourier ModifiedExpansion with just the first mode .
Definition: BasisType.h:58
Nektar::MultiRegions::ExpList3DHomogeneous1DSharedPtr
boost::shared_ptr< ExpList3DHomogeneous1D > ExpList3DHomogeneous1DSharedPtr
Shared pointer to an ExpList3DHomogeneous1D object.
Definition: ExpList3DHomogeneous1D.h:50
Nektar::FieldUtils::Field::m_fielddef
vector< LibUtilities::FieldDefinitionsSharedPtr > m_fielddef
Definition: Field.hpp:90

Member Data Documentation

bool Nektar::FieldUtils::Field::m_addNormals

Definition at line 110 of file Field.hpp.

vector<unsigned int> Nektar::FieldUtils::Field::m_bndRegionsToWrite

Definition at line 108 of file Field.hpp.

LibUtilities::CommSharedPtr Nektar::FieldUtils::Field::m_comm

Definition at line 101 of file Field.hpp.

vector<vector<double> > Nektar::FieldUtils::Field::m_data

Definition at line 91 of file Field.hpp.

bool Nektar::FieldUtils::Field::m_declareExpansionAsContField

Definition at line 94 of file Field.hpp.

bool Nektar::FieldUtils::Field::m_declareExpansionAsDisContField

Definition at line 95 of file Field.hpp.

vector<MultiRegions::ExpListSharedPtr> Nektar::FieldUtils::Field::m_exp

Definition at line 92 of file Field.hpp.

vector<LibUtilities::FieldDefinitionsSharedPtr> Nektar::FieldUtils::Field::m_fielddef

Definition at line 90 of file Field.hpp.

LibUtilities::FieldMetaDataMap Nektar::FieldUtils::Field::m_fieldMetaDataMap

Definition at line 118 of file Field.hpp.

LibUtilities::PtsFieldSharedPtr Nektar::FieldUtils::Field::m_fieldPts

Definition at line 114 of file Field.hpp.

map<string, LibUtilities::FieldIOSharedPtr> Nektar::FieldUtils::Field::m_fld
private

Map to store FieldIO instances. Key is the reader type, value is the FieldIO object.

Definition at line 737 of file Field.hpp.

bool Nektar::FieldUtils::Field::m_fldToBnd

Definition at line 109 of file Field.hpp.

SpatialDomains::MeshGraphSharedPtr Nektar::FieldUtils::Field::m_graph

Definition at line 103 of file Field.hpp.

map<string, vector<string> > Nektar::FieldUtils::Field::m_inputfiles

Definition at line 105 of file Field.hpp.

MultiRegions::AssemblyMapCGSharedPtr Nektar::FieldUtils::Field::m_locToGlobalMap

Definition at line 116 of file Field.hpp.

LibUtilities::PtsIOSharedPtr Nektar::FieldUtils::Field::m_ptsIO

Definition at line 104 of file Field.hpp.

bool Nektar::FieldUtils::Field::m_requireBoundaryExpansion

Definition at line 97 of file Field.hpp.

LibUtilities::SessionReaderSharedPtr Nektar::FieldUtils::Field::m_session

Definition at line 102 of file Field.hpp.

bool Nektar::FieldUtils::Field::m_setUpEquiSpacedFields

Definition at line 112 of file Field.hpp.

bool Nektar::FieldUtils::Field::m_useFFT

Definition at line 99 of file Field.hpp.

bool Nektar::FieldUtils::Field::m_verbose

Definition at line 89 of file Field.hpp.

bool Nektar::FieldUtils::Field::m_writeBndFld

Definition at line 107 of file Field.hpp.

Generated on Wed Mar 8 2017 12:29:40 for Nektar++ by   doxygen 1.8.8