I have used the SCALAPACK routines PDSYEV and PZHEEV for diagonalization of symmetric and Hermitian matrices respectively. Using these routines to diagonalize a square matrix in a square process grid (for ex.- 2 x 2, 4 x 4) I have got correct eigenvalues and eigenvectors. The problem arises when I choose a rectangular process grid (for ex.- 2 x 4, 4 x 12 ) I get incorrect eigenvalues and eigenvectors. I wanted to know whether the SCALAPACK routines PDSYEV and PZHEEV can be applied to rectangular process grid or not.
Link Copied
Hi Pranay,
Symmetric eigensolver drivers in ScaLAPACK support rectangular process grid. Can you please provide a reproducer of the problem?
Regards,
Konstantin
For more complete information about compiler optimizations, see our Optimization Notice.