How work with updated Lagrangian meshes

Leaving this here as a future to-do to
1) Try to create an example/tutorial for solving something with an updated Lagrangian mesh description
2) In the process, identify and implement the (if any) new infrastructure that is required.

(I think it would require all requests from FEValues to also include the derivatives wrt. the cell coordinates)