mcmc: use cholesky decomp rather than SVD for "matrix sqrt"