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