[Commits] r25307 - in /fsf/glibc-2_19-branch/libc: ChangeLog NEWS sysdeps/ieee754/dbl-64/s_sin.c

[eglibc@xxxxxxxxxx]

Author: eglibc
Date: Fri Feb 28 00:02:01 2014
New Revision: 25307

Log:
Import glibc-2.19 for 2014-02-28

Modified:
    fsf/glibc-2_19-branch/libc/ChangeLog
    fsf/glibc-2_19-branch/libc/NEWS
    fsf/glibc-2_19-branch/libc/sysdeps/ieee754/dbl-64/s_sin.c

Modified: fsf/glibc-2_19-branch/libc/ChangeLog
==============================================================================
--- fsf/glibc-2_19-branch/libc/ChangeLog (original)
+++ fsf/glibc-2_19-branch/libc/ChangeLog Fri Feb 28 00:02:01 2014
@@ -1,3 +1,10 @@
+2014-02-27  Siddhesh Poyarekar  <siddhesh@xxxxxxxxxx>
+
+	[BZ #16623]
+	* sysdeps/ieee754/dbl-64/s_sin.c (__sin): Preserve sign of A
+	and DA.
+	(__cos): Likewise.
+
 2014-02-12  Dylan Alex Simon  <dylan@xxxxxxxxx>
 
 	[BZ #16545]

Modified: fsf/glibc-2_19-branch/libc/NEWS
==============================================================================
--- fsf/glibc-2_19-branch/libc/NEWS (original)
+++ fsf/glibc-2_19-branch/libc/NEWS Fri Feb 28 00:02:01 2014
@@ -9,7 +9,7 @@
 
 * The following bugs are resolved with this release:
 
-  16545.
+  16545, 16623.
 
 Version 2.19
 

Modified: fsf/glibc-2_19-branch/libc/sysdeps/ieee754/dbl-64/s_sin.c
==============================================================================
--- fsf/glibc-2_19-branch/libc/sysdeps/ieee754/dbl-64/s_sin.c (original)
+++ fsf/glibc-2_19-branch/libc/sysdeps/ieee754/dbl-64/s_sin.c Fri Feb 28 00:02:01 2014
@@ -447,19 +447,21 @@
 	    }
 	  else
 	    {
+	      double t;
 	      if (a > 0)
 		{
 		  m = 1;
+		  t = a;
 		  db = da;
 		}
 	      else
 		{
 		  m = 0;
-		  a = -a;
+		  t = -a;
 		  db = -da;
 		}
-	      u.x = big + a;
-	      y = a - (u.x - big);
+	      u.x = big + t;
+	      y = t - (u.x - big);
 	      res = do_sin (u, y, db, &cor);
 	      cor = (cor > 0) ? 1.035 * cor + eps : 1.035 * cor - eps;
 	      retval = ((res == res + cor) ? ((m) ? res : -res)
@@ -671,19 +673,21 @@
 	    }
 	  else
 	    {
+	      double t;
 	      if (a > 0)
 		{
 		  m = 1;
+		  t = a;
 		  db = da;
 		}
 	      else
 		{
 		  m = 0;
-		  a = -a;
+		  t = -a;
 		  db = -da;
 		}
-	      u.x = big + a;
-	      y = a - (u.x - big);
+	      u.x = big + t;
+	      y = t - (u.x - big);
 	      res = do_sin (u, y, db, &cor);
 	      cor = (cor > 0) ? 1.035 * cor + eps : 1.035 * cor - eps;
 	      retval = ((res == res + cor) ? ((m) ? res : -res)

_______________________________________________
Commits mailing list
Commits@xxxxxxxxxx
http://eglibc.org/cgi-bin/mailman/listinfo/commits