1 /*- 2 * Copyright (C) 2007 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 nan(), nanf(), and nanl(). We also test that strtod("nan(...)") 29 * and sscanf("nan(...)", ...) work identically. 30 */ 31 32 #include <sys/cdefs.h> 33 __FBSDID("$FreeBSD$"); 34 35 #include <sys/param.h> 36 #include <fenv.h> 37 #include <float.h> 38 #include <locale.h> 39 #include <math.h> 40 #include <stdio.h> 41 #include <stdlib.h> 42 #include <string.h> 43 44 #include "test-utils.h" 45 46 static void 47 testnan(const char *nan_format) 48 { 49 char nan_str[128]; 50 char *end; 51 long double ald[4]; 52 double ad[4]; 53 float af[4]; 54 unsigned i; 55 56 snprintf(nan_str, sizeof(nan_str), "nan(%s)", nan_format); 57 for (i = 0; i < nitems(ad); i++) { 58 /* 59 * x86 has an 80-bit long double stored in 96 bits, 60 * so we need to initialize the memory for the memcmp() 61 * checks below to work. 62 */ 63 bzero(&af[i], sizeof(float)); 64 bzero(&ad[i], sizeof(double)); 65 bzero(&ald[i], sizeof(long double)); 66 } 67 68 af[0] = nanf(nan_format); 69 ATF_REQUIRE(isnan(af[0])); 70 af[1] = strtof(nan_str, &end); 71 ATF_REQUIRE(end == nan_str + strlen(nan_str)); 72 ATF_REQUIRE(sscanf(nan_str, "%e", &af[2]) == 1); 73 ATF_REQUIRE(memcmp(&af[0], &af[1], sizeof(float)) == 0); 74 ATF_REQUIRE(memcmp(&af[1], &af[2], sizeof(float)) == 0); 75 if (*nan_format == '\0') { 76 /* nanf("") == strtof("nan") */ 77 af[3] = strtof("nan", NULL); 78 ATF_REQUIRE(memcmp(&af[2], &af[3], sizeof(float)) == 0); 79 } 80 81 ad[0] = nan(nan_format); 82 ATF_REQUIRE(isnan(ad[0])); 83 ad[1] = strtod(nan_str, &end); 84 ATF_REQUIRE(end == nan_str + strlen(nan_str)); 85 ATF_REQUIRE(sscanf(nan_str, "%le", &ad[2]) == 1); 86 ATF_REQUIRE(memcmp(&ad[0], &ad[1], sizeof(double)) == 0); 87 ATF_REQUIRE(memcmp(&ad[1], &ad[2], sizeof(double)) == 0); 88 if (*nan_format == '\0') { 89 /* nan("") == strtod("nan") */ 90 ad[3] = strtod("nan", NULL); 91 ATF_REQUIRE(memcmp(&ad[2], &ad[3], sizeof(double)) == 0); 92 } 93 94 ald[0] = nanl(nan_format); 95 ATF_REQUIRE(isnan(ald[0])); 96 ald[1] = strtold(nan_str, &end); 97 ATF_REQUIRE(end == nan_str + strlen(nan_str)); 98 ATF_REQUIRE(sscanf(nan_str, "%Le", &ald[2]) == 1); 99 ATF_REQUIRE(memcmp(&ald[0], &ald[1], sizeof(long double)) == 0); 100 ATF_REQUIRE(memcmp(&ald[1], &ald[2], sizeof(long double)) == 0); 101 if (*nan_format == '\0') { 102 /* nanl("") == strtold("nan") */ 103 ald[3] = strtold("nan", NULL); 104 ATF_REQUIRE(memcmp(&ald[2], &ald[3], sizeof(long double)) == 0); 105 } 106 } 107 108 ATF_TC_WITHOUT_HEAD(nan); 109 ATF_TC_BODY(nan, tc) 110 { 111 /* Die if a signalling NaN is returned */ 112 feenableexcept(FE_INVALID); 113 114 testnan("0x1234"); 115 testnan(""); 116 } 117 118 ATF_TP_ADD_TCS(tp) 119 { 120 ATF_TP_ADD_TC(tp, nan); 121 122 return (atf_no_error()); 123 } 124