[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