From ea8aa16f88b6d4662838c4d8e754967b1b20b2f2 Mon Sep 17 00:00:00 2001 From: Tortar <68152031+Tortar@users.noreply.github.com> Date: Sat, 27 Apr 2024 03:50:26 +0200 Subject: [PATCH] Avoid aliasing problems in sortedrandrange --- src/SortedSamplingMulti.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/SortedSamplingMulti.jl b/src/SortedSamplingMulti.jl index 2ccc8c1..d7814db 100644 --- a/src/SortedSamplingMulti.jl +++ b/src/SortedSamplingMulti.jl @@ -62,7 +62,7 @@ end function sortedrandrange(rng, range, n) exp_rands = randexp(rng, n) - sorted_rands = accumulate!(+, exp_rands, exp_rands) + sorted_rands = cumsum(exp_rands) a, b = range.start, range.stop range_size = b-a+1 cum_step = (sorted_rands[end] + randexp(rng)) / range_size