diff options
| author | Rich Felker <dalias@aerifal.cx> | 2012-11-15 08:37:58 -0500 |
|---|---|---|
| committer | Rich Felker <dalias@aerifal.cx> | 2012-11-15 08:37:58 -0500 |
| commit | 8bb181622222f2ee3462c8b021bcae4fcdbbd37a (patch) | |
| tree | 166a1be834c8cd5a4c2521d41d7ee506bd8c8e51 /src/math/logbf.c | |
| parent | update copyright file for recent contributions (diff) | |
| parent | fenv: return FE_TONEAREST in dummy fegetround (diff) | |
| download | grovel-8bb181622222f2ee3462c8b021bcae4fcdbbd37a.tar.gz grovel-8bb181622222f2ee3462c8b021bcae4fcdbbd37a.tar.xz | |
Merge remote-tracking branch 'nsz/math'
Diffstat (limited to 'src/math/logbf.c')
| -rw-r--r-- | src/math/logbf.c | 14 |
1 files changed, 6 insertions, 8 deletions
diff --git a/src/math/logbf.c b/src/math/logbf.c index 934827f8..950d3569 100644 --- a/src/math/logbf.c +++ b/src/math/logbf.c @@ -1,12 +1,10 @@ -#include <limits.h> #include "libm.h" -float logbf(float x) { - int i = ilogbf(x); - - if (i == FP_ILOGB0) - return -1.0f/fabsf(x); - if (i == FP_ILOGBNAN || i == INT_MAX) +float logbf(float x) +{ + if (!isfinite(x)) return x * x; - return i; + if (x == 0) + return -1/(x+0); + return ilogbf(x); } |