X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=tests%2Ftest-frexpl.c;h=5a5eae52c44cf94de1b32de38f62321d6ec2be69;hb=fba20d159fe191b9d67688738b7c31e769171aac;hp=f98fcf109db3d1fd12d98aa9a569dd0cd7e8139c;hpb=80689324582ff25f245c1ff52a401a59072c1928;p=gnulib.git diff --git a/tests/test-frexpl.c b/tests/test-frexpl.c index f98fcf109..5a5eae52c 100644 --- a/tests/test-frexpl.c +++ b/tests/test-frexpl.c @@ -1,10 +1,10 @@ /* Test of splitting a 'long double' into fraction and mantissa. - Copyright (C) 2007 Free Software Foundation, Inc. + Copyright (C) 2007-2013 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,160 +20,59 @@ #include -#include -#include +#include "signature.h" +SIGNATURE_CHECK (frexpl, long double, (long double, int *)); -#define ASSERT(expr) if (!(expr)) abort (); +#include -static long double -my_ldexp (long double x, int d) -{ - for (; d > 0; d--) - x *= 2.0L; - for (; d < 0; d++) - x *= 0.5L; - return x; -} +#include "fpucw.h" +#include "isnanl-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 long double +#define VOLATILE +#define ISNAN isnanl +#define INFINITY Infinityl () +#define NAN NaNl () +#define L_(literal) literal##L +#define MINUS_ZERO minus_zerol +#define MAX_EXP LDBL_MAX_EXP +#define MIN_EXP LDBL_MIN_EXP +/* On MIPS IRIX machines, LDBL_MIN_EXP is -1021, but the smallest reliable + exponent for 'long double' is -964. Similarly, on PowerPC machines, + LDBL_MIN_EXP is -1021, but the smallest reliable exponent for 'long double' + is -968. For exponents below that, the precision may be truncated to the + precision used for 'double'. */ +#ifdef __sgi +# define MIN_NORMAL_EXP (LDBL_MIN_EXP + 57) +#elif defined __ppc || defined __ppc__ || defined __powerpc || defined __powerpc__ +# define MIN_NORMAL_EXP (LDBL_MIN_EXP + 53) +#else +# define MIN_NORMAL_EXP LDBL_MIN_EXP +#endif +#define FREXP frexpl +#define RANDOM randoml +#include "test-frexp.h" int main () { - int i; - long double x; - - { /* NaN. */ - int exp = -9999; - long double mantissa; - x = 0.0L / 0.0L; - mantissa = frexpl (x, &exp); - ASSERT (mantissa != mantissa); - } - - { /* Positive infinity. */ - int exp = -9999; - long double mantissa; - x = 1.0L / 0.0L; - mantissa = frexpl (x, &exp); - ASSERT (mantissa == x); - } - - { /* Negative infinity. */ - int exp = -9999; - long double mantissa; - x = -1.0L / 0.0L; - mantissa = frexpl (x, &exp); - ASSERT (mantissa == x); - } - - { /* Positive zero. */ - int exp = -9999; - long double mantissa; - x = 0.0L; - mantissa = frexpl (x, &exp); - ASSERT (exp == 0); - ASSERT (mantissa == x); - } - - { /* Negative zero. */ - int exp = -9999; - long double mantissa; - x = -0.0L; - mantissa = frexpl (x, &exp); - ASSERT (exp == 0); - ASSERT (mantissa == x); - } - - for (i = 1, x = 1.0L; i <= LDBL_MAX_EXP; i++, x *= 2.0L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.5L); - } - for (i = 1, x = 1.0L; i >= LDBL_MIN_EXP; i--, x *= 0.5L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.5L); - } - for (; i >= LDBL_MIN_EXP - 100 && x > 0.0L; i--, x *= 0.5L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.5L); - } - - for (i = 1, x = -1.0L; i <= LDBL_MAX_EXP; i++, x *= 2.0L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == -0.5L); - } - for (i = 1, x = -1.0L; i >= LDBL_MIN_EXP; i--, x *= 0.5L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == -0.5L); - } - for (; i >= LDBL_MIN_EXP - 100 && x < 0.0L; i--, x *= 0.5L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == -0.5L); - } - - for (i = 1, x = 1.01L; i <= LDBL_MAX_EXP; i++, x *= 2.0L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.505L); - } - for (i = 1, x = 1.01L; i >= LDBL_MIN_EXP; i--, x *= 0.5L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.505L); - } - for (; i >= LDBL_MIN_EXP - 100 && x > 0.0L; i--, x *= 0.5L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa >= 0.5L); - ASSERT (mantissa < 1.0L); - ASSERT (mantissa == my_ldexp (x, - exp)); - } - - for (i = 1, x = 1.73205L; i <= LDBL_MAX_EXP; i++, x *= 2.0L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.866025L); - } - for (i = 1, x = 1.73205L; i >= LDBL_MIN_EXP; i--, x *= 0.5L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i); - ASSERT (mantissa == 0.866025L); - } - for (; i >= LDBL_MIN_EXP - 100 && x > 0.0L; i--, x *= 0.5L) - { - int exp = -9999; - long double mantissa = frexpl (x, &exp); - ASSERT (exp == i || exp == i + 1); - ASSERT (mantissa >= 0.5L); - ASSERT (mantissa < 1.0L); - ASSERT (mantissa == my_ldexp (x, - exp)); - } + DECL_LONG_DOUBLE_ROUNDING + + BEGIN_LONG_DOUBLE_ROUNDING (); + + test_function (); return 0; }