/* Quad-precision floating point argument reduction.
- Copyright (C) 1999, 2007, 2009-2011 Free Software Foundation, Inc.
+ Copyright (C) 1999, 2007, 2009-2014 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Jakub Jelinek <jj@ultra.linux.cz>
#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[] = {
0xa2f983, 0x6e4e44, 0x1529fc, 0x2757d1, 0xf534dd, 0xc0db62,
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));