about summary refs log tree commit diff
path: root/src/math/logbl.c
diff options
context:
space:
mode:
authorSzabolcs Nagy <nsz@port70.net>2013-09-02 23:35:46 +0000
committerSzabolcs Nagy <nsz@port70.net>2013-09-05 11:30:07 +0000
commit98be442ee8a2b8b7e0802b604e384d5a2c43282e (patch)
treefe20ecbce7be9e457fd678a4e4cea6b0fe32e692 /src/math/logbl.c
parent4cec31fc23b28271a2ed0b290e65699dffac3fee (diff)
downloadmusl-98be442ee8a2b8b7e0802b604e384d5a2c43282e.tar.gz
musl-98be442ee8a2b8b7e0802b604e384d5a2c43282e.tar.xz
musl-98be442ee8a2b8b7e0802b604e384d5a2c43282e.zip
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
Diffstat (limited to 'src/math/logbl.c')
-rw-r--r--src/math/logbl.c4
1 files changed, 2 insertions, 2 deletions
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