We describe a Markov chain Monte Carlo method to approximately simulate high dimensional Gaussian vectors whose covariance matrix is easy to compute. The standard Monte Carlo method is based on the Cholesky decomposition, which is performed in cubic time and quadratic storage cost in the...