46{
47 int i;
49
50 if (argc != 3)
51 {
52 fprintf(stderr,
53 "Usage: ./ExtractCriticalLayer meshfile fieldfile \n");
54 exit(1);
55 }
56
57
58
60 LibUtilities::SessionReader::CreateInstance(argc, argv);
61
62
63
64
66 SpatialDomains::MeshGraph::Read(vSession);
67
68
69
70
72
74 graphShPt);
75
76
77
78
79 string fieldfile(argv[argc - 1]);
80 vector<LibUtilities::FieldDefinitionsSharedPtr> fielddef;
81 vector<vector<NekDouble>> fielddata;
83
84
85
86
87 string streak_field("w");
88 for (unsigned int i = 0; i < fielddata.size(); ++i)
89 {
90 streak->ExtractDataToCoeffs(fielddef[i], fielddata[i], streak_field,
91 streak->UpdateCoeffs());
92 }
93
94
95 int npts;
96 vSession->LoadParameter("NumCriticalLayerPts", npts, 30);
99
101 vSession->LoadParameter("WidthOfLayers", trans, 0.1);
102
104
105 cout << "# x_c y_c" << endl;
106 for (i = 0; i < npts; ++i)
107 {
108 fprintf(stdout, "%12.10lf %12.10lf \n", x_c[i], y_c[i]);
109
110 }
111}
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
void Import(const std::string &infilename, std::vector< FieldDefinitionsSharedPtr > &fielddefs, std::vector< std::vector< NekDouble > > &fielddata, FieldMetaDataMap &fieldinfomap, const Array< OneD, int > &ElementIDs)
This function allows for data to be imported from an FLD file when a session and/or communicator is n...
std::shared_ptr< SessionReader > SessionReaderSharedPtr
std::shared_ptr< ExpList > ExpListSharedPtr
Shared pointer to an ExpList object.
std::shared_ptr< MeshGraph > MeshGraphSharedPtr