math: fix logb(-0.0) in downward rounding mode

use -1/(x*x) instead of -1/(x+0) to return -inf, -0+0 is -0 in
downward rounding mode

diff --git a/src/math/logb.c b/src/math/logb.c
index 624425a..7f8bdfa 100644 (file)
--- a/src/math/logb.c
+++ b/src/math/logb.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
 
 /*
 special cases:
@@ -12,6 +12,6 @@ double logb(double x)
        if (!isfinite(x))
                return x * x;
        if (x == 0)
-               return -1/(x+0);
+               return -1/(x*x);
        return ilogb(x);
 }

diff --git a/src/math/logbf.c b/src/math/logbf.c
index 950d356..a0a0b5e 100644 (file)
--- a/src/math/logbf.c
+++ b/src/math/logbf.c
@@ -1,10 +1,10 @@
-#include "libm.h"
+#include <math.h>
 
 float logbf(float x)
 {
        if (!isfinite(x))
                return x * x;
        if (x == 0)
-               return -1/(x+0);
+               return -1/(x*x);
        return ilogbf(x);
 }

diff --git a/src/math/logbl.c b/src/math/logbl.c
index f385074..962973a 100644 (file)
--- a/src/math/logbl.c
+++ b/src/math/logbl.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
 long double logbl(long double x)
 {
@@ -10,7 +10,7 @@ long double logbl(long double x)
        if (!isfinite(x))
                return x * x;
        if (x == 0)
-               return -1/(x+0);
+               return -1/(x*x);
        return ilogbl(x);
 }
 #endif