Implementation of Latent Dirichlet Allocation where the topic distrubutions are multivariate Gaussians. This implementation uses Cholesky decomposition of covariance matrices to perform an order of magnitude times faster posterior inference using Gibbs sampling.