Fix custom interpolant

This MR fixes the new interpolation introduced in !917 (merged). I also added some testing this time, so hopefully, it will be more robust.

I've verified that the marginalized likelihood example runs now.