def METHOD_NAME(self):
    for ndim in [1, 2, 3]:
        for dtype in [np.float32, np.complex64]:
            with self.subTest(ndim=ndim, dtype=dtype):
                shape = [2, 20] + [1] * (ndim - 1)
                coord = np.random.random([10, ndim])
                input = np.random.random(shape).astype(dtype=dtype)
                output_cpu = interp.interpolate(
                    input, coord, kernel="kaiser_bessel"
                )
                output_gpu = interp.interpolate(
                    cp.array(input),
                    cp.array(coord),
                    kernel="kaiser_bessel",
                ).get()
                np.testing.assert_allclose(
                    output_cpu, output_gpu, atol=1e-5
                )