Assume that we have a well-defined real-valued functionf[x,y,z] whose (inverse) Fourier transform is impossible to solve for symbolically and numerically using NIntegrate. I would like to sample the function and use (inverse) DFT (i.e. (Inverse)Fourier) and then use Interpolation to get a function I can further use. Doing this for 1d functions is a bit easier, but I don’t really understand how to do it in more dimensions. Any help is appreciated, thanks!