an efficient C++ finite element environment
|
|
Go to the documentation of this file.
31 int main (
int argc,
char**argv) {
34 cerr <<
"usage: " << argv[0] <<
" <geo> <Re> <err> <n_adapt>" << endl;
39 Float Re = (argc > 2) ? atof(argv[2]) : 100;
40 options.
err = (argc > 3) ? atof(argv[3]) : 1e-2;
41 size_t n_adapt = (argc > 4) ? atoi(argv[4]) : 5;
46 space Qh (omega,
"P1");
50 for (
size_t i = 0;
true; i++) {
51 size_t max_iter = 1000;
56 <<
catchmark(
"delta_t") << delta_t << endl
60 if (i >= n_adapt)
break;
62 omega =
adapt (ch, options);
63 o.open (omega.name(),
"geo");
67 Qh =
space (omega,
"P1");
int main(int argc, char **argv)
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
see the catchmark page for the full documentation
see the field page for the full documentation
field_basic< Float > field
see the field page for the full documentation
adapt_option: see the adapt page for the full documentation
see the space page for the full documentation
rheolef - reference manual
int navier_stokes_solve(Float Re, Float delta_t, field l0h, field &uh, field &ph, size_t &max_iter, Float &tol, odiststream *p_derr=0)
space_basic< Float > space
see the environment page for the full documentation
The Navier-Stokes equations – adaptive mesh criterion.
This file is part of Rheolef.
The Navier-Stokes equations with the method of characteristics – solver function.
odiststream derr(cerr)
see the diststream page for the full documentation
odiststream: see the diststream page for the full documentation
see the Float page for the full documentation
The driven cavity benchmark: 2D boundary conditions.
static space velocity_space(const geo &omega, string approx)
static field velocity_field(const space &Xh, Float alpha=1)
field navier_stokes_criterion(Float Re, const field &uh)
see the geo page for the full documentation