/*
Hauptteil des Codes.
Die Eingangsdaten werden in Blöcke von je 1024 Samples unterteilt und auf eine 1024-Punkte FFT gegeben.
*/
while(z<size-SAMPLE)
	{		

		for(k=0;k<SAMPLE;k++)
		{
			real[k] = daten[k+z];				// daten -> Eingangsdaten aus wav
		}

		zeros(imag,0,size);				


		// Forward FFT
		FFT(1,FFTEXPO,real,imag);					// imag -> Array für Imaginärteil, mit Nullen initialisiert

		// Kartesisch -> Polar
		abscplxarr(abscpl, real, imag);
		atancplxarr(atancpl, real, imag);

		// Polar -> Kartesisch
		rueckreal(real, abscpl, atancpl);
		rueckimag(imag, abscpl, atancpl);
		
		// Reverse FFT
		FFT(-1,FFTEXPO,real,imag);


		// Werte in Ergebnisarray
		for(i=0;i<SAMPLE;i++)		
			aus[i+z] = real[i];

		z+=SAMPLE;

	}





/////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// Kartesisch -> Polar
void abscplxarr(double *x, double *re, double *im)
{
	int i;
	for(i=0;i<SAMPLE;i++)	
	{
		x[i] = sqrt((re[i] * re[i]) + (im[i] * im[i]));
	}
}

void atancplxarr(double *x, double *re, double *im)
{
	int i = 0;
	for(i=0;i<SAMPLE;i++)
	{
		x[i] = atan2( im[i], re[i] );
	}
}

// Polar -> Kartesisch -------------------------------------------------------------------------
void rueckreal(double *x, double *re, double *im)
{
	int i;
	for(i=0;i<SAMPLE;i++)
	{
		x[i] = 0,5 * re[i] * cos(im[i]);
	}
}

void rueckimag(double *x, double *re, double *im)
{
	int i;
	for(i=0;i<SAMPLE;i++)
	{
		x[i] = 0,5 * re[i] * sin(im[i]);
	}
}


/*
Funktion für Nullen
*/
void zeros(double *x, int start, int end)
{
	double * ende = x+end;
	double * zaehl = x+start;
	while(zaehl < ende)
		*(zaehl++) = 0.0;
}