EXPM1MINUS evaluation stubs
PiperOrigin-RevId: 344736247
diff --git a/src/math/expm1minus-scalar-rr2-lut16-p4.c b/src/math/expm1minus-scalar-rr2-lut16-p4.c
index dcc6c81..11a9834 100644
--- a/src/math/expm1minus-scalar-rr2-lut16-p4.c
+++ b/src/math/expm1minus-scalar-rr2-lut16-p4.c
@@ -29,9 +29,9 @@
const uint32_t vindex_mask = UINT32_C(0xF);
// The largest x for which expm1f(x) is saturated at -1.0f.
const float vsat_cutoff = -0x1.154246p+4f;
- // Last 11 bits are zeroes
- const float vminus_ln2_hi = -0x1.62E000p-1f;
- const float vminus_ln2_lo = -0x1.0BFBE8p-15f;
+ // Last 9 bits are zeroes
+ const float vminus_ln2_hi = -0x1.62E400p-1f;
+ const float vminus_ln2_lo = -0x1.7F7D1Cp-20f;
// Coefficient of polynomial approximation
// exp(t) - 1 ~ t * (1 + t * (c2 + t * (c3 + t * c4)))
// on [-log(2)/32, log(2)/32]
@@ -88,7 +88,7 @@
vp *= vt;
// Reconstruct the exp(x) - 1 value:
- // exp(x) - 1 = s * (1 + t * (1 + t * (c2 + t * c3))) - 1
+ // exp(x) - 1 = s * (1 + t * (1 + t * (c2 + t * (c3 + t * c4)))) - 1
// = (s - 1) + s * (t + t * p)
// = ((t * s) + (t * s) * p) + (s - 1)
vt *= vs;