[SciPy-dev] test_lstsq failure

[Nils Wagner]

Arnd Baecker wrote:
> Hi,
>
> on a 64 Bit opteron I get with current svn:
>
> ======================================================================
> FAIL: check_random_complex_overdet
> (scipy.linalg.tests.test_basic.test_lstsq)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
>   File
> "/home/abaecker/BUILDS3/BuildDir/inst_numpy/lib/python2.4/site-packages/scipy/linalg/tests/test_basic.py",
> line 413, in check_random_complex_overdet
>     assert_array_almost_equal(x,direct_lstsq(a,b),3)
>   File
> "/home/abaecker/BUILDS3/BuildDir/inst_numpy/lib/python2.4/site-packages/numpy/testing/utils.py",
> line 233, in assert_array_almost_equal
>     assert cond,\
> AssertionError:
> Arrays are not almost equal (mismatch 77.7777777778%):
>         Array 1: [[-0.0137+0.0173j  0.0037-0.0173j -0.0114+0.0119j]
>  [ 0.0029-0.0356j  0.0086-0.034j   0.033 -0.0879j]
>  [ 0.0041-0.0097j ...
>         Array 2: [[-0.016 +0.0162j  0.003 -0.0171j -0.0148+0.009j ]
>  [-0.0017-0.0405j  0.003 -0.036j   0.0256-0.0977j]
>  [ 0.0038-0.0112j ...
>
> Best, Arnd
>
> _______________________________________________
> Scipy-dev mailing list
> Scipy-dev at scipy.net
> http://www.scipy.net/mailman/listinfo/scipy-dev
>   
r1934 | cookedm | 2006-06-08 05:12:30 +0200 (Thu, 08 Jun 2006) | 4 lines

Fix scipy.linalg.tests.test_basic.test_lstsq.check_random_complex_overdet
to use # of decimals, not a tolerance, as argument to
assert_array_almost_equal.
This exposes the bug again.