LAPACK Archives

[Lapack] One more bug in LAPACK: real workspace size in (C/Z)GELSD is un

Thank you very much Alexander,
I applied your patch in our repository. The bug is bug0059.
It is included in the 3.2.2 release.
Do not hesitate to contact us if you spot other bugs.
Best Regards,
Julie and Julien

On Jun 23, 2010, at 10:03 PM, Kobotov, Alexander V wrote:

Hello,
 
ZGELSD crashes on large matrices (eq M=N=NHRS>=88), if real workspace 
allocated exactly as reported on workspace query.
Root cause is wrong estimation of real workspace size in ZLALSA (on path 
zgelsd->zlalsd->zlalsa). It assumes that RWORK array length should be max( N, 
(SMLSZ+1)*NRHS*3 ), where SMLSIZ=25. But in fact the function calls zlals0 
which requires RWORK to be ( K*(1+NRHS) + 2*NRHS ), where K could be as large 
as number of columns in matrix A.
The same is actual for CGELSD.
Example of diff which fixes the problem is attached.
 
W.B.R.
Alexander
 
--------------------------------------------------------------------
Closed Joint Stock Company Intel A/O
Registered legal address: Krylatsky Hills Business Park, 
17 Krylatskaya Str., Bldg 4, Moscow 121614, 
Russian Federation

This e-mail and any attachments may contain confidential material for
the sole use of the intended recipient(s). Any review or distribution
by others is strictly prohibited. If you are not the intended
recipient, please contact the sender and delete all copies.
<diff.txt>_______________________________________________
Lapack mailing list
Lapack@Domain.Removed
http://lists.eecs.utk.edu/mailman/listinfo/lapack

**********************************************
Julie Langou; Research Associate in Computer Science
Innovative Computing Laboratory;
University of Tennessee from Denver, Colorado ;-)
julie@Domain.Removed; http://www.cs.utk.edu/~julie/







-------------- next part --------------
An HTML attachment was scrubbed...
URL: 
http://lists.eecs.utk.edu/mailman/private/lapack/attachments/20100624/7d77af53/attachment.html
 

<Prev in Thread] Current Thread [Next in Thread>


For additional information you may use the LAPACK/ScaLAPACK Forum.
Or one of the mailing lists, or