rheolef  6.5
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 }
20