mirrors/cosmopolitan
1
0
Fork 0
mirror of https://github.com/jart/cosmopolitan.git synced 2025-01-31 03:27:39 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

Import optimized routines changes to exp10

Browse source
This commit is contained in:
Justine Tunney 2024-08-15 18:33:40 -07:00
parent 2045e87b7c
commit 3fd275f59f
No known key found for this signature in database
GPG key ID: BE714B4575D6E328
1 changed files with 3 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
libc/tinymath/exp10.c
View file

@ -43,7 +43,7 @@ special_case (uint64_t sbits, double_t tmp, uint64_t ki)
{ {
double_t scale, y; double_t scale, y;
if (ki - (1ull << 16) < 0x80000000) if ((ki & 0x80000000) == 0)
{ {
/* The exponent of scale might have overflowed by 1. */ /* The exponent of scale might have overflowed by 1. */
sbits -= 1ull << 52; sbits -= 1ull << 52;
@ -109,14 +109,14 @@ exp10 (double x)
/* Reduce x: z = x * N / log10(2), k = round(z). */ /* Reduce x: z = x * N / log10(2), k = round(z). */
double_t z = __exp_data.invlog10_2N * x; double_t z = __exp_data.invlog10_2N * x;
double_t kd; double_t kd;
int64_t ki; uint64_t ki;
#if TOINT_INTRINSICS #if TOINT_INTRINSICS
kd = roundtoint (z); kd = roundtoint (z);
ki = converttoint (z); ki = converttoint (z);
#else #else
kd = eval_as_double (z + Shift); kd = eval_as_double (z + Shift);
ki = asuint64 (kd);
kd -= Shift; kd -= Shift;
ki = kd;
#endif #endif
/* r = x - k * log10(2), r in [-0.5, 0.5]. */ /* r = x - k * log10(2), r in [-0.5, 0.5]. */