1 diff --git a/math_private.h b/math_private.h
7 * Emulate FreeBSD internal double types.
8 * Adapted from https://github.com/freebsd/freebsd-src/search?q=__double_t
11 typedef double __double_t;
12 typedef __double_t double_t;
13 +typedef float __float_t;
16 * The original fdlibm code used statements like:
17 * n0 = ((*(int*)&one)>>29)^1; * index of high word *
18 * ix0 = *(n0+(int*)&x); * high word of x *
19 * ix1 = *((1-n0)+(int*)&x); * low word of x *
20 * to dig two 32 bit words out of the 64 bit IEEE floating point
21 * value. That is non-ANSI, and, moreover, the gcc instruction
22 diff --git a/e_rem_pio2f.cpp b/e_rem_pio2f.cpp
25 @@ -49,17 +49,17 @@ int
32 /* 33+53 bit pi is good enough for medium size */
33 if(ix<0x4dc90fdb) { /* |x| ~< 2^28*(pi/2), medium size */
34 - fn = rnint((float_t)x*invpio2);
35 + fn = rnint((float)x*invpio2);
43 * all other (large) arguments
44 diff --git a/s_scalbnf.cpp b/s_scalbnf.cpp
48 * for all contributors to musl.
52 #include "math_private.h"
54 float scalbnf(float x, int n)
56 union {float f; uint32_t i;} u;