[Numpy-discussion] Bug in built in SVD?

Nadav Horesh nadavh at visionsense.com
Thu Aug 12 00:15:07 EDT 2004

Sorry it is a gcc3.4.1 optimization bug. Either lowering to optimization to -O1 or stepping back to gcc3.3.4 (with all optimization guns on) fix it.


-----Original Message-----
From:	Nadav Horesh
Sent:	Wed 11-Aug-04 15:48
To:	numpy-discussion
Subject:	[Numpy-discussion] Bug in built in SVD?
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?


SF.Net email is sponsored by Shop4tech.com-Lowest price on Blank Media
100pk Sonic DVD-R 4x for only $29 -100pk Sonic DVD+R for only $33
Save 50% off Retail on Ink & Toner - Free Shipping and Free Gift.
Numpy-discussion mailing list
Numpy-discussion at lists.sourceforge.net

More information about the NumPy-Discussion mailing list