Add elasticity diamond solver skeleton

This defines a simple skeleton to ease implementation of elasticity
solver.

It is linked to the pugs' language. Called for instance by
`elasticity(mesh, bc_list, lambda, mu, f);`