OLD | NEW |
1 #include "libm.h" | 1 #include "libm.h" |
2 | 2 |
3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 | 3 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 |
4 long double frexpl(long double x, int *e) | 4 long double frexpl(long double x, int* e) { |
5 { | 5 return frexp(x, e); |
6 » return frexp(x, e); | |
7 } | 6 } |
8 #elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 | 7 #elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 |
9 long double frexpl(long double x, int *e) | 8 long double frexpl(long double x, int* e) { |
10 { | 9 union ldshape u = {x}; |
11 » union ldshape u = {x}; | 10 int ee = u.i.se & 0x7fff; |
12 » int ee = u.i.se & 0x7fff; | |
13 | 11 |
14 » if (!ee) { | 12 if (!ee) { |
15 » » if (x) { | 13 if (x) { |
16 » » » x = frexpl(x*0x1p120, e); | 14 x = frexpl(x * 0x1p120, e); |
17 » » » *e -= 120; | 15 *e -= 120; |
18 » » } else *e = 0; | 16 } else |
19 » » return x; | 17 *e = 0; |
20 » } else if (ee == 0x7fff) { | 18 return x; |
21 » » return x; | 19 } else if (ee == 0x7fff) { |
22 » } | 20 return x; |
| 21 } |
23 | 22 |
24 » *e = ee - 0x3ffe; | 23 *e = ee - 0x3ffe; |
25 » u.i.se &= 0x8000; | 24 u.i.se &= 0x8000; |
26 » u.i.se |= 0x3ffe; | 25 u.i.se |= 0x3ffe; |
27 » return u.f; | 26 return u.f; |
28 } | 27 } |
29 #endif | 28 #endif |
OLD | NEW |