wl1ball_runtime_benchmark.jl

#w1ball_runtime_benchmark.jl

using BenchmarkTools, Random, Distributions
BenchmarkTools.DEFAULT_PARAMETERS.gcsample = true
BenchmarkTools.DEFAULT_PARAMETERS.samples = 100
BenchmarkTools.DEFAULT_PARAMETERS.seconds = 20

include("weighted_simplex_and_ball/wl1ball.jl")

println("You are using ", nthreads(), " threads for parallel computing")
println("Warning: following experiments are for 80 threads!")

function wsortscan_test()
    Random.seed!(12345); res = @benchmark wball_sortscan_s($(rand(Normal(0,1), 10^8)), $(rand(10^8)), 1)
    println("serial: ", median(res))
    for i = 1:10
        Random.seed!(12345); res = @benchmark wball_sortscan_p($(rand(Normal(0,1), 10^8)), $(rand(10^8)), 1, )
        println(i*8, " : ", median(res))
    end
end

function wmichelot_test()
    Random.seed!(12345); res = @benchmark wball_michelot_s($(rand(Normal(0,1), 10^8)), $(rand(10^8)), 1)
    println("serial: ", median(res))
    for i = 1:10
        Random.seed!(12345); res = @benchmark wball_michelot_p($(rand(Normal(0,1), 10^8)), $(rand(10^8)), 1, $(i*8))
        println(i*8, " : ", median(res))
    end
end

function wcondat_test()
    Random.seed!(12345); res = @benchmark wball_condat_s($(rand(Normal(0,1), 10^8)), $(rand(10^8)), 1)
    println("serial: ", median(res))
    for i = 1:10
        Random.seed!(12345); res = @benchmark wball_condat_p($(rand(Normal(0,1), 10^8)), $(rand(10^8)), 1, $(i*8))
        println(i*8, " : ", median(res))
    end
end


wsortscan_test()
wmichelot_test()
wcondat_test()