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");
12
field b = m*field(Wh,1);
13
field lh = riesz(Wh, f(d));
14
csr<Float> A = {{ a.uu(), b.u()},
16
vec<Float> B = { lh.u(), 0 };
18
vec<Float> U = sa.solve (B);
20
uh.set_u() = U [range(0,uh.u().size())];