Add StrainInterpolateXi routine

Currently only implemented for finite elasticity,
this calculates the Green-Lagrange strain tensor at a given xi location.