[SciPy-Dev] LAPACK problems with *gges wrappers

Denis Laxalde denis at laxalde.org
Thu May 3 14:08:13 EDT 2012


Skipper Seabold a écrit :
> Thanks. Can you check again with pushed fix?

All test pass now.

> 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.



More information about the SciPy-Dev mailing list