The algorithm uses univariate normal and chi-square random values and then generates multivariate t values via a spectral decomposition.This function is used by rcompt.
References
Aitchison J. (1986). The statistical analysis of compositional data. Chapman & Hall.