Interpolation seems to work, but still problem with Python reference implementation