[SciPy-dev] Status on ubuntu jaunty 64bits.

Xavier Gnata xavier.gnata@gmail....
Tue Mar 17 18:03:33 CDT 2009


Nils Wagner wrote:
> On Tue, 17 Mar 2009 21:21:46 +0100
>   Xavier Gnata <xavier.gnata@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