mirror/kissfft
1
0
Fork 0
mirror of https://github.com/mborgerding/kissfft.git synced 2025-06-04 01:28:23 -04:00
Code Issues Packages Projects Releases Wiki Activity

about to try to split up k into two loops

Browse Source
This commit is contained in:
Mark Borgerding 2003-10-10 21:03:50 +00:00
parent 66b0646c9c
commit 93de2a9410
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
fft.py
View File

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