Correction to bug120

Fix workspace definition in comment following bug report from Rasmus Munk Larsen, Google Knowledge & Research
See forum https://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4587

See revision 1349, and 649 for history

diff --git a/SRC/dgesdd.f b/SRC/dgesdd.f
index 6fc7533..886d9b7 100644 (file)
--- a/SRC/dgesdd.f
+++ b/SRC/dgesdd.f
@@ -175,7 +175,7 @@
 *>            LWORK >= 3*min(M,N) + 
 *>                     max(max(M,N),5*min(M,N)*min(M,N)+4*min(M,N)).
 *>          If JOBZ = 'S' or 'A'
-*>            LWORK >= min(M,N)*(6+4*min(M,N))+max(M,N)
+*>            LWORK >= min(M,N)*(7+4*min(M,N))
 *>          For good performance, LWORK should generally be larger.
 *>          If LWORK = -1 but other input arguments are legal, WORK(1)
 *>          returns the optimal LWORK.

diff --git a/SRC/sgesdd.f b/SRC/sgesdd.f
index 746fd3b..6dccb36 100644 (file)
--- a/SRC/sgesdd.f
+++ b/SRC/sgesdd.f
@@ -175,7 +175,7 @@
 *>            LWORK >= 3*min(M,N) + 
 *>                     max(max(M,N),5*min(M,N)*min(M,N)+4*min(M,N)).
 *>          If JOBZ = 'S' or 'A'
-*>            LWORK >= min(M,N)*(6+4*min(M,N))+max(M,N)
+*>            LWORK >= min(M,N)*(7+4*min(M,N))
 *>          For good performance, LWORK should generally be larger.
 *>          If LWORK = -1 but other input arguments are legal, WORK(1)
 *>          returns the optimal LWORK.