• check if covariance matrix is pd in rmvnorm (by Fritz Leisch)