add subnormal strtod tests

diff --git a/src/stdlib/strtod.c b/src/stdlib/strtod.c
index a4a8b53..ff85c6a 100644 (file)
--- a/src/stdlib/strtod.c
+++ b/src/stdlib/strtod.c
@@ -14,7 +14,7 @@
        ((r) = (f)) == (x) || \
        (error("%s failed (" m ")\n", #f, r, x, r-x), 0) )
 
-void test_strtod(void) {
+void test_strtod_simple(void) {
        int i;
        double d, d2;
        char buf[1000];
@@ -28,3 +28,104 @@ void test_strtod(void) {
        TEST(d, strtod("0x1p4", 0), 16.0, "hex float %a != %a");
        TEST(d, strtod("0x1.1p4", 0), 17.0, "hex float %a != %a");
 }
+
+#define length(x) (sizeof(x) / sizeof(*(x)))
+
+/* TODO: float exceptions, rounding mode, endptr check */
+
+static struct {
+       char *s;
+       long double f;
+} tl[] = {
+       {"12.345", 12.345L},
+       {"1.2345e1", 12.345L},
+       // 2^-16445 * 0.5 - eps
+       {".1822599765941237301264202966809709908199525407846781671860490243514185844316698e-4950", 0},
+       // 2^-16445 * 0.5 + eps
+       {".1822599765941237301264202966809709908199525407846781671860490243514185844316699e-4950", 0x1p-16445L},
+       // 2^-16445 * 1.5 - eps
+       {".5467799297823711903792608900429129724598576223540345015581470730542557532950096e-4950", 0x1p-16445L},
+       // 2^-16445 * 1.5 + eps
+       {".5467799297823711903792608900429129724598576223540345015581470730542557532950097e-4950", 0x1p-16444L},
+       // 2^-16382 + 2^-16446 - eps
+       {".3362103143112093506444937793915876332724499641527442230928779770593420866576777e-4931", 0x1p-16382L},
+       // 2^-16382 + 2^-16446 + eps
+       {".3362103143112093506444937793915876332724499641527442230928779770593420866576778e-4931", 0x1.0000000000000002p-16382L},
+};
+static struct {
+       char *s;
+       double f;
+} t[] = {
+//     {"-.00000", -0.0},
+       {"1e+1000000", INFINITY},
+       {"1e-1000000", 0},
+       // 2^-1074 * 0.5 - eps
+       {".2470328229206232720882843964341106861825299013071623822127928412503377536351043e-323", 0},
+       // 2^-1074 * 0.5 + eps
+       {".2470328229206232720882843964341106861825299013071623822127928412503377536351044e-323", 0x1p-1074},
+       // 2^-1074 * 1.5 - eps
+       {".7410984687618698162648531893023320585475897039214871466383785237510132609053131e-323", 0x1p-1074},
+       // 2^-1074 * 1.5 + eps
+       {".7410984687618698162648531893023320585475897039214871466383785237510132609053132e-323", 0x1p-1073},
+       // 2^-1022 + 2^-1075 - eps
+       {".2225073858507201630123055637955676152503612414573018013083228724049586647606759e-307", 0x1p-1022},
+       // 2^-1022 + 2^-1075 + eps
+       {".2225073858507201630123055637955676152503612414573018013083228724049586647606760e-307", 0x1.0000000000001p-1022},
+};
+static struct {
+       char *s;
+       float f;
+} tf[] = {
+       // 2^-149 * 0.5 - eps
+       {".7006492321624085354618647916449580656401309709382578858785341419448955413429303e-45", 0},
+       // 2^-149 * 0.5 + eps
+       {".7006492321624085354618647916449580656401309709382578858785341419448955413429304e-45", 0x1p-149},
+       // 2^-149 * 0.5 - eps
+       {".2101947696487225606385594374934874196920392912814773657635602425834686624028790e-44", 0x1p-149},
+       // 2^-149 * 0.5 + eps
+       {".2101947696487225606385594374934874196920392912814773657635602425834686624028791e-44", 0x1p-148},
+       // 2^-126 + 2^-150 - eps
+       {".1175494420887210724209590083408724842314472120785184615334540294131831453944281e-37", 0x1p-126},
+       // 2^-126 + 2^-150 + eps
+       {".1175494420887210724209590083408724842314472120785184615334540294131831453944282e-37", 0x1.000002p-126},
+};
+
+
+void test_strtold()
+{
+       int i;
+       long double x;
+       char *p;
+
+       for (i = 0; i < length(tl); i++) {
+               x = strtold(tl[i].s, &p);
+               if (x != tl[i].f)
+                       error("strtold(\"%s\") want %La got %La\n", tl[i].s, tl[i].f, x);
+       }
+}
+
+void test_strtod()
+{
+       int i;
+       double x;
+       char *p;
+
+       for (i = 0; i < length(t); i++) {
+               x = strtod(t[i].s, &p);
+               if (x != t[i].f)
+                       error("strtod(\"%s\") want %a got %a\n", t[i].s, t[i].f, x);
+       }
+}
+
+void test_strtof()
+{
+       int i;
+       float x;
+       char *p;
+
+       for (i = 0; i < length(tf); i++) {
+               x = strtof(tf[i].s, &p);
+               if (x != tf[i].f)
+                       error("strtof(\"%s\") want %a got %a\n", tf[i].s, tf[i].f, x);
+       }
+}