Reproject using biquadratic and bicubic interpolation returns Nans

I am working on reprojecting SIN images to an HPX projection. While using nearest-neighbor and bilinear interpolation generates images successfully, these methods aren't optimal for conserving flux. I prefer to use biquadratic or bicubic interpolation, but unfortunately, all output images turn out as NaNs when I apply these methods.

And there is no error message. There was similar issue before, but not sure if the issue was resolved.

Reproduce the SIN-to-HPX reprojection described in the issue using biquadratic and bicubic interpolation, then compare it with nearest-neighbor and bilinear results. Inspect the interpolation path and existing handling for NaN outputs; done means the higher-order methods produce valid output images or report a clear error.

Written by the indexing model from the issue text.