2
using namespace rheolef;
5
int main(int argc, char**argv) {
6
environment rheolef(argc, argv);
8
space Wh (gamma, argv[2]);
9
size_t d = gamma.dimension();
10
form m (Wh, Wh, "mass");
11
form a (Wh, Wh, "grad_grad");
13
field lh = riesz(Wh, f(d));
16
uh.set_u() = sa.solve (lh.u() - a.ub()*uh.b());