Hermitian Matrices [message #10971] |
Thu, 29 January 1998 00:00 |
Evilio del Rio
Messages: 17 Registered: December 1997
|
Junior Member |
|
|
Hello,
I am trying to get the eigenvalues and eigenvectors of an
Hermitian matrix. An hermitian matrix has the property of being
its own transpose-conjugate:
A = CONJ(TRANSPOSE(A)) => A[i,j] = CONJ(A[j,i])
In fact, the real-symmetric matrices are just a special case of hermitian
and finding eigenvalues and vectors should be as easy as for real-symm.,
but all the routines in the IDL standard library ask you for real ones
(symmetric or not, depending on provedure).
Has anybody implemented an equivalent of EIGENQL (or even
ELMHES/HQR/EIGENVEC set)? If not, does anybody know the procedure to
reduce a Hermitian matrix to tridiagonal form (with complex arithmetics)?
Many thanks,
____________________________________________________________ ________
Evilio Jose del Rio Silvan Institut de Ciencies del Mar
E-mail: edelrio@icm.csic.es URL: http://www.ieec.fcr.es/~evilio/
"Anywhere you choose,/ Anyway, you're gonna lose"- Mike Oldfield
|
|
|