This becomes obsolete because the volume variables are now only computed
twice per dof in each newton iteration, a necessary update and one time
with a deflected solution to compute the derivatives correctly.