Fix ranks when calling ksl code

tt/eigb/setup.py

@@ -25,6 +25,7 @@ def configuration(parent_package='', top_path=None):
         'primme',
         sources=[join(PRIMME_DIR, '*.c')] + [join(PRIMME_DIR, '*.f')],
         include_dirs=['.'],
+        #extra_compile_args=[' -O0 -g -fcheck=all,no-array-temps -fbacktrace -ffpe-trap=zero,overflow,underflow'],
     )
     config.add_extension(
         'tt_eigb',
@@ -39,6 +40,7 @@ def configuration(parent_package='', top_path=None):
             'mytt',
             'print_lib',
         ],
+        #extra_compile_args=[' -O0 -g -fcheck=all,no-array-temps -fbacktrace -ffpe-trap=zero,overflow,underflow'],
     )
 
     return config

tt/ksl/ksl.py

@@ -60,6 +60,9 @@ def ksl(A, y0, tau, verb=1, scheme='symm', space=8, rmax=2000):
         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
@@ -167,6 +170,9 @@ 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

tt/ksl/setup.py

@@ -37,6 +37,8 @@ def configuration(parent_package='', top_path=None):
     config.add_library(
         'expokit',
         sources=expokit_src,
+        f2py_options=[' --debug-capi']
+        #extra_compile_args=[' -O0 -g -fcheck=all,no-array-temps -fbacktrace -ffpe-trap=zero,overflow,underflow'],
     )
     config.add_extension(
         'dyn_tt',
@@ -51,6 +53,8 @@ def configuration(parent_package='', top_path=None):
             'expokit',
             'mytt',
         ],
+        # extra_compile_args=[' -O0 -g -fcheck=all,no-array-temps -fbacktrace -ffpe-trap=zero,overflow,underflow'],
+        # f2py_options=[' --debug-capi'],
     )
 
     return config