Rheolef  7.1
an efficient C++ finite element environment
stokes_taylor_dg.cc
Go to the documentation of this file.
1 #include "rheolef.h"
26 using namespace rheolef;
27 using namespace std;
28 #include "taylor.h"
29 #include "stokes_dirichlet_dg.icc"
30 int main(int argc, char**argv) {
31  environment rheolef (argc, argv);
32  geo omega (argv[1]);
33  space Xh (omega, argv[2], "vector");
34  space Qh (omega, argv[2]);
35  form a, b, c, mp;
36  field lh, kh;
37  stokes_dirichlet_dg (Xh, Qh, a, b, c, mp, lh, kh);
38  field uh (Xh, 0), ph (Qh, 0);
39  problem_mixed stokes (a, b, c);
40  stokes.set_metric (mp);
41  stokes.solve (lh, kh, uh, ph);
42  dout << catchmark("u") << uh
43  << catchmark("p") << ph;
44 }
field lh(Float epsilon, Float t, const test &v)
see the field page for the full documentation
see the form page for the full documentation
see the geo page for the full documentation
see the problem_mixed 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:115
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:430
see the space page for the full documentation
This file is part of Rheolef.
rheolef - reference manual
The Stokes problem with Dirichlet boundary condition by the discontinuous Galerkin method – solver fu...
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())
int main(int argc, char **argv)
The Taylor benchmark – right-hand-side and boundary condition.