[SciPy-dev] [SciPy-user] scipy.linalg.eig() returns transposed eigenvector matrix
Fernando Perez
Fernando.Perez at colorado.edu
Mon Nov 14 13:56:18 EST 2005
Travis Oliphant wrote:
> Pearu Peterson wrote:
>>This is a matter of definition. scipy.linalg.eig and
>>scipy.basic.linalg.eig return correct results according to their
>>documentation. Just scipy.linalg.eig assumes that eigenvectors are
>>returned column-wise, i.e.
>>
>> a * vr[:,i] = w[i] * vr[:,i]
>>
>>holds. While scipy.basic.linalg.eig, that is copied from Numeric, assumes
>>that eigenvectors are returned row-wise, i.e
>>
>>
>> a * vr[i] = w[i] * vr[i]
>>
>>holds.
>>
>>
>
> Thanks for the clarification, Pearu. I'm glad things are actually
> working as advertised.
If I may suggest, I think these two should be unified, though. It will be
seriously disconcerting for new users to find that
scipy.basic.linalg.eig = transpose(scipy.linalg.eig)
Just my 2c.
Cheers,
f
More information about the SciPy-Dev
mailing list