[SciPy-Dev] Strange results by eigs
Fri Oct 21 12:50:13 CDT 2011
Ah, thanks Nils. I missed that detail when I wrote the wrapper.
It's all a bit confusing, because each wrapper calls one of six arpack
routines, depending on the type of the input.
Looking at the code, I think everything should work as expected with
hermitian input for M, but the docs need to be updated, as Pauli suggested.
The bit on M in both eigs() and eigsh() should specify that if the
matrix A is complex hermitian, then M may be complex hermitian as well.
It still must be real symmetric in any other case.
Nils Wagner wrote:
> On Fri, 21 Oct 2011 10:12:38 -0700
> Jacob VanderPlas <email@example.com>
>>>> The docstring of eigs is misleading in that context.
>>>> M must represent a real symmetric matrix.
>>> Yep, it should say hermitian instead of real sym. here.
>> Actually, the piece of the documentation quoted here has
>> to do with the second matrix supplied for the general
>> eigenvalue problem. This matrix (M) must be real
>> symmetric, as desribed in the ARPACK function
> The Hermitian case is here
> SciPy-Dev mailing list
More information about the SciPy-Dev