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