//定点
/*Source: Numerical Recipe In C
* Modified by IIT Madras
* The following code computes 32point FFT
* The bit reversal positions are precomputed and stored in the header files.
* The Twiddle Factor is decomposed into following form:
* exp(ia + ib) = exp(ia)*exp(ib) = exp(ia) + Z*exp(ia)
* Where Z = exp(ib) - 1 = -(1 - cos(b)) + isin(b) = -2*sin^2(b/2) + isin(b)
* Z is precomputed and stored in header file.
* */
#include "32fft.h"
#include "32fftvalues.h"
using namespace std;
void
