diff options
Diffstat (limited to 'src/math/fabs.c')
| -rw-r--r-- | src/math/fabs.c | 8 |
1 files changed, 7 insertions, 1 deletions
diff --git a/src/math/fabs.c b/src/math/fabs.c index 9664b702..5fc42750 100644 --- a/src/math/fabs.c +++ b/src/math/fabs.c @@ -6,6 +6,12 @@ /** absolute value **/ TYPE TGFN(fabs)(TYPE x) { + switch (fpclassify(x)) { + case FP_ZERO: return 0.0; + case FP_INFINITE: return INFINITY; + default: break; + } + if (0) { errno = ERANGE; /* The result cannot be represented */ /* RETURN_FAILURE(CONSTANT(HUGE_VAL), A range error occurred); */ @@ -13,7 +19,7 @@ TYPE TGFN(fabs)(TYPE x) } /* RETURN_SUCCESS(ABS(ARGUMENT(x))); */ - return x < 0.0 ? -x : x; + return signbit(x) ? -x : x; } /*** |