void HPRS::draw_mu_given_zero_mean_sufficient_statistics() { Vector mu = rmvn_suf_mt(rng(), xtx_, xtu_); model_->data_parent_model()->set_mu(mu); }
Vector rmvn_suf(const SpdMatrix &Ivar, const Vector &IvarMu) { return rmvn_suf_mt(GlobalRng::rng, Ivar, IvarMu); }