1 /* $NetBSD: t_strtod.c,v 1.30 2012/05/15 18:46:20 alnsn Exp $ */ 2 3 /*- 4 * Copyright (c) 2011 The NetBSD Foundation, Inc. 5 * All rights reserved. 6 * 7 * This code is derived from software contributed to The NetBSD Foundation 8 * by Jukka Ruohonen. 9 * 10 * Redistribution and use in source and binary forms, with or without 11 * modification, are permitted provided that the following conditions 12 * are met: 13 * 1. Redistributions of source code must retain the above copyright 14 * notice, this list of conditions and the following disclaimer. 15 * 2. Redistributions in binary form must reproduce the above copyright 16 * notice, this list of conditions and the following disclaimer in the 17 * documentation and/or other materials provided with the distribution. 18 * 19 * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS 20 * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED 21 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR 22 * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS 23 * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 24 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF 25 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 26 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN 27 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) 28 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 29 * POSSIBILITY OF SUCH DAMAGE. 30 */ 31 32 /* Public domain, Otto Moerbeek <otto@drijf.net>, 2006. */ 33 34 #include <sys/cdefs.h> 35 __RCSID("$NetBSD: t_strtod.c,v 1.30 2012/05/15 18:46:20 alnsn Exp $"); 36 37 #include <errno.h> 38 #include <math.h> 39 #include <stdio.h> 40 #include <stdlib.h> 41 #include <string.h> 42 43 #include <atf-c.h> 44 #include <atf-c/config.h> 45 46 #if defined(__i386__) || defined(__amd64__) || defined(__sparc__) 47 #include <fenv.h> 48 #endif 49 50 #if !defined(__vax__) 51 static const char * const inf_strings[] = 52 { "Inf", "INF", "-Inf", "-INF", "Infinity", "+Infinity", 53 "INFINITY", "-INFINITY", "InFiNiTy", "+InFiNiTy" }; 54 const char *nan_string = "NaN(x)y"; 55 #endif 56 57 ATF_TC(strtod_basic); 58 ATF_TC_HEAD(strtod_basic, tc) 59 { 60 atf_tc_set_md_var(tc, "descr", "A basic test of strtod(3)"); 61 } 62 63 ATF_TC_BODY(strtod_basic, tc) 64 { 65 static const size_t n = 1024 * 1000; 66 67 for (size_t i = 1; i < n; i = i + 1024) { 68 char buf[512]; 69 (void)snprintf(buf, sizeof(buf), "%zu.%zu", i, i + 1); 70 71 errno = 0; 72 double d = strtod(buf, NULL); 73 74 ATF_REQUIRE(d > 0.0); 75 ATF_REQUIRE(errno == 0); 76 } 77 } 78 79 ATF_TC(strtod_hex); 80 ATF_TC_HEAD(strtod_hex, tc) 81 { 82 atf_tc_set_md_var(tc, "descr", "A strtod(3) with hexadecimals"); 83 } 84 85 #ifdef __vax__ 86 #define SMALL_NUM 1.0e-38 87 #else 88 #define SMALL_NUM 1.0e-40 89 #endif 90 91 ATF_TC_BODY(strtod_hex, tc) 92 { 93 const char *str; 94 char *end; 95 volatile double d; 96 97 str = "-0x0"; 98 d = strtod(str, &end); /* -0.0 */ 99 100 ATF_REQUIRE(end == str + 4); 101 ATF_REQUIRE(signbit(d) != 0); 102 ATF_REQUIRE(fabs(d) < SMALL_NUM); 103 104 str = "-0x"; 105 d = strtod(str, &end); /* -0.0 */ 106 107 ATF_REQUIRE(end == str + 2); 108 ATF_REQUIRE(signbit(d) != 0); 109 ATF_REQUIRE(fabs(d) < SMALL_NUM); 110 } 111 112 ATF_TC(strtod_inf); 113 ATF_TC_HEAD(strtod_inf, tc) 114 { 115 atf_tc_set_md_var(tc, "descr", "A strtod(3) with INF (PR lib/33262)"); 116 } 117 118 ATF_TC_BODY(strtod_inf, tc) 119 { 120 #ifndef __vax__ 121 for (size_t i = 0; i < __arraycount(inf_strings); i++) { 122 volatile double d = strtod(inf_strings[i], NULL); 123 ATF_REQUIRE(isinf(d) != 0); 124 } 125 #else 126 atf_tc_skip("vax not supported"); 127 #endif 128 } 129 130 ATF_TC(strtof_inf); 131 ATF_TC_HEAD(strtof_inf, tc) 132 { 133 atf_tc_set_md_var(tc, "descr", "A strtof(3) with INF (PR lib/33262)"); 134 } 135 136 ATF_TC_BODY(strtof_inf, tc) 137 { 138 #ifndef __vax__ 139 for (size_t i = 0; i < __arraycount(inf_strings); i++) { 140 volatile float f = strtof(inf_strings[i], NULL); 141 ATF_REQUIRE(isinf(f) != 0); 142 } 143 #else 144 atf_tc_skip("vax not supported"); 145 #endif 146 } 147 148 ATF_TC(strtold_inf); 149 ATF_TC_HEAD(strtold_inf, tc) 150 { 151 atf_tc_set_md_var(tc, "descr", "A strtold(3) with INF (PR lib/33262)"); 152 } 153 154 ATF_TC_BODY(strtold_inf, tc) 155 { 156 #ifndef __vax__ 157 # ifdef __HAVE_LONG_DOUBLE 158 if (system("cpuctl identify 0 | grep -q QEMU") == 0) 159 atf_tc_expect_fail("PR misc/44767"); 160 161 for (size_t i = 0; i < __arraycount(inf_strings); i++) { 162 volatile long double ld = strtold(inf_strings[i], NULL); 163 ATF_REQUIRE(isinf(ld) != 0); 164 } 165 # else 166 atf_tc_skip("Requires long double support"); 167 # endif 168 #else 169 atf_tc_skip("vax not supported"); 170 #endif 171 } 172 173 ATF_TC(strtod_nan); 174 ATF_TC_HEAD(strtod_nan, tc) 175 { 176 atf_tc_set_md_var(tc, "descr", "A strtod(3) with NaN"); 177 } 178 179 ATF_TC_BODY(strtod_nan, tc) 180 { 181 #ifndef __vax__ 182 char *end; 183 184 volatile double d = strtod(nan_string, &end); 185 ATF_REQUIRE(isnan(d) != 0); 186 ATF_REQUIRE(strcmp(end, "y") == 0); 187 #else 188 atf_tc_skip("vax not supported"); 189 #endif 190 } 191 192 ATF_TC(strtof_nan); 193 ATF_TC_HEAD(strtof_nan, tc) 194 { 195 atf_tc_set_md_var(tc, "descr", "A strtof(3) with NaN"); 196 } 197 198 ATF_TC_BODY(strtof_nan, tc) 199 { 200 #ifndef __vax__ 201 char *end; 202 203 volatile float f = strtof(nan_string, &end); 204 ATF_REQUIRE(isnanf(f) != 0); 205 ATF_REQUIRE(strcmp(end, "y") == 0); 206 #else 207 atf_tc_skip("vax not supported"); 208 #endif 209 } 210 211 ATF_TC(strtold_nan); 212 ATF_TC_HEAD(strtold_nan, tc) 213 { 214 atf_tc_set_md_var(tc, "descr", "A strtold(3) with NaN (PR lib/45020)"); 215 } 216 217 ATF_TC_BODY(strtold_nan, tc) 218 { 219 #ifndef __vax__ 220 # ifdef __HAVE_LONG_DOUBLE 221 222 char *end; 223 224 if (system("cpuctl identify 0 | grep -q QEMU") == 0) 225 atf_tc_expect_fail("PR misc/44767"); 226 227 volatile long double ld = strtold(nan_string, &end); 228 ATF_REQUIRE(isnan(ld) != 0); 229 ATF_REQUIRE(__isnanl(ld) != 0); 230 ATF_REQUIRE(strcmp(end, "y") == 0); 231 # else 232 atf_tc_skip("Requires long double support"); 233 # endif 234 #else 235 atf_tc_skip("vax not supported"); 236 #endif 237 } 238 239 ATF_TC(strtod_round); 240 ATF_TC_HEAD(strtod_round, tc) 241 { 242 atf_tc_set_md_var(tc, "descr", "Test rouding in strtod(3)"); 243 } 244 245 ATF_TC_BODY(strtod_round, tc) 246 { 247 #if defined(__i386__) || defined(__amd64__) || defined(__sparc__) 248 249 /* 250 * Test that strtod(3) honors the current rounding mode. 251 * The used value is somewhere near 1 + DBL_EPSILON + FLT_EPSILON. 252 */ 253 const char *val = 254 "1.00000011920928977282585492503130808472633361816406"; 255 256 (void)fesetround(FE_UPWARD); 257 258 volatile double d1 = strtod(val, NULL); 259 260 (void)fesetround(FE_DOWNWARD); 261 262 volatile double d2 = strtod(val, NULL); 263 264 if (fabs(d1 - d2) > 0.0) 265 return; 266 else { 267 atf_tc_expect_fail("PR misc/44767"); 268 atf_tc_fail("strtod(3) did not honor fesetround(3)"); 269 } 270 #else 271 atf_tc_skip("Requires one of i386, amd64 or sparc"); 272 #endif 273 } 274 275 ATF_TC(strtod_underflow); 276 ATF_TC_HEAD(strtod_underflow, tc) 277 { 278 atf_tc_set_md_var(tc, "descr", "Test underflow in strtod(3)"); 279 } 280 281 ATF_TC_BODY(strtod_underflow, tc) 282 { 283 284 const char *tmp = 285 "0.0000000000000000000000000000000000000000000000000000" 286 "000000000000000000000000000000000000000000000000000000" 287 "000000000000000000000000000000000000000000000000000000" 288 "000000000000000000000000000000000000000000000000000000" 289 "000000000000000000000000000000000000000000000000000000" 290 "000000000000000000000000000000000000000000000000000000" 291 "000000000000000000000000000000000000000000000000000000" 292 "000000000000000002"; 293 294 errno = 0; 295 volatile double d = strtod(tmp, NULL); 296 297 if (d != 0 || errno != ERANGE) 298 atf_tc_fail("strtod(3) did not detect underflow"); 299 } 300 301 /* 302 * Bug found by Geza Herman. 303 * See 304 * http://www.exploringbinary.com/a-bug-in-the-bigcomp-function-of-david-gays-strtod/ 305 */ 306 ATF_TC(strtod_gherman_bug); 307 ATF_TC_HEAD(strtod_gherman_bug, tc) 308 { 309 atf_tc_set_md_var(tc, "descr", "Test a bug found by Geza Herman"); 310 } 311 312 ATF_TC_BODY(strtod_gherman_bug, tc) 313 { 314 315 const char *str = 316 "1.8254370818746402660437411213933955878019332885742187"; 317 318 errno = 0; 319 volatile double d = strtod(str, NULL); 320 321 ATF_CHECK(d == 0x1.d34fd8378ea83p+0); 322 } 323 324 ATF_TP_ADD_TCS(tp) 325 { 326 327 ATF_TP_ADD_TC(tp, strtod_basic); 328 ATF_TP_ADD_TC(tp, strtod_hex); 329 ATF_TP_ADD_TC(tp, strtod_inf); 330 ATF_TP_ADD_TC(tp, strtof_inf); 331 ATF_TP_ADD_TC(tp, strtold_inf); 332 ATF_TP_ADD_TC(tp, strtod_nan); 333 ATF_TP_ADD_TC(tp, strtof_nan); 334 ATF_TP_ADD_TC(tp, strtold_nan); 335 ATF_TP_ADD_TC(tp, strtod_round); 336 ATF_TP_ADD_TC(tp, strtod_underflow); 337 ATF_TP_ADD_TC(tp, strtod_gherman_bug); 338 339 return atf_no_error(); 340 } 341