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.
