Hi, so it is done, thanks again, Pearu, for your hint. could someone try on the fresh SVN version (4114): $ cd <your_path_to_scipy_sources>/scipy/sparse/linalg/eigen/lobpcg $ python lobpcg.py the output should end with: [ 4. 5. 6.] r.