The Navier-Stokes equations on the driven cavity benchmark with the method of characteristics
int main (
int argc,
char**argv) {
if (argc < 2) {
cerr << "usage: " << argv[0] << " <geo> <Re> <err> <n_adapt>" << endl;
exit (1);
}
Float Re = (argc > 2) ? atof(argv[2]) : 100;
options.
err = (argc > 3) ? atof(argv[3]) : 1e-2;
size_t n_adapt = (argc > 4) ? atoi(argv[4]) : 5;
for (size_t i = 0; true; i++) {
size_t max_iter = 1000;
o.close();
if (i >= n_adapt) break;
omega =
adapt (ch, options);
o.open (omega.name(), "geo");
o << omega;
o.close();
Qh =
space (omega,
"P1");
}
}