Some simple parallel FFT benchmarks in Julia

fftbench.jl

using BenchmarkTools

x = rand(Complex{Float64}, 2048,2048);
y = similar(x)

times = map(1:8) do i
    FFTW.set_num_threads(i)
    p = plan_fft(x, flags=FFTW.MEASURE)
    median((@benchmark A_mul_B!($y,$p,$x)).times)/1e6
end

println(times)