Vector rmvt(const Vector &mu, const SpdMatrix &Sigma, double nu){ return rmvt_mt(GlobalRng::rng, mu, Sigma, nu);}
Vec rmvt(const Vec &mu, const Spd &Sigma, double nu){ return rmvt_mt(GlobalRng::rng, mu, Sigma, nu);}