[SciPy-dev] sparse eigensolvers
Thu Oct 11 08:55:17 CDT 2007
> There is an ARPACK wrapper in the scipy sandbox. Right now the Python
> interface only handles the standard eigenvalue problem Ax=wx
> (symmetric and nonsymmetric) but the f2py wrapper should be able to
> handle the generalized eigenproblem and all of the various shifted modes.
> There are some design decisions to be made regarding the interface
> and some more code to be written to accomplish that.
> There is a description of the current state and some ideas for
> an interface at
Oh, I misread you are interested in eigensolvers, sorry about that. :)
Yes, I don't know about any opensource systematic approach, in the
same spirit as petsc4py. And I need this as well, for my master
thesis. There is slepc4py:
which is free, but slepc itself is not free. Arpack is unfortunately
too archaic (doesn't perform well on my problems), but I have a very
good experience with blzpack and pysparse, also primme is
However, I'd be interested in some kind of unified interface in
Python, the same way like one can use petsc4py for Ax=b problems.
There are many libraries for eigensolvers and linear solvers though,
so I am not sure if they all should be in SciPy? Maybe it's better to
create a separater project for a unified python interface to
eigensolvers, in a similar way like petsc4py is a unified python
interface for almost all opensource linear solvers available?
More information about the Scipy-dev