[Scipy-tickets] [SciPy] #623: scipy.special.jv incorrect for certain argument(s)

SciPy scipy-tickets@scipy....
Mon Jan 19 13:03:35 CST 2009


#623: scipy.special.jv incorrect for certain argument(s)
-------------------------------+--------------------------------------------
 Reporter:  wa03               |        Owner:  cookedm 
     Type:  defect             |       Status:  assigned
 Priority:  normal             |    Milestone:  0.7.0   
Component:  scipy.special      |      Version:          
 Severity:  major              |   Resolution:          
 Keywords:  special.jv bessel  |  
-------------------------------+--------------------------------------------
Comment (by pv):

 Something else is rotten in `jv` than the value at (3,4):
 {{{
 >>> sc.jv(301, 1296.0681)
 -0,018079906032848522
 >>> sc.jv(301, 1296.0682)
 -0,025101462068653005
 }}}
 The asymptotic expansion seems to suffer from problems here, and the
 output from `jv(301,x)` seems to be completely wrong beyond x >= 1300.

 Specfun's JYNB may do a better job here, at least `jn_zeros(301,300)`
 (which calls JYNB down the line) computes the zeros in this region
 correctly.

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


More information about the Scipy-tickets mailing list