Matrix variate Dirichlet distribution