rheolef  6.6
navier_stokes_criterion.cc
Go to the documentation of this file.
1 // zcat square-4.field.gz | ./navier_stokes_criterion | geo -
2 #include "rheolef.h"
3 using namespace rheolef;
4 using namespace std;
5 #include "navier_stokes_criterion.icc"
6 int main (int argc, char**argv) {
7  environment rheolef (argc, argv);
8  adapt_option_type options;
9  options.err = (argc > 1) ? atof(argv[1]) : 1e-4;
10  options.hmin = (argc > 2) ? atof(argv[2]) : 0.004;
11  options.additional = "-RelError";
12  options.hmax = 0.1;
13  Float Re;
14  field uh;
15  din >> catchmark("Re") >> Re
16  >> catchmark("u") >> uh;
17  field ch = navier_stokes_criterion (Re, uh);
18  dout << adapt (ch, options);
19 }
field - piecewise polynomial finite element field
Definition: field.h:225
std::string additional
Definition: adapt.h:112
STL namespace.
irheostream, orheostream - large data streams
Definition: compiler.h:7
double Float
Definition: compiler.h:177
geo_basic< T, M > adapt(const field_basic< T, M > &uh, const adapt_option_type &opts)
Definition: adapt.cc:123
idiststream din(cin)
catchmark - iostream manipulator
Definition: catchmark.h:30
adapt - mesh adaptation
Definition: adapt.h:95
odiststream dout(cout)
Definition: diststream.h:317