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.
1.26.4
3.10.14 | packaged by conda-forge | (main, Mar 20 2024, 12:45:18) [GCC 12.3.0]
[{'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'}]
Would like to understand whether there is a bug in the calculation of the Nyquist frequency component of 2D FFT.
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:
Error message:
Python and NumPy Versions:
Runtime Environment:
Context for the issue:
Would like to understand whether there is a bug in the calculation of the Nyquist frequency component of 2D FFT.