Spe math fixes (#379)

* Wrong Printf version requirement printing

* Set null to inf when exiting from thread

* Fix on pthread_get_specific if called from main thread

* Fixed critical (and minor) bugs in pthread library. At moment all test program seems to work correctly

* Some fixes for libresolv. Added DNS cache to requests

* Some fixes for libresolv. Added DNS cache to requests

* stdlib: fix multiple bugs and correctness issues

* stdio: fix buffer overwrite, error handling, and overflow issues

* posix: fix critical bugs in IPC, signals, mmap, and statvfs

* socket: fix multiple bugs found

* Fix multiple bugs across clib4 library

* math: fixed 8 bugs (2 general + 6 SPE-specific)

Author: afxgroup

library/c.lib_rev.h

@@ -2,7 +2,7 @@
 #define REVISION		1
 #define SUBREVISION		0
 
-#define DATE			"06.03.2026"
+#define DATE			"09.03.2026"
 #define VERS			"clib4.library 2.1"
-#define VSTRING			"clib4.library 2.1 (06.03.2026)\r\n"
-#define VERSTAG			"\0$VER: clib4.library 2.1-ab40f98 (06.03.2026)"
+#define VSTRING			"clib4.library 2.1 (09.03.2026)\r\n"
+#define VERSTAG			"\0$VER: clib4.library 2.1-ab40f98 (09.03.2026)"

library/math/e_sqrt.c

@@ -11,8 +11,6 @@ static	const float	one	= 1.0, tiny = 1.0e-30;
 double
 __ieee754_sqrt(double x) {
     double z;
-
-#ifndef __SPE__
     int32_t sign = (int32_t) 0x80000000;
     int32_t ix0, s0, q, m, t, i;
     uint32_t r, t1, s1, ix1, q1;
@@ -109,6 +107,5 @@ __ieee754_sqrt(double x) {
     if ((q & 1) == 1) ix1 |= sign;
     ix0 += (m << 20);
     INSERT_WORDS(z, ix0, ix1);
-#endif
     return z;
 }

library/math/fpsetmask.c

@@ -19,10 +19,10 @@ fpsetmask(fp_except_t mask) {
     return (old);
 #else
     uint32_t fpscr;
-	fp_rnd_t old;
+	fp_except_t old;
 
 	__asm__ __volatile("mfspr %0, %1" : "=r"(fpscr) : "K"(SPR_SPEFSCR));
-	old = (fp_rnd_t)((fpscr >> 3) & 0x1f);
+	old = (fp_except_t)((fpscr >> 3) & 0x1f);
 	fpscr = (fpscr & 0xffffff07) | (mask << 3);
 	__asm__ __volatile("mtspr %1,%0" :: "r"(fpscr), "K"(SPR_SPEFSCR));
 	return (old);

library/math/s_fpclassify.c

@@ -71,7 +71,7 @@ __fpclassify_double(double d) {
         /* Both exponent and fraction are zero -> zero */
         result = FP_ZERO;
     }
-    else if ((x.raw[0] & 0x7fff0000) == 0) {
+    else if ((x.raw[0] & 0x7ff00000) == 0) {
         SHOWMSG("subnormal");
 
         /* Exponent = 0 -> subnormal (IEEE 754) */

library/math/s_isinf.c

@@ -20,11 +20,10 @@ __isinf(double d) {
 #else
     union ieee_double x;
 
-	/* Exponent = 2047 and fraction = 0.0 -> infinity */
-	x.raw[0] = 0x7ff00000;
-	x.raw[1] = 0x00000000;
+	x.value = d;
 
-	return(x.value);
+	/* Exponent = 2047 and fraction = 0.0 -> infinity */
+	return (((x.raw[0] & 0x7fffffff) == 0x7ff00000) && (x.raw[1] == 0));
 #endif
 }
 
@@ -39,7 +38,7 @@ __isinff(float f) {
     int32_t ix;
     GET_FLOAT_WORD(ix, f);
     ix &= 0x7fffffff;
-    return FLT_UWORD_IS_NAN(ix);
+    return (ix == 0x7f800000);
 #endif
 }
 

library/math/s_isnan.c

@@ -16,11 +16,10 @@ __isnan(double d) {
 #else
     union ieee_double x;
 
-	/* Exponent = 2047 and fraction != 0.0; this must be a quiet nan. */
-	x.raw[0] = 0x7ff80000;
-	x.raw[1] = 0x00000001;
+	x.value = d;
 
-	return x.value;
+	/* Exponent = 2047 and fraction != 0.0 -> NaN */
+	return (((x.raw[0] & 0x7ff00000) == 0x7ff00000) && ((x.raw[0] & 0x000fffff) != 0 || (x.raw[1] != 0)));
 #endif
 }
 

library/math/s_ldexp.c

@@ -15,7 +15,7 @@ ldexp(double x, int exp) {
     } else {
         result = scalbn(x, exp);
 
-        if (isinf(result) || (result < DBL_MIN || result > -DBL_MIN))
+        if (isinf(result) || (result < DBL_MIN && result > -DBL_MIN))
             __set_errno(ERANGE);
     }
 

library/math/s_ldexpf.c

@@ -15,7 +15,7 @@ ldexpf(float x, int exp) {
     } else {
         result = scalbnf(x, exp);
 
-        if (isinf(result) || (result < FLT_MIN || result > -FLT_MIN))
+        if (isinf(result) || (result < FLT_MIN && result > -FLT_MIN))
             __set_errno(ERANGE);
     }