[Scipy-tickets] [SciPy] #1334: Spurious "divide by zero" warnings from orthogonal polynomials

SciPy Trac scipy-tickets@scipy....
Fri Nov 26 22:16:33 CST 2010


#1334: Spurious "divide by zero" warnings from orthogonal polynomials
------------------------------+---------------------------------------------
 Reporter:  warren.weckesser  |       Owner:  pv   
     Type:  defect            |      Status:  new  
 Priority:  normal            |   Milestone:  0.9.0
Component:  scipy.special     |     Version:  0.8.0
 Keywords:  orthogonal        |  
------------------------------+---------------------------------------------
 The following script:
 {{{
 from scipy.special.orthogonal import legendre, laguerre, gegenbauer

 print "\ncalling legendre(1)(0)"
 y = legendre(1)(0)

 print "\ncalling laguerre(1)(1)"
 y = laguerre(1)(1)

 print "\ncalling gegenbauer(1,1)(0)"
 y = gegenbauer(1,1)(0)
 }}}
 produces the following output:
 {{{
 calling legendre(1)(0)
 Warning: divide by zero encountered in hyp2f1

 calling laguerre(1)(1)
 Warning: divide by zero encountered in hyp1f1

 calling gegenbauer(1,1)(0)
 Warning: divide by zero encountered in hyp2f1
 }}}
 (using r6945).  Evaluating a polynomial should not produce such a warning.

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


More information about the Scipy-tickets mailing list