#include <stdio.h>
#include <stdlib.h>
-#include "isnand.h"
+#include "isnand-nolibm.h"
+#include "nan.h"
+
+/* Avoid some warnings from "gcc -Wshadow".
+ This file doesn't use the exp() function. */
+#undef exp
+#define exp exponent
#define ASSERT(expr) \
do \
if (!(expr)) \
{ \
fprintf (stderr, "%s:%d: assertion failed\n", __FILE__, __LINE__); \
+ fflush (stderr); \
abort (); \
} \
} \
while (0)
-/* The Compaq (ex-DEC) C 6.4 compiler chokes on the expression 0.0 / 0.0. */
-#ifdef __DECC
-static double
-NaN ()
-{
- static double zero = 0.0;
- return zero / zero;
-}
-#else
-# define NaN() (0.0 / 0.0)
-#endif
+/* 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;
static double
my_ldexp (double x, int d)
{ /* NaN. */
int exp = -9999;
double mantissa;
- x = NaN ();
+ x = NaNd ();
mantissa = frexp (x, &exp);
ASSERT (isnand (mantissa));
}
{ /* Negative zero. */
int exp = -9999;
double mantissa;
- x = -0.0;
+ x = -zero;
mantissa = frexp (x, &exp);
ASSERT (exp == 0);
ASSERT (mantissa == x);