Sample random values from the multivariate normal distribution
sample_mvnorm(mu, sigma)
mu
: mean vector
sigma
: covariance matrix
Samples multivariate normal variables by multiplying univariate random normal variables by the cholesky decomposition of the covariance matrix.
If mu is length 1 then just uses rnorm instead.
Useful links