nsz Git - libc-test/blob - 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         // 2^16384 - 2^16319 - eps
  55         {"118973149535723176505351158982948.86679662540046955672e4900", 0x1.fffffffffffffffep16383L},
  56         // 2^16384 - 2^16319 + eps
  57         {"118973149535723176505351158982948.86679662540046955673e4900", INFINITY},
  58 };
  59 static struct {
  60         char *s;
  61         double f;
  62 } t[] = {
  63         {"0", 0.0},
  64         {"00.00", 0.0},
  65         {"-.00000", -0.0},
  66         {"1e+1000000", INFINITY},
  67         {"1e-1000000", 0},
  68         // 2^-1074 * 0.5 - eps
  69         {".2470328229206232720882843964341106861825299013071623822127928412503377536351043e-323", 0},
  70         // 2^-1074 * 0.5 + eps
  71         {".2470328229206232720882843964341106861825299013071623822127928412503377536351044e-323", 0x1p-1074},
  72         // 2^-1074 * 1.5 - eps
  73         {".7410984687618698162648531893023320585475897039214871466383785237510132609053131e-323", 0x1p-1074},
  74         // 2^-1074 * 1.5 + eps
  75         {".7410984687618698162648531893023320585475897039214871466383785237510132609053132e-323", 0x1p-1073},
  76         // 2^-1022 + 2^-1075 - eps
  77         {".2225073858507201630123055637955676152503612414573018013083228724049586647606759e-307", 0x1p-1022},
  78         // 2^-1022 + 2^-1075 + eps
  79         {".2225073858507201630123055637955676152503612414573018013083228724049586647606760e-307", 0x1.0000000000001p-1022},
  80         // 2^1024 - 2^970 - eps
  81         {"17976931348623158079372897140530341507993413271003782693617377898044"
  82         "49682927647509466490179775872070963302864166928879109465555478519404"
  83         "02630657488671505820681908902000708383676273854845817711531764475730"
  84         "27006985557136695962284291481986083493647529271907416844436551070434"
  85         "2711559699508093042880177904174497791.999999999999999999999999999999", 0x1.fffffffffffffp1023},
  86         // 2^1024 - 2^970
  87         {"17976931348623158079372897140530341507993413271003782693617377898044"
  88         "49682927647509466490179775872070963302864166928879109465555478519404"
  89         "02630657488671505820681908902000708383676273854845817711531764475730"
  90         "27006985557136695962284291481986083493647529271907416844436551070434"
  91         "2711559699508093042880177904174497792", INFINITY},
  92         // some random numbers
  93         {".5961860348131807091861002266453941950428e00", 0.59618603481318067}, // 0x1.313f4bc3b584cp-1
  94         {"1.815013169218038729887460898733526957442e-1", 0.18150131692180388}, // 0x1.73b6f662e1712p-3
  95         {"42.07082357534453600681618685682257590772e-2", 0.42070823575344535}, // 0x1.aece23c6e028dp-2
  96         {"665.4686306516261456328973225579833470816e-3", 0.66546863065162609}, // 0x1.54b84dea53453p-1
  97         {"6101.852922970868621786690495485449831753e-4", 0.61018529229708685}, // 0x1.386a34e5d516bp-1
  98         {"76966.95208236968077849464348875471158549e-5", 0.76966952082369677}, // 0x1.8a121f9954dfap-1
  99         {"250506.5322228682496132604807222923702304e-6", 0.25050653222286823}, // 0x1.0084c8cd538c2p-2
 100         {"2740037.230228005325852424697698331177377e-7", 0.27400372302280052}, // 0x1.18946e9575ef4p-2
 101         {"20723093.50049742645941529268715428324490e-8", 0.20723093500497428}, // 0x1.a868b14486e4dp-3
 102         {"0.7900280238081604956226011047460238748912e1", 7.9002802380816046}, // 0x1.f99e3100f2eaep+2
 103         {"0.9822860653737296848190558448760465863597e2", 98.228606537372968}, // 0x1.88ea17d506accp+6
 104         {"0.7468949723190370809405570560160405324869e3", 746.89497231903704}, // 0x1.75728e73f48b7p+9
 105         {"0.1630268320282728475980459844271031751665e4", 1630.2683202827284}, // 0x1.97912c28d5cbp+10
 106         {"0.4637168629719170695109918769645492022088e5", 46371.686297191707}, // 0x1.6a475f6258737p+15
 107         {"0.6537805944497711554209461686415872067523e6", 653780.59444977110}, // 0x1.3f3a9305bb86cp+19
 108         {"0.2346324356502437045212230713960457676531e6", 234632.43565024371}, // 0x1.ca4437c3631eap+17
 109         {"0.9709481716420048341897258980454298205278e8", 97094817.164200485}, // 0x1.7263284a8242cp+26
 110         {"0.4996908522051874110779982354932499499602e9", 499690852.20518744}, // 0x1.dc8ad6434872ap+28
 111 };
 112 static struct {
 113         char *s;
 114         float f;
 115 } tf[] = {
 116         // 2^-149 * 0.5 - eps
 117         {".7006492321624085354618647916449580656401309709382578858785341419448955413429303e-45", 0},
 118         // 2^-149 * 0.5 + eps
 119         {".7006492321624085354618647916449580656401309709382578858785341419448955413429304e-45", 0x1p-149},
 120         // 2^-149 * 0.5 - eps
 121         {".2101947696487225606385594374934874196920392912814773657635602425834686624028790e-44", 0x1p-149},
 122         // 2^-149 * 0.5 + eps
 123         {".2101947696487225606385594374934874196920392912814773657635602425834686624028791e-44", 0x1p-148},
 124         // 2^-126 + 2^-150 - eps
 125         {".1175494420887210724209590083408724842314472120785184615334540294131831453944281e-37", 0x1p-126},
 126         // 2^-126 + 2^-150 + eps
 127         {".1175494420887210724209590083408724842314472120785184615334540294131831453944282e-37", 0x1.000002p-126},
 128         // 2^128 - 2^103 - eps
 129         {"340282356779733661637539395458142568447.9999999999999999999", 0x1.fffffep127},
 130         // 2^128 - 2^103
 131         {"340282356779733661637539395458142568448", INFINITY},
 132 };
 133
 134
 135 void test_strtold()
 136 {
 137         int i;
 138         long double x;
 139         char *p;
 140
 141         for (i = 0; i < length(tl); i++) {
 142                 x = strtold(tl[i].s, &p);
 143                 if (x != tl[i].f)
 144                         error("strtold(\"%s\") want %La got %La\n", tl[i].s, tl[i].f, x);
 145         }
 146 }
 147
 148 void test_strtod()
 149 {
 150         int i;
 151         double x;
 152         char *p;
 153
 154         for (i = 0; i < length(t); i++) {
 155                 x = strtod(t[i].s, &p);
 156                 if (x != t[i].f)
 157                         error("strtod(\"%s\") want %a got %a\n", t[i].s, t[i].f, x);
 158         }
 159 }
 160
 161 void test_strtof()
 162 {
 163         int i;
 164         float x;
 165         char *p;
 166
 167         for (i = 0; i < length(tf); i++) {
 168                 x = strtof(tf[i].s, &p);
 169                 if (x != tf[i].f)
 170                         error("strtof(\"%s\") want %a got %a\n", tf[i].s, tf[i].f, x);
 171         }
 172 }
 173
 174 void bench_strtod(int N)
 175 {
 176         volatile double y;
 177         char *p;
 178         int i;
 179
 180         for (i = 0; i < N; i++)
 181                 y = strtod("2740037.230228005325852424697698331177377e-7", &p);
 182 }
 183
 184
 185 void bench_strtold_big(int N)
 186 {
 187         volatile long double y;
 188         char *p;
 189         int i;
 190
 191         for (i = 0; i < N; i++)
 192                 y = strtold("118973149535723176505351158982948.86679662540046955672e4900", &p);
 193 }
 194