Correcting Bug111 reported by Kevin Squire on LAPACK Mailing - August 15th 2013
authorjulie <julielangou@users.noreply.github.com>
Wed, 18 Mar 2015 05:09:27 +0000 (05:09 +0000)
committerjulie <julielangou@users.noreply.github.com>
Wed, 18 Mar 2015 05:09:27 +0000 (05:09 +0000)
commit7819b90287ba881b8cddde78d45f891407ffcc48
treec333ec47b1c861a708e71c0d47f11d246c0fed70
parent52508f6632c094c4431c6fe5275ab9d9ea1c3842
Correcting Bug111 reported by Kevin Squire on LAPACK Mailing - August 15th 2013

From Kevin: 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.

JL: bug confirmed and reproduced - correction in comments
If JOBZ = N, LRWORK >= 7*min(M,N)
SRC/cgesdd.f
SRC/zgesdd.f