interp_points = (n_samples, n_sample_dims) - CONFIRMED
data = always N-dimensional, where N = cube.ndim - CONFIRMED
result = not entirely sure

@pp-mo can you confirm the last one for me while I continue through the review?  Iv spent ~15min on it but its not immediately obvious.