[Numpy-discussion] Bug in built in SVD?

[Nadav Horesh]

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.

   Nadav.

-----Original Message-----
From:	Nadav Horesh
Sent:	Wed 11-Aug-04 15:48
To:	numpy-discussion
Cc:	
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

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.



-------------------------------------------------------
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.
http://www.shop4tech.com/z/Inkjet_Cartridges/9_108_r285
_______________________________________________
Numpy-discussion mailing list
Numpy-discussion at lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/numpy-discussion