[SciPy-user] Failures in symeig.test.test()

[Nils Wagner]

Hi all,

Finally I have build scipy from cvs using ATLAS.

For the external package symeig I have used FORTRAN LAPACK/BLAS
(see http://mdp-toolkit.sourceforge.net/symeig.html)

by setting the environment variables

setenv ATLAS None
setenv LAPACK None
setenv BLAS None
setenv LAPACK_SRC /var/tmp/src/lapack
setenv BLAS_SRC /var/tmp/src/blas

 >>> import symeig
 >>> symeig.test.test()
Random Seed:  (1268049219, 2102953867)
testComplex (symeig.test.test_symeig.SymeigTestCase) ... ERROR
testComplexGeneralized (symeig.test.test_symeig.SymeigTestCase) ... ERROR
testIntegerMatrix (symeig.test.test_symeig.SymeigTestCase) ... ok
testNonContiguousMatrix (symeig.test.test_symeig.SymeigTestCase) ... ok
testOverwriteBug (symeig.test.test_symeig.SymeigTestCase) ... ok
testReal (symeig.test.test_symeig.SymeigTestCase) ... ERROR
testRealGeneralized (symeig.test.test_symeig.SymeigTestCase) ... ERROR
testTypecodeConversion (symeig.test.test_symeig.SymeigTestCase) ... ok

======================================================================
ERROR: testComplex (symeig.test.test_symeig.SymeigTestCase)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/usr/lib/python2.3/site-packages/symeig/test/test_symeig.py", 
line 126, in testComplex
    self.eigenproblem(DIM,'D',0,'off',range)
  File "/usr/lib/python2.3/site-packages/symeig/test/test_symeig.py", 
line 53, in eigenproblem
    assert_array_almost_equal(diag, w, DIGITS[typecode])
  File "/usr/lib/python2.3/site-packages/symeig/test/testing_tools.py", 
line 17, in assert_array_almost_equal_diff
    maxdiff = max(scipy.ravel(abs(x-y)))/\
ZeroDivisionError: float division

======================================================================
ERROR: testComplexGeneralized (symeig.test.test_symeig.SymeigTestCase)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/usr/lib/python2.3/site-packages/symeig/test/test_symeig.py", 
line 141, in testComplexGeneralized
    self.geneigenproblem(DIM,'D',0,'off',range)
  File "/usr/lib/python2.3/site-packages/symeig/test/test_symeig.py", 
line 68, in geneigenproblem
    assert_array_almost_equal(diag1, w, DIGITS[typecode])
  File "/usr/lib/python2.3/site-packages/symeig/test/testing_tools.py", 
line 17, in assert_array_almost_equal_diff
    maxdiff = max(scipy.ravel(abs(x-y)))/\
ZeroDivisionError: float division

======================================================================
ERROR: testReal (symeig.test.test_symeig.SymeigTestCase)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/usr/lib/python2.3/site-packages/symeig/test/test_symeig.py", 
line 99, in testReal
    self.eigenproblem(DIM,'d',0,'off',range)
  File "/usr/lib/python2.3/site-packages/symeig/test/test_symeig.py", 
line 53, in eigenproblem
    assert_array_almost_equal(diag, w, DIGITS[typecode])
  File "/usr/lib/python2.3/site-packages/symeig/test/testing_tools.py", 
line 17, in assert_array_almost_equal_diff
    maxdiff = max(scipy.ravel(abs(x-y)))/\
ZeroDivisionError: float division

======================================================================
ERROR: testRealGeneralized (symeig.test.test_symeig.SymeigTestCase)
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/usr/lib/python2.3/site-packages/symeig/test/test_symeig.py", 
line 114, in testRealGeneralized
    self.geneigenproblem(DIM,'d',0,'off',range)
  File "/usr/lib/python2.3/site-packages/symeig/test/test_symeig.py", 
line 68, in geneigenproblem
    assert_array_almost_equal(diag1, w, DIGITS[typecode])
  File "/usr/lib/python2.3/site-packages/symeig/test/testing_tools.py", 
line 17, in assert_array_almost_equal_diff
    maxdiff = max(scipy.ravel(abs(x-y)))/\
ZeroDivisionError: float division

----------------------------------------------------------------------
Ran 8 tests in 1.127s

FAILED (errors=4)

The former segmentation fault vanishes if  I  use FORTRAN LAPACK/BLAS 
instead of ATLAS.

Can someone reproduce these failures ?

Nils