Wrong lwmin in pdstedc

[bE554357]

Wrong lwmin in pdstedc

[bE554357]

UPD.

Looks like, this will lead to similar bug in p(z|c)heevd function, since minimal size of workspace for real numbers is about 2*np*nq for it; a few lines from pzheevd.f:

*  LRWORK  (local input) INTEGER
*          Size of RWORK array.
*          LRWORK >= 1 + 9*N + 3*NP*NQ,
.....
      INDZ = INDE + N
      INDRWORK = INDZ + NP0*NQ0
      LLRWORK = LRWORK - INDRWORK + 1
....
      CALL PDSTEDC( 'I', N, W, RWORK( INDE+OFFSET ), RWORK( INDZ ), IZ,
     $              JZ, DESCRZ, RWORK( INDRWORK ), LLRWORK, IWORK,
     $              LIWORK, IINFO )