Calculation of volume change (Linear Elasticity)

Use a projection, as explained by @kamensky in: Project(x) works, interpolate(x) does not - #2 by kamensky