double logp() const {
double ans(0);
const uint mu_size = mu_.size();
const uint tau_size = tau_.size();
for(size_t i = 0; i < MCMCStochastic<double,ArmaT>::size(); i++) {
ans += normal_logp(MCMCStochastic<double,ArmaT>::value_[i], mu_[i % mu_size], tau_[i % tau_size]);
}
return ans;
}
inline double calc() const {
return normal_logp(x_,mu_,tau_);
}