Hi,
I need to calculate all the eigenvalues of a large (double) complex matrix. I have a program to do this using LAPACK with a call to zgehrd, followed by a call to zlahqr. The problem is, when I try to emulate the program in SacLAPACK using pzgehrd and pzlahqr, the eigenvalues returned do not match.
I've tried by-passing pzgehrd (by reading in the output of zgehrd instead) and so I know that the problem is isolated to pzlahqr. I've looked at the "current notes / restrictions" section of pzlahqr and tried to satisfy all those conditions, but still no joy.
Does anyone know of any other unrecorded problems with pzlahqr, or something I might be doing wrong?
Thanks,
Ian