From 98be442ee8a2b8b7e0802b604e384d5a2c43282e Mon Sep 17 00:00:00 2001
From: Szabolcs Nagy <nsz@port70.net>
Date: Mon, 2 Sep 2013 23:35:46 +0000
Subject: 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
---
 src/math/logbl.c | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

(limited to 'src/math/logbl.c')

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
-- 
cgit v1.2.1