fix: use rand! in rand for MvGaussianMeanScalePrecision

ReactiveBayes · Oct 22, 2024 · fc1103d · fc1103d
1 parent 98f2343
commit fc1103d
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/src/distributions/normal_family/mv_normal_mean_scale_precision.jl b/src/distributions/normal_family/mv_normal_mean_scale_precision.jl
@@ -183,8 +183,9 @@ function BayesBase.prod(
 end
 
 function BayesBase.rand(rng::AbstractRNG, dist::MvGaussianMeanScalePrecision{T}) where {T}
-    μ, γ = mean(dist), scale(dist)
-    return μ + 1 / γ .* randn(rng, T, length(μ))
+    μ, γ = params(dist)
+    d = length(μ)
+    return rand!(rng, dist, Vector{T}(undef, d))
 end
 
 function BayesBase.rand(rng::AbstractRNG, dist::MvGaussianMeanScalePrecision{T}, size::Int64) where {T}