diff options
Diffstat (limited to 'src/math/log10.c')
| -rw-r--r-- | src/math/log10.c | 13 |
1 files changed, 13 insertions, 0 deletions
diff --git a/src/math/log10.c b/src/math/log10.c index cea1d45b..e504df44 100644 --- a/src/math/log10.c +++ b/src/math/log10.c @@ -2,11 +2,24 @@ #include <math.h> #include "_tgmath.h" #include "errno.h" +#include "fenv.h" /** base-10 logarithm **/ TYPE TGFN(log10)(TYPE x) { + switch (fpclassify(x)) { + case FP_ZERO: feraiseexcept(FE_INVALID); return - INFINITY; + case FP_INFINITE: if (!signbit(x)) { return x; } break; + default: break; + } + + if (x == 1.0) { + return 0.0; + } + if (x < 0) { + feraiseexcept(FE_INVALID); + return NAN; errno = EDOM; /* ARGUMENT(x) is negative */ return TGHUGE; } |
