Prefer https to http for gnu.org and fsf.org URLs

[thirdparty/glibc.git] / sysdeps / ieee754 / dbl-64 / s_round.c

diff --git a/sysdeps/ieee754/dbl-64/s_round.c b/sysdeps/ieee754/dbl-64/s_round.c
index 770f8e63863d37655a07112c76823d2fb5d481b9..1793e575e9ae1b122b066f6e38a0ee2b8542e75d 100644 (file)
--- a/sysdeps/ieee754/dbl-64/s_round.c
+++ b/sysdeps/ieee754/dbl-64/s_round.c
@@ -1,5 +1,5 @@
 /* Round double to integer away from zero.
-   Copyright (C) 1997-2015 Free Software Foundation, Inc.
+   Copyright (C) 1997-2019 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
@@ -15,21 +15,20 @@
 
    You should have received a copy of the GNU Lesser General Public
    License along with the GNU C Library; if not, see
-   <http://www.gnu.org/licenses/>.  */
+   <https://www.gnu.org/licenses/>.  */
 
+#define NO_MATH_REDIRECT
 #include <math.h>
 
 #include <math_private.h>
-
-
-static const double huge = 1.0e300;
+#include <libm-alias-double.h>
 
 
 double
 __round (double x)
 {
   int32_t i0, j0;
-  u_int32_t i1;
+  uint32_t i1;
 
   EXTRACT_WORDS (i0, i1, x);
   j0 = ((i0 >> 20) & 0x7ff) - 0x3ff;
@@ -37,8 +36,6 @@ __round (double x)
     {
       if (j0 < 0)
        {
-         math_force_eval (huge + x);
-
          i0 &= 0x80000000;
          if (j0 == -1)
            i0 |= 0x3ff00000;
@@ -46,13 +43,11 @@ __round (double x)
        }
       else
        {
-         u_int32_t i = 0x000fffff >> j0;
+         uint32_t i = 0x000fffff >> j0;
          if (((i0 & i) | i1) == 0)
            /* X is integral.  */
            return x;
-         math_force_eval (huge + x);
 
-         /* Raise inexact if x != 0.  */
          i0 += 0x00080000 >> j0;
          i0 &= ~i;
          i1 = 0;
@@ -68,15 +63,12 @@ __round (double x)
     }
   else
     {
-      u_int32_t i = 0xffffffff >> (j0 - 20);
+      uint32_t i = 0xffffffff >> (j0 - 20);
       if ((i1 & i) == 0)
        /* X is integral.  */
        return x;
 
-      math_force_eval (huge + x);
-
-      /* Raise inexact if x != 0.  */
-      u_int32_t j = i1 + (1 << (51 - j0));
+      uint32_t j = i1 + (1 << (51 - j0));
       if (j < i1)
        i0 += 1;
       i1 = j;
@@ -86,8 +78,4 @@ __round (double x)
   INSERT_WORDS (x, i0, i1);
   return x;
 }
-weak_alias (__round, round)
-#ifdef NO_LONG_DOUBLE
-strong_alias (__round, __roundl)
-weak_alias (__round, roundl)
-#endif
+libm_alias_double (__round, round)