Merge pull request #96 from oseledets/hotfix/ksl

Revert fix on ranks growth in KSL
oseledets · Dec 14, 2022 · a50d5e0 · a50d5e0
2 parents 21509a3 + 9b36c44
commit a50d5e0
Showing 1 changed file with 0 additions and 6 deletions.
diff --git a/tt/ksl/ksl.py b/tt/ksl/ksl.py
@@ -62,9 +62,6 @@ def ksl(A, y0, tau, verb=1, scheme='symm', space=8, rmax=2000, use_normest=1):
         0.0
     """
 
-    y0 = y0.round(1e-14)    # This will fix ranks
-                            # to be no more than maximal reasonable.
-                            # Fortran part doesn't  handle excessive ranks
     ry = y0.r.copy()
     if scheme is 'symm':
         tp = 2
@@ -178,9 +175,6 @@ def diag_ksl(A, y0, tau, verb=1, scheme='symm', space=8, rmax=2000):
         >>> print tt.dot(y1, y0) / (y1.norm() * y0.norm()) - 1 #Eigenvectors should not change
         0.0
     """
-    y0 = y0.round(1e-14)    # This will fix ranks
-                            # to be no more than maximal reasonable.
-                            # Fortran part doesn't  handle excessive ranks
     ry = y0.r.copy()
     if scheme is 'symm':
         tp = 2