search

jakevdp / nfft

Lightweight non-uniform Fast Fourier Transform in Python
MIT License
201 stars 29 forks source link

Divide by N to make round trip work (nfft and nfft_adjoint) inverses. #18

Closed cottrell closed 2 years ago

jakevdp commented 2 years ago

Hi - thanks for looking into this. I don't think we'll want to make this change, because the current implementation does precisely what is documented in the package README, which follow the convention in the referenced paper: https://www-user.tu-chemnitz.de/~potts/paper/nfft3.pdf

There may be good arguments for using a different normalization, but I'd prefer to stick to those conventions here.

cottrell commented 2 years ago

Ah if that is convention then leave it. Wish I had noticed that. adjoint is not an inverse in this case.