Function for sampling Dirichlet random variables (base-2 normalized via log-sum-exp for stability)
Useful links