| 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 |