BMNSamples(Theta, numSamples, burnIn, skip)numSamples times p where p is the number of rows of Theta.
BMNSamples generates numSamples by using Gibbs sampling. When using Gibbs sampling, it is necessary to discard the initial samples, which is controlled by the parameter burnIn. In order for the drawn samples to be independent, samples in-between also have to be discarded, which is controlled by skip.
BMNPseudo, BMNExact