[SciPy-Dev] LAPACK problems with *gges wrappers

[Skipper Seabold]

On Thu, May 3, 2012 at 2:08 PM, Denis Laxalde <denis at laxalde.org> wrote:
> Skipper Seabold a écrit :
>> Thanks. Can you check again with pushed fix?
>
> All test pass now.

Thanks again.

>
>> I think the documentation is wrong and the actual minimum value for
>> lwork should be max(8*n, 6*n+16) and not the 8*n+16 in the
>> documentation.
>>
>> If anyone cares to double check my reading of the source.
>>
>> http://www.netlib.no/netlib/lapack/double/dgges.f
>
> That's indeed what the code checks for. However 8*n+16 > max(8*n,
> 6*n+16) so the documentation indicates a conservative value.

Sure. I guess I read the docs as saying that LWORK _should be_ >=
8*n+16, so that's what I also checked in the wrappers. However, this
is not always the case given an LWORK query.

Skipper