Dear lapack team,
I tried to compute the left-eigenvector of the matrix A=[[0,1]; [-1,0]]
using the function zgeev and I got some results that seem wrong.
I get the eigenvalue eig=I and the associated left-eigenvector vect =
[-0.7071 * I, 0.7071].
However, this is not a left eigenvector:
vect * A = [-0.701, -0.701]
eig * vect = [0.701, 0.701*I]
The correct left-eigenvector for I should be [0.701, -0.701*I]
Also the other eigenvalue, eigenvector pair that is returned by zgeev
seems wrong (eig= -I, vect=[0.707*I, 0.707])
I bumped into this problem using the lapacke c interface, but I
reproduced it also with fortran (I used the example code provided online
I also tried the examples using the svn version of the library, but I
get the same results.
I attach a tar.gz archive with the source code and the test case to
reproduce the issue.
As I am new to lapack I could have missed some details and I could have
misused the library.
PS: I also tried to use the forum, but I cannot post the bug report
since the website always mark it as spam.
-------------- next part --------------
A non-text attachment was scrubbed...
Size: 17795 bytes
Desc: not available