LAPACK Archives

[Lapack] incorrect RWORK size for ZGESDD

I'm one of the developers of Julia.  We were following the recommendation
for RWORK size calculation from the ZGESDD header here (, and were finding that when
INFO='N', the size was insufficient, and generally led to a segmentation

The current definition says that "If JOBZ = 'N', LRWORK >= 5*min(M,N)", but
we found that the minimum size needed to be at least 7*min(M,N).  This was
through experimentation and comparison with code from Numpy and Octave, so
you may wish to verify this through the code.

Relevant Julia links:

      - fix:
      - fix:<>


-------------- next part --------------
An HTML attachment was scrubbed...

<Prev in Thread] Current Thread [Next in Thread>

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