Add _gibbs_sample() and unit test.