Hi,
I'm new to Lapack libraries and Intel MKL.
I have to solve an eigensystem problem, with 3x3 complex matrixes, and I don't know if I'm calling Lapack routines in the correct order. I use these routines from C. Here is the list:
First I call
zgebal
with job ='P'.
Then I call
zgehdr.
Then i copy the Hessenberg matrix to a temporary matrix, and the temporary one I put as input to:
unghr
to get the Q matrix.
After that:
zhseqr
with
job = 'S' and compz = 'V'; matrix inputs are the Hessemberg matrix, and as vl the Q matrix.
Then I call
zgebak
with the same Q matrix (modified by hseqr, so eigenvectors are in).
Finally
ztrevc
is called.
Also, I wish to know if somebody have some example of how to use this functions.
Thank you.
