def test_log_likelihood_scikits_exact(self):
tau = 2 ** (-11.35)
kappa = 2 ** (-13.1)
o = OzonePosterior(prior=None, logdet_method="scikits",
solve_method="scikits")
self.assertAlmostEqual(o.log_likelihood(tau, kappa),
- 9.336375798558606e+05,
delta=0.001)
def precompute_likelihood_estimates(self, tau, kappa):
logging.debug("Entering")
estimates = asarray([OzonePosterior.log_likelihood(self, tau, kappa) for _ in range(self.num_estimates)])
logging.debug("Leaving")
return estimates
