Yang Gao's implementation of Jimmy Fould's SCVB0 algorithm, coded in C++ with Eigen matrix package, burn-in for each doc is modified after Michael Hankin's vectorized method