Mercurial > hg > octave-nkf > gnulib-hg

diff lib/trigl.c @ 9923:810b08d769f8

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
fix typos in mathl 2008-04-15 Paolo Bonzini <bonzini@gnu.org> * lib/tanl.c (kernel_tanl): Rename flag to invert, initialize it to 0. * lib/trigl.c (ieee754_rem_pio2l): Fix range checks.
author Paolo Bonzini <bonzini@gnu.org>
date Tue, 15 Apr 2008 15:25:38 +0200
parents bbbbbf4cd1c5
children e8d2c6fc33ad
line wrap: on
line diff
--- a/lib/trigl.c
+++ b/lib/trigl.c
@@ -205,7 +205,7 @@
   int exp, n;
 
   if (x >= -0.78539816339744830961566084581987572104929234984377
-      && x < 0.78539816339744830961566084581987572104929234984377)
+      && x <= 0.78539816339744830961566084581987572104929234984377)
     /* x in <-pi/4, pi/4> */
     {
       y[0] = x;
@@ -213,9 +213,7 @@
       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;
@@ -223,7 +221,8 @@
 	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;