From ca7f1b4863603c483550030dcbadcb748a7b2702 Mon Sep 17 00:00:00 2001
From: hackbard <hackbard>
Date: Fri, 8 Oct 2004 10:06:09 +0000
Subject: [PATCH] improved "clever way" method - seems to work now .. strange
 it didnt before!

---
 fourier/fourier.c | 14 +++++++-------
 1 file changed, 7 insertions(+), 7 deletions(-)

diff --git a/fourier/fourier.c b/fourier/fourier.c
index 9f810ec..89f42c2 100644
--- a/fourier/fourier.c
+++ b/fourier/fourier.c
@@ -220,22 +220,22 @@ int fourier_dft_2d(t_fourier *fourier) {
   }
   //  dft on index 2 of 'index 1 transformed data'
   for(x=0;x<X;x++) {
-    off_f=0;
+    off_f=x;
     for(v=0;v<Y;v++) {
       //dprintf(fourier->outfd,"[fourier] (u=%d,v=%d)\n",v,x);
-      off_r=0;
+      off_r=x;
       for(y=0;y<Y;y++) {
         arg=-2.0*M_PI*v*y/Y;
-        fourier->ftdata[off_f+x].r+=(cos(arg)*data[off_r+x].r-sin(arg)*data[off_r+x].i);
-        fourier->ftdata[off_f+x].r+=(sin(arg)*data[off_r+x].r+cos(arg)*data[off_r+x].i);
+        fourier->ftdata[off_f].r+=(cos(arg)*data[off_r].r-sin(arg)*data[off_r].i);
+        fourier->ftdata[off_f].r+=(sin(arg)*data[off_r].r+cos(arg)*data[off_r].i);
         off_r+=X;
       }
-      fourier->ftdata[off_f+x].r/=Y;
-      fourier->ftdata[off_f+x].i/=Y;
+      fourier->ftdata[off_f].r/=Y;
+      fourier->ftdata[off_f].i/=Y;
       off_f+=X;
     }
   }
-
+ 
   free(data);
 
   return F_SUCCESS;
-- 
2.20.1