X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=tests%2Ftest-floorf2.c;h=2d15a47be033ddbf7b784be89ab10207d003de43;hb=46f5f314f34a08c9305758482d7d2fdb0e999d09;hp=d668487e52f0d98286f9a154903ae271ec400608;hpb=441aa3044f43e5572f58c354f01e6bc070acd5c7;p=gnulib.git diff --git a/tests/test-floorf2.c b/tests/test-floorf2.c index d668487e5..2d15a47be 100644 --- a/tests/test-floorf2.c +++ b/tests/test-floorf2.c @@ -1,5 +1,5 @@ /* Test of rounding towards negative infinity. - Copyright (C) 2007-2008 Free Software Foundation, Inc. + Copyright (C) 2007-2014 Free Software Foundation, Inc. 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 @@ -16,6 +16,9 @@ /* Written by Bruno Haible , 2007. */ +/* When this test fails on some platform, build it together with the gnulib + module 'fprintf-posix' for optimal debugging output. */ + #include #include @@ -24,21 +27,15 @@ #include #include #include -#include #include "isnanf-nolibm.h" +#include "macros.h" -#define ASSERT(expr) \ - do \ - { \ - if (!(expr)) \ - { \ - fprintf (stderr, "%s:%d: assertion failed\n", __FILE__, __LINE__); \ - fflush (stderr); \ - abort (); \ - } \ - } \ - while (0) +/* MSVC with option -fp:strict refuses to compile constant initializers that + contain floating-point operations. Pacify this compiler. */ +#ifdef _MSC_VER +# pragma fenv_access (off) +#endif /* The reference implementation, taken from lib/floor.c. */ @@ -72,8 +69,12 @@ floorf_reference (DOUBLE x) if (z > L_(0.0)) { + /* For 0 < x < 1, return +0.0 even if the current rounding mode is + FE_DOWNWARD. */ + if (z < L_(1.0)) + z = L_(0.0); /* Avoid rounding errors for values near 2^k, where k >= MANT_DIG-1. */ - if (z < TWO_MANT_DIG) + else if (z < TWO_MANT_DIG) { /* Round to the next integer (nearest or up or down, doesn't matter). */ z += TWO_MANT_DIG; @@ -85,6 +86,9 @@ floorf_reference (DOUBLE x) } else if (z < L_(0.0)) { + /* Work around ICC's desire to optimize denormal floats to 0. */ + if (z > -FLT_MIN) + return L_(-1.0); /* Avoid rounding errors for values near -2^k, where k >= MANT_DIG-1. */ if (z > - TWO_MANT_DIG) { @@ -131,8 +135,10 @@ check (float x) return 0; else { +#if GNULIB_TEST_FPRINTF_POSIX fprintf (stderr, "floorf %g(%a) = %g(%a) or %g(%a)?\n", x, x, reference, reference, result, result); +#endif return 1; } }