Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
elasticity_criterion.icc
Go to the documentation of this file.
1
26 string grad_approx = "P" + to_string(uh.get_space().degree()-1) + "d";
27 space Xh (uh.get_geo(), grad_approx);
28 if (grad_approx == "P0d") return lazy_interpolate (Xh, norm(uh));
29 space T0h (uh.get_geo(), grad_approx);
30 size_t d = uh.get_geo().dimension();
31 tensor I = tensor::eye (d);
32 return lazy_interpolate (T0h, sqrt(2*norm2(D(uh)) + lambda*sqr(div(uh))));
33}
see the Float page for the full documentation
see the field page for the full documentation
see the space page for the full documentation
see the tensor page for the full documentation
field elasticity_criterion(Float lambda, const field &uh)