Yang Gao's implementation of Jimmy Fould's SCVB0 algorithm, burn-in for each doc is modified after Michael Hankin's vectorized method