random.multivariate_normal() now uses cholesky decomposition