def test_log_det_exact_toy_large_scikits(self):
n = 1e6
d = abs(randn(n))
Q = spdiags(d, 0, n, n)
self.assertAlmostEqual(OzonePosterior.log_det_scikits(Q), sum(log(d)),
delta=1e-5)
def test_log_det_ozone_scikits_exact(self):
o = OzonePosterior()
kappa = 2 ** (-13.1)
Q = o.create_Q_matrix(kappa)
self.assertAlmostEqual(OzonePosterior.log_det_scikits(Q),
2.317769370813052e+06,
delta=0.001)
