def test_g1k_multiple_equals_g1k(self): N = 10 X = np.random.randn(N) me = GaussianQuadraticTest(self.grad_log_normal) G1K = me.g1k_multiple(X) for i in range(N): for j in range(N): g1k = me.g1k(X[i], X[j]) assert_almost_equal(G1K[i, j], g1k)
def test_g1k_multiple_equals_g1k(self): N = 10 X = np.random.randn(N) me = GaussianQuadraticTest(self.grad_log_normal) G1K = me.g1k_multiple(X) for i in range(N): for j in range(N): g1k = me.g1k(X[i], X[j]) assert_almost_equal(G1K[i, j], g1k)