Allen,
Probably it needs to work in n-dimensions, like the existing
np.fft.fftshift function does, with an optional axis=tuple parameter. I
recall that fftshift is just an array indexing trick? It would be helpful
to see what's faster, two fftshifts and a edge padding or your
inter-padding. Prob
I'd like to get some feedback on my [pull request](
https://github.com/numpy/numpy/pull/7593).
This pull request adds a function `ifftpad` which pads a spectrum by
inserting zeros where the highest frequencies would be. This is necessary
because the padding that `ifft` does simply inserts zeros at
