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 <assert.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 void 45 testnan(const char *nan_format) 46 { 47 char nan_str[128]; 48 char *end; 49 long double ald[4]; 50 double ad[4]; 51 float af[4]; 52 int i; 53 54 snprintf(nan_str, sizeof(nan_str), "nan(%s)", nan_format); 55 for (i = 0; i < 4; i++) { 56 /* 57 * x86 has an 80-bit long double stored in 96 bits, 58 * so we need to initialize the memory for the memcmp() 59 * checks below to work. 60 */ 61 bzero(&af[i], sizeof(float)); 62 bzero(&ad[i], sizeof(double)); 63 bzero(&ald[i], sizeof(long double)); 64 65 } 66 67 af[0] = nanf(nan_format); 68 assert(isnan(af[0])); 69 af[1] = strtof(nan_str, &end); 70 assert(end == nan_str + strlen(nan_str)); 71 assert(sscanf(nan_str, "%e", &af[2]) == 1); 72 assert(memcmp(&af[0], &af[1], sizeof(float)) == 0); 73 assert(memcmp(&af[1], &af[2], sizeof(float)) == 0); 74 if (*nan_format == '\0') { 75 /* nanf("") == strtof("nan") */ 76 af[3] = strtof("nan", NULL); 77 assert(memcmp(&af[2], &af[3], sizeof(float)) == 0); 78 } 79 80 ad[0] = nan(nan_format); 81 assert(isnan(ad[0])); 82 ad[1] = strtod(nan_str, &end); 83 assert(end == nan_str + strlen(nan_str)); 84 assert(sscanf(nan_str, "%le", &ad[2]) == 1); 85 assert(memcmp(&ad[0], &ad[1], sizeof(double)) == 0); 86 assert(memcmp(&ad[1], &ad[2], sizeof(double)) == 0); 87 if (*nan_format == '\0') { 88 /* nan("") == strtod("nan") */ 89 ad[3] = strtod("nan", NULL); 90 assert(memcmp(&ad[2], &ad[3], sizeof(double)) == 0); 91 } 92 93 ald[0] = nanl(nan_format); 94 assert(isnan(ald[0])); 95 ald[1] = strtold(nan_str, &end); 96 assert(end == nan_str + strlen(nan_str)); 97 assert(sscanf(nan_str, "%Le", &ald[2]) == 1); 98 assert(memcmp(&ald[0], &ald[1], sizeof(long double)) == 0); 99 assert(memcmp(&ald[1], &ald[2], sizeof(long double)) == 0); 100 if (*nan_format == '\0') { 101 /* nanl("") == strtold("nan") */ 102 ald[3] = strtold("nan", NULL); 103 assert(memcmp(&ald[2], &ald[3], sizeof(long double)) == 0); 104 } 105 } 106 107 int 108 main(int argc, char *argv[]) 109 { 110 111 printf("1..1\n"); 112 113 /* Die if a signalling NaN is returned */ 114 feenableexcept(FE_INVALID); 115 116 testnan("0x1234"); 117 testnan(""); 118 119 printf("ok 1 - nan\n"); 120 121 return (0); 122 } 123