CORE-MATH ea7a421: log2.c, sinh.c, tanpi.c (on 81a899d's nine)

Upstream 10-12: sinh restores int64_t for ml (81a899d had narrowed it to int),
a log2 patch, and a tanpi change. Old (040ee48, main) against new in all four
rounding modes: the same results and flags on 686,469,600 calls on x86-64
(cfarm421) and aarch64 (cfarm424), with every input of the ten functions'
worst-case files and its negative (1,922,435), control differs. 81a899d against
ea7a421 is identical on the same calls: the narrowed ml cost time, not results.
All carried files equal ea7a421 (harness/cm-diff/carried-identical.py in
openpocl). Supersedes prep/cm-81a899d. For the next release.

Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>

Author: anun333

README.md

@@ -129,8 +129,8 @@ The scalar functions, their tables, and the error analyses the vector paths
 rely on are [CORE-MATH](https://core-math.gitlabpages.inria.fr/)'s, by Alexei
 Sibidanov, Paul Zimmermann, Tom Hubrecht and others. Their files are
 included unmodified under their own MIT license and copyright notices. All
-165 are byte-identical to CORE-MATH's master branch at `81a899d`
-(2026-10-11). The
+165 are byte-identical to CORE-MATH's master branch at `ea7a421`
+(2026-10-12). The
 `crmvec-*-tab.h` headers copy their tables. Everything else is under the MIT
 license in `LICENSE`.
 

log2.c

@@ -302,9 +302,10 @@ static double __attribute__((noinline)) as_log2_refine(double x, double a){
   double v2, v0 = fasttwosum(L[0], sh, &v2), v1 = fasttwosum(v2, sl, &v2);
   t.f = v1;
   if(__builtin_expect(!(t.u&(~(u64)0>>12)), 0)){
-    // even t and v2 are always of the same sign here, leave it aligned with other functions
+    // even if t and v2 are always of the same sign here, leave it aligned with other functions
     b64u64_u w = {.f = v2};
-    int d = ((int)((t.u^w.u)>>63))*2 + 1;
+    int d = 1 - 2*(int)((t.u^w.u)>>63);
+    // d=1 if t and v2 are of the same sign, d=-1 otherwise
     t.u += d;
     v1 = t.f;
   }

sinh.c

@@ -352,7 +352,7 @@ double cr_sinh(double x){
       // b64u64_u uh = {.f = th};
       b64u64_u ul = {.f = tl};
       // int64_t eh = (uh.u>>52)&0x7ff, el = (ul.u>>52)&0x7ff;
-      int ml = (ul.u + 8)&(~0ul>>12);
+      int64_t ml = (ul.u + 8)&(~0ul>>12);
       th += tl;
       th *= 2;
       th *= sp.f;

tanpi.c

@@ -142,7 +142,7 @@ static __attribute__((noinline)) double as_tanpi_database(double x, double f){
   uint64_t p = 0;
   double ax = __builtin_fabs(x);
   int32_t e = (aix>>52)-1022;
-  if(e>=0){ // take into account periodic structure of tan
+  if(e>=0){ // |x|>=0.5, take into account periodic structure of tan
     b64u64_u t = {.u = aix&~(~(uint64_t)0>>(12+e))};
     ax -= t.f;
     p = ((aix|(uint64_t)1<<52)>>(52-e))&1;