I have a 30 by 30 covariance matrix that is symmetric (passes with both the isSymmetric function as well as is.symmetric.matrix from matrixcalc library). However, I can’t sample from it using rmvnorm(), as I get the following error: