请问：急！！！！如何调用傅立叶反变化函数IFFT()!!!!!!!!!!!!!!!!!!

[orangeliu]

void IFFT(complex<double>*FD,complex<double>*TD,int r)
{
long count;
int i;
complex<double>*X;
count=1<<r;
X=new complex<double>[count];
memcpy(X,FD,sizeof(complex<double>)*count);
for(i=0;i<count;i++)
{
X[i]=complex<double>(X[i].real(),-X[i].imag());
}
FFT(X,TD,r);
for(i=0;i<count;i++)
{
TD[i]=complex<double>(TD[i].real()/count,-TD[i].imag()/count);
}
delete X;
}

[LionEagle]

to orangeliu () :
FD是输入，TD是输出

delete X;不对，应该用 delete []X;

memcpy(X,FD,sizeof(complex<double>)*count);没必要，直接这样：

for( i=0; i<count; i++ )
{
X[i] = complex<double>( FD[i].real(),-FD[i].imag() );
}

to zlma(zlma):
这个算法没什么问题，也不慢，借用FFT而已

[zlma]

这种算法也能用？
怎么样也要搞个快速算法吧。。