[Scipy-tickets] [SciPy] #640: special package numerical errors

SciPy scipy-tickets@scipy....
Tue Apr 15 17:51:04 CDT 2008


#640: special package numerical errors
---------------------------+------------------------------------------------
 Reporter:  veplaini       |        Owner:  cookedm 
     Type:  defect         |       Status:  assigned
 Priority:  normal         |    Milestone:  0.7     
Component:  scipy.special  |      Version:          
 Severity:  normal         |   Resolution:          
 Keywords:                 |  
---------------------------+------------------------------------------------
Comment (by veplaini):

 Thanks very much for your quick action.

 I can report that
 1) iv & ivp seem to work as expected for complex arguments:

 from scipy.special import *
 from pylab import *

 x = arange(600, 800,1)
 for i in range(0,100):
  s = []
  for j in x:
   s.append(ivp(i,complex(j,0))/(2*exp(j)))
  plot(x,s,linewidth=0.5)
 show()

 You can see it fails around 700, which is what
 I expect (same as in matlab)

 2) Using jv(n, complex(0, k)) seems to work for
 some orders but not all

 Thanks again for your help. I was *so glad* when
 I discovered scipy, because I was very reluctant
 to use matlab/octave, in my view, lousy programming
 languages. Python is so much better and with scipy
 I can do everything I need.

 Victor Lazzarini

-- 
Ticket URL: <http://scipy.org/scipy/scipy/ticket/640#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