[Scipy-tickets] [SciPy] #1589: special.kv goes from 0 to nan for large numbers

SciPy Trac scipy-tickets@scipy....
Mon Jan 23 23:24:36 CST 2012


#1589: special.kv goes from 0 to nan for large numbers
---------------------------+------------------------------------------------
 Reporter:  johannct       |       Owner:  pv         
     Type:  defect         |      Status:  new        
 Priority:  normal         |   Milestone:  Unscheduled
Component:  scipy.special  |     Version:  0.10.0     
 Keywords:  bessel         |  
---------------------------+------------------------------------------------
 In [1]: import scipy

 In [2]: scipy.__version__
 Out[2]: '0.10.0'

 In [3]: import scipy.special as sp

 In [4]: sp.kv(1./3.,1.08e+09)
 Out[4]: nan

 In [5]: sp.kv(1./3.,1.07e+09)
 Out[5]: 0.0

 The scipy implementation apparently comes from Boost :
 http://www.boost.org/doc/libs/1_44_0/libs/math/doc/sf_and_dist/html/math_toolkit/special/bessel/mbessel.html

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


More information about the Scipy-tickets mailing list