diff options
| author | Jason Grout <jason.grout@drake.edu> | 2011-04-26 10:18:50 -0500 |
|---|---|---|
| committer | Pauli Virtanen <pav@iki.fi> | 2011-04-26 19:56:57 +0200 |
| commit | cc0c0cbde7fac047bdbebd33db8508fd07665f6d (patch) | |
| tree | 9b1593d8530c90f0ef32de912afa9f80a55429b9 /numpy/linalg | |
| parent | 9752349d0e6f21298b99d2175e9a7846c0589c91 (diff) | |
| download | numpy-cc0c0cbde7fac047bdbebd33db8508fd07665f6d.tar.gz | |
BUG: linalg: Update the rwork size in the call to zgesdd to reflect a documentation bugfix in Lapack SVN revision 729.
See http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=2&t=1779
Diffstat (limited to 'numpy/linalg')
| -rw-r--r-- | numpy/linalg/linalg.py | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/numpy/linalg/linalg.py b/numpy/linalg/linalg.py index 44396f075..36602e43e 100644 --- a/numpy/linalg/linalg.py +++ b/numpy/linalg/linalg.py @@ -1297,7 +1297,8 @@ def svd(a, full_matrices=1, compute_uv=1): iwork = zeros((8*min(m, n),), fortran_int) if isComplexType(t): lapack_routine = lapack_lite.zgesdd - rwork = zeros((5*min(m, n)*min(m, n) + 5*min(m, n),), real_t) + lrwork = min(m,n)*max(5*min(m,n)+7, 2*max(m,n)+2*min(m,n)+1) + rwork = zeros((lrwork,), real_t) lwork = 1 work = zeros((lwork,), t) results = lapack_routine(option, m, n, a, m, s, u, m, vt, nvt, |