[Numpy-discussion] Bug in built in SVD?

[Nadav Horesh]

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

example:

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

  Nadav.