Rheolef  7.1
an efficient C++ finite element environment
stokes_taylor_dg.cc

The Stokes problem for the Taylor benchmark by the discontinuous Galerkin method

#include "rheolef.h"
using namespace rheolef;
using namespace std;
#include "taylor.h"
int main(int argc, char**argv) {
environment rheolef (argc, argv);
geo omega (argv[1]);
space Xh (omega, argv[2], "vector");
space Qh (omega, argv[2]);
form a, b, c, mp;
field lh, kh;
stokes_dirichlet_dg (Xh, Qh, a, b, c, mp, lh, kh);
field uh (Xh, 0), ph (Qh, 0);
problem_mixed stokes (a, b, c);
stokes.set_metric (mp);
stokes.solve (lh, kh, uh, ph);
dout << catchmark("u") << uh
<< catchmark("p") << ph;
}
form
see the form page for the full documentation
rheolef::catchmark
see the catchmark page for the full documentation
Definition: catchmark.h:67
taylor.h
The Taylor benchmark – right-hand-side and boundary condition.
mkgeo_ball.b
int b
Definition: mkgeo_ball.sh:152
field
see the field page for the full documentation
problem_mixed
see the problem_mixed page for the full documentation
main
int main(int argc, char **argv)
Definition: stokes_taylor_dg.cc:30
space
see the space page for the full documentation
rheolef.h
rheolef - reference manual
stokes_dirichlet_dg.icc
The Stokes problem with Dirichlet boundary condition by the discontinuous Galerkin method – solver fu...
a
Definition: diffusion_isotropic.h:25
rheolef::environment
see the environment page for the full documentation
Definition: environment.h:115
lh
field lh(Float epsilon, Float t, const test &v)
Definition: burgers_diffusion_operators.icc:25
stokes_dirichlet_dg
void stokes_dirichlet_dg(const space &Xh, const space &Qh, form &a, form &b, form &c, form &mp, field &lh, field &kh, integrate_option iopt=integrate_option())
Definition: stokes_dirichlet_dg.icc:25
rheolef
This file is part of Rheolef.
Definition: compiler_eigen.h:37
mkgeo_ball.a
int a
Definition: mkgeo_ball.sh:151
rheolef::dout
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:430
rheolef::std
Definition: vec_expr_v2.h:402
geo
see the geo page for the full documentation
mkgeo_ball.c
int c
Definition: mkgeo_ball.sh:153