Python ProductOperator Examples

Example #1

File: test_matfuncs.py Project: YautongNg/gengis

 def test_product_operator(self):
     random.seed(1234)
     n = 5
     k = 2
     nsamples = 10
     for i in range(nsamples):
         A = np.random.randn(n, n)
         B = np.random.randn(n, n)
         C = np.random.randn(n, n)
         D = np.random.randn(n, k)
         op = ProductOperator(A, B, C)
         assert_allclose(op.matmat(D), A.dot(B).dot(C).dot(D))
         assert_allclose(op.T.matmat(D), (A.dot(B).dot(C)).T.dot(D))

Example #2

File: test_matfuncs.py Project: NeedAName/scipy

 def test_product_operator(self):
     random.seed(1234)
     n = 5
     k = 2
     nsamples = 10
     for i in range(nsamples):
         A = np.random.randn(n, n)
         B = np.random.randn(n, n)
         C = np.random.randn(n, n)
         D = np.random.randn(n, k)
         op = ProductOperator(A, B, C)
         assert_allclose(op.matmat(D), A.dot(B).dot(C).dot(D))
         assert_allclose(op.T.matmat(D), (A.dot(B).dot(C)).T.dot(D))