[SciPy-dev] lobpcg: removed symeig dependence

Nathan Bell wnbell at gmail.com
Tue Apr 8 12:28:30 EDT 2008


On Tue, Apr 8, 2008 at 8:58 AM, Robert Cimrman <cimrman3 at ntc.zcu.cz> wrote:
> 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.]

Works for me!

Thanks for tackling this issue.  I was worried that lobpcg wouldn't be
ready for the next release.

Robert, could we switch the order of the A and BlockVectorX
parameters?  Having A first matches the convention used by other
iterative solvers in SciPy.

Also, we should provide a short description of what the method
computes in terms of A, B, lambda, etc.  A small example would be nice
too (since the interface is a little different than, e.g. ARPACK).

Thanks again!

-- 
Nathan Bell wnbell at gmail.com
http://graphics.cs.uiuc.edu/~wnbell/



More information about the SciPy-Dev mailing list