def test_g2k_multiple_equals_g2k(self):
     N = 10
     X = np.random.randn(N)
     me = GaussianQuadraticTest(self.grad_log_normal)
     G2K = me.g2k_multiple(X)
      
     for i in range(N):
         for j in range(N):
             g2k = me.g2k(X[i], X[j])
             assert_almost_equal(G2K[i, j], g2k)
    def test_g2k_multiple_equals_g2k(self):
        N = 10
        X = np.random.randn(N)
        me = GaussianQuadraticTest(self.grad_log_normal)
        G2K = me.g2k_multiple(X)

        for i in range(N):
            for j in range(N):
                g2k = me.g2k(X[i], X[j])
                assert_almost_equal(G2K[i, j], g2k)