gettimeofday: port recent C++ fix to Emacs

[gnulib.git] / lib / trigl.c

diff --git a/lib/trigl.c b/lib/trigl.c
index 1a8a791..418ded9 100644 (file)
--- a/lib/trigl.c
+++ b/lib/trigl.c
@@ -1,5 +1,5 @@
 /* Quad-precision floating point argument reduction.
-   Copyright (C) 1999, 2007, 2009, 2010 Free Software Foundation, Inc.
+   Copyright (C) 1999, 2007, 2009-2013 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Jakub Jelinek <jj@ultra.linux.cz>
 
@@ -19,9 +19,13 @@
 #include <config.h>
 
 /* Specification.  */
-#include <math.h>
+#include "trigl.h"
 
 #include <float.h>
+#include <math.h>
+
+/* Code based on glibc/sysdeps/ieee754/ldbl-128/e_rem_pio2l.c
+   and           glibc/sysdeps/ieee754/dbl-64/k_rem_pio2.c.  */
 
 /* Table of constants for 2/pi, 5628 hexadecimal digits of 2/pi */
 static const int two_over_pi[] = {
@@ -455,7 +459,7 @@ kernel_rem_pio2 (double *x, double *y, int e0, int nx, int prec,
 
   jz = jk;
 recompute:
-  /* distill q[] into iq[] reversingly */
+  /* distill q[] into iq[] in reverse order */
   for (i = 0, j = jz, z = q[jz]; j > 0; i++, j--)
     {
       fw = (double) ((int) (twon24 * z));