Nektar++
Functions
utilities/ExtractCriticalLayer.cpp File Reference
#include <cstdio>
#include <cstdlib>
#include <MultiRegions/ExpList.h>
#include <MultiRegions/ExpList2D.h>
#include "ExtractCriticalLayerFunctions.h"
Include dependency graph for utilities/ExtractCriticalLayer.cpp:

Go to the source code of this file.

Functions

int main (int argc, char *argv[])
 

Function Documentation

int main ( int  argc,
char *  argv[] 
)

Definition at line 44 of file utilities/ExtractCriticalLayer.cpp.

References Computestreakpositions(), Nektar::LibUtilities::Import(), and npts.

45 {
46  int i;
47  NekDouble cr = 0;
48 
49  if(argc !=3)
50  {
51  fprintf(stderr,"Usage: ./ExtractCriticalLayer meshfile fieldfile \n");
52  exit(1);
53  }
54 
55  //------------------------------------------------------------
56  // Create Session file.
58  = LibUtilities::SessionReader::CreateInstance(argc, argv);
59  //-----------------------------------------------------------
60 
61  //-------------------------------------------------------------
62  // Read in mesh from input file
63  SpatialDomains::MeshGraphSharedPtr graphShPt = SpatialDomains::MeshGraph::Read(vSession);
64  //------------------------------------------------------------
65 
66  //-------------------------------------------------------------
67  // Define Streak Expansion
69 
71  ::AllocateSharedPtr(vSession,graphShPt);
72  //---------------------------------------------------------------
73 
74  //----------------------------------------------
75  // Import field file.
76  string fieldfile(argv[argc-1]);
77  vector<LibUtilities::FieldDefinitionsSharedPtr> fielddef;
78  vector<vector<NekDouble> > fielddata;
79  LibUtilities::Import(fieldfile,fielddef,fielddata);
80  //----------------------------------------------
81 
82  //----------------------------------------------
83  // Copy data from field file
84  string streak_field("w");
85  for(unsigned int i = 0; i < fielddata.size(); ++i)
86  {
87  streak->ExtractDataToCoeffs(fielddef [i],
88  fielddata[i],
89  streak_field,
90  streak->UpdateCoeffs());
91  }
92  //----------------------------------------------
93 
94  int npts;
95  vSession->LoadParameter("NumCriticalLayerPts",npts,30);
96  Array<OneD, NekDouble> x_c(npts);
97  Array<OneD, NekDouble> y_c(npts);
98 
99  NekDouble trans;
100  vSession->LoadParameter("WidthOfLayers",trans,0.1);
101 
102  Computestreakpositions(streak,x_c, y_c,cr,trans);
103 
104  cout << "# x_c y_c" << endl;
105  for(i = 0; i < npts; ++i)
106  {
107  fprintf(stdout,"%12.10lf %12.10lf \n",x_c[i],y_c[i]);
108  //cout << x_c[i] << " " << y_c[i] << endl;
109  }
110 
111 }
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
Definition: MeshPartition.h:50
void Import(const std::string &infilename, std::vector< FieldDefinitionsSharedPtr > &fielddefs, std::vector< std::vector< NekDouble > > &fielddata, FieldMetaDataMap &fieldinfomap, const Array< OneD, int > ElementiDs)
Imports an FLD file.
Definition: FieldIO.cpp:106
void Computestreakpositions(MultiRegions::ExpListSharedPtr &streak, Array< OneD, NekDouble > &xc, Array< OneD, NekDouble > &yc, NekDouble cr, NekDouble trans)
boost::shared_ptr< ExpList > ExpListSharedPtr
Shared pointer to an ExpList object.
Definition: ExpList.h:1340
static std::string npts
Definition: InputFld.cpp:43
double NekDouble
boost::shared_ptr< MeshGraph > MeshGraphSharedPtr
Definition: MeshGraph.h:432