rheolef  6.5
dirichlet2-ball.cc
Go to the documentation of this file.
1 #include "rheolef.h"
2 using namespace rheolef;
3 using namespace std;
4 #include "sinusrad_laplace.icc"
5 int main(int argc, char**argv) {
6  environment rheolef(argc, argv);
7  geo omega (argv[1]);
8  size_t d = omega.dimension();
9  space Xh (omega, argv[2]);
10  Xh.block ("boundary");
11  trial u (Xh); test v (Xh);
12  form a = integrate (dot(grad(u),grad(v)));
13  field lh = integrate (f(d)*v);
14  field uh (Xh);
15  uh ["boundary"] = 0;
16  solver sa (a.uu());
17  uh.set_u() = sa.solve (lh.u() - a.ub()*uh.b());
18  dout << uh;
19 }
20