[Scipy-tickets] [SciPy] #1317: Arpack wrappers should tolerate odd # of complex eigenvalues
SciPy Trac
scipy-tickets@scipy....
Sat Oct 23 09:51:42 CDT 2010
#1317: Arpack wrappers should tolerate odd # of complex eigenvalues
---------------------------------+------------------------------------------
Reporter: pv | Owner: wnbell
Type: defect | Status: new
Priority: normal | Milestone: 0.9.0
Component: scipy.sparse.linalg | Version: 0.8.0
Keywords: |
---------------------------------+------------------------------------------
Remaining item from #1313:
III) This had been the hardest issue to find the reason for since I've
never done anything with ARPACK. What I've got was an InvalidIndexError?
here: arpack.py:264 z[:,i] = zr[:,i] + 1.0j * zr[:,i+1] called
(indirectly) from test_nonsymmetric_modes. The reason is the following:
I've got three eigenvalues in d with imaginary part nonzero. Why do I get
it and you don't: If you don't give arpack a starting vector, it'll
generate a random one (THIS SHOULD BE DOCUMENTED!) and this influences -
which eigenvalues are returned by ARPACK - in which order they're
returned. I've written two pure fortran test programs resembling python's
ARPACK usage if you want to verify this (also happens on my Ubuntu Lucid
Lynx x86_64 box). One should think about fixing that piece of code there
in arpack.py (taking into account the possibility of an odd number of
imaginary eigenvalues). But the testsuite will still fail due to the fact
that the set of eigenvalues returned is not fixed. The solution is to fix
them by giving starting vectors in the testsuite (see
04_arpack_test_always_set_v0.diff)
--
Ticket URL: <http://projects.scipy.org/scipy/ticket/1317>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.
More information about the Scipy-tickets
mailing list