diff options
| author | Joseph Myers <joseph@codesourcery.com> | 2013-12-03 16:25:18 +0000 |
|---|---|---|
| committer | Joseph Myers <joseph@codesourcery.com> | 2013-12-03 16:25:18 +0000 |
| commit | 34e16df5a1a46e128edb9eb44a09ac5762957136 (patch) | |
| tree | d8751230a2270d89531af1f7c12868dae50f534f /sysdeps/ieee754/dbl-64/s_erf.c | |
| parent | d8e2dbe3e380729a1552d546da582b02202dde0a (diff) | |
| download | glibc-34e16df5a1a46e128edb9eb44a09ac5762957136.tar.gz | |
Fix erfc errno setting on underflow (bug 6786).
Diffstat (limited to 'sysdeps/ieee754/dbl-64/s_erf.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/s_erf.c | 17 |
1 files changed, 15 insertions, 2 deletions
diff --git a/sysdeps/ieee754/dbl-64/s_erf.c b/sysdeps/ieee754/dbl-64/s_erf.c index aab4ed593e..3f37397224 100644 --- a/sysdeps/ieee754/dbl-64/s_erf.c +++ b/sysdeps/ieee754/dbl-64/s_erf.c @@ -112,6 +112,8 @@ static char rcsid[] = "$NetBSD: s_erf.c,v 1.8 1995/05/10 20:47:05 jtc Exp $"; */ +#include <errno.h> +#include <float.h> #include <math.h> #include <math_private.h> @@ -391,14 +393,25 @@ __erfc (double x) r = __ieee754_exp (-z * z - 0.5625) * __ieee754_exp ((z - x) * (z + x) + R / S); if (hx > 0) - return r / x; + { +#if FLT_EVAL_METHOD != 0 + volatile +#endif + double ret = r / x; + if (ret == 0) + __set_errno (ERANGE); + return ret; + } else return two - r / x; } else { if (hx > 0) - return tiny * tiny; + { + __set_errno (ERANGE); + return tiny * tiny; + } else return two - tiny; } |