File fix-test-sg03ad.patch of Package python-slycot

Overview Repositories Revisions Requests Users Attributes Meta

File fix-test-sg03ad.patch of Package python-slycot (Revision 3c0d3aa1821ee646523276e1ba0d7230)

Currently displaying revision 3c0d3aa1821ee646523276e1ba0d7230 , Show latest

diff --git a/slycot/tests/test_sg03ad.py b/slycot/tests/test_sg03ad.py
index 7b498d6..35d23e6 100644
--- a/slycot/tests/test_sg03ad.py
+++ b/slycot/tests/test_sg03ad.py
@@ -8,52 +8,73 @@ import unittest
 from slycot import synthesis
 import numpy as np
 
-from numpy.testing import assert_raises, assert_almost_equal
+from numpy.testing import assert_almost_equal
 
 # test cases from
 # http://www.qucosa.de/fileadmin/data/qucosa/documents/4168/data/b002.pdf
 
+
 class test_sg03ad(unittest.TestCase):
 
-    def test_sg03ad_a(self):
+    def test_sg03ad_1c(self):
         # Example 1
         n = 100
-        Xref = np.ones((n,n)) 
+        Xref = np.ones((n, n))
         U = np.tril(Xref)
         for t in range(0, 50, 10):
-            A = 2.0**(-t) - np.eye(n) + np.diag(range(1,n+1)) + U.T
-            E = np.eye(n) + 2**(-t)*U
+            A = (2**(-t) - 1) * np.eye(n) + np.diag(np.arange(1., n+1.)) + U.T
+            E = np.eye(n) + 2**(-t) * U
             Y = A.T.dot(Xref).dot(E) + E.T.dot(Xref).dot(A)
-            Q = np.zeros((n,n))
-            Z = np.zeros((n,n))
+            Q = np.zeros((n, n))
+            Z = np.zeros((n, n))
             A, E, Q, Z, X, scale, sep, ferr, alphar, alphai, beta = \
                 synthesis.sg03ad('C', 'B', 'N', 'N', 'L', n, A, E, Q, Z, Y)
-            assert_almost_equal(Xref, X)
-        
+            try:
+                assert_almost_equal(X, Xref, decimal=5)
+            except AssertionError as e:
+                print("Failed on t={}".format(t))
+                raise e
+
+    def test_sg03ad_1d(self):
+        # Example 1
+        n = 100
+        Xref = np.ones((n, n))
+        U = np.tril(Xref)
+        for t in range(0, 50, 10):
+            A = 2**(-t) * np.eye(n) + np.diag(np.arange(1., n+1.)) + U.T
+            E = np.eye(n) + 2**(-t) * U
+            Y = A.T.dot(Xref).dot(A) - E.T.dot(Xref).dot(E)
+            Q = np.zeros((n, n))
+            Z = np.zeros((n, n))
+            A, E, Q, Z, X, scale, sep, ferr, alphar, alphai, beta = \
+                synthesis.sg03ad('D', 'B', 'N', 'N', 'L', n, A, E, Q, Z, Y)
+            try:
+                assert_almost_equal(X, Xref, decimal=5)
+            except AssertionError as e:
+                print("Failed on t={}".format(t))
+                raise e
+
     def test_sg03ad_3(self):
         n = 3
         A = np.array([[3.0, 1.0, 1.0],
                       [1.0, 3.0, 0.0],
                       [1.0, 0.0, 2.0]])
         E = np.array([[1.0, 3.0, 0.0],
-                       [3.0, 2.0, 1.0],
-                       [1.0, 0.0, 1.0]])
+                      [3.0, 2.0, 1.0],
+                      [1.0, 0.0, 1.0]])
         Y = np.array([[64.0, 73.0, 28.0],
-                       [73.0, 70.0, 25.0],
-                       [28.0, 25.0, 18.0]])
+                      [73.0, 70.0, 25.0],
+                      [28.0, 25.0, 18.0]])
         Xref = np.array([[-2.0000, -1.0000, 0.0000],
-                          [-1.0000, -3.0000, -1.0000],
-                          [0.0000, -1.0000, -3.0000]])
-        Q = np.zeros((3,3))
-        Z = np.zeros((3,3))
+                         [-1.0000, -3.0000, -1.0000],
+                         [0.0000, -1.0000, -3.0000]])
+        Q = np.zeros((3, 3))
+        Z = np.zeros((3, 3))
         A, E, Q, Z, X, scale, sep, ferr, alphar, alphai, beta = \
             synthesis.sg03ad('C', 'B', 'N', 'N', 'L', n, A, E, Q, Z, -Y)
-        #print(A, E, Q, Z, X, scale, sep)
+        # print(A, E, Q, Z, X, scale, sep)
         assert_almost_equal(X, Xref)
 
-def suite():
-   return unittest.TestLoader().loadTestsFromTestCase(TestConvert)
-
 
 if __name__ == "__main__":
     unittest.main()

Places

File fix-test-sg03ad.patch of Package python-slycot (Revision 3c0d3aa1821ee646523276e1ba0d7230)

Places