[SciPy-dev] Status on ubuntu jaunty 64bits.
Xavier Gnata
xavier.gnata at gmail.com
Tue Mar 17 19:03:33 EDT 2009
Nils Wagner wrote:
> On Tue, 17 Mar 2009 21:21:46 +0100
> Xavier Gnata <xavier.gnata at gmail.com> wrote:
>
>> Hi,
>>
>> Here it is:
>>
>> ======================================================================
>>
>>
>> ERROR: test_implicit
>> (test_odr.TestODR)
>>
>> ----------------------------------------------------------------------
>>
>>
>> Traceback (most recent call
>> last):
>>
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/odr/tests/test_odr.py",
>> line 88, in
>> test_implicit
>>
>>
>>
>> out =
>> implicit_odr.run()
>>
>>
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/odr/odrpack.py",
>> line 1055, in run
>> self.output = Output(apply(odr, args,
>> kwds))
>> TypeError: y must be a sequence or integer (if model is
>> implicit)
>>
>> ======================================================================
>> FAIL: test_random_real (test_basic.TestSingleIFFT)
>>
>> ----------------------------------------------------------------------
>> Traceback (most recent call last):
>>
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/fftpack/tests/test_basic.py",
>> line 206, in
>> test_random_real
>>
>>
>>
>> assert_array_almost_equal (y2,
>> x)
>> File
>> "/usr/lib/python2.6/dist-packages/numpy/testing/utils.py",
>> line
>> 311, in
>> assert_array_almost_equal
>>
>>
>>
>> header='Arrays are not almost
>> equal')
>> File
>> "/usr/lib/python2.6/dist-packages/numpy/testing/utils.py",
>> line
>> 296, in
>> assert_array_compare
>>
>>
>>
>> raise
>> AssertionError(msg)
>>
>>
>> AssertionError:
>>
>>
>>
>> Arrays are not almost
>> equal
>>
>>
>>
>> (mismatch 0.900900900901%)
>> x: array([ 0.89560729 -4.65661287e-09j, 0.87991965
>> +7.21774995e-09j,
>> 0.44631395 -2.04890966e-08j, 0.71974921
>> +4.46598869e-11j,
>> 0.20776373 +1.41855736e-08j, 0.83089650
>> -1.69798398e-09j,...
>> y: array([ 0.89560729, 0.87991947, 0.44631401,
>> 0.71974903, 0.20776364,
>> 0.83089662, 0.86079419, 0.93193549,
>> 0.20852582, 0.51215041,
>> 0.91066802, 0.99397069, 0.74227983,
>> 0.67712617, 0.244197 ,...
>>
>> ======================================================================
>> FAIL: test_iv_cephes_vs_amos (test_basic.TestBessel)
>>
>> ----------------------------------------------------------------------
>> Traceback (most recent call last):
>>
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 1653, in
>> test_iv_cephes_vs_amos
>>
>>
>> self.check_cephes_vs_amos(iv, iv, rtol=1e-8,
>> atol=1e-305)
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 1642, in
>> check_cephes_vs_amos
>>
>>
>> assert_tol_equal(c1, c2, err_msg=(v, z), rtol=rtol,
>> atol=atol)
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 38, in
>> assert_tol_equal
>>
>>
>>
>> verbose=verbose,
>> header=header)
>>
>> File
>> "/usr/lib/python2.6/dist-packages/numpy/testing/utils.py",
>> line
>> 296, in
>> assert_array_compare
>>
>>
>>
>> raise
>> AssertionError(msg)
>>
>>
>> AssertionError:
>>
>>
>>
>> Not equal to tolerance rtol=1e-08,
>> atol=1e-305
>> (-120,
>> -11)
>>
>>
>>
>> (mismatch
>> 100.0%)
>>
>>
>>
>> x:
>> array(1.3384173609003782e-110)
>>
>>
>> y:
>> array((1.3384173859242368e-110+0j))
>>
>>
>>
>> ======================================================================
>> FAIL: test_yn_zeros (test_basic.TestBessel)
>>
>> ----------------------------------------------------------------------
>> Traceback (most recent call last):
>>
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 1598, in
>> test_yn_zeros
>>
>>
>>
>> 488.98055964441374646],
>> rtol=1e-19)
>>
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 38, in
>> assert_tol_equal
>>
>>
>>
>> verbose=verbose,
>> header=header)
>>
>> File
>> "/usr/lib/python2.6/dist-packages/numpy/testing/utils.py",
>> line
>> 296, in
>> assert_array_compare
>>
>>
>>
>> raise
>> AssertionError(msg)
>>
>>
>> AssertionError:
>>
>>
>>
>> Not equal to tolerance rtol=1e-19,
>> atol=0
>>
>> (mismatch 100.0%)
>> x: array([ 450.136, 463.057, 472.807, 481.274,
>> 488.981])
>> y: array([ 450.136, 463.057, 472.807, 481.274,
>> 488.981])
>>
>> ======================================================================
>> FAIL: test_ynp_zeros (test_basic.TestBessel)
>>
>> ----------------------------------------------------------------------
>> Traceback (most recent call last):
>>
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 1604, in
>> test_ynp_zeros
>>
>>
>>
>> assert_tol_equal(yvp(443, ao), 0,
>> atol=1e-15)
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 38, in
>> assert_tol_equal
>>
>>
>>
>> verbose=verbose,
>> header=header)
>>
>> File
>> "/usr/lib/python2.6/dist-packages/numpy/testing/utils.py",
>> line
>> 296, in
>> assert_array_compare
>>
>>
>>
>> raise
>> AssertionError(msg)
>>
>>
>> AssertionError:
>>
>>
>>
>> Not equal to tolerance rtol=1e-07,
>> atol=1e-15
>>
>> (mismatch 100.0%)
>> x: array([ 1.239e-10, -8.119e-16, 3.608e-16,
>> 5.898e-16, 1.226e-15])
>> y: array(0)
>>
>>
>>
>> ======================================================================
>> FAIL: test_yv_cephes_vs_amos (test_basic.TestBessel)
>>
>> ----------------------------------------------------------------------
>> Traceback (most recent call last):
>>
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 1650, in
>> test_yv_cephes_vs_amos
>>
>>
>> self.check_cephes_vs_amos(yv, yn, rtol=1e-11,
>> atol=1e-305)
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 1640, in
>> check_cephes_vs_amos
>>
>>
>> assert c2.imag != 0, (v,
>> z)
>>
>> AssertionError: (301,
>> 1.0)
>>
>>
>>
>> ======================================================================
>> FAIL: test_pbdv (test_basic.TestCephes)
>>
>> ----------------------------------------------------------------------
>> Traceback (most recent call last):
>>
>> File
>> "/usr/local/lib/python2.6/dist-packages/scipy/special/tests/test_basic.py",
>> line 370, in test_pbdv
>> assert_equal(cephes.pbdv(1,0),(0.0,0.0))
>> File
>> "/usr/lib/python2.6/dist-packages/numpy/testing/utils.py",
>> line
>> 176, in assert_equal
>> assert_equal(actual[k], desired[k], 'item=%r\n%s' %
>> (k,err_msg),
>> verbose)
>> File
>> "/usr/lib/python2.6/dist-packages/numpy/testing/utils.py",
>> line
>> 183, in assert_equal
>> raise AssertionError(msg)
>> AssertionError:
>> Items are not equal:
>> item=1
>>
>> ACTUAL: 1.0
>> DESIRED: 0.0
>>
>> ----------------------------------------------------------------------
>>
>> The last error on pbdv is a quite old one.
>> Can someone reproduce this error?
>>
>
> Yes I can.
>
> ======================================================================
> FAIL: test_pbdv (test_basic.TestCephes)
> ----------------------------------------------------------------------
> Traceback (most recent call last):
> File
> "/home/nwagner/local/lib64/python2.6/site-packages/scipy/special/tests/test_basic.py",
> line 370, in test_pbdv
> assert_equal(cephes.pbdv(1,0),(0.0,0.0))
> File
> "/home/nwagner/local/lib64/python2.6/site-packages/numpy/testing/utils.py",
> line 183, in assert_equal
> assert_equal(actual[k], desired[k], 'item=%r\n%s' %
> (k,err_msg), verbose)
> File
> "/home/nwagner/local/lib64/python2.6/site-packages/numpy/testing/utils.py",
> line 190, in assert_equal
> raise AssertionError(msg)
> AssertionError:
> Items are not equal:
> item=1
>
> ACTUAL: 1.0
> DESIRED: 0.0
>
> ----------------------------------------------------------------------
> Ran 3542 tests in 80.817s
>
> FAILED (KNOWNFAIL=2, SKIP=17, errors=1, failures=6)
> <nose.result.TextTestResult run=3542 errors=1 failures=6>
>
> Cheers,
> Nils
> _______________________________________________
>
Ok.
According to mathematica:
ParabolicCylinderD[1., 0] : Same result with scipy. Ok.
ParabolicCylinderD[2., 0]=-1 : Same result with scipy. Ok.
D[ParabolicCylinderD[a, x], x]=1/2 x ParabolicCylinderD[a, x] -
ParabolicCylinderD[1 + a, x]
As a result, with a=1 and x=0, we have:
1/2 * 1 * 1 - (-1) = 0
It means that:
assert_equal(cephes.pbdv(1,0),(0.0,0.0)) is wrong.
It should be :
assert_equal(cephes.pbdv(1,0),(0.0,1.0))
It also looks like there is anyhow another problem:
I cannot reproduce the value of the derivative for "large" value of a and x:
pbdv(.2,.3)=(0.90436167932402323, 0.067568701127784084) : Ok.
but
pbdv(2,3)=(0.84319379649491466, 1.4228895315851684) : according to mathematica, it should be (0.84319379649491466,-0.632395)
Cheers,
Xavier
More information about the SciPy-Dev
mailing list