arma::rowvec rmvlx(arma::mat R, arma::rowvec Z){ Rcpp::RNGScope scope; int m = R.n_rows; double e = Rf_rexp(1); arma::rowvec ans(m); ans = sqrt(e) * rmvnormx(R, Z); return ans; }
Type rexp(Type rate) { return Rf_rexp(asDouble(rate)); }