我们需要改变/重新实现在GSL标准DFT实现,这是 int
FUNCTION(gsl_dft_complex,transform) (const BASE data[],
const size_t stride, const size_t n,
BASE result[],
const gsl_fft_directio
我试图实现3D DFT,但我遇到了一些麻烦。我相信我应该做的只是连续3次一维DFT,每个方向一个。假设一维DFT是正确的,你可以看到什么是错,此代码: def dft3d(self, real3d, img3d, nx, ny, nz, dir):
#Transform depth
for i in range(nx):
for j in range(ny):