OLD | NEW |
1 #define _GNU_SOURCE | 1 #define _GNU_SOURCE |
2 #include <float.h> | 2 #include <float.h> |
3 #include <math.h> | 3 #include <math.h> |
4 #include "libc.h" | 4 #include "libc.h" |
5 #include "libm.h" | 5 #include "libm.h" |
6 | 6 |
7 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 | 7 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024 |
8 long double exp10l(long double x) | 8 long double exp10l(long double x) { |
9 { | 9 return exp10(x); |
10 » return exp10(x); | |
11 } | 10 } |
12 #elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 | 11 #elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384 |
13 long double exp10l(long double x) | 12 long double exp10l(long double x) { |
14 { | 13 static const long double p10[] = { |
15 » static const long double p10[] = { | 14 1e-15L, 1e-14L, 1e-13L, 1e-12L, 1e-11L, 1e-10L, 1e-9L, 1e-8L, |
16 » » 1e-15L, 1e-14L, 1e-13L, 1e-12L, 1e-11L, 1e-10L, | 15 1e-7L, 1e-6L, 1e-5L, 1e-4L, 1e-3L, 1e-2L, 1e-1L, 1, |
17 » » 1e-9L, 1e-8L, 1e-7L, 1e-6L, 1e-5L, 1e-4L, 1e-3L, 1e-2L, 1e-1L, | 16 1e1, 1e2, 1e3, 1e4, 1e5, 1e6, 1e7, 1e8, |
18 » » 1, 1e1, 1e2, 1e3, 1e4, 1e5, 1e6, 1e7, 1e8, 1e9, | 17 1e9, 1e10, 1e11, 1e12, 1e13, 1e14, 1e15}; |
19 » » 1e10, 1e11, 1e12, 1e13, 1e14, 1e15 | 18 long double n, y = modfl(x, &n); |
20 » }; | 19 union ldshape u = {n}; |
21 » long double n, y = modfl(x, &n); | 20 /* fabsl(n) < 16 without raising invalid on nan */ |
22 » union ldshape u = {n}; | 21 if ((u.i.se & 0x7fff) < 0x3fff + 4) { |
23 » /* fabsl(n) < 16 without raising invalid on nan */ | 22 if (!y) |
24 » if ((u.i.se & 0x7fff) < 0x3fff+4) { | 23 return p10[(int)n + 15]; |
25 » » if (!y) return p10[(int)n+15]; | 24 y = exp2l(3.32192809488736234787031942948939L * y); |
26 » » y = exp2l(3.32192809488736234787031942948939L * y); | 25 return y * p10[(int)n + 15]; |
27 » » return y * p10[(int)n+15]; | 26 } |
28 » } | 27 return powl(10.0, x); |
29 » return powl(10.0, x); | |
30 } | 28 } |
31 #endif | 29 #endif |
32 | 30 |
33 weak_alias(exp10l, pow10l); | 31 weak_alias(exp10l, pow10l); |
OLD | NEW |