BUG: 2D fft doesn't reconstruct function at Nyquist rate #27413
Description
Describe the issue:
I have a unit test here for a 2D FFT that only passes above the Nyquist rate, when it should pass at the Nyquist rate.
The interpolation code is here. I'm certain there are no implementation issues. Could someone please review this?
I am creating the issue here rather than with JAX and scipy because the issue is common to all three and I think those libraries just use numpy's implementation.
Reproduce the code example:
See link above.Error message:
See link above.Python and NumPy Versions:
1.26.4
3.10.14 | packaged by conda-forge | (main, Mar 20 2024, 12:45:18) [GCC 12.3.0]
Runtime Environment:
[{'numpy_version': '1.26.4',
'python': '3.10.14 | packaged by conda-forge | (main, Mar 20 2024, 12:45:18) '
'[GCC 12.3.0]',
'uname': uname_result(system='Linux', node='fedora', release='6.10.9-200.fc40.x86_64', version='#1 SMP PREEMPT_DYNAMIC Sun Sep 8 17:23:55 UTC 2024', machine='x86_64')},
{'simd_extensions': {'baseline': ['SSE', 'SSE2', 'SSE3'],
'found': ['SSSE3',
'SSE41',
'POPCNT',
'SSE42',
'AVX',
'F16C',
'FMA3',
'AVX2'],
'not_found': ['AVX512F',
'AVX512CD',
'AVX512_KNL',
'AVX512_KNM',
'AVX512_SKX',
'AVX512_CLX',
'AVX512_CNL',
'AVX512_ICL',
'AVX512_SPR']}},
{'architecture': 'Haswell',
'filepath': '/home/kaya/miniconda3/envs/desc-env/lib/libopenblasp-r0.3.27.so',
'internal_api': 'openblas',
'num_threads': 12,
'prefix': 'libopenblas',
'threading_layer': 'pthreads',
'user_api': 'blas',
'version': '0.3.27'}]
Context for the issue:
Would like to understand whether there is a bug in the calculation of the Nyquist frequency component of 2D FFT.