Generating Random Numbers from Multivariate Normal Distributions with Non-Positive Definite Covariance Matrices in R

Here's a more detailed and formatted version of the response:

Fixing Performance Issues with RcppArmadillo: A Solution for pmvnorm_cpp Function