add subnormal strtod tests
[libc-test] / src / stdlib / strtod.c
1 #include <stdio.h>
2 #include <stdlib.h>
3 #include <string.h>
4 #include <math.h>
5 #include "test.h"
6
7 /* r = place to store result
8  * f = function call to test (or any expression)
9  * x = expected result
10  * m = message to print on failure (with formats for r & x)
11 **/
12
13 #define TEST(r, f, x, m) ( \
14         ((r) = (f)) == (x) || \
15         (error("%s failed (" m ")\n", #f, r, x, r-x), 0) )
16
17 void test_strtod_simple(void) {
18         int i;
19         double d, d2;
20         char buf[1000];
21
22         for (i=0; i<100; i++) {
23                 d = sin(i);
24                 snprintf(buf, sizeof buf, "%.300f", d);
25                 TEST(d2, strtod(buf, 0), d, "round trip fail %a != %a (%a)");
26         }
27
28         TEST(d, strtod("0x1p4", 0), 16.0, "hex float %a != %a");
29         TEST(d, strtod("0x1.1p4", 0), 17.0, "hex float %a != %a");
30 }
31
32 #define length(x) (sizeof(x) / sizeof(*(x)))
33
34 /* TODO: float exceptions, rounding mode, endptr check */
35
36 static struct {
37         char *s;
38         long double f;
39 } tl[] = {
40         {"12.345", 12.345L},
41         {"1.2345e1", 12.345L},
42         // 2^-16445 * 0.5 - eps
43         {".1822599765941237301264202966809709908199525407846781671860490243514185844316698e-4950", 0},
44         // 2^-16445 * 0.5 + eps
45         {".1822599765941237301264202966809709908199525407846781671860490243514185844316699e-4950", 0x1p-16445L},
46         // 2^-16445 * 1.5 - eps
47         {".5467799297823711903792608900429129724598576223540345015581470730542557532950096e-4950", 0x1p-16445L},
48         // 2^-16445 * 1.5 + eps
49         {".5467799297823711903792608900429129724598576223540345015581470730542557532950097e-4950", 0x1p-16444L},
50         // 2^-16382 + 2^-16446 - eps
51         {".3362103143112093506444937793915876332724499641527442230928779770593420866576777e-4931", 0x1p-16382L},
52         // 2^-16382 + 2^-16446 + eps
53         {".3362103143112093506444937793915876332724499641527442230928779770593420866576778e-4931", 0x1.0000000000000002p-16382L},
54 };
55 static struct {
56         char *s;
57         double f;
58 } t[] = {
59 //      {"-.00000", -0.0},
60         {"1e+1000000", INFINITY},
61         {"1e-1000000", 0},
62         // 2^-1074 * 0.5 - eps
63         {".2470328229206232720882843964341106861825299013071623822127928412503377536351043e-323", 0},
64         // 2^-1074 * 0.5 + eps
65         {".2470328229206232720882843964341106861825299013071623822127928412503377536351044e-323", 0x1p-1074},
66         // 2^-1074 * 1.5 - eps
67         {".7410984687618698162648531893023320585475897039214871466383785237510132609053131e-323", 0x1p-1074},
68         // 2^-1074 * 1.5 + eps
69         {".7410984687618698162648531893023320585475897039214871466383785237510132609053132e-323", 0x1p-1073},
70         // 2^-1022 + 2^-1075 - eps
71         {".2225073858507201630123055637955676152503612414573018013083228724049586647606759e-307", 0x1p-1022},
72         // 2^-1022 + 2^-1075 + eps
73         {".2225073858507201630123055637955676152503612414573018013083228724049586647606760e-307", 0x1.0000000000001p-1022},
74 };
75 static struct {
76         char *s;
77         float f;
78 } tf[] = {
79         // 2^-149 * 0.5 - eps
80         {".7006492321624085354618647916449580656401309709382578858785341419448955413429303e-45", 0},
81         // 2^-149 * 0.5 + eps
82         {".7006492321624085354618647916449580656401309709382578858785341419448955413429304e-45", 0x1p-149},
83         // 2^-149 * 0.5 - eps
84         {".2101947696487225606385594374934874196920392912814773657635602425834686624028790e-44", 0x1p-149},
85         // 2^-149 * 0.5 + eps
86         {".2101947696487225606385594374934874196920392912814773657635602425834686624028791e-44", 0x1p-148},
87         // 2^-126 + 2^-150 - eps
88         {".1175494420887210724209590083408724842314472120785184615334540294131831453944281e-37", 0x1p-126},
89         // 2^-126 + 2^-150 + eps
90         {".1175494420887210724209590083408724842314472120785184615334540294131831453944282e-37", 0x1.000002p-126},
91 };
92
93
94 void test_strtold()
95 {
96         int i;
97         long double x;
98         char *p;
99
100         for (i = 0; i < length(tl); i++) {
101                 x = strtold(tl[i].s, &p);
102                 if (x != tl[i].f)
103                         error("strtold(\"%s\") want %La got %La\n", tl[i].s, tl[i].f, x);
104         }
105 }
106
107 void test_strtod()
108 {
109         int i;
110         double x;
111         char *p;
112
113         for (i = 0; i < length(t); i++) {
114                 x = strtod(t[i].s, &p);
115                 if (x != t[i].f)
116                         error("strtod(\"%s\") want %a got %a\n", t[i].s, t[i].f, x);
117         }
118 }
119
120 void test_strtof()
121 {
122         int i;
123         float x;
124         char *p;
125
126         for (i = 0; i < length(tf); i++) {
127                 x = strtof(tf[i].s, &p);
128                 if (x != tf[i].f)
129                         error("strtof(\"%s\") want %a got %a\n", tf[i].s, tf[i].f, x);
130         }
131 }