[SciPy-dev] [SciPy-user] scipy.linalg.eig() returns transposed eigenvector matrix

Fernando Perez Fernando.Perez at colorado.edu
Mon Nov 14 12:56:18 CST 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