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/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 }