[Numpy-discussion] svd error checking vs. speed

[Dave Hirschfeld]

Sturla Molden <sturla.molden <at> gmail.com> writes:

> 
> <josef.pktd <at> gmail.com> wrote:
> 
> > I use official numpy release for development, Windows, 32bit python,
> > i.e. MingW 3.5 and whatever old ATLAS the release includes.
> > 
> > a constant 13% cpu usage is 1/8 th of my 8 virtual cores.
> 
> Based on this and Alex' message it seems the offender is the f2c generated
> lapack_lite library.
> 
> So do we do with lapack_lite? Should we patch it?
> 
> Sturla
> 

Even if lapack_lite always performed the isfinite check and threw a python 
error if False, it would be much better than either hanging or segfaulting and 
people who care about the isfinite cost probably would be linking to a fast 
lapack anyway.

-Dave