I would like to merge the branch interpolation_with_projection_into_element where the interpolation function will automatically include projection into element call when a point reference coordinate is computed provided that extrapolation function is not active.