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 <fenv.h> 32 #include <float.h> 33 #include <locale.h> 34 #include <math.h> 35 #include <stdio.h> 36 #include <stdlib.h> 37 #include <string.h> 38 39 #include <atf-c.h> 40 41 #define eq(type, a, b) _eq(type##_EPSILON, (a), (b)) 42 static int 43 _eq(long double epsilon, long double a, long double b) 44 { 45 long double delta; 46 47 delta = fabsl(a - b); 48 return (delta <= epsilon); 49 } 50 51 ATF_TC_WITHOUT_HEAD(normalized_numbers); 52 ATF_TC_BODY(normalized_numbers, tc) 53 { 54 char buf[128]; 55 long double ld = 0.0; 56 double d = 0.0; 57 float f = 0.0; 58 59 buf[0] = '\0'; 60 ATF_REQUIRE(setlocale(LC_NUMERIC, "")); 61 62 sscanf("3.141592", "%e", &f); 63 ATF_REQUIRE(eq(FLT, f, 3.141592)); 64 65 sscanf("3.141592653589793", "%lf", &d); 66 ATF_REQUIRE(eq(DBL, d, 3.141592653589793)); 67 68 sscanf("1.234568e+06", "%E", &f); 69 ATF_REQUIRE(eq(FLT, f, 1.234568e+06)); 70 71 sscanf("-1.234568e6", "%lF", &d); 72 ATF_REQUIRE(eq(DBL, d, -1.234568e6)); 73 74 sscanf("+1.234568e-52", "%LG", &ld); 75 ATF_REQUIRE(eq(LDBL, ld, 1.234568e-52L)); 76 77 sscanf("0.1", "%la", &d); 78 ATF_REQUIRE(eq(DBL, d, 0.1)); 79 80 sscanf("00.2", "%lA", &d); 81 ATF_REQUIRE(eq(DBL, d, 0.2)); 82 83 sscanf("123456", "%5le%s", &d, buf); 84 ATF_REQUIRE(eq(DBL, d, 12345.)); 85 ATF_REQUIRE(strcmp(buf, "6") == 0); 86 87 sscanf("1.0Q", "%*5le%s", buf); 88 ATF_REQUIRE(strcmp(buf, "Q") == 0); 89 90 sscanf("-1.23e", "%e%s", &f, buf); 91 ATF_REQUIRE(eq(FLT, f, -1.23)); 92 ATF_REQUIRE(strcmp(buf, "e") == 0); 93 94 sscanf("1.25e+", "%le%s", &d, buf); 95 ATF_REQUIRE(eq(DBL, d, 1.25)); 96 ATF_REQUIRE(strcmp(buf, "e+") == 0); 97 98 sscanf("1.23E4E5", "%le%s", &d, buf); 99 ATF_REQUIRE(eq(DBL, d, 1.23e4)); 100 ATF_REQUIRE(strcmp(buf, "E5") == 0); 101 102 sscanf("12e6", "%le", &d); 103 ATF_REQUIRE(eq(DBL, d, 12e6)); 104 105 sscanf("1.a", "%le%s", &d, buf); 106 ATF_REQUIRE(eq(DBL, d, 1.0)); 107 ATF_REQUIRE(strcmp(buf, "a") == 0); 108 109 sscanf(".0p4", "%le%s", &d, buf); 110 ATF_REQUIRE(eq(DBL, d, 0.0)); 111 ATF_REQUIRE(strcmp(buf, "p4") == 0); 112 113 d = 0.25; 114 ATF_REQUIRE(sscanf(".", "%le", &d) == 0); 115 ATF_REQUIRE(d == 0.25); 116 117 sscanf("0x08", "%le", &d); 118 ATF_REQUIRE(d == 0x8p0); 119 120 sscanf("0x90a.bcdefP+09a", "%le%s", &d, buf); 121 ATF_REQUIRE(d == 0x90a.bcdefp+09); 122 ATF_REQUIRE(strcmp(buf, "a") == 0); 123 124 #if (LDBL_MANT_DIG > DBL_MANT_DIG) && !defined(__i386__) 125 sscanf("3.14159265358979323846", "%Lg", &ld); 126 ATF_REQUIRE(eq(LDBL, ld, 3.14159265358979323846L)); 127 128 sscanf(" 0X.0123456789abcdefffp-3g", "%Le%s", &ld, buf); 129 ATF_REQUIRE(ld == 0x0.0123456789abcdefffp-3L); 130 ATF_REQUIRE(strcmp(buf, "g") == 0); 131 #endif 132 133 sscanf("0xg", "%le%s", &d, buf); 134 ATF_REQUIRE(d == 0.0); 135 ATF_REQUIRE(strcmp(buf, "xg") == 0); 136 137 ATF_REQUIRE(setlocale(LC_NUMERIC, "ru_RU.ISO8859-5")); /* decimalpoint==, */ 138 139 sscanf("1.23", "%le%s", &d, buf); 140 ATF_REQUIRE(d == 1.0); 141 ATF_REQUIRE(strcmp(buf, ".23") == 0); 142 143 sscanf("1,23", "%le", &d); 144 ATF_REQUIRE(d == 1.23); 145 146 ATF_REQUIRE(setlocale(LC_NUMERIC, "")); 147 } 148 149 ATF_TC_WITHOUT_HEAD(infinities_and_nans); 150 ATF_TC_BODY(infinities_and_nans, tc) 151 { 152 char buf[128]; 153 long double ld = 0.0; 154 double d = 0.0; 155 float f = 0.0; 156 157 ATF_REQUIRE(setlocale(LC_NUMERIC, "C")); 158 159 sscanf("-Inf", "%le", &d); 160 ATF_REQUIRE(d < 0.0 && isinf(d)); 161 162 sscanf("iNfInItY and beyond", "%le%s", &d, buf); 163 ATF_REQUIRE(d > 0.0 && isinf(d)); 164 ATF_REQUIRE(strcmp(buf, " and beyond")); 165 166 sscanf("NaN", "%le", &d); 167 ATF_REQUIRE(isnan(d)); 168 169 sscanf("NAN(123Y", "%le%s", &d, buf); 170 ATF_REQUIRE(isnan(d)); 171 ATF_REQUIRE(strcmp(buf, "(123Y") == 0); 172 173 sscanf("nan(f00f)plugh", "%le%s", &d, buf); 174 ATF_REQUIRE(isnan(d)); 175 ATF_REQUIRE(strcmp(buf, "plugh") == 0); 176 177 sscanf("-nan", "%le", &d); 178 ATF_REQUIRE(isnan(d)); 179 180 /* Only quiet NaNs should be returned. */ 181 sscanf("NaN", "%e", &f); 182 sscanf("nan", "%le", &d); 183 sscanf("nan", "%Le", &ld); 184 feclearexcept(FE_ALL_EXCEPT); 185 ATF_REQUIRE(f != f); 186 ATF_REQUIRE(d != d); 187 ATF_REQUIRE(ld != ld); 188 ATF_REQUIRE(fetestexcept(FE_INVALID) == 0); 189 sscanf("nan(1234)", "%e", &f); 190 sscanf("nan(1234)", "%le", &d); 191 sscanf("nan(1234)", "%Le", &ld); 192 feclearexcept(FE_ALL_EXCEPT); 193 ATF_REQUIRE(f != f); 194 ATF_REQUIRE(d != d); 195 ATF_REQUIRE(ld != ld); 196 /* POSIX says we should only generate quiet NaNs. */ 197 ATF_REQUIRE(fetestexcept(FE_INVALID) == 0); 198 } 199 200 ATF_TC_WITHOUT_HEAD(rounding_tests); 201 ATF_TC_BODY(rounding_tests, tc) 202 { 203 long double ld = 0.0; 204 double d = 0.0; 205 206 ATF_REQUIRE(setlocale(LC_NUMERIC, "C")); 207 208 fesetround(FE_DOWNWARD); 209 210 sscanf("1.999999999999999999999999999999999", "%le", &d); 211 ATF_REQUIRE(d < 2.0); 212 sscanf("0x1.ffffffffffffffp0", "%le", &d); 213 ATF_REQUIRE(d < 2.0); 214 sscanf("1.999999999999999999999999999999999", "%Le", &ld); 215 ATF_REQUIRE(ld < 2.0); 216 217 sscanf("1.0571892669084007", "%le", &d); 218 ATF_REQUIRE(d == 0x1.0ea3f4af0dc59p0); 219 sscanf("-1.0571892669084007", "%le", &d); 220 ATF_REQUIRE(d == -0x1.0ea3f4af0dc5ap0); 221 sscanf("1.0571892669084010", "%le", &d); 222 ATF_REQUIRE(d == 0x1.0ea3f4af0dc5ap0); 223 224 sscanf("0x1.23p-5000", "%le", &d); 225 ATF_REQUIRE(d == 0.0); 226 227 sscanf("0x1.2345678p-1050", "%le", &d); 228 ATF_REQUIRE(d == 0x1.234567p-1050); 229 230 fesetround(FE_UPWARD); 231 232 sscanf("1.0571892669084007", "%le", &d); 233 ATF_REQUIRE(d == 0x1.0ea3f4af0dc5ap0); 234 sscanf("-1.0571892669084007", "%le", &d); 235 ATF_REQUIRE(d == -0x1.0ea3f4af0dc59p0); 236 sscanf("1.0571892669084010", "%le", &d); 237 ATF_REQUIRE(d == 0x1.0ea3f4af0dc5bp0); 238 239 sscanf("0x1.23p-5000", "%le", &d); 240 ATF_REQUIRE(d == 0x1p-1074); 241 242 sscanf("0x1.2345678p-1050", "%le", &d); 243 ATF_REQUIRE(d == 0x1.234568p-1050); 244 245 fesetround(FE_TOWARDZERO); 246 247 sscanf("1.0571892669084007", "%le", &d); 248 ATF_REQUIRE(d == 0x1.0ea3f4af0dc59p0); 249 sscanf("-1.0571892669084007", "%le", &d); 250 ATF_REQUIRE(d == -0x1.0ea3f4af0dc59p0); 251 sscanf("1.0571892669084010", "%le", &d); 252 ATF_REQUIRE(d == 0x1.0ea3f4af0dc5ap0); 253 254 sscanf("0x1.23p-5000", "%le", &d); 255 ATF_REQUIRE(d == 0.0); 256 257 sscanf("0x1.2345678p-1050", "%le", &d); 258 ATF_REQUIRE(d == 0x1.234567p-1050); 259 260 fesetround(FE_TONEAREST); 261 262 /* 1.0571892669084007 is slightly closer to 0x1.0ea3f4af0dc59p0 */ 263 sscanf("1.0571892669084007", "%le", &d); 264 ATF_REQUIRE(d == 0x1.0ea3f4af0dc59p0); 265 sscanf("-1.0571892669084007", "%le", &d); 266 ATF_REQUIRE(d == -0x1.0ea3f4af0dc59p0); 267 sscanf("1.0571892669084010", "%le", &d); 268 ATF_REQUIRE(d == 0x1.0ea3f4af0dc5bp0); 269 270 /* strtod() should round small numbers to 0. */ 271 sscanf("0x1.23p-5000", "%le", &d); 272 ATF_REQUIRE(d == 0.0); 273 274 /* Extra digits in a denormal shouldn't break anything. */ 275 sscanf("0x1.2345678p-1050", "%le", &d); 276 ATF_REQUIRE(d == 0x1.234568p-1050); 277 } 278 279 ATF_TC_WITHOUT_HEAD(strtod); 280 ATF_TC_BODY(strtod, tc) 281 { 282 char *endp; 283 284 ATF_REQUIRE(setlocale(LC_NUMERIC, "C")); 285 286 ATF_REQUIRE(strtod("0xy", &endp) == 0); 287 ATF_REQUIRE(strcmp("xy", endp) == 0); 288 289 /* This used to cause an infinite loop and round the wrong way. */ 290 fesetround(FE_DOWNWARD); 291 ATF_REQUIRE(strtof("3.5e38", &endp) == FLT_MAX); 292 ATF_REQUIRE(strtod("2e308", &endp) == DBL_MAX); 293 fesetround(FE_UPWARD); 294 ATF_REQUIRE(strtof("3.5e38", &endp) == INFINITY); 295 ATF_REQUIRE(strtod("2e308", &endp) == INFINITY); 296 fesetround(FE_TOWARDZERO); 297 ATF_REQUIRE(strtof("3.5e38", &endp) == FLT_MAX); 298 ATF_REQUIRE(strtod("2e308", &endp) == DBL_MAX); 299 fesetround(FE_TONEAREST); 300 ATF_REQUIRE(strtof("3.5e38", &endp) == INFINITY); 301 ATF_REQUIRE(strtod("2e308", &endp) == INFINITY); 302 } 303 304 ATF_TP_ADD_TCS(tp) 305 { 306 307 ATF_TP_ADD_TC(tp, normalized_numbers); 308 ATF_TP_ADD_TC(tp, infinities_and_nans); 309 ATF_TP_ADD_TC(tp, rounding_tests); 310 ATF_TP_ADD_TC(tp, strtod); 311 312 return (atf_no_error()); 313 } 314