FFT klappt nicht richtig

[drops]

Hallo,
ich hänge hier schon etwas länger dran, aber irgendwie bkomme ich nicht ganz zu rande. Ich habe aus dem Numerical Recipes in C eine Funktion für eine FFT rausgesucht, aber irgendwie bekomme ich die nicht zum laufen. Vielleicht könnt ihr mir ja weiterhelfen. Vielleicht habt ihr ja eine funktionierende. Sie muß ja nicht schnell sein, sie muß nur klappen.

void Wave::FFT(void)
{
	waveIn[0]=0;
	float data[NUMPTS * 2+2];

	for (int Schleife = 1; Schleife <= NUMPTS; Schleife++)
	{
		data[2 * Schleife-1] = waveIn[Schleife];
		data[2*Schleife] = 0;
	}
	unsigned long nn = NUMPTS;
	int isign = 1;
	/*Replaces data[1..2*nn] by its discrete Fourier transform, if isign is input as 1; or replaces
	data[1..2*nn] by nn times its inverse discrete Fourier transform, if isign is input as -1.
	data is a complex array of length nn or, equivalently, a real array of length 2*nn. nn MUST
	be an integer power of 2 (this is not checked for!).*/

	unsigned long n,mmax,m,j,istep,i;
	double wtemp,wr,wpr,wpi,wi,theta;// Double precision for the trigonometric recurrences. 
	float tempr,tempi;
	n=nn << 1;
	j=1;
	for (i=1;i<n;i+=2)//This is the bit-reversal section of the routine.
	{
		if (j > i)
		{
			SWAP(data[j],data[i]);// Exchange the two complex numbers.
			SWAP(data[j+1],data[i+1]);
		}
		m=nn;
		while (m >= 2 && j > m)
		{
			j -= m;
			m >>= 1;
		}
		j += m;
	}
	//Here begins the Danielson-Lanczos section of the routine.
	mmax=2;
	while (n > mmax)//Outer loop executed log2 nn times.
	{
		istep=mmax << 1;
		theta=isign*(6.28318530717959/mmax);// Initialize the trigonometric recurrence.
		wtemp=sin(0.5*theta);
		wpr = -2.0*wtemp*wtemp;
		wpi=sin(theta);
		wr=1.0;
		wi=0.0;
		for (m=1;m<mmax;m+=2) // Here are the two nested inner loops.
		{
			for (i=m;i<=n;i+=istep)
			{
				j=i+mmax; //This is the Danielson-Lanczos formula:
				tempr=wr*data[j]-wi*data[j+1];
				tempi=wr*data[j+1]+wi*data[j];
				data[j]=data[i]-tempr;
				data[j+1]=data[i+1]-tempi;
				data[i] += tempr;
				data[i+1] += tempi;
			}
			wr=(wtemp=wr)*wpr-wi*wpi+wr; //Trigonometric recurrence.
			wi=wi*wpr+wtemp*wpi+wi;
		}
		mmax=istep;
	}
}

[Marc++us]

hier habe ich eine drin... allerdings noch in Pascal:

http://www.baeckmann.de/DFT.pdf