diff options
Diffstat (limited to 'sysdeps/ieee754/dbl-64/e_gamma_r.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/e_gamma_r.c | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/sysdeps/ieee754/dbl-64/e_gamma_r.c b/sysdeps/ieee754/dbl-64/e_gamma_r.c index adeb61a248..d1acaa6a12 100644 --- a/sysdeps/ieee754/dbl-64/e_gamma_r.c +++ b/sysdeps/ieee754/dbl-64/e_gamma_r.c @@ -194,6 +194,11 @@ __ieee754_gamma_r (double x, int *signgamp) double tret = M_PI / (-x * sinpix * gamma_positive (-x, &exp2_adj)); ret = __scalbn (tret, -exp2_adj); + if (ret < DBL_MIN) + { + double force_underflow = ret * ret; + math_force_eval (force_underflow); + } } } } |