[SciPy-user] Sparse eigenvalues/eigenvectors problem

Nils Wagner nwagner@iam.uni-stuttgart...
Fri Jun 13 11:54:13 CDT 2008


On Fri, 13 Jun 2008 15:57:54 +0100
  "James Philbin" <philbinj@gmail.com> wrote:
> Hmm... Has something changed in eigen_symmetric? Now, 
>when I run
> eigen_symmetric(S,k=10,which='LA') i get the following 
>eigenvalues
> returned: [ 1.  1.  1.  1.  1.  1.  1.  1.  1.  1.], 
>which is good!
> Curious as to what was changed.
> 
> Thanks,
> James
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user

Is your matrix positive definite ?
Did you try symeig ?
http://mdp-toolkit.sourceforge.net/symeig.html

from scipy import *
from pylab import spy, show
from symeig import symeig
A = io.mmread('S.mtx.gz')
w,v = symeig(A.todense(),range=(1,25))
spy(A.todense())
show()

Nils


More information about the SciPy-user mailing list