[Numpy-discussion] Test failures r7300

[lukshuntim at gmail.com]

David Cournapeau wrote:
> On Sat, Aug 8, 2009 at 10:33 PM, <lukshuntim at gmail.com> wrote:
>> David Cournapeau wrote:
>>> On Sat, Aug 8, 2009 at 9:38 PM, <lukshuntim at gmail.com> wrote:
>>>> Hi,
>>>>
>>>> I got 16 test failures after building r7300 from svn on debian/sid/i386.
>>>> Seems all related to complex linear algebra modules.
>>> Are you using atlas ? (numpy.show_config() output)
>> Yes, it's libatlas-sse2 3.6.0-24 debian/sid package.
> 
> I wonder if debian atlas package has the same problem as on recent Ubuntu.

[snipped]

> What I would try is first install libatlas-base (or whatever it is
> called on sid), i.e. the non sse version, and compare test output with
> both sse2/nosse (e.g. using LD_LIBRARY_PATH to point to /usr/lib so
> that the nosse is loaded, you can check using ldd which one is loaded
> by ld).

Just to clarify, you mean doing a "ldd lapack_lite.so" to check which
blas and lapack is used at runtime. Right?

I also removed site.cfg when building. With no atlas, and with both
libatlas-base and libatlas-sse, the complex linear algebra errors went
away and I got only 1 error:

<quote>
FAIL: Test bug in reduceat with structured arrays copied for speed.
----------------------------------------------------------------------
Traceback (most recent call last):
  File "/usr/lib/pymodules/python2.5/nose/case.py", line 183, in runTest
    self.test(*self.arg)
  File
"/var/opt/py/lib/python2.5/site-packages/numpy/core/tests/test_umath.py",
line 818, in test_reduceat
    assert_array_almost_equal(h1, h2)
  File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py",
line 726, in assert_array_almost_equal
    header='Arrays are not almost equal')
  File "/var/opt/py/lib/python2.5/site-packages/numpy/testing/utils.py",
line 571, in assert_array_compare
    raise AssertionError(msg)
AssertionError:
Arrays are not almost equal

(mismatch 100.0%)
 x: array([ -5.17543592e+11,  -5.17543592e+11,  -5.17543592e+11,
        -5.17543592e+11], dtype=float32)
 y: array([  700.,   800.,  1000.,  7500.], dtype=float32)
</quote>

So it appears that it's the sse2 variant that is causing the problem.

Regards,
ST
--