/* Test of splitting a double into fraction and mantissa.
- Copyright (C) 2007 Free Software Foundation, Inc.
+ Copyright (C) 2007-2010 Free Software Foundation, Inc.
- This program is free software; you can redistribute it and/or modify
+ This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
- the Free Software Foundation; either version 2, or (at your option)
- any later version.
+ the Free Software Foundation; either version 3 of the License, or
+ (at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
- along with this program; if not, write to the Free Software Foundation,
- Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. */
+ along with this program. If not, see <http://www.gnu.org/licenses/>. */
/* Written by Bruno Haible <bruno@clisp.org>, 2007. */
#include <math.h>
+#include "signature.h"
+SIGNATURE_CHECK (frexp, double, (double, int *));
+
#include <float.h>
-#include <stdlib.h>
-#define ASSERT(expr) if (!(expr)) abort ();
+#include "isnand-nolibm.h"
+#include "nan.h"
+#include "macros.h"
+
+/* Avoid some warnings from "gcc -Wshadow".
+ This file doesn't use the exp() function. */
+#undef exp
+#define exp exponent
+
+/* 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 = 0.0 / 0.0;
+ x = NaNd ();
mantissa = frexp (x, &exp);
- ASSERT (mantissa != mantissa);
+ ASSERT (isnand (mantissa));
}
{ /* Positive infinity. */
mantissa = frexp (x, &exp);
ASSERT (exp == 0);
ASSERT (mantissa == x);
+ ASSERT (!signbit (mantissa));
}
{ /* Negative zero. */
int exp = -9999;
double mantissa;
- x = -0.0;
+ x = -zero;
mantissa = frexp (x, &exp);
ASSERT (exp == 0);
ASSERT (mantissa == x);
+ ASSERT (signbit (mantissa));
}
for (i = 1, x = 1.0; i <= DBL_MAX_EXP; i++, x *= 2.0)