Improve non-matching interpolation for non-affine grids.

[jorgensd]

For non-affine grids, the convex hull search often prove insufficient to get the correct colliding cell, for instance for non-matching interpolation (ref: https://fenicsproject.discourse.group/t/second-order-meshes-interpolation-error/19544/9?u=dokken).
This PR splits the non-affine (local) collision detection into two steps:

1. Check if point is within convex hull of cell.
2. If within convex hull, use the non-affine pullback to get point in reference space.

The logic for affine geometries is unchanged.
Requires #4126 as it showed up as a major performance bottleneck for this problem.