Efficient Simulation of High Dimensional Gaussian Vectors
Efficient Simulation of High Dimensional Gaussian Vectors
We describe a Markov chain Monte Carlo method to approximately simulate a centered d-dimensional Gaussian vector X with given covariance matrix.The standard Monte Carlo method is based on the Cholesky decomposition, which takes cubic time and has quadratic storage cost in d.In contrast, the storage cost of our algorithm is …