[Scipy-tickets] [SciPy] #1750: hyp0f1(v, z) is nan for many (all?) v combined with many (all?) negative float z

SciPy Trac scipy-tickets@scipy....
Tue Oct 30 16:50:53 CDT 2012


#1750: hyp0f1(v, z) is nan for many (all?) v combined with many (all?) negative
float z
------------------------------+---------------------------------------------
 Reporter:  unknown (google)  |       Owner:  pv          
     Type:  defect            |      Status:  needs_review
 Priority:  normal            |   Milestone:  Unscheduled 
Component:  scipy.special     |     Version:  devel       
 Keywords:  hyp0f1            |  
------------------------------+---------------------------------------------
Changes (by rgommers):

  * priority:  high => normal
  * status:  new => needs_review
  * version:  0.10.1 => devel


Old description:

> >>> scipy.special.hyp0f1(3.0, numpy.linspace(-3, 3, 21))
> array([        nan,         nan,         nan,         nan,         nan,
>                nan,         nan,         nan,         nan,         nan,
>         1.        ,  1.10382595,  1.21561526,  1.33585293,  1.46504844,
>         1.60373685,  1.7524799 ,  1.9118671 ,  2.082517  ,  2.26507831,
>         2.46023128])
>

>
> >>> scipy.special.hyp0f1(0.5, numpy.linspace(-3, 3, 21))
> array([         nan,          nan,          nan,          nan,
>                 nan,          nan,          nan,          nan,
>                 nan,          nan,   1.        ,   1.66245212,
>          2.4600452 ,   3.40913698,   4.52749411,   5.83438641,
>          7.35068624,   9.09897326,  11.10364481,  13.39103224,
> 15.98952331])
>

> You can work around this by setting adding 0j to the z argument to force
> the typechecking in special/basic.py to use more general code.
>
> I'm setting the priority to 'high' to match
> http://projects.scipy.org/scipy/ticket/1034

New description:

 >>> scipy.special.hyp0f1(3.0, numpy.linspace(-3, 3, 21))
 array([        nan,         nan,         nan,         nan,         nan,
                nan,         nan,         nan,         nan,         nan,
         1.        ,  1.10382595,  1.21561526,  1.33585293,  1.46504844,
         1.60373685,  1.7524799 ,  1.9118671 ,  2.082517  ,  2.26507831,
         2.46023128])



 >>> scipy.special.hyp0f1(0.5, numpy.linspace(-3, 3, 21))
 array([         nan,          nan,          nan,          nan,
                 nan,          nan,          nan,          nan,
                 nan,          nan,   1.        ,   1.66245212,
          2.4600452 ,   3.40913698,   4.52749411,   5.83438641,
          7.35068624,   9.09897326,  11.10364481,  13.39103224,
 15.98952331])


 You can work around this by setting adding 0j to the z argument to force
 the typechecking in special/basic.py to use more general code.

 I'm setting the priority to 'high' to match
 http://projects.scipy.org/scipy/ticket/1034

--

Comment:

 Thanks, fixed by https://github.com/scipy/scipy/pull/345

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


More information about the Scipy-tickets mailing list