2
using namespace rheolef;
4
#include "cosinusrad_laplace.icc"
5
int main(int argc, char**argv) {
6
environment rheolef(argc, argv);
8
size_t d = omega.dimension();
9
space Xh (omega, argv[2]);
10
Xh.block ("boundary");
11
form a (Xh, Xh, "grad_grad");
12
field lh = riesz (Xh, f(d));
14
space Wh (omega["boundary"], argv[2]);
15
uh ["boundary"] = interpolate(Wh, g(d));
17
uh.set_u() = sa.solve (lh.u() - a.ub()*uh.b());