The algorithm is straightforward, for each vector, independent gamma values are generated and then divided by their total sum.
References
Ng Kai Wang, Guo-Liang Tian and Man-Lai Tang (2011). Dirichlet and related distributions: Theory, methods and applications. John Wiley & Sons.
Aitchison J. (1986). The statistical analysis of compositional data. Chapman & Hall.