fix #3966

the size of the RWORK array in zgesdd was wrong.
JuliaLang · Aug 7, 2013 · 484a9f0 · ViralBShah · Aug 7, 2013 · JeffBezanson
1 parent 24189fc
commit 484a9f0
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/base/linalg/lapack.jl b/base/linalg/lapack.jl
@@ -858,7 +858,8 @@ for (geev, gesvd, gesdd, ggsvd, elty, relty) in
             S      = Array($relty, minmn)
             cmplx  = iseltype(A,Complex)
             if cmplx
-                rwork = Array($relty, job == 'N' ? 7*minmn : 5*minmn*minmn + 5*minmn)
+                rwork = Array($relty, job == 'N' ? 5*minmn :
+                              minmn*max(5*minmn+7, 2*max(m,n)+2*minmn+1))
             end
             iwork  = Array(BlasInt, 8*minmn)
             info   = Array(BlasInt, 1)