Esempio n. 1
0
 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;
 }
Esempio n. 2
0
 inline double calc() const {
   return normal_logp(x_,mu_,tau_);
 }