/usr/share/doc/rheolef-doc/examples/elasticity_criterion.icc is in rheolef-doc 6.5-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 | field elasticity_criterion (Float lambda, const field& uh) {
string grad_approx = (uh.get_approx() == "P2") ? "P1d" : "P0";
space Xh (uh.get_geo(), grad_approx);
if (grad_approx == "P0") return interpolate (Xh, norm(uh));
space T0h (uh.get_geo(), grad_approx);
size_t d = uh.get_geo().dimension();
tensor I = tensor::eye (d);
return interpolate (T0h, sqrt(2*norm2(D(uh)) + lambda*sqr(div(uh))));
}
|