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$"); 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)) 45 static 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 54 ATF_TC_WITHOUT_HEAD(normalized_numbers); 55 ATF_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 153 ATF_TC_WITHOUT_HEAD(infinities_and_nans); 154 ATF_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 205 ATF_TC_WITHOUT_HEAD(rounding_tests); 206 ATF_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 287 ATF_TC_WITHOUT_HEAD(strtod); 288 ATF_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 316 ATF_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