rheolef  6.5
robin2.cc
Go to the documentation of this file.
1 #include "rheolef.h"
2 using namespace rheolef;
3 using namespace std;
4 #include "sinusprod_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  trial u (Xh); test v (Xh);
11  form a = integrate (dot(grad(u),grad(v))) + integrate ("boundary", u*v);
12  field lh = integrate (f(d)*v) + integrate ("boundary", g(d)*v);
13  field uh (Xh);
14  solver sa (a.uu());
15  uh.set_u() = sa.solve (lh.u() - a.ub()*uh.b());
16  dout << uh;
17 }
18