rand!(a, Multinomial(100, p)) StackOverflowError bug

[cossio]

using Random, Distributions
p = rand(10)
p = p / sum(p)
a = zeros(10)
rand!(a, Multinomial(100, p))

Results in an error,

ERROR: StackOverflowError:
Stacktrace:
 [1] rand!(rng::TaskLocalRNG, A::Vector{Float64}, X::Distributions.MultinomialSampler{Float64}) (repeats 79984 times)
   @ Random ~/.julia/juliaup/julia-1.11.2+0.x64.linux.gnu/share/julia/stdlib/v1.11/Random/src/Random.jl:268

[devmotion]

The arguments of rand! are incorrectly ordered in your example: https://juliastats.org/Distributions.jl/stable/multivariate/#Random.rand!-Tuple{AbstractRNG,%20MultivariateDistribution,%20AbstractArray}

julia> rand!(Multinomial(100, p), a)
10-element Vector{Float64}:
  3.0
 13.0
 16.0
  5.0
 12.0
 17.0
  7.0
 24.0
  1.0
  2.0

[cossio]

I missed that. Isn't this order different from the one at Base?
https://docs.julialang.org/en/v1/stdlib/Random/#Random.rand!

[devmotion]

Yes: #1316