lib/msun/tests/nan_test.c

   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 #include <sys/param.h>
  34 #include <fenv.h>
  35 #include <float.h>
  36 #include <locale.h>
  37 #include <math.h>
  38 #include <stdio.h>
  39 #include <stdlib.h>
  40 #include <string.h>
  41
  42 #include "test-utils.h"
  43
  44 static 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         unsigned i;
  53
  54         snprintf(nan_str, sizeof(nan_str), "nan(%s)", nan_format);
  55         for (i = 0; i < nitems(ad); 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         af[0] = nanf(nan_format);
  67         ATF_REQUIRE(isnan(af[0]));
  68         af[1] = strtof(nan_str, &end);
  69         ATF_REQUIRE(end == nan_str + strlen(nan_str));
  70         ATF_REQUIRE(sscanf(nan_str, "%e", &af[2]) == 1);
  71         ATF_REQUIRE(memcmp(&af[0], &af[1], sizeof(float)) == 0);
  72         ATF_REQUIRE(memcmp(&af[1], &af[2], sizeof(float)) == 0);
  73         if (*nan_format == '\0') {
  74                 /* nanf("") == strtof("nan") */
  75                 af[3] = strtof("nan", NULL);
  76                 ATF_REQUIRE(memcmp(&af[2], &af[3], sizeof(float)) == 0);
  77         }
  78
  79         ad[0] = nan(nan_format);
  80         ATF_REQUIRE(isnan(ad[0]));
  81         ad[1] = strtod(nan_str, &end);
  82         ATF_REQUIRE(end == nan_str + strlen(nan_str));
  83         ATF_REQUIRE(sscanf(nan_str, "%le", &ad[2]) == 1);
  84         ATF_REQUIRE(memcmp(&ad[0], &ad[1], sizeof(double)) == 0);
  85         ATF_REQUIRE(memcmp(&ad[1], &ad[2], sizeof(double)) == 0);
  86         if (*nan_format == '\0') {
  87                 /* nan("") == strtod("nan") */
  88                 ad[3] = strtod("nan", NULL);
  89                 ATF_REQUIRE(memcmp(&ad[2], &ad[3], sizeof(double)) == 0);
  90         }
  91
  92         ald[0] = nanl(nan_format);
  93         ATF_REQUIRE(isnan(ald[0]));
  94         ald[1] = strtold(nan_str, &end);
  95         ATF_REQUIRE(end == nan_str + strlen(nan_str));
  96         ATF_REQUIRE(sscanf(nan_str, "%Le", &ald[2]) == 1);
  97         ATF_REQUIRE(memcmp(&ald[0], &ald[1], sizeof(long double)) == 0);
  98         ATF_REQUIRE(memcmp(&ald[1], &ald[2], sizeof(long double)) == 0);
  99         if (*nan_format == '\0') {
 100                 /* nanl("") == strtold("nan") */
 101                 ald[3] = strtold("nan", NULL);
 102                 ATF_REQUIRE(memcmp(&ald[2], &ald[3], sizeof(long double)) == 0);
 103         }
 104 }
 105
 106 ATF_TC_WITHOUT_HEAD(nan);
 107 ATF_TC_BODY(nan, tc)
 108 {
 109         /* Die if a signalling NaN is returned */
 110         feenableexcept(FE_INVALID);
 111
 112         testnan("0x1234");
 113         testnan("");
 114 }
 115
 116 ATF_TP_ADD_TCS(tp)
 117 {
 118         ATF_TP_ADD_TC(tp, nan);
 119
 120         return (atf_no_error());
 121 }