|
nsz |
1be84a |
#include <stdio.h>
|
|
nsz |
1be84a |
#include <stdlib.h>
|
|
nsz |
1be84a |
#include <string.h>
|
|
nsz |
1be84a |
#include <math.h>
|
|
nsz |
1be84a |
#include "test.h"
|
|
nsz |
1be84a |
|
|
nsz |
1be84a |
/* r = place to store result
|
|
nsz |
1be84a |
* f = function call to test (or any expression)
|
|
nsz |
1be84a |
* x = expected result
|
|
nsz |
1be84a |
* m = message to print on failure (with formats for r & x)
|
|
nsz |
1be84a |
**/
|
|
nsz |
1be84a |
|
|
nsz |
1be84a |
#define TEST(r, f, x, m) ( \
|
|
nsz |
1be84a |
((r) = (f)) == (x) || \
|
|
nsz |
1be84a |
(error("%s failed (" m ")\n", #f, r, x, r-x), 0) )
|
|
nsz |
1be84a |
|
|
nsz |
5bd394 |
void test_strtod_simple(void) {
|
|
nsz |
1be84a |
int i;
|
|
nsz |
1be84a |
double d, d2;
|
|
nsz |
1be84a |
char buf[1000];
|
|
nsz |
1be84a |
|
|
nsz |
1be84a |
for (i=0; i<100; i++) {
|
|
nsz |
1be84a |
d = sin(i);
|
|
nsz |
1be84a |
snprintf(buf, sizeof buf, "%.300f", d);
|
|
nsz |
1be84a |
TEST(d2, strtod(buf, 0), d, "round trip fail %a != %a (%a)");
|
|
nsz |
1be84a |
}
|
|
nsz |
1be84a |
|
|
nsz |
1be84a |
TEST(d, strtod("0x1p4", 0), 16.0, "hex float %a != %a");
|
|
nsz |
1be84a |
TEST(d, strtod("0x1.1p4", 0), 17.0, "hex float %a != %a");
|
|
nsz |
1be84a |
}
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
#define length(x) (sizeof(x) / sizeof(*(x)))
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
/* TODO: float exceptions, rounding mode, endptr check */
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
static struct {
|
|
nsz |
5bd394 |
char *s;
|
|
nsz |
5bd394 |
long double f;
|
|
nsz |
5bd394 |
} tl[] = {
|
|
nsz |
5bd394 |
{"12.345", 12.345L},
|
|
nsz |
5bd394 |
{"1.2345e1", 12.345L},
|
|
nsz |
5bd394 |
// 2^-16445 * 0.5 - eps
|
|
nsz |
5bd394 |
{".1822599765941237301264202966809709908199525407846781671860490243514185844316698e-4950", 0},
|
|
nsz |
5bd394 |
// 2^-16445 * 0.5 + eps
|
|
nsz |
5bd394 |
{".1822599765941237301264202966809709908199525407846781671860490243514185844316699e-4950", 0x1p-16445L},
|
|
nsz |
5bd394 |
// 2^-16445 * 1.5 - eps
|
|
nsz |
5bd394 |
{".5467799297823711903792608900429129724598576223540345015581470730542557532950096e-4950", 0x1p-16445L},
|
|
nsz |
5bd394 |
// 2^-16445 * 1.5 + eps
|
|
nsz |
5bd394 |
{".5467799297823711903792608900429129724598576223540345015581470730542557532950097e-4950", 0x1p-16444L},
|
|
nsz |
5bd394 |
// 2^-16382 + 2^-16446 - eps
|
|
nsz |
5bd394 |
{".3362103143112093506444937793915876332724499641527442230928779770593420866576777e-4931", 0x1p-16382L},
|
|
nsz |
5bd394 |
// 2^-16382 + 2^-16446 + eps
|
|
nsz |
5bd394 |
{".3362103143112093506444937793915876332724499641527442230928779770593420866576778e-4931", 0x1.0000000000000002p-16382L},
|
|
nsz |
279834 |
// 2^16384 - 2^16319 - eps
|
|
nsz |
279834 |
{"118973149535723176505351158982948.86679662540046955672e4900", 0x1.fffffffffffffffep16383L},
|
|
nsz |
279834 |
// 2^16384 - 2^16319 + eps
|
|
nsz |
279834 |
{"118973149535723176505351158982948.86679662540046955673e4900", INFINITY},
|
|
nsz |
5bd394 |
};
|
|
nsz |
5bd394 |
static struct {
|
|
nsz |
5bd394 |
char *s;
|
|
nsz |
5bd394 |
double f;
|
|
nsz |
5bd394 |
} t[] = {
|
|
nsz |
f37d87 |
{"0", 0.0},
|
|
nsz |
f37d87 |
{"00.00", 0.0},
|
|
nsz |
f37d87 |
{"-.00000", -0.0},
|
|
nsz |
5bd394 |
{"1e+1000000", INFINITY},
|
|
nsz |
5bd394 |
{"1e-1000000", 0},
|
|
nsz |
5bd394 |
// 2^-1074 * 0.5 - eps
|
|
nsz |
5bd394 |
{".2470328229206232720882843964341106861825299013071623822127928412503377536351043e-323", 0},
|
|
nsz |
5bd394 |
// 2^-1074 * 0.5 + eps
|
|
nsz |
5bd394 |
{".2470328229206232720882843964341106861825299013071623822127928412503377536351044e-323", 0x1p-1074},
|
|
nsz |
5bd394 |
// 2^-1074 * 1.5 - eps
|
|
nsz |
5bd394 |
{".7410984687618698162648531893023320585475897039214871466383785237510132609053131e-323", 0x1p-1074},
|
|
nsz |
5bd394 |
// 2^-1074 * 1.5 + eps
|
|
nsz |
5bd394 |
{".7410984687618698162648531893023320585475897039214871466383785237510132609053132e-323", 0x1p-1073},
|
|
nsz |
5bd394 |
// 2^-1022 + 2^-1075 - eps
|
|
nsz |
5bd394 |
{".2225073858507201630123055637955676152503612414573018013083228724049586647606759e-307", 0x1p-1022},
|
|
nsz |
5bd394 |
// 2^-1022 + 2^-1075 + eps
|
|
nsz |
5bd394 |
{".2225073858507201630123055637955676152503612414573018013083228724049586647606760e-307", 0x1.0000000000001p-1022},
|
|
nsz |
279834 |
// 2^1024 - 2^970 - eps
|
|
nsz |
4f47ab |
{"17976931348623158079372897140530341507993413271003782693617377898044"
|
|
nsz |
4f47ab |
"49682927647509466490179775872070963302864166928879109465555478519404"
|
|
nsz |
4f47ab |
"02630657488671505820681908902000708383676273854845817711531764475730"
|
|
nsz |
4f47ab |
"27006985557136695962284291481986083493647529271907416844436551070434"
|
|
nsz |
279834 |
"2711559699508093042880177904174497791.999999999999999999999999999999", 0x1.fffffffffffffp1023},
|
|
nsz |
4f47ab |
// 2^1024 - 2^970
|
|
nsz |
4f47ab |
{"17976931348623158079372897140530341507993413271003782693617377898044"
|
|
nsz |
4f47ab |
"49682927647509466490179775872070963302864166928879109465555478519404"
|
|
nsz |
4f47ab |
"02630657488671505820681908902000708383676273854845817711531764475730"
|
|
nsz |
4f47ab |
"27006985557136695962284291481986083493647529271907416844436551070434"
|
|
nsz |
4f47ab |
"2711559699508093042880177904174497792", INFINITY},
|
|
nsz |
5bd394 |
};
|
|
nsz |
5bd394 |
static struct {
|
|
nsz |
5bd394 |
char *s;
|
|
nsz |
5bd394 |
float f;
|
|
nsz |
5bd394 |
} tf[] = {
|
|
nsz |
5bd394 |
// 2^-149 * 0.5 - eps
|
|
nsz |
5bd394 |
{".7006492321624085354618647916449580656401309709382578858785341419448955413429303e-45", 0},
|
|
nsz |
5bd394 |
// 2^-149 * 0.5 + eps
|
|
nsz |
5bd394 |
{".7006492321624085354618647916449580656401309709382578858785341419448955413429304e-45", 0x1p-149},
|
|
nsz |
5bd394 |
// 2^-149 * 0.5 - eps
|
|
nsz |
5bd394 |
{".2101947696487225606385594374934874196920392912814773657635602425834686624028790e-44", 0x1p-149},
|
|
nsz |
5bd394 |
// 2^-149 * 0.5 + eps
|
|
nsz |
5bd394 |
{".2101947696487225606385594374934874196920392912814773657635602425834686624028791e-44", 0x1p-148},
|
|
nsz |
5bd394 |
// 2^-126 + 2^-150 - eps
|
|
nsz |
5bd394 |
{".1175494420887210724209590083408724842314472120785184615334540294131831453944281e-37", 0x1p-126},
|
|
nsz |
5bd394 |
// 2^-126 + 2^-150 + eps
|
|
nsz |
5bd394 |
{".1175494420887210724209590083408724842314472120785184615334540294131831453944282e-37", 0x1.000002p-126},
|
|
nsz |
279834 |
// 2^128 - 2^103 - eps
|
|
nsz |
279834 |
{"340282356779733661637539395458142568447.9999999999999999999", 0x1.fffffep127},
|
|
nsz |
4f47ab |
// 2^128 - 2^103
|
|
nsz |
4f47ab |
{"340282356779733661637539395458142568448", INFINITY},
|
|
nsz |
5bd394 |
};
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
void test_strtold()
|
|
nsz |
5bd394 |
{
|
|
nsz |
5bd394 |
int i;
|
|
nsz |
5bd394 |
long double x;
|
|
nsz |
5bd394 |
char *p;
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
for (i = 0; i < length(tl); i++) {
|
|
nsz |
5bd394 |
x = strtold(tl[i].s, &p);
|
|
nsz |
5bd394 |
if (x != tl[i].f)
|
|
nsz |
5bd394 |
error("strtold(\"%s\") want %La got %La\n", tl[i].s, tl[i].f, x);
|
|
nsz |
5bd394 |
}
|
|
nsz |
5bd394 |
}
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
void test_strtod()
|
|
nsz |
5bd394 |
{
|
|
nsz |
5bd394 |
int i;
|
|
nsz |
5bd394 |
double x;
|
|
nsz |
5bd394 |
char *p;
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
for (i = 0; i < length(t); i++) {
|
|
nsz |
5bd394 |
x = strtod(t[i].s, &p);
|
|
nsz |
5bd394 |
if (x != t[i].f)
|
|
nsz |
5bd394 |
error("strtod(\"%s\") want %a got %a\n", t[i].s, t[i].f, x);
|
|
nsz |
5bd394 |
}
|
|
nsz |
5bd394 |
}
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
void test_strtof()
|
|
nsz |
5bd394 |
{
|
|
nsz |
5bd394 |
int i;
|
|
nsz |
5bd394 |
float x;
|
|
nsz |
5bd394 |
char *p;
|
|
nsz |
5bd394 |
|
|
nsz |
5bd394 |
for (i = 0; i < length(tf); i++) {
|
|
nsz |
5bd394 |
x = strtof(tf[i].s, &p);
|
|
nsz |
5bd394 |
if (x != tf[i].f)
|
|
nsz |
5bd394 |
error("strtof(\"%s\") want %a got %a\n", tf[i].s, tf[i].f, x);
|
|
nsz |
5bd394 |
}
|
|
nsz |
5bd394 |
}
|