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
|