我們需要改變/重新實現在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):