[Scipy-tickets] [SciPy] #1114: scipy.special.sph_jn produces NaN result for large z

SciPy Trac scipy-tickets@scipy....
Tue Feb 16 02:58:02 CST 2010


#1114: scipy.special.sph_jn produces NaN result for large z
---------------------------+------------------------------------------------
 Reporter:  juya           |       Owner:  pv   
     Type:  defect         |      Status:  new  
 Priority:  normal         |   Milestone:  0.8.0
Component:  scipy.special  |     Version:  0.7.0
 Keywords:  sph_jn, NaN    |  
---------------------------+------------------------------------------------
 An example:

 from scipy.special import sph_jn
 print sph_jn(2,10000)

 The result is:
 (array([ NaN,  NaN,  NaN]), array([ -9.52124807e-05,              NaN,
 NaN]))

 with the analytical expression of j_2(z) = (3/z**3 - 1/z)*sinz -
 3/z**2*cos(z), it should not give NaN result.

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


More information about the Scipy-tickets mailing list