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 (
http://netlib.org/lapack/complex16/zgesdd.f), and were finding that when
INFO='N', the size was insufficient, and generally led to a segmentation
fault.

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:
<https://github.com/JuliaLang/julia/commit/484a9f0892>

   - https://github.com/JuliaLang/julia/issues/3966
      - fix: https://github.com/JuliaLang/julia/commit/484a9f0892
   - https://github.com/JuliaLang/julia/issues/4016
      - fix: 
https://github.com/JuliaLang/julia/commit/9f6058804c<https://github.com/JuliaLang/julia/commit/9f6058804ca2b>

Thanks,

   Kevin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: 
<http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20130813/a51ca8c7/attachment.html>

<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