diff options
| author | Szabolcs Nagy <nsz@port70.net> | 2013-09-02 23:35:46 +0000 |
|---|---|---|
| committer | Szabolcs Nagy <nsz@port70.net> | 2013-09-05 11:30:07 +0000 |
| commit | 98be442ee8a2b8b7e0802b604e384d5a2c43282e (patch) | |
| tree | fe20ecbce7be9e457fd678a4e4cea6b0fe32e692 /src/math/logbl.c | |
| parent | 4cec31fc23b28271a2ed0b290e65699dffac3fee (diff) | |
| download | musl-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.c | 4 |
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 |