LAPACK Archives

[Lapack] IWORK size autodetection in DGELSD

Dear Tamas,
Thank you for your bug report.
This already has been corrected in our SVN repository.
See our Errata page: http://www.netlib.org/lapack/Errata/ (bug0038)
Sincerely,
Julie Langou
On Jun 8, 2010, at 9:09 AM, Tamas Papp wrote:

Dear LAPACK Team,

I am writing a Common Lisp wrapper for LAPACK.  It appears that for
all the xGELSD routines _but_ DGELSD, the first element of the IWORK
array can be used for autodetecting the workspace size, eg in ZGELSD
[1]:

129: *  IWORK   (workspace) INTEGER array, dimension (MAX(1,LIWORK))
130: *          LIWORK >= max(1, 3*MINMN*NLVL + 11*MINMN),
131: *          where MINMN = MIN( M,N ).
132: *          On exit, if INFO = 0, IWORK(1) returns the minimum LIWORK.

But for DGELSD, it only says that:

116: *  IWORK   (workspace) INTEGER array, dimension (MAX(1,LIWORK))
117: *          LIWORK >= 3 * MINMN * NLVL + 11 * MINMN,
118: *          where MINMN = MIN( M,N ).

with no reference to IWORK(1).  And it appears that IWORK(1) is
sometimes 0 after the call.  I am wondering if this discrepancy is
intentional, and if not, whether it could be fixed.  Please let me
know if I misunderstood something.

Thanks,

Tamas Papp

[1] http://www.netlib.org/lapack/explore-html/zgelsd.f.html
[2] http://www.netlib.org/lapack/explore-html/dgelsd.f.html
_______________________________________________
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/20100609/6374af1f/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