diff options
Diffstat (limited to 'sysdeps/libm-ieee754/s_cbrtl.c')
| -rw-r--r-- | sysdeps/libm-ieee754/s_cbrtl.c | 6 |
1 files changed, 4 insertions, 2 deletions
diff --git a/sysdeps/libm-ieee754/s_cbrtl.c b/sysdeps/libm-ieee754/s_cbrtl.c index b3a53a3..1d021b7 100644 --- a/sysdeps/libm-ieee754/s_cbrtl.c +++ b/sysdeps/libm-ieee754/s_cbrtl.c @@ -48,8 +48,10 @@ __cbrtl (long double x) xm = __frexpl (fabs (x), &xe); /* If X is not finite or is null return it (with raising exceptions - if necessary. */ - if (xe == 0) + if necessary. + Note: *Our* version of `frexp' sets XE to zero if the argument is + Inf or NaN. This is not portable but faster. */ + if (xe == 0 && fpclassify (x) <= FP_ZERO) return x + x; u = (0.338058687610520237 |