Multivariate Normal distribution sampler (Rcpp version), derived using the eigendecomposition of the covariance matrix Sigma. The function utilizes the arma random normal generator
.mvrnorm_arma(n, Mu, Xmat, eigen = TRUE)
an n
sample from a multivariate Normal distribution
sample size
mean vector. Will set the dimension
covariance matrix, of same dimension as Mu
(and square matrix).
No sanity check is performed to validate that the matrix is symmetric, so use at own risk