Mercurial > hg > octave-kai > gnulib-hg

diff lib/sinl.c @ 8563:c3d334708928

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Better support of signalling NaNs.
author Bruno Haible <bruno@clisp.org>
date Mon, 26 Mar 2007 22:26:05 +0000
parents c115624bf99a
children 678640901735
line wrap: on
line diff
--- a/lib/sinl.c
+++ b/lib/sinl.c
@@ -52,6 +52,7 @@
 #include "trigl.h"
 #include "trigl.c"
 #include "sincosl.c"
+#include "isnanl.h"
 
 long double
 sinl (long double x)
@@ -59,13 +60,17 @@
   long double y[2], z = 0.0L;
   int n;
 
+  /* sinl(NaN) is NaN */
+  if (isnanl (x))
+    return x;
+
   /* |x| ~< pi/4 */
   if (x >= -0.7853981633974483096156608458198757210492 &&
       x <= 0.7853981633974483096156608458198757210492)
     return kernel_sinl (x, z, 0);
 
-    /* sinl(Inf or NaN) is NaN, sinl(0) is 0 */
-  else if (x + x == x || x != x)
+    /* sinl(Inf) is NaN, sinl(0) is 0 */
+  else if (x + x == x)
     return x - x;		/* NaN */
 
   /* argument reduction needed */