[Scipy-tickets] [SciPy] #1442: expi not reliable on real axis (windows, mkl)

SciPy Trac scipy-tickets@scipy....
Fri May 20 23:24:35 CDT 2011


#1442: expi not reliable on real axis (windows, mkl)
------------------------------+---------------------------------------------
 Reporter:  charris           |       Owner:  somebody
     Type:  defect            |      Status:  new     
 Priority:  normal            |   Milestone:          
Component:  Other             |     Version:  0.9.0   
 Keywords:  windows mkl expi  |  
------------------------------+---------------------------------------------
 This is from Christolph Gohlke;


 {{{
 FAIL: test_mpmath.test_expi_complex
 ----------------------------------------------------------------------
 Traceback (most recent call last):
   File "C:\Python32\lib\site-packages\nose\case.py", line 188, in runTest
     self.test(*self.arg)
   File "C:\Python32\lib\site-packages\numpy\testing\decorators.py",
 line 147, in  skipper_func
     return f(*args, **kwargs)
   File
 "C:\Python32\lib\site-packages\scipy\special\tests\test_mpmath.py", line
 48, in test_expi_complex
     FuncData(sc.expi, dataset, 0, 1).check()
   File
 "C:\Python32\lib\site-packages\scipy\special\tests\testutils.py", line
 224, in check
     assert_(False, "\n".join(msg))
   File "C:\Python32\lib\site-packages\numpy\testing\utils.py", line 34,
 in assert_
     raise AssertionError(msg)
 AssertionError:
 Max |adiff|: 3.86856e+25
 Max |rdiff|: 0.304157
 Bad results for the following points (in output 0):
                     (1e-99+0j) =>
 (-227.37870854150898+6.283185307179586j) !=
     (-227.37870854150898+0j)  (rdiff            0.02763312953742352)
     (1.668100537200083e-88+0j) =>
 }}}

 The problem seems to be the use of the wrong branch in some complex
 logarithm evaluations. See the
 [http://article.gmane.org/gmane.comp.python.scientific.devel/15181 thread]
 on scipy-devel.

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1442>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list