ScaLAPACK Archives

[Scalapack] question on scalapack use

Dear Scalapack gurus,

I would like to compute all eigenvalues of a large nonsymmetric matrix (perhaps 
1 million x 1 million). The matrix is dense but block off-diagonal:

          0      A
M =                          where A and B are full.
          B      0

As I understand, in Scalapack I have to first use PDGEHRD to reduce the matrix 
to upper Hessenberg form and then PDLAHQR to compute the eigenvalues.

But these subroutines cannot take advantage of the block structure, can they? 
Is there another possibility how to compute the eigenvalues efficiently?

Thank you for your answers, Tomas Dohnal.
-------------- next part --------------
An HTML attachment was scrubbed...

<Prev in Thread] Current Thread [Next in Thread>
  • [Scalapack] question on scalapack use, tomas dohnal <=

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