From 93de2a94100bcfde891acfc7cab564dc85bc9998 Mon Sep 17 00:00:00 2001
From: Mark Borgerding <mark@borgerding.net>
Date: Fri, 10 Oct 2003 21:03:50 +0000
Subject: [PATCH] about to try to split up k into two loops

---
 fft.py | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/fft.py b/fft.py
index 5fc4ae2..f5633a8 100644
--- a/fft.py
+++ b/fft.py
@@ -24,15 +24,16 @@ def fft(f):
 
     Fout=[0]*n
     for k in range(n):
-        val = 0
-        for q in range(p):
+        km = (k%m)
+        val = Fm[ km ]
+        for q in range(1,p):
             t = e ** ( j*2*pi*k*q/n )
-            val += Fm[ q*m + (k%m) ] * t
+            val += Fm[ q*m + km ] * t
         Fout[k] = val
 
     return Fout
 
-def test(f=range(256),ntimes=10):
+def test(f=range(1024),ntimes=10):
     import time
     t0 = time.time()
     for i in range(ntimes):