Internal function used for generating truncated multivariate normal data
rtmvnorm (N, mu, sigma, A, b,..., rhoThr=NULL, maxSample=NULL)
Useful links