FFTW在某些循环中返回不同的结果。

对于具有相同大小和类型的输入，FFTW 应该始终返回相同的输出。如果在某些循环中观察到不同的结果，则可能是由于以下问题之一：

1. 数据类型不匹配。确保输入数组与 FFTW 算法所需的数据类型匹配。例如，如果使用 double 类型的 FFTW 算法，则输入数组也应为 double 类型。
2. 内存对齐问题。确保输入数组在内存中正确对齐。使用 FFTW 提供的“ aligned_malloc”函数来分配内存，以确保内存对齐。
3. FFTW 模式未设置正确。在使用 FFTW 进行多线程计算时，确保正确设置 FFTW 模式，以便 FFTW 可以正确识别要处理的数据和线程数。例如，如果要通过 FFTW 处理 1D 输入数组并启用 4 个线程，则应使用“FFTW_ESTIMATE | FFTW_PRESERVE_INPUT | FFTW_DESTROY_INPUT”标志设置 FFTW 模式。

示例代码：

double *in, *out; const int N = 1024;

in = (double*) fftw_aligned_malloc(sizeof(double)N); out = (double) fftw_aligned_malloc(sizeof(double)*N);

// ... fill in[] with data ...

fftw_complex out_complex; out_complex = (fftw_complex) fftw_aligned_malloc(sizeof(fftw_complex)*N);

fftw_plan plan = fftw_plan_dft_r2c_1d(N, in, out_complex, FFTW_ESTIMATE); fftw_execute(plan);

// ... do something with out_complex ...

fftw_destroy_plan(plan); fftw_aligned_free(in); fftw_aligned_free(out); fftw_aligned_free(out_complex);