[SciPy-dev] float exception, fblas cgemv - newscipy

[Arnd Baecker]

Hi,

on another debian sarge based 32Bit machine (where I made a couple of
links so that ATLAS is found) I get the following

check_default_beta_y (scipy.linalg.fblas.test_fblas.test_cgemv) ... zsh:
8527 floating point exception

The shortest example to reproduce this is

gdb
file /usr/bin/python
run

from scipy import *
from scipy.basic.random import normal
a = normal(0.,1.,(3,3))
x = arange(shape(a)[0])
print a
print x
linalg.fblas.cgemv(1,a,x)

Giving

>>> print a
[[-0.70685647  0.65494695  0.9179007 ]
 [-0.30862233  0.58137533  2.01614468]
 [ 0.18585164 -0.8134587   1.23613689]]
>>> print x
[0 1 2]
>>> linalg.fblas.cgemv(1,a,x)

Program received signal SIGFPE, Arithmetic exception.
[Switching to Thread 16384 (LWP 26791)]
0x404167a0 in ATL_cgemvC_a1_x1_bXi0_y1 () from /usr/lib/sse2/libatlas.so.3

Does anyone else see this?
Any ideas what to do with this?

Best, Arnd


System details:

  ipython
  import scipy
  scipy.__core_version__
  scipy.__scipy_version__
  scipy.__scipy_config__.show()

gives


Python 2.3.5 (#2, Sep  4 2005, 22:01:42)

In [1]:import scipy
Importing io to scipy
Importing fftpack to scipy
Importing special to scipy
Importing cluster to scipy
Importing sparse to scipy
Importing utils to scipy
Importing interpolate to scipy
Importing integrate to scipy
Importing signal to scipy
Importing optimize to scipy
Importing linalg to scipy
Importing stats to scipy
In [2]:scipy.__core_version__
Out[2]:'0.4.3.1376'
In [3]:scipy.__scipy_version__
Out[3]:'0.4.2_1400'
In [4]:scipy.__scipy_config__.show()
dfftw_info:
  NOT AVAILABLE

blas_opt_info:
    libraries = ['f77blas', 'cblas', 'atlas']
    library_dirs = ['/usr/lib/sse2']
    define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')]
    language = c

djbfft_info:
  NOT AVAILABLE

atlas_blas_threads_info:
  NOT AVAILABLE

lapack_opt_info:
    libraries = ['lapack', 'f77blas', 'cblas', 'atlas']
    library_dirs = ['/usr/lib/atlas', '/usr/lib/']
    define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')]
    language = f77

atlas_info:
    libraries = ['lapack', 'f77blas', 'cblas', 'atlas']
    library_dirs = ['/usr/lib/atlas', '/usr/lib/']
    language = f77
    define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')]

fftw_info:
  NOT AVAILABLE

atlas_blas_info:
    libraries = ['f77blas', 'cblas', 'atlas']
    library_dirs = ['/usr/lib/sse2']
    language = c
    define_macros = [('ATLAS_INFO', '"\\"3.6.0\\""')]

atlas_threads_info:
  NOT AVAILABLE




Atlas:

ls -l /usr/lib/atlas
total 17168
lrwxrwxrwx  1 root root      12 Oct 25 00:32 libblas.so -> libblas.so.3
lrwxrwxrwx  1 root root      14 Mar 23  2005 libblas.so.2 ->
libblas.so.2.3
-rw-r--r--  1 root root 3051888 Mar  4  2005 libblas.so.2.3
lrwxrwxrwx  1 root root      14 Nov 18  2004 libblas.so.3 ->
libblas.so.3.0
-rw-r--r--  1 root root 3443424 Oct 27  2004 libblas.so.3.0
lrwxrwxrwx  1 root root      14 Oct 25 00:32 liblapack.so ->
liblapack.so.3
lrwxrwxrwx  1 root root      16 Mar 23  2005 liblapack.so.2 ->
liblapack.so.2.3
-rw-r--r--  1 root root 5500688 Mar  4  2005 liblapack.so.2.3
lrwxrwxrwx  1 root root      16 Nov 18  2004 liblapack.so.3 ->
liblapack.so.3.0
-rw-r--r--  1 root root 5537840 Oct 27  2004 liblapack.so.3.0
drwxr-xr-x  2 root root    4096 Nov 18  2004 sse2


ls -l /usr/lib/  | grep atlas
drwxr-xr-x   3 root root     4096 Oct 25 00:32 atlas
lrwxrwxrwx   1 root root       16 Mar 22  2004 libatlas.so ->
sse2/libatlas.so
lrwxrwxrwx   1 root root       15 Mar 23  2005 libatlas.so.2 ->
libatlas.so.2.3
-rw-r--r--   1 root root  2835272 Mar  4  2005 libatlas.so.2.3
lrwxrwxrwx   1 root root       15 Nov 18  2004 libatlas.so.3 ->
libatlas.so.3.0
-rw-r--r--   1 root root  3234288 Oct 27  2004 libatlas.so.3.0
lrwxrwxrwx   1 root root       34 Mar 22  2004 liblapack.so ->
/usr/lib/atlas/sse2/liblapack.so.3
lrwxrwxrwx   1 root root       25 Mar 22  2004 liblapack_atlas.so ->
sse2/liblapack_atlas.so.3
lrwxrwxrwx   1 root root       22 Mar 23  2005 liblapack_atlas.so.2 ->
liblapack_atlas.so.2.3
-rw-r--r--   1 root root    60568 Mar  4  2005 liblapack_atlas.so.2.3
lrwxrwxrwx   1 root root       22 Nov 18  2004 liblapack_atlas.so.3 ->
liblapack_atlas.so.3.0
-rw-r--r--   1 root root   131224 Oct 27  2004 liblapack_atlas.so.3.0