r9488 - /1.3/test_suite/unit_tests/_maths_fns/test_kronecker_product.py -- September 09, 2009 - 16:04

Author: bugman
Date: Wed Sep  9 16:04:40 2009
New Revision: 9488

URL: http://svn.gna.org/viewcvs/relax?rev=9488&view=rev
Log:
Added a unit test to check if all the 
maths_fns.kronecker_product.transpose_*() functions revert.

Calling the function twice should return the matrix to the original state.


Modified:
    1.3/test_suite/unit_tests/_maths_fns/test_kronecker_product.py

Modified: 1.3/test_suite/unit_tests/_maths_fns/test_kronecker_product.py
URL: 
http://svn.gna.org/viewcvs/relax/1.3/test_suite/unit_tests/_maths_fns/test_kronecker_product.py?rev=9488&r1=9487&r2=9488&view=diff
==============================================================================
--- 1.3/test_suite/unit_tests/_maths_fns/test_kronecker_product.py (original)
+++ 1.3/test_suite/unit_tests/_maths_fns/test_kronecker_product.py Wed Sep  9 
16:04:40 2009
@@ -282,3 +282,25 @@
             for j in range(9):
                 print(("i = %2s, j = %2s, daeg[i,j] = %s" % (i, j, daegT[i, 
j])))
                 self.assertEqual(self.daeg[i, j], daegT[i, j])
+
+
+    def test_transpose_reversions(self):
+        """Check that the transposes revert back to the original matrix."""
+
+        # Make a copy of the frame order matrix.
+        daeg_orig = self.to_numpy(self.daeg_str)
+
+        # List of transpose functions.
+        Tij = [transpose_12, transpose_13, transpose_14, transpose_23, 
transpose_24, transpose_34]
+
+        # Check the transpose reversions.
+        for transpose in Tij:
+            # Transpose twice.
+            transpose(self.daeg)
+            transpose(self.daeg)
+
+            # Check.
+            for i in range(9):
+                for j in range(9):
+                    print(("i = %2s, j = %2s, daeg[i,j] = %s" % (i, j, 
daeg_orig[i, j])))
+                    self.assertEqual(self.daeg[i, j], daeg_orig[i, j])