Add a single gibbs step function to LDA and use it in the tests