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)