benchmark.jl

# export JULIA_NUM_THREADS=6
# ./julia -p 6
# BLAS.set_num_threads(6)
flps=Array{Float64}(0)
iterations=10
for n in 1:iterations
  one_run=peakflops(3000,parallel=true) / 1e9;
  push!(flps, one_run); 
  println("run no. $(n) (of $(iterations)): ~ $(round(one_run, 2)) Gflops");
end; println("minimum: $(round(minimum(flps),2)), maximum: $(round(maximum(flps),2)), mean: $(round(mean(flps),2)), median: $(round(median(flps),2)) with std: $(round(std(flps),2))")