Rheolef  7.2
an efficient C++ finite element environment
navier_stokes_criterion.cc
Go to the documentation of this file.
1 // zcat square-4.field.gz | ./navier_stokes_criterion | geo -
22 #include "rheolef.h"
23 using namespace rheolef;
24 using namespace std;
26 int main (int argc, char**argv) {
27  environment rheolef (argc, argv);
28  adapt_option options;
29  options.err = (argc > 1) ? atof(argv[1]) : 1e-4;
30  options.hmin = (argc > 2) ? atof(argv[2]) : 0.004;
31  options.additional = "-RelError";
32  options.hmax = 0.1;
33  Float Re;
34  field uh;
35  din >> catchmark("Re") >> Re
36  >> catchmark("u") >> uh;
37  field ch = navier_stokes_criterion (Re, uh);
38  dout << adapt (ch, options);
39 }
see the Float page for the full documentation
see the field page for the full documentation
see the catchmark page for the full documentation
Definition: catchmark.h:67
see the environment page for the full documentation
Definition: environment.h:121
idiststream din(cin)
see the diststream page for the full documentation
Definition: diststream.h:464
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:467
This file is part of Rheolef.
geo_basic< T, M > adapt(const field_basic< T, M > &uh, const adapt_option &opts)
adapt(uh,opts): see the adapt page for the full documentation
Definition: adapt.cc:179
int main(int argc, char **argv)
The Navier-Stokes equations – adaptive mesh criterion.
field navier_stokes_criterion(Float Re, const field &uh)
rheolef - reference manual
adapt_option: see the adapt page for the full documentation
Definition: adapt.h:147
std::string additional
Definition: adapt.h:164