[Scipy-tickets] [SciPy] #1757: new upstream of complex erf code (special.wofz) - performance, accuracy improvements

SciPy Trac scipy-tickets@scipy....
Tue Oct 30 19:15:31 CDT 2012


#1757: new upstream of complex erf code (special.wofz) - performance, accuracy
improvements
---------------------------+------------------------------------------------
 Reporter:  stevenj        |       Owner:  pv         
     Type:  defect         |      Status:  new        
 Priority:  normal         |   Milestone:  Unscheduled
Component:  scipy.special  |     Version:  devel      
 Keywords:                 |  
---------------------------+------------------------------------------------

Comment(by stevenj):

 Note that the latest upstream version also incorporates special-case code
 for the real-z axis, where the imaginary part of w is Dawson's integral.
 In my benchmarks (g++ on x86-64), it seems a bit faster than the CEPHES
 library code (scipy/special/cephes/dawsn.c) that SciPy is currently using
 for the special.dawsn function, while attaining similar accuracy.  So, if
 you wanted, you could probably replace the special.dawsn implementation
 with something like:

 {{{
 double Dawson(double x) {
   return 0.88622692545275801364908374167057 // sqrt(pi)/2
     * ImFaddeeva_w(x);
 }
 }}}

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


More information about the Scipy-tickets mailing list