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/logb.c | |
| parent | 4cec31fc23b28271a2ed0b290e65699dffac3fee (diff) | |
| download | musl-98be442ee8a2b8b7e0802b604e384d5a2c43282e.tar.gz | |
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/logb.c')
| -rw-r--r-- | src/math/logb.c | 4 |
1 files changed, 2 insertions, 2 deletions
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); } |