Avoid a compilation warning.

[gnulib.git] / lib / trigl.c

diff --git a/lib/trigl.c b/lib/trigl.c
index ed9fd51..6f9dcf5 100644 (file)
--- a/lib/trigl.c
+++ b/lib/trigl.c
@@ -3,19 +3,18 @@
    This file is part of the GNU C Library.
    Contributed by Jakub Jelinek <jj@ultra.linux.cz>
 
-   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
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    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.  */
+   You should have received a copy of the GNU General Public License
+   along with this program.  If not, see <http://www.gnu.org/licenses/>.  */
 
 #include <config.h>
 
@@ -206,7 +205,7 @@ ieee754_rem_pio2l (long double x, long double *y)
   int exp, n;
 
   if (x >= -0.78539816339744830961566084581987572104929234984377
-      && x < 0.78539816339744830961566084581987572104929234984377)
+      && x <= 0.78539816339744830961566084581987572104929234984377)
     /* x in <-pi/4, pi/4> */
     {
       y[0] = x;
@@ -214,9 +213,7 @@ ieee754_rem_pio2l (long double x, long double *y)
       return 0;
     }
 
-  if (x >= 2.35619449019234492884698253745962716314787704953131
-      && x < 2.35619449019234492884698253745962716314787704953131)
-    if (x > 0)
+  if (x > 0 && x < 2.35619449019234492884698253745962716314787704953131)
       {
        /* 113 + 93 bit PI is ok */
        z = x - PI_2_1;
@@ -224,7 +221,8 @@ ieee754_rem_pio2l (long double x, long double *y)
        y[1] = (z - y[0]) - PI_2_1t;
        return 1;
       }
-    else
+
+  if (x < 0 && x > -2.35619449019234492884698253745962716314787704953131)
       {
        /* 113 + 93 bit PI is ok */
        z = x + PI_2_1;
@@ -233,7 +231,7 @@ ieee754_rem_pio2l (long double x, long double *y)
        return -1;
       }
 
-  if (x + x == x || x != x)    /* x is +=oo or NaN */
+  if (x + x == x)      /* x is ±oo */
     {
       y[0] = x - x;
       y[1] = y[0];
@@ -423,7 +421,7 @@ static const double PIo2[] = {
 static const double zero = 0.0, one = 1.0, two24 = 1.67772160000000000000e+07, /* 0x41700000, 0x00000000 */
   twon24 = 5.96046447753906250000e-08; /* 0x3E700000, 0x00000000 */
 
-int
+static int
 kernel_rem_pio2 (double *x, double *y, int e0, int nx, int prec,
                 const int *ipio2)
 {