changeset 19651:520ee3fec8eb stable

update Faddeeva.cc to latest upstream version
author John W. Eaton <jwe@octave.org>
date Fri, 23 Jan 2015 15:20:00 -0500
parents 7005044f3080
children 75a671fcdd73
files liboctave/cruft/Faddeeva/Faddeeva.cc
diffstat 1 files changed, 2 insertions(+), 1 deletions(-) [+]
line wrap: on
line diff
--- a/liboctave/cruft/Faddeeva/Faddeeva.cc	Wed Jan 21 20:32:40 2015 -0800
+++ b/liboctave/cruft/Faddeeva/Faddeeva.cc	Fri Jan 23 15:20:00 2015 -0500
@@ -139,6 +139,7 @@
                        slight accuracy improvements to erf and dawson
                        functions near the origin.  Use gnulib functions
                        if GNULIB_NAMESPACE is defined.
+     18 December 2012: Slight tweaks (remove recomputation of x*x in Dawson)
 */
 
 /////////////////////////////////////////////////////////////////////////
@@ -584,7 +585,7 @@
     }
     else {
       double D = spi2 * FADDEEVA(w_im)(x);
-      double x2 = x*x, y2 = y*y;
+      double y2 = y*y;
       return C
         (D + y2 * (D + x - 2*D*x2)
          + y2*y2 * (D * (0.5 - x2 * (2 - 0.66666666666666666667*x2))