From: Szabolcs Nagy <nsz@port70.net>
Date: Mon, 2 Sep 2013 23:35:46 +0000 (+0000)
Subject: math: fix logb(-0.0) in downward rounding mode
X-Git-Tag: v0.9.14~42
X-Git-Url: https://git.librecmc.org/?a=commitdiff_plain;h=98be442ee8a2b8b7e0802b604e384d5a2c43282e;p=oweals%2Fmusl.git

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 624425a8..7f8bdfae 100644
--- 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 950d3569..a0a0b5ed 100644
--- 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 f3850745..962973a7 100644
--- 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