10
10
form a (Xh, Xh, "grad_grad");
11
11
form ab (Xh, Xh, "mass", omega["boundary"]);
13
field lh = riesz(Xh, f(d)) + riesz(Xh, "boundary", g(d));
13
field lh = riesz(Xh, f(d)) + riesz(Xh, g(d), "boundary");
16
16
uh.set_u() = sa.solve (lh.u() - a.ub()*uh.b());