X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=tests%2Ftest-frexp.c;h=f7a10f5bc1f5f24cbd2bc2a40c527b21617e8005;hb=1276a2c5f24c0c932426aca9c899fa524d2443f2;hp=9ce989b06db40540d786b5636ca73c4477abe3e7;hpb=75d146361fe980e1e1bbb9639edc6fbed5345e82;p=gnulib.git diff --git a/tests/test-frexp.c b/tests/test-frexp.c index 9ce989b06..f7a10f5bc 100644 --- a/tests/test-frexp.c +++ b/tests/test-frexp.c @@ -1,10 +1,10 @@ /* Test of splitting a double into fraction and mantissa. - Copyright (C) 2007 Free Software Foundation, Inc. + Copyright (C) 2007-2014 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 @@ -12,8 +12,7 @@ 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 . */ /* Written by Bruno Haible , 2007. */ @@ -21,191 +20,48 @@ #include -#include -#include -#include - -#include "isnan.h" - -#define ASSERT(expr) \ - do \ - { \ - if (!(expr)) \ - { \ - fprintf (stderr, "%s:%d: assertion failed\n", __FILE__, __LINE__); \ - abort (); \ - } \ - } \ - while (0) +#include "signature.h" +SIGNATURE_CHECK (frexp, double, (double, int *)); -/* 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 +#include -static double -my_ldexp (double x, int d) -{ - for (; d > 0; d--) - x *= 2.0; - for (; d < 0; d++) - x *= 0.5; - return x; -} +#include "isnand-nolibm.h" +#include "minus-zero.h" +#include "infinity.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 + +#undef INFINITY +#undef NAN + +#define DOUBLE double +/* The use of 'volatile' guarantees that excess precision bits are dropped + when dealing with denormalized numbers. It is necessary on x86 systems + where double-floats are not IEEE compliant by default, to avoid that the + results become platform and compiler option dependent. 'volatile' is a + portable alternative to gcc's -ffloat-store option. */ +#define VOLATILE volatile +#define ISNAN isnand +#define INFINITY Infinityd () +#define NAN NaNd () +#define L_(literal) literal +#define MINUS_ZERO minus_zerod +#define MAX_EXP DBL_MAX_EXP +#define MIN_EXP DBL_MIN_EXP +#define MIN_NORMAL_EXP DBL_MIN_EXP +#define FREXP frexp +#define RANDOM randomd +#include "test-frexp.h" int main () { - int i; - /* The use of 'volatile' guarantees that excess precision bits are dropped - when dealing with denormalized numbers. It is necessary on x86 systems - where double-floats are not IEEE compliant by default, to avoid that the - results become platform and compiler option dependent. 'volatile' is a - portable alternative to gcc's -ffloat-store option. */ - volatile double x; - - { /* NaN. */ - int exp = -9999; - double mantissa; - x = NaN (); - mantissa = frexp (x, &exp); - ASSERT (isnan (mantissa)); - } - - { /* Positive infinity. */ - int exp = -9999; - double mantissa; - x = 1.0 / 0.0; - mantissa = frexp (x, &exp); - ASSERT (mantissa == x); - } - - { /* Negative infinity. */ - int exp = -9999; - double mantissa; - x = -1.0 / 0.0; - mantissa = frexp (x, &exp); - ASSERT (mantissa == x); - } - - { /* Positive zero. */ - int exp = -9999; - double mantissa; - x = 0.0; - mantissa = frexp (x, &exp); - ASSERT (exp == 0); - ASSERT (mantissa == x); - ASSERT (!signbit (mantissa)); - } - - { /* Negative zero. */ - int exp = -9999; - double mantissa; - x = -0.0; - 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) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.5); - } - for (i = 1, x = 1.0; i >= DBL_MIN_EXP; i--, x *= 0.5) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.5); - } - for (; i >= DBL_MIN_EXP - 100 && x > 0.0; i--, x *= 0.5) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.5); - } - - for (i = 1, x = -1.0; i <= DBL_MAX_EXP; i++, x *= 2.0) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == -0.5); - } - for (i = 1, x = -1.0; i >= DBL_MIN_EXP; i--, x *= 0.5) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == -0.5); - } - for (; i >= DBL_MIN_EXP - 100 && x < 0.0; i--, x *= 0.5) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == -0.5); - } - - for (i = 1, x = 1.01; i <= DBL_MAX_EXP; i++, x *= 2.0) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.505); - } - for (i = 1, x = 1.01; i >= DBL_MIN_EXP; i--, x *= 0.5) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.505); - } - for (; i >= DBL_MIN_EXP - 100 && x > 0.0; i--, x *= 0.5) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa >= 0.5); - ASSERT (mantissa < 1.0); - ASSERT (mantissa == my_ldexp (x, - exp)); - } - - for (i = 1, x = 1.73205; i <= DBL_MAX_EXP; i++, x *= 2.0) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.866025); - } - for (i = 1, x = 1.73205; i >= DBL_MIN_EXP; i--, x *= 0.5) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.866025); - } - for (; i >= DBL_MIN_EXP - 100 && x > 0.0; i--, x *= 0.5) - { - int exp = -9999; - double mantissa = frexp (x, &exp); - ASSERT (exp == i || exp == i + 1); - ASSERT (mantissa >= 0.5); - ASSERT (mantissa < 1.0); - ASSERT (mantissa == my_ldexp (x, - exp)); - } + test_function (); return 0; }