Rapid Posterior Exploration in Bayesian Non-negative Matrix Factorization