Add a new unit test suite for the dunshire.matrices module.

[dunshire.git] / dunshire / matrices.py

diff --git a/dunshire/matrices.py b/dunshire/matrices.py
index 2d4bb17c98a9187e2cacaa0094ef9c0cfe3c4600..bcf83778d62752436f7894b5a6fbad4cc3e1012e 100644 (file)
--- a/dunshire/matrices.py
+++ b/dunshire/matrices.py
@@ -140,7 +140,10 @@ def eigenvalues(symmat):
 
     domain_dim = symmat.size[0]
     eigs = matrix(0, (domain_dim, 1), tc='d')
-    syevr(symmat, eigs)
+
+    # Create a copy of ``symmat`` here because ``syevr`` clobbers it.
+    dummy = matrix(symmat, symmat.size)
+    syevr(dummy, eigs)
     return list(eigs)
 
 
@@ -435,7 +438,7 @@ def condition_number(mat):
 
     """
     num_eigs = min(mat.size)
-    eigs = matrix(0, (num_eigs,1), tc='d')
+    eigs = matrix(0, (num_eigs, 1), tc='d')
     gesdd(mat, eigs)
 
     if len(eigs) > 0: