benchmark.jl updated. run it to benchmark the proposed algorithm against the current Julia implementation