Use gibbs_step! in the run_lda script