Accelerated Gibbs sampling of normal distributions using matrix splittings and polynomials

Colin Fox, Albert Parker

Standard Gibbs sampling applied to a multivariate normal distribution with a specified precision matrix is equivalent in fundamental ways to the Gauss-Seidel iterative solution of…