不用cuda的cufft库实现fft的代码

[sabin541]

我想实现fft，但不希望用cufft库，也就是自己写的意思，下面是一篇论文中提供的核函数。不过怎么调用核函数的部分没有找到，有谁知道吗？或者和下面完全不同的思路的也可以，高手指教！
static __global__ void FFT( Complex * DataIn, Complex *DataOut, const unsigned int N)
{
extern __shared__ Complex sdata[];
const unsigned int tid_in_block = threadIdx.x;
if(tid_in_block<N)
{
sdata[tid_in_block] = DataIn[tid_in_block];
sdata[tid_in_block + N/2] = DataIn[tid_in_block + N/2];
__syncthreads();
if(tid_in_block < N/2){
int p,q;
Complex Wn,Xp,XqWn;
float stage = 0.0;
for(int Ns=1; Ns<N; Ns=Ns*2)
{
p = tid_in_block/Ns*Ns*2 + tid_in_block%Ns;
q = p + Ns;
Wn = tex2D(texRef,tid_in_block,stage++);
XqWn = ComplexMul(sdata[q], Wn);
Xp = sdata[p];
sdata[p] = ComplexAdd(Xp, XqWn);
sdata[q] = ComplexSub(Xp, XqWn);
__syncthreads();
}
DataOut[p] = sdata[p];
DataOut[q] = sdata[q];
}
}
}

[sabin541]

自己回复一下，上面上论文上的ｃｕｄａ代码到底也没搞明白，当点数是２的幂次方的时候，我用快速傅里叶变化来做的，然后进行ｃｕｄａ优化，当点数不是２的幂次方的时候，我用ＤＦＴ的变形ｍｉｘｄｆｔ算法，然后进行ｃｕｄａ优化，具体代码就不贴了。总之最后速度还可以。

[NickelCao]

他那个论文的纹理存储器来实现Wn的按表查找，以及 p q的值的 那段 应该是用来定位这个thread要计算的2个点。这2部分搞不大懂。你是在他基础上改的么，能提供一点代码参考下么？谢谢！