LAPACK Archives

[Lapack] A fault in lapack?



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).

I added my program with my data file.
The program is written in C.

For calling the programm from a console if the data file 100eV513.txt and the
programm in the same folder:

./Eig.x  100eV513.txt non  // for the calculation of the eigvalues and
./Eig.x  100eV513.txt trans // for calulating the eigvalues of the transposed
matrix



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 ?

Greetings
Mike
-------------- next part --------------
A non-text attachment was scrubbed...
Name: A.tex
Type: text/x-tex
Size: 25279 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/A-0001.bin
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Atrans.tex
Type: text/x-tex
Size: 25775 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/Atrans-0001.bin
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Eig.x
Type: application/octet-stream
Size: 79868 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/Eig-0001.obj
 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: 100eV513.txt
Url: 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/100eV513-0001.txt
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Hom_DGL_Solver.c
Type: text/x-c++src
Size: 6102 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/Hom_DGL_Solver-0002.bin
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Hom_DGL_Solver.h
Type: text/x-chdr
Size: 487 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/Hom_DGL_Solver-0003.bin
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: complex_funktion.h
Type: text/x-chdr
Size: 579 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/complex_funktion-0001.bin
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: lapack.h
Type: text/x-chdr
Size: 1016 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/lapack-0001.bin
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: matrix.h
Type: text/x-chdr
Size: 147 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/matrix-0001.bin
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: qss_main.c
Type: text/x-c++src
Size: 1165 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/qss_main-0001.bin
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: read_data.c
Type: text/x-c++src
Size: 2442 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/read_data-0002.bin
 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: read_data.h
Type: text/x-chdr
Size: 121 bytes
Desc: not available
Url : 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100728/3871cc3c/read_data-0003.bin
 

<Prev in Thread] Current Thread [Next in Thread>
  • [Lapack] A fault in lapack?, mmuecke at smail.uni-koeln.de <=


For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or