comparison data/fft/FFTapi.cpp @ 227:e919a2b97c2a

* Fix in-house FFT provision

 #include "FFTapi.h"
 
 #ifndef HAVE_FFTW3F
 
 #include <cmath>
+#include <iostream>
 
 void
 fft(unsigned int n, bool inverse, double *ri, double *ii, double *ro, double *io)
 {
     if (!ri || !ro || !io) return;

 	    bits = i;
 	    break;
 	}
     }
 
-    static unsigned int tableSize = 0;
-    static int *table = 0;
-
-    if (tableSize != n) {
-
-	delete[] table;
-
-	table = new int[n];
-
-	for (i = 0; i < n; ++i) {
+    int *table = new int[n];
+
+    for (i = 0; i < n; ++i) {
 	
-	    m = i;
-
-	    for (j = k = 0; j < bits; ++j) {
-		k = (k << 1) | (m & 1);
-		m >>= 1;
-	    }
-
-	    table[i] = k;
-	}
-
-	tableSize = n;
+        m = i;
+        
+        for (j = k = 0; j < bits; ++j) {
+            k = (k << 1) | (m & 1);
+            m >>= 1;
+        }
+        
+        table[i] = k;
     }
 
     if (ii) {
 	for (i = 0; i < n; ++i) {
 	    ro[table[i]] = ri[i];

 	}
 
 	blockEnd = blockSize;
     }
 
+/* fftw doesn't normalise, so nor will we
+
     if (inverse) {
 
 	double denom = (double)n;
 
 	for (i = 0; i < n; i++) {
 	    ro[i] /= denom;
 	    io[i] /= denom;
 	}
     }
+*/
+    delete[] table;
 }
 
 struct fftf_plan_ {
     int size;
     int inverse;

 }
 
 void
 fftf_execute(const fftf_plan p)
 {
-    //!!! doesn't work yet for inverse transforms
-
     float *real = p->real;
     fftf_complex *cplx = p->cplx;
     int n = p->size;
-    int inv = p->inverse;
+    int forward = !p->inverse;
 
     double *ri = new double[n];
     double *ro = new double[n];
     double *io = new double[n];
 
     double *ii = 0;
-    if (inv) ii = new double[n];
-
-    if (!inv) {
+    if (!forward) ii = new double[n];
+
+    if (forward) {
         for (int i = 0; i < n; ++i) {
             ri[i] = real[i];
         }
     } else {
         for (int i = 0; i < n/2+1; ++i) {

                 ii[n-i] = -ii[i];
             }
         }
     }
 
-    fft(n, inv, ri, ii, ro, io);
-
-    if (!inv) {
+    fft(n, !forward, ri, ii, ro, io);
+
+    if (forward) {
         for (int i = 0; i < n/2+1; ++i) {
             cplx[i][0] = ro[i];
             cplx[i][1] = io[i];
         }
     } else {

     }
 
     delete[] ri;
     delete[] ro;
     delete[] io;
-    delete[] ii;
+    if (ii) delete[] ii;
 }
 
 #endif