added option to estimate pi by variational Bayes approach (need testing)
add function to compute KL divergence between dirichlet distns