From 4e2567ecd3643cd80af11f1e939f21bf6898b4f4 Mon Sep 17 00:00:00 2001
From: simon
Date: Fri, 24 Oct 2008 15:16:42 +0000
Subject: added some comment and clean some code

---
 src/signalProcessing/ifft/fft842.c | 52 +++++++-------------------------------
 1 file changed, 9 insertions(+), 43 deletions(-)

(limited to 'src/signalProcessing/ifft/fft842.c')

diff --git a/src/signalProcessing/ifft/fft842.c b/src/signalProcessing/ifft/fft842.c
index 27682ae6..2742463f 100644
--- a/src/signalProcessing/ifft/fft842.c
+++ b/src/signalProcessing/ifft/fft842.c
@@ -46,28 +46,20 @@ void fft842 (doubleComplex* b, int size , int in)
   fn = nthpo;
 
 
-  if(in==FORWARD)
-    /* take conjugate  */
-    for(i=0;i< size ;i++) {
-	b[i] = DoubleComplex ( zreals( b[i]) , - zimags (b[i]));
-
-     /* b[i].im *= -1.0;*/
-    }
-
   if(in==INVERSE)
     /*scramble inputs*/
-        
+
     for(i=0,j=size/2;j<size;i++,j++)
       {
-	temp = DoubleComplex ( zreals ( b[j] ) , zimags( b[j] ));
-	b[j] = DoubleComplex ( zreals ( b[i] ) , zimags( b[i] ));
-	b[i] = DoubleComplex ( zreals ( temp ) , zimags( temp ));
+        temp = DoubleComplex ( zreals ( b[j] ) , zimags( b[j] ));
+        b[j] = DoubleComplex ( zreals ( b[i] ) , zimags( b[i] ));
+        b[i] = DoubleComplex ( zreals ( temp ) , zimags( temp ));
 
-	/*
-	r = b[j].re; fi = b[j].im;
-	b[j].re = b[i].re; b[j].im = b[i].im;
-	b[i].re = r; b[i].im = fi;
-	*/
+        /*
+        r = b[j].re; fi = b[j].im;
+        b[j].re = b[i].re; b[j].im = b[i].im;
+        b[i].re = r; b[i].im = fi;
+        */
       }
 
   n8pow = n2pow/3;
@@ -158,42 +150,16 @@ void fft842 (doubleComplex* b, int size , int in)
       ij++;
     }
 
-  if(in==FORWARD)
-    {    /* take conjugates & unscramble outputs */
-        for(i=0,j=size/2;j<size;i++,j++)
-        {
-        temp = DoubleComplex ( zreals ( b[j] ) ,- zimags( b[j] ));
-        b[j] = DoubleComplex ( zreals ( b[i] ) ,- zimags( b[i] ));
-        b[i] = DoubleComplex ( zreals ( temp ) , zimags( temp ));
 
 
-        /*  r = b[j].re; fi = b[j].im;
-          b[j].re = b[i].re; b[j].im = -b[i].im;
-          b[i].re = r; b[i].im = -fi;
-        */
-        }
-    }
-
   if(in==INVERSE) /* scale outputs */
     {
         for(i=0;i<nthpo;i++)
           {
-            printf ( "plopl\n");  
             b[i] =  DoubleComplex ( zreals( b[i] )/fn , zimags(b[i])/fn);
             fn *= -1 ;
           }
     }
 
-/*
-       for ( i = 0 ; i < size /2 ; i++)
-        {
-          temp =  b[i] ;
-          b[i] =  b[i+(size/2)];
-          b[i+(size/2)]= temp ;
-
-
 
-        }
-
-*/
 }
-- 
cgit