| /** |
| * This file has no copyright assigned and is placed in the Public Domain. |
| * This file is part of the mingw-w64 runtime package. |
| * No warranty is given; refer to the file DISCLAIMER.PD within this package. |
| */ |
| /* |
| * Written by J.T. Conklin <jtc@netbsd.org>. |
| * Changes for long double by Ulrich Drepper <drepper@cygnus.com> |
| * Public domain. |
| */ |
| |
| #include <math.h> |
| |
| float |
| logbf (float x) |
| { |
| #ifdef __x86_64__ |
| int v; |
| __mingw_flt_type_t hlp; |
| |
| hlp.x = x; |
| v = hlp.val & 0x7fffffff; /* high |x| */ |
| if (!v) |
| return (float)-1.0 / fabsf (x); |
| if (v >= 0x7f800000) |
| return x * x; |
| if ((v >>= 23) == 0) /* IEEE 754 logb */ |
| return -126.0; |
| return (float) (v - 127); |
| #else |
| float res = 0.0F; |
| asm ("fxtract\n\t" |
| "fstp %%st" : "=t" (res) : "0" (x)); |
| return res; |
| #endif |
| } |