Avoid aliasing problems in sortedrandrange

JuliaDynamics · Apr 27, 2024 · ea8aa16 · ea8aa16
1 parent c26a496
commit ea8aa16
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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