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

[David Martin]

Robert Kern wrote:
> On Mon, Jul 13, 2009 at 16:43, David Martin<frostedcheerios@gmail.com> wrote:
>   
>> Has anyone else gotten these scipy.test() errors when running on RH3 and
>> RH4 (both 64bit)? The following tests pass on RH5 64-bit and RH3,4,5
>> 32-bit, so this seems to be a 64-bit RH issue prior to RH5. It looks
>> like eigvals isn't returning the correct values:
>>     
>
> Try to run the ATLAS tests for the version you built scipy against, if
> you can. Incorrect values from eigvals() usually mean that the
> underlying library is having problems, not scipy.
>
>   
Unfortunately I only had the .so libraries to link against so I couldn't 
run the ATLAS sanity tests against them. I built a new set of ATLAS 
libraries on RH3 64-bit and ran "make check" after the build, which 
showed all tests passing.

DONE BUILDING TESTERS, RUNNING:
SCOPING FOR FAILURES IN BIN TESTS:
fgrep -e fault -e FAULT -e error -e ERROR -e fail -e FAIL \
        bin/sanity.out
8 cases: 8 passed, 0 skipped, 0 failed
4 cases: 4 passed, 0 skipped, 0 failed
8 cases: 8 passed, 0 skipped, 0 failed
4 cases: 4 passed, 0 skipped, 0 failed
8 cases: 8 passed, 0 skipped, 0 failed
4 cases: 4 passed, 0 skipped, 0 failed
8 cases: 8 passed, 0 skipped, 0 failed
4 cases: 4 passed, 0 skipped, 0 failed
DONE
SCOPING FOR FAILURES IN CBLAS TESTS:
fgrep -e fault -e FAULT -e error -e ERROR -e fail -e FAIL \
        interfaces/blas/C/testing/sanity.out | \
                fgrep -v PASSED
DONE
SCOPING FOR FAILURES IN F77BLAS TESTS:
fgrep -e fault -e FAULT -e error -e ERROR -e fail -e FAIL \
        interfaces/blas/F77/testing/sanity.out | \
                fgrep -v PASSED
DONE

I rebuilt Numpy and Scipy against this new ATLAS build and re-ran both 
the Numpy and Scipy tests. The Numpy tests all pass, but most of the 
Scipy failures persist (apparently the new ATLAS libraries fixed the 
test_speigs.TestEigs issue):


======================================================================
FAIL: test_simple (test_decomp.TestEig)
----------------------------------------------------------------------
Traceback (most recent call last):
  File 
"/home/builder/master-1.0.4/lib/python2.5/site-packages/scipy/linalg/tests/test_decomp.py", 
line 145, in test_simple
    assert_array_almost_equal(w,exact_w)
  File 
"/home/builder/recipes/install/numpy-1.3.0-1.egg/numpy/testing/utils.py", 
line 537, in assert_array_almost_equal
    header='Arrays are not almost equal')
  File 
"/home/builder/recipes/install/numpy-1.3.0-1.egg/numpy/testing/utils.py", 
line 395, in assert_array_compare
    raise AssertionError(msg)
AssertionError:
Arrays are not almost equal
 
(mismatch 66.6666666667%)
 x: array([ 9.43719064+0.j, -0.11536526+0.j, -0.32182538+0.j])
 y: array([ 9.32182538,  0.        , -0.32182538])
 
======================================================================
FAIL: test_simple (test_decomp.TestEigVals)
----------------------------------------------------------------------
Traceback (most recent call last):
  File 
"/home/builder/master-1.0.4/lib/python2.5/site-packages/scipy/linalg/tests/test_decomp.py", 
line 114, in test_simple
    assert_array_almost_equal(w,exact_w)
  File 
"/home/builder/recipes/install/numpy-1.3.0-1.egg/numpy/testing/utils.py", 
line 537, in assert_array_almost_equal
    header='Arrays are not almost equal')
  File 
"/home/builder/recipes/install/numpy-1.3.0-1.egg/numpy/testing/utils.py", 
line 395, in assert_array_compare
    raise AssertionError(msg)
AssertionError:
Arrays are not almost equal
 
(mismatch 66.6666666667%)
 x: array([ 9.43719064+0.j, -0.11536526+0.j, -0.32182538+0.j])
 y: array([ 9.32182538,  0.        , -0.32182538])
 
======================================================================
FAIL: test_simple_tr (test_decomp.TestEigVals)
----------------------------------------------------------------------
Traceback (most recent call last):
  File 
"/home/builder/master-1.0.4/lib/python2.5/site-packages/scipy/linalg/tests/test_decomp.py", 
line 122, in test_simple_tr
    assert_array_almost_equal(w,exact_w)
  File 
"/home/builder/recipes/install/numpy-1.3.0-1.egg/numpy/testing/utils.py", 
line 537, in assert_array_almost_equal
    header='Arrays are not almost equal')
  File 
"/home/builder/recipes/install/numpy-1.3.0-1.egg/numpy/testing/utils.py", 
line 395, in assert_array_compare
    raise AssertionError(msg)
AssertionError:
Arrays are not almost equal
 
(mismatch 66.6666666667%)
 x: array([ 9.43719064+0.j, -0.11536526+0.j, -0.32182538+0.j])
 y: array([ 9.32182538,  0.        , -0.32182538])
 
----------------------------------------------------------------------


-David