octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[PATCH] Add workaround for broken xGELSD workspace queries.


From: John W. Eaton
Subject: [PATCH] Add workaround for broken xGELSD workspace queries.
Date: Tue, 12 Feb 2008 23:49:34 -0500

On 12-Feb-2008, Jason Riedy wrote:

| The xGELSD workspace queries are inconsistant with the logic inside
| those routines in the short and fat case (m > n).  For that case,
| calculate a workspace size large enough to enable the efficient
| algorithm.  This may allocate more memory than necessary for small
| cases.  The calculated alternate size is at least the documented
| minimum workspace requirement, so this amount always will work even
| if the routines change.
| ---
|    Sorry it's been so long.  I don't want to describe the LAPACK
|    developer dynamic, but it's severely busted.  ergh.  I think
|    this is a minimal change, and the ILAENV declaration should
|    work on every Fortran compiler I know. -- Jason

Just to avoid possible incompatible C<->Fortran calling conventions
for Fortran funtions, we have always wrapped them in Fortran
subroutines (see the files in libcruft/lapack-xtra).  I did that for
ilaenv and checked in your change.

Thanks,

jwe


reply via email to

[Prev in Thread] Current Thread [Next in Thread]