import random from numpy import * import time n=10000 #x=[1,1,1,1,1] x=random.random(n) k=matrix(arange(n)) l=matrix(arange(n)) F = exp( - 2*pi*1j*k.T*l/n ) start=time.time() fft.fft(x) end=time.time() print(end-start) start=time.time() F*reshape(x,(n,1)) end=time.time() print(end-start)