我目前正在使用 FFTW2(西方最快的傅立叶变换)库,在用 Fortran 编写成功的例程后,我将转向 C。但是,当我尝试输入要转换的数据(即sin(x) 的值)。目前,我的代码是:
#include <stdio.h>
#include <fftw.h>
#include <math.h>
#include <complex.h>
void compute_fft(){
int i;
const int n[8];
fftw_complex in[8];
fftwnd_plan p;
in[0]->re = 0;
in[1]->re = (sqrt(2)/2);
in[2]->re = 1;
in[3]->re = (sqrt(2)/2);
in[4]->re = 0;
in[5]->re = -(sqrt(2)/2);
in[6]->re = -1;
in[7]->re = -(sqrt(2)/2);
for(i = 0; i < 8; i++){
(in[i])->im = 0;
}
p = fftwnd_create_plan(8, n, FFTW_FORWARD, FFTW_ESIMATE | FFTW_IN_PLACE);
fftwnd_one(p, &in[0], NULL);
fftwnd_destroy_plan(p);
printf("Sin\n");
for(i = 0; i < 8; i++){
printf("%d\n", n[i]);
}
}
我一直在用这个http://fftw.org/fftw2_doc/用于文档/教程目的的链接,目前我的错误是“a->a 的类型参数无效(具有 afftw_complexa)”,并且 -> 和 fftw_complex 两侧的“a”字符上方有克拉。我使用的模式与我在 Fortran 中编写此内容时的模式基本相同,而且我似乎正在遵循教程,实际上只是这里的这个作业,我搞乱了语法。根据记录,我使用的编译器是 nvcc,如果它有区别的话(我之前尝试过 gcc)。如果这里有人以前用过 C 中的 FFTW2,你能帮助纠正我的错误吗?谢谢!
最佳答案
这可能是因为你的数组“in”是 fftw_complex 的数组,所以不要使用 in[0]->re = 0 你应该使用它作为 in[0].re = 0 。除非 fftw_complex 再次出现类型定义为某个数组。
fftw_complex in[8];
in[0].re = 0;
in[1].re = (sqrt(2)/2);
in[2].re = 1;
in[3].re = (sqrt(2)/2);
in[4].re = 0;
in[5].re = -(sqrt(2)/2);
in[6].re = -1;
in[7].re = -(sqrt(2)/2);
关于c - C 中的 FFTW2 例程问题,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/5958766/