[Numpy-discussion] Bug in built in SVD?

Nadav Horesh nadavh at visionsense.com
Wed Aug 11 04:49:02 CDT 2004

When numarray is compiled with the buit-in linear-algebra routine the svd routine runs into an infinite loop.

Platform2: python 2.3.4 and 2.4a2; numarray 1.1; gcc 3.4.1; RH9 and gentoo 1.4


>> from numarray.linear_algebra.mlab import *
>>> svd(arange(20.0, shape = (5,4)))

... infine loop ...


svd runs OK when compiled againts external LAPACK.

BTW: Why numarray used dgesdd routine for SVD, and not the svd routine that has the ATLAS equivalent?


More information about the Numpy-discussion mailing list