From d1253e0bacf2770d1b077b0332d0b58142d96a50 Mon Sep 17 00:00:00 2001 From: Miles Lucas Date: Thu, 6 Aug 2020 01:49:01 -0500 Subject: [PATCH] Update matrixnormal.jl --- src/matrix/matrixnormal.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/matrix/matrixnormal.jl b/src/matrix/matrixnormal.jl index 75845b1e61..d34fdf6c79 100644 --- a/src/matrix/matrixnormal.jl +++ b/src/matrix/matrixnormal.jl @@ -125,7 +125,7 @@ function _rand!(rng::AbstractRNG, d::MatrixNormal, Y::AbstractMatrix) X = randn(rng, n, p) A = cholesky(d.U).L B = cholesky(d.V).U - Y .= d.M + A * X * B + Y .= d.M .+ A * X * B end # ----------------------------------------------------------------------------- @@ -154,7 +154,7 @@ end function _rand_params(::Type{MatrixNormal}, elty, n::Int, p::Int) M = randn(elty, n, p) - U = (X = 2rand(elty, n, n) .- 1; X * X') - V = (Y = 2rand(elty, p, p) .- 1; Y * Y') + U = (X = 2 .* rand(elty, n, n) .- 1; X * X') + V = (Y = 2 .* rand(elty, p, p) .- 1; Y * Y') return M, U, V end