[SciPy-dev] scipy.linalg.eigvals error on 64bit RH3 & RH4

[David Cournapeau]

David Martin wrote:
> Python 2.5.4 (r254:67916, May 20 2009, 22:57:32) [GCC 3.2.3 20030502 
> (Red Hat Linux 3.2.3-56)] on linux2
> Type "help", "copyright", "credits" or "license" for more information.
>  >>> from scipy.linalg import atlas_version
>  >>> atlas_version
> <module 'scipy.linalg.atlas_version' from 
> '/home/tester/master/lib/python2.5/site-packages/scipy/linalg/atlas_version.so'>
>  >>> atlas_version.version()
> ATLAS version 3.8.3 built by builder on Mon Jul 13 20:09:54 CDT 2009:
>    UNAME    : Linux centos-3-x86_64-build 2.4.21-47.EL #1 SMP Tue Aug 1 
> 07:58:43 EDT 2006 x86_64 x86_64 x86_64 GNU/Linux
>    INSTFLG  : -1 0 -a 1
>    ARCHDEFS : -DATL_OS_Linux -DATL_ARCH_Core2 -DATL_CPUMHZ=2492 
> -DATL_SSE3 -DATL_SSE2 -DATL_SSE1 -DATL_USE64BITS -DATL_GAS_x8664
>    F2CDEFS  : -DAdd__ -DF77_INTEGER=int -DStringSunStyle
>    CACHEEDGE: 262144
>    F77      : g77, version GNU Fortran (GCC 3.2.3 20030502 (Red Hat 
> Linux 3.2.3-56)) 3.2.3 20030502 (Red Hat Linux 3.2.3-56)
>    F77FLAGS : -O -fPIC -g -m64
>    SMC      : gcc, version gcc (GCC) 3.2.3 20030502 (Red Hat Linux 3.2.3-56)
>    SMCFLAGS : -fomit-frame-pointer -mfpmath=sse -msse3 -O2 -fPIC -g -m64
>    SKC      : gcc, version gcc (GCC) 3.2.3 20030502 (Red Hat Linux 3.2.3-56)
>    SKCFLAGS : -fomit-frame-pointer -mfpmath=sse -msse3 -O2 -fPIC -g -m64
>   

There are some problems with gcc 3.2/3.3 and using sse as the math unit
(-mfpmath=sse), and I would not be surprised if g77 inherits them as
well. Maybe Red Hat backported fixes, though, you would have to see this
with your OS provider. One way to confirm this would be to compile your
own gcc/gfortran, and recompile everything (lapack, atlas, numpy and
scipy) with it.

cheers,

David