Bug in dgeevx_ ?

[Mikeekim]

Hi,

I calculated the eigenvalues from a matrix A and of her transposed matrix A^t.
In the first case i got complex eigenvalues and in the second real
eigenvalues.

But i would expect that both have the same eigenvalues, because

det(A-s*Id)=det((A-s*Id)^t)=det(A^t -s*Id).

The program is written in C.
I compiled it by
gcc -Wall -g -llapack -lblas -lgfortran qss_main.c read_data.c Hom_DGL_Solver.c
-o Eig.x

Did I make a mistake or lapack ?

[cottrell]

Did I make a mistake or lapack ?

Who knows? The former is more likely, the latter not impossible.
Show us your calls to dgeevx if you want an answer.

Allin Cottrell