Workspace for DGESDD in LAPACK 3.4.2

Open discussion regarding features, bugs, issues, vendors, etc.

Workspace for DGESDD in LAPACK 3.4.2

Postby simon2212 » Mon Oct 08, 2012 4:36 pm

The size of the workspace for DGESDD for JOBZ = 'S' or 'A' has been changed in LAPACK 3.4.2, and is set by LWORK. The current description of LWORK reads:

Code: Select all
LWORK >= 3*min(M,N) +
         max(max(M,N),4*min(M,N)*min(M,N)+3*min(M,N)+max(M,N)).


In the second line, both arguments of max contain max(M,N), and so the the function is redundant - the second argument will always be realised.
So the description could be simplified to:

Code: Select all
LWORK >= 3*min(M,N) +
         4*min(M,N)*min(M,N)+3*min(M,N)+max(M,N)


or further simplified to

Code: Select all
LWORK >= min(M,N)*(6+4*min(M,N))+max(M,N)
simon2212
 
Posts: 1
Joined: Mon Oct 08, 2012 4:22 pm

Re: Workspace for DGESDD in LAPACK 3.4.2

Postby admin » Mon Oct 08, 2012 5:13 pm

Indeed Simon,
Thank you this has been corrected in our repository.
Julie
admin
Site Admin
 
Posts: 468
Joined: Wed Dec 08, 2004 7:07 pm


Return to User Discussion

Who is online

Users browsing this forum: Yahoo [Bot] and 2 guests