From b41b7aa47263069794edfa57a4a0dc8ed598def7 Mon Sep 17 00:00:00 2001 From: julie Date: Mon, 22 Sep 2014 15:38:39 +0000 Subject: [PATCH] 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 --- SRC/dgesdd.f | 2 +- SRC/sgesdd.f | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/SRC/dgesdd.f b/SRC/dgesdd.f index 6fc7533..886d9b7 100644 --- 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 --- 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. -- 2.7.4