[Scipy-tickets] [SciPy] #1802: Infinite time taken by scipy hypergeom function for some particular values

SciPy Trac scipy-tickets@scipy....
Thu Jan 3 10:31:57 CST 2013


#1802: Infinite time taken by scipy hypergeom function for some particular values
-------------------------+--------------------------------------------------
 Reporter:  imsc         |       Owner:  somebody    
     Type:  defect       |      Status:  needs_review
 Priority:  normal       |   Milestone:  Unscheduled 
Component:  scipy.stats  |     Version:  0.11.0      
 Keywords:               |  
-------------------------+--------------------------------------------------

Comment(by josefpktd):

 There is still a problem with the other degenerate case

 {{{
 >>> spd.hypergeom.cdf(0, 32335, 0, 895)  #ok
 1.0
 >>> spd.hypergeom.pmf(0, 32335, 0, 895)
 nan
 >>> spd.hypergeom.pmf(1, 32335, 0, 895)
 nan
 >>> spd.hypergeom.cdf(-1, 32335, 0, 895)
 nan
 >>> spd.hypergeom.pmf(-1, 32335, 0, 895)
 nan
 >>> spd.hypergeom.ppf(0.02, 32335, 0, 895)
 nan
 >>> spd.hypergeom.ppf(0.0, 32335, 0, 895)
 -1.0
 >>> spd.hypergeom.ppf(1.0, 32335, 0, 895)
 nan
 >>> spd.hypergeom.a
 0
 >>> spd.hypergeom.b
 array(0)
 }}}

 problem was also in older scipy, (spd is scipy.stats.distributions)

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


More information about the Scipy-tickets mailing list