summaryrefslogtreecommitdiff
path: root/numpy/linalg/linalg.py
diff options
context:
space:
mode:
authorTravis Oliphant <oliphant@enthought.com>2006-11-02 04:35:04 +0000
committerTravis Oliphant <oliphant@enthought.com>2006-11-02 04:35:04 +0000
commit39242cbff4ad887be3d0d46dd348f769d6729d82 (patch)
treedc2eb581310ce44f71daa0b3ec16c514e3c38692 /numpy/linalg/linalg.py
parenta39c1cd2e4ec42adfdd548da86ae36a9241a58da (diff)
downloadnumpy-39242cbff4ad887be3d0d46dd348f769d6729d82.tar.gz
Fix problem with dgesdd where the optimal work size is not computed correctly.
Diffstat (limited to 'numpy/linalg/linalg.py')
-rw-r--r--numpy/linalg/linalg.py9
1 files changed, 0 insertions, 9 deletions
diff --git a/numpy/linalg/linalg.py b/numpy/linalg/linalg.py
index b2db82064..5b26f356c 100644
--- a/numpy/linalg/linalg.py
+++ b/numpy/linalg/linalg.py
@@ -590,15 +590,6 @@ def svd(a, full_matrices=1, compute_uv=1):
results = lapack_routine(option, m, n, a, m, s, u, m, vt, nvt,
work, -1, iwork, 0)
lwork = int(work[0])
- if option == 'N' and lwork==1:
- # there seems to be a bug in dgesdd of lapack
- # (NNemec, 060310)
- # returning the wrong lwork size for option == 'N'
- results = lapack_routine('A', m, n, a, m, s, u, m, vt, n,
- work, -1, iwork, 0)
- lwork = int(work[0])
- assert lwork > 1
-
work = zeros((lwork,), t)
results = lapack_routine(option, m, n, a, m, s, u, m, vt, nvt,
work, lwork, iwork, 0)