rdir function

Random variate generation for Dirichlet distribution on SdSd