/*
- * Copyright (C) 2008 Free Software Foundation
+ * Copyright (C) 2008, 2009, 2010 Free Software Foundation, Inc.
* Written by Eric Blake
*
* This program is free software: you can redistribute it and/or modify
#include <stdlib.h>
+#include "signature.h"
+SIGNATURE_CHECK (strtod, double, (char const *, char **));
+
#include <errno.h>
#include <float.h>
#include <math.h>
-#include <stdio.h>
#include <string.h>
-#define ASSERT(expr) \
- do \
- { \
- if (!(expr)) \
- { \
- fprintf (stderr, "%s:%d: assertion failed\n", __FILE__, __LINE__); \
- /* FIXME abort ();*/status = 1; \
- } \
- } \
- while (0)
+#include "isnand-nolibm.h"
+#include "macros.h"
+
+/* Avoid requiring -lm just for fabs. */
+#define FABS(d) ((d) < 0.0 ? -(d) : (d))
+
+/* HP cc on HP-UX 10.20 has a bug with the constant expression -0.0.
+ So we use -zero instead. */
+double zero = 0.0;
int
-main ()
+main (void)
{
int status = 0;
/* Subject sequence empty or invalid. */
double result;
errno = 0;
result = strtod (input, &ptr);
- ASSERT (result == 0.5);
+ /* FIXME - gnulib's version is rather inaccurate. It would be
+ nice to guarantee an exact result, but for now, we settle for a
+ 1-ulp error. */
+ ASSERT (FABS (result - 0.5) < DBL_EPSILON);
ASSERT (ptr == input + 2);
ASSERT (errno == 0);
}
double result;
errno = 0;
result = strtod (input, &ptr);
- ASSERT (result == 0.5);
+ /* FIXME - gnulib's version is rather inaccurate. It would be
+ nice to guarantee an exact result, but for now, we settle for a
+ 1-ulp error. */
+ ASSERT (FABS (result - 0.5) < DBL_EPSILON);
ASSERT (ptr == input + 4);
ASSERT (errno == 0);
}
errno = 0;
result = strtod (input, &ptr);
ASSERT (result == 0.0);
- ASSERT (!!signbit (result) == !!signbit (-0.0)); /* IRIX 6.5, OSF/1 4.0 */
+ ASSERT (!!signbit (result) == !!signbit (-zero)); /* IRIX 6.5, OSF/1 4.0 */
ASSERT (ptr == input + 2);
ASSERT (errno == 0);
}
errno = 0;
result = strtod (input, &ptr);
ASSERT (result == 0.0);
- ASSERT (!!signbit (result) == !!signbit (-0.0)); /* MacOS X 10.3, FreeBSD 6.2, IRIX 6.5, OSF/1 4.0 */
+ ASSERT (!!signbit (result) == !!signbit (-zero)); /* MacOS X 10.3, FreeBSD 6.2, IRIX 6.5, OSF/1 4.0 */
ASSERT (ptr == input + 2); /* glibc-2.3.6, MacOS X 10.3, FreeBSD 6.2 */
ASSERT (errno == 0);
}
double result;
errno = 0;
result = strtod (input, &ptr);
- ASSERT (0.0 <= result && result <= FLT_MIN);
+ ASSERT (0.0 <= result && result <= DBL_MIN);
ASSERT (!signbit (result));
ASSERT (ptr == input + 9);
ASSERT (errno == ERANGE);
double result;
errno = 0;
result = strtod (input, &ptr);
- ASSERT (-FLT_MIN <= result && result <= 0.0);
+ ASSERT (-DBL_MIN <= result && result <= 0.0);
#if 0
/* FIXME - this is glibc bug 5995; POSIX allows returning positive
0 on negative underflow, even though quality of implementation
demands preserving the sign. Disable this test until fixed
glibc is more prevalent. */
- ASSERT (!!signbit (result) == !!signbit (-0.0)); /* glibc-2.3.6, mingw */
+ ASSERT (!!signbit (result) == !!signbit (-zero)); /* glibc-2.3.6, mingw */
#endif
ASSERT (ptr == input + 10);
ASSERT (errno == ERANGE);
result1 = strtod (input, &ptr1);
result2 = strtod (input + 1, &ptr2);
#if 1 /* All known CPUs support NaNs. */
- ASSERT (isnan (result1)); /* OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
- ASSERT (isnan (result2)); /* OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (isnand (result1)); /* OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (isnand (result2)); /* OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
# if 0
/* Sign bits of NaN is a portability sticking point, not worth
worrying about. */
result1 = strtod (input, &ptr1);
result2 = strtod (input + 1, &ptr2);
#if 1 /* All known CPUs support NaNs. */
- ASSERT (isnan (result1)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
- ASSERT (isnan (result2)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (isnand (result1)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (isnand (result2)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
ASSERT (!!signbit (result1) == !!signbit (result2));
ASSERT (ptr1 == input + 4); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, Solaris 2.5.1, mingw */
ASSERT (ptr2 == input + 4); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, Solaris 2.5.1, mingw */
result1 = strtod (input, &ptr1);
result2 = strtod (input + 1, &ptr2);
#if 1 /* All known CPUs support NaNs. */
- ASSERT (isnan (result1)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
- ASSERT (isnan (result2)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (isnand (result1)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (isnand (result2)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
# if 0
/* Sign bits of NaN is a portability sticking point, not worth
worrying about. */
errno = 0;
result = strtod (input, &ptr);
#if 1 /* All known CPUs support NaNs. */
- ASSERT (isnan (result)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (isnand (result)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
ASSERT (ptr == input + 6); /* glibc-2.3.6, MacOS X 10.3, FreeBSD 6.2, OpenBSD 4.0, AIX 5.1, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
ASSERT (errno == 0);
#else
result1 = strtod (input, &ptr1);
result2 = strtod (input + 1, &ptr2);
#if 1 /* All known CPUs support NaNs. */
- ASSERT (isnan (result1)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
- ASSERT (isnan (result2)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (isnand (result1)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (isnand (result2)); /* OpenBSD 4.0, HP-UX 11.11, IRIX 6.5, OSF/1 5.1, mingw */
# if 0
/* Sign bits of NaN is a portability sticking point, not worth
worrying about. */
char *input = malloc (m + 1);
if (input)
{
- char *ptr;
- double result;
- memset (input, '\t', m - 1);
- input[m - 1] = '1';
- input[m] = '\0';
- errno = 0;
- result = strtod (input, &ptr);
- ASSERT (result == 1.0);
- ASSERT (ptr == input + m);
- ASSERT (errno == 0);
+ char *ptr;
+ double result;
+ memset (input, '\t', m - 1);
+ input[m - 1] = '1';
+ input[m] = '\0';
+ errno = 0;
+ result = strtod (input, &ptr);
+ ASSERT (result == 1.0);
+ ASSERT (ptr == input + m);
+ ASSERT (errno == 0);
}
free (input);
}
char *input = malloc (m + 1);
if (input)
{
- char *ptr;
- double result;
- memset (input, '0', m - 1);
- input[m - 1] = '1';
- input[m] = '\0';
- errno = 0;
- result = strtod (input, &ptr);
- ASSERT (result == 1.0);
- ASSERT (ptr == input + m);
- ASSERT (errno == 0);
+ char *ptr;
+ double result;
+ memset (input, '0', m - 1);
+ input[m - 1] = '1';
+ input[m] = '\0';
+ errno = 0;
+ result = strtod (input, &ptr);
+ ASSERT (result == 1.0);
+ ASSERT (ptr == input + m);
+ ASSERT (errno == 0);
}
free (input);
}
char *input = malloc (m + 1);
if (input)
{
- char *ptr;
- double result;
- input[0] = '.';
- memset (input + 1, '0', m - 10);
- input[m - 9] = '1';
- input[m - 8] = 'e';
- input[m - 7] = '+';
- input[m - 6] = '9';
- input[m - 5] = '9';
- input[m - 4] = '9';
- input[m - 3] = '9';
- input[m - 2] = '9';
- input[m - 1] = '1';
- input[m] = '\0';
- errno = 0;
- result = strtod (input, &ptr);
- ASSERT (result == 1.0); /* MacOS X 10.3, FreeBSD 6.2, NetBSD 3.0, OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
- ASSERT (ptr == input + m); /* OSF/1 5.1 */
- ASSERT (errno == 0); /* MacOS X 10.3, FreeBSD 6.2, NetBSD 3.0, OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
+ char *ptr;
+ double result;
+ input[0] = '.';
+ memset (input + 1, '0', m - 10);
+ input[m - 9] = '1';
+ input[m - 8] = 'e';
+ input[m - 7] = '+';
+ input[m - 6] = '9';
+ input[m - 5] = '9';
+ input[m - 4] = '9';
+ input[m - 3] = '9';
+ input[m - 2] = '9';
+ input[m - 1] = '1';
+ input[m] = '\0';
+ errno = 0;
+ result = strtod (input, &ptr);
+ ASSERT (result == 1.0); /* MacOS X 10.3, FreeBSD 6.2, NetBSD 3.0, OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (ptr == input + m); /* OSF/1 5.1 */
+ ASSERT (errno == 0); /* MacOS X 10.3, FreeBSD 6.2, NetBSD 3.0, OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
}
free (input);
}
char *input = malloc (m + 1);
if (input)
{
- char *ptr;
- double result;
- input[0] = '1';
- memset (input + 1, '0', m - 9);
- input[m - 8] = 'e';
- input[m - 7] = '-';
- input[m - 6] = '9';
- input[m - 5] = '9';
- input[m - 4] = '9';
- input[m - 3] = '9';
- input[m - 2] = '9';
- input[m - 1] = '1';
- input[m] = '\0';
- errno = 0;
- result = strtod (input, &ptr);
- ASSERT (result == 1.0); /* MacOS X 10.3, FreeBSD 6.2, NetBSD 3.0, OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
- ASSERT (ptr == input + m);
- ASSERT (errno == 0); /* MacOS X 10.3, FreeBSD 6.2, NetBSD 3.0, OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
+ char *ptr;
+ double result;
+ input[0] = '1';
+ memset (input + 1, '0', m - 9);
+ input[m - 8] = 'e';
+ input[m - 7] = '-';
+ input[m - 6] = '9';
+ input[m - 5] = '9';
+ input[m - 4] = '9';
+ input[m - 3] = '9';
+ input[m - 2] = '9';
+ input[m - 1] = '1';
+ input[m] = '\0';
+ errno = 0;
+ result = strtod (input, &ptr);
+ ASSERT (result == 1.0); /* MacOS X 10.3, FreeBSD 6.2, NetBSD 3.0, OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
+ ASSERT (ptr == input + m);
+ ASSERT (errno == 0); /* MacOS X 10.3, FreeBSD 6.2, NetBSD 3.0, OpenBSD 4.0, IRIX 6.5, OSF/1 5.1, mingw */
}
free (input);
}
char *input = malloc (m + 1);
if (input)
{
- char *ptr;
- double result;
- input[0] = '-';
- input[1] = '0';
- input[2] = 'e';
- input[3] = '1';
- memset (input + 4, '0', m - 3);
- input[m] = '\0';
- errno = 0;
- result = strtod (input, &ptr);
- ASSERT (result == 0.0);
- ASSERT (!!signbit (result) == !!signbit (-0.0)); /* IRIX 6.5, OSF/1 4.0 */
- ASSERT (ptr == input + m);
- ASSERT (errno == 0);
+ char *ptr;
+ double result;
+ input[0] = '-';
+ input[1] = '0';
+ input[2] = 'e';
+ input[3] = '1';
+ memset (input + 4, '0', m - 3);
+ input[m] = '\0';
+ errno = 0;
+ result = strtod (input, &ptr);
+ ASSERT (result == 0.0);
+ ASSERT (!!signbit (result) == !!signbit (-zero)); /* IRIX 6.5, OSF/1 4.0 */
+ ASSERT (ptr == input + m);
+ ASSERT (errno == 0);
}
free (input);
}
/* Rounding. */
/* TODO - is it worth some tests of rounding for typical IEEE corner
cases, such as .5 ULP rounding up to the smallest denormal and
- not causing underflow, or FLT_MIN - .5 ULP not causing an
+ not causing underflow, or DBL_MIN - .5 ULP not causing an
infinite loop? */
return status;
projects / gnulib.git / blobdiff