1
field elasticity_solve(const space& Xh, Float lambda) {
2
size_t d = Xh.get_geo().dimension();
5
field lh = riesz (Xh, f);
6
form m (Xh, Xh, "mass");
7
form a1 (Xh, Xh, "div_div");
8
form a2 (Xh, Xh, "2D_D");
9
form a = lambda*a1 + a2;
12
uh.set_u() = sa.solve (lh.u() - a.ub()*uh.b());