scanfloat_test.c revision 290537 
1/*-
2 * Copyright (C) 2003, 2005 David Schultz <das@FreeBSD.org>
3 * All rights reserved.
4 *
5 * Redistribution and use in source and binary forms, with or without
6 * modification, are permitted provided that the following conditions
7 * are met:
8 * 1. Redistributions of source code must retain the above copyright
9 *    notice, this list of conditions and the following disclaimer.
10 * 2. Redistributions in binary form must reproduce the above copyright
11 *    notice, this list of conditions and the following disclaimer in the
12 *    documentation and/or other materials provided with the distribution.
13 *
14 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
15 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
16 * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
17 * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
18 * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
19 * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
20 * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
21 * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
22 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
23 * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
24 * SUCH DAMAGE.
25 */
26
27/*
28 * Test for scanf() floating point formats.
29 */
30
31#include <sys/cdefs.h>
32__FBSDID("$FreeBSD: head/lib/libc/tests/stdio/scanfloat_test.c 290537 2015-11-08 06:37:50Z ngie $");
33
34#include <fenv.h>
35#include <float.h>
36#include <locale.h>
37#include <math.h>
38#include <stdio.h>
39#include <stdlib.h>
40#include <string.h>
41
42#include <atf-c.h>
43
44#define	eq(type, a, b)	_eq(type##_EPSILON, (a), (b))
45static int
46_eq(long double epsilon, long double a, long double b)
47{
48	long double delta;
49
50	delta = fabsl(a - b);
51	return (delta <= epsilon);
52}
53
54ATF_TC_WITHOUT_HEAD(normalized_numbers);
55ATF_TC_BODY(normalized_numbers, tc)
56{
57	char buf[128];
58	long double ld = 0.0;
59	double d = 0.0;
60	float f = 0.0;
61	char *endp;
62
63	buf[0] = '\0';
64	ATF_REQUIRE(setlocale(LC_NUMERIC, ""));
65
66	sscanf("3.141592", "%e", &f);
67	ATF_REQUIRE(eq(FLT, f, 3.141592));
68
69	sscanf("3.141592653589793", "%lf", &d);
70	ATF_REQUIRE(eq(DBL, d, 3.141592653589793));
71
72	sscanf("1.234568e+06", "%E", &f);
73	ATF_REQUIRE(eq(FLT, f, 1.234568e+06));
74
75	sscanf("-1.234568e6", "%lF", &d);
76	ATF_REQUIRE(eq(DBL, d, -1.234568e6));
77
78	sscanf("+1.234568e-52", "%LG", &ld);
79	ATF_REQUIRE(eq(LDBL, ld, 1.234568e-52L));
80
81	sscanf("0.1", "%la", &d);
82	ATF_REQUIRE(eq(DBL, d, 0.1));
83
84	sscanf("00.2", "%lA", &d);
85	ATF_REQUIRE(eq(DBL, d, 0.2));
86
87	sscanf("123456", "%5le%s", &d, buf);
88	ATF_REQUIRE(eq(DBL, d, 12345.));
89	ATF_REQUIRE(strcmp(buf, "6") == 0);
90
91	sscanf("1.0Q", "%*5le%s", buf);
92	ATF_REQUIRE(strcmp(buf, "Q") == 0);
93
94	sscanf("-1.23e", "%e%s", &f, buf);
95	ATF_REQUIRE(eq(FLT, f, -1.23));
96	ATF_REQUIRE(strcmp(buf, "e") == 0);
97
98	sscanf("1.25e+", "%le%s", &d, buf);
99	ATF_REQUIRE(eq(DBL, d, 1.25));
100	ATF_REQUIRE(strcmp(buf, "e+") == 0);
101
102	sscanf("1.23E4E5", "%le%s", &d, buf);
103	ATF_REQUIRE(eq(DBL, d, 1.23e4));
104	ATF_REQUIRE(strcmp(buf, "E5") == 0);
105
106	sscanf("12e6", "%le", &d);
107	ATF_REQUIRE(eq(DBL, d, 12e6));
108
109	sscanf("1.a", "%le%s", &d, buf);
110	ATF_REQUIRE(eq(DBL, d, 1.0));
111	ATF_REQUIRE(strcmp(buf, "a") == 0);
112
113	sscanf(".0p4", "%le%s", &d, buf);
114	ATF_REQUIRE(eq(DBL, d, 0.0));
115	ATF_REQUIRE(strcmp(buf, "p4") == 0);
116
117	d = 0.25;
118	ATF_REQUIRE(sscanf(".", "%le", &d) == 0);
119	ATF_REQUIRE(d == 0.25);
120
121	sscanf("0x08", "%le", &d);
122	ATF_REQUIRE(d == 0x8p0);
123
124	sscanf("0x90a.bcdefP+09a", "%le%s", &d, buf);
125	ATF_REQUIRE(d == 0x90a.bcdefp+09);
126	ATF_REQUIRE(strcmp(buf, "a") == 0);
127
128#if (LDBL_MANT_DIG > DBL_MANT_DIG) && !defined(__i386__)
129	sscanf("3.14159265358979323846", "%Lg", &ld);
130	ATF_REQUIRE(eq(LDBL, ld, 3.14159265358979323846L));
131
132	sscanf("  0X.0123456789abcdefffp-3g", "%Le%s", &ld, buf);
133	ATF_REQUIRE(ld == 0x0.0123456789abcdefffp-3L);
134	ATF_REQUIRE(strcmp(buf, "g") == 0);
135#endif
136
137	sscanf("0xg", "%le%s", &d, buf);
138	ATF_REQUIRE(d == 0.0);
139	ATF_REQUIRE(strcmp(buf, "xg") == 0);
140
141	ATF_REQUIRE(setlocale(LC_NUMERIC, "ru_RU.ISO8859-5")); /* decimalpoint==, */
142
143	sscanf("1.23", "%le%s", &d, buf);
144	ATF_REQUIRE(d == 1.0);
145	ATF_REQUIRE(strcmp(buf, ".23") == 0);
146
147	sscanf("1,23", "%le", &d);
148	ATF_REQUIRE(d == 1.23);
149
150	ATF_REQUIRE(setlocale(LC_NUMERIC, ""));
151}
152
153ATF_TC_WITHOUT_HEAD(infinities_and_nans);
154ATF_TC_BODY(infinities_and_nans, tc)
155{
156	char buf[128];
157	long double ld = 0.0;
158	double d = 0.0;
159	float f = 0.0;
160	char *endp;
161
162	ATF_REQUIRE(setlocale(LC_NUMERIC, "C"));
163
164	sscanf("-Inf", "%le", &d);
165	ATF_REQUIRE(d < 0.0 && isinf(d));
166
167	sscanf("iNfInItY and beyond", "%le%s", &d, buf);
168	ATF_REQUIRE(d > 0.0 && isinf(d));
169	ATF_REQUIRE(strcmp(buf, " and beyond"));
170
171	sscanf("NaN", "%le", &d);
172	ATF_REQUIRE(isnan(d));
173
174	sscanf("NAN(123Y", "%le%s", &d, buf);
175	ATF_REQUIRE(isnan(d));
176	ATF_REQUIRE(strcmp(buf, "(123Y") == 0);
177
178	sscanf("nan(f00f)plugh", "%le%s", &d, buf);
179	ATF_REQUIRE(isnan(d));
180	ATF_REQUIRE(strcmp(buf, "plugh") == 0);
181
182	sscanf("-nan", "%le", &d);
183	ATF_REQUIRE(isnan(d));
184
185	/* Only quiet NaNs should be returned. */
186	sscanf("NaN", "%e", &f);
187	sscanf("nan", "%le", &d);
188	sscanf("nan", "%Le", &ld);
189	feclearexcept(FE_ALL_EXCEPT);
190	ATF_REQUIRE(f != f);
191	ATF_REQUIRE(d != d);
192	ATF_REQUIRE(ld != ld);
193	ATF_REQUIRE(fetestexcept(FE_INVALID) == 0);
194	sscanf("nan(1234)", "%e", &f);
195	sscanf("nan(1234)", "%le", &d);
196	sscanf("nan(1234)", "%Le", &ld);
197	feclearexcept(FE_ALL_EXCEPT);
198	ATF_REQUIRE(f != f);
199	ATF_REQUIRE(d != d);
200	ATF_REQUIRE(ld != ld);
201	/* POSIX says we should only generate quiet NaNs. */
202	ATF_REQUIRE(fetestexcept(FE_INVALID) == 0);
203}
204
205ATF_TC_WITHOUT_HEAD(rounding_tests);
206ATF_TC_BODY(rounding_tests, tc)
207{
208	char buf[128];
209	long double ld = 0.0;
210	double d = 0.0;
211	float f = 0.0;
212	char *endp;
213
214	ATF_REQUIRE(setlocale(LC_NUMERIC, "C"));
215
216	fesetround(FE_DOWNWARD);
217
218	sscanf("1.999999999999999999999999999999999", "%le", &d);
219	ATF_REQUIRE(d < 2.0);
220	sscanf("0x1.ffffffffffffffp0", "%le", &d);
221	ATF_REQUIRE(d < 2.0);
222	sscanf("1.999999999999999999999999999999999", "%Le", &ld);
223	ATF_REQUIRE(ld < 2.0);
224
225	sscanf("1.0571892669084007", "%le", &d);
226	ATF_REQUIRE(d == 0x1.0ea3f4af0dc59p0);
227	sscanf("-1.0571892669084007", "%le", &d);
228	ATF_REQUIRE(d == -0x1.0ea3f4af0dc5ap0);
229	sscanf("1.0571892669084010", "%le", &d);
230	ATF_REQUIRE(d == 0x1.0ea3f4af0dc5ap0);
231
232	sscanf("0x1.23p-5000", "%le", &d);
233	ATF_REQUIRE(d == 0.0);
234
235	sscanf("0x1.2345678p-1050", "%le", &d);
236	ATF_REQUIRE(d == 0x1.234567p-1050);
237
238	fesetround(FE_UPWARD);
239
240	sscanf("1.0571892669084007", "%le", &d);
241	ATF_REQUIRE(d == 0x1.0ea3f4af0dc5ap0);
242	sscanf("-1.0571892669084007", "%le", &d);
243	ATF_REQUIRE(d == -0x1.0ea3f4af0dc59p0);
244	sscanf("1.0571892669084010", "%le", &d);
245	ATF_REQUIRE(d == 0x1.0ea3f4af0dc5bp0);
246
247	sscanf("0x1.23p-5000", "%le", &d);
248	ATF_REQUIRE(d == 0x1p-1074);
249
250	sscanf("0x1.2345678p-1050", "%le", &d);
251	ATF_REQUIRE(d == 0x1.234568p-1050);
252
253	fesetround(FE_TOWARDZERO);
254
255	sscanf("1.0571892669084007", "%le", &d);
256	ATF_REQUIRE(d == 0x1.0ea3f4af0dc59p0);
257	sscanf("-1.0571892669084007", "%le", &d);
258	ATF_REQUIRE(d == -0x1.0ea3f4af0dc59p0);
259	sscanf("1.0571892669084010", "%le", &d);
260	ATF_REQUIRE(d == 0x1.0ea3f4af0dc5ap0);
261
262	sscanf("0x1.23p-5000", "%le", &d);
263	ATF_REQUIRE(d == 0.0);
264
265	sscanf("0x1.2345678p-1050", "%le", &d);
266	ATF_REQUIRE(d == 0x1.234567p-1050);
267
268	fesetround(FE_TONEAREST);
269
270	/* 1.0571892669084007 is slightly closer to 0x1.0ea3f4af0dc59p0 */
271	sscanf("1.0571892669084007", "%le", &d);
272	ATF_REQUIRE(d == 0x1.0ea3f4af0dc59p0);
273	sscanf("-1.0571892669084007", "%le", &d);
274	ATF_REQUIRE(d == -0x1.0ea3f4af0dc59p0);
275	sscanf("1.0571892669084010", "%le", &d);
276	ATF_REQUIRE(d == 0x1.0ea3f4af0dc5bp0);
277
278	/* strtod() should round small numbers to 0. */
279	sscanf("0x1.23p-5000", "%le", &d);
280	ATF_REQUIRE(d == 0.0);
281
282	/* Extra digits in a denormal shouldn't break anything. */
283	sscanf("0x1.2345678p-1050", "%le", &d);
284	ATF_REQUIRE(d == 0x1.234568p-1050);
285}
286
287ATF_TC_WITHOUT_HEAD(strtod);
288ATF_TC_BODY(strtod, tc)
289{
290	char buf[128];
291	long double ld = 0.0;
292	double d = 0.0;
293	float f = 0.0;
294	char *endp;
295
296	ATF_REQUIRE(setlocale(LC_NUMERIC, "C"));
297
298	ATF_REQUIRE(strtod("0xy", &endp) == 0);
299	ATF_REQUIRE(strcmp("xy", endp) == 0);
300
301	/* This used to cause an infinite loop and round the wrong way. */
302	fesetround(FE_DOWNWARD);
303	ATF_REQUIRE(strtof("3.5e38", &endp) == FLT_MAX);
304	ATF_REQUIRE(strtod("2e308", &endp) == DBL_MAX);
305	fesetround(FE_UPWARD);
306	ATF_REQUIRE(strtof("3.5e38", &endp) == INFINITY);
307	ATF_REQUIRE(strtod("2e308", &endp) == INFINITY);
308	fesetround(FE_TOWARDZERO);
309	ATF_REQUIRE(strtof("3.5e38", &endp) == FLT_MAX);
310	ATF_REQUIRE(strtod("2e308", &endp) == DBL_MAX);
311	fesetround(FE_TONEAREST);
312	ATF_REQUIRE(strtof("3.5e38", &endp) == INFINITY);
313	ATF_REQUIRE(strtod("2e308", &endp) == INFINITY);
314}
315
316ATF_TP_ADD_TCS(tp)
317{
318
319	ATF_TP_ADD_TC(tp, normalized_numbers);
320	ATF_TP_ADD_TC(tp, infinities_and_nans);
321	ATF_TP_ADD_TC(tp, rounding_tests);
322	ATF_TP_ADD_TC(tp, strtod);
323
324	return (atf_no_error());
325}
326