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)