about summary refs log tree commit diff
path: root/sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2013-05-09 21:30:08 +0000
committerJoseph Myers <joseph@codesourcery.com>2013-05-09 21:30:08 +0000
commited41ffefc3f947f14d565ea8d239ff2d31f6a7fe (patch)
tree6a83a53b6d441694683731cf4c3568a37581e388 /sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c
parentd0213cd0b629cae95ed9607fcf26793a3a0d686c (diff)
downloadglibc-ed41ffefc3f947f14d565ea8d239ff2d31f6a7fe.tar.gz
glibc-ed41ffefc3f947f14d565ea8d239ff2d31f6a7fe.tar.xz
glibc-ed41ffefc3f947f14d565ea8d239ff2d31f6a7fe.zip
Fix ldbl-128ibm cos range reduction near pi/2 (bug 15359).
Diffstat (limited to 'sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c')
-rw-r--r--sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c10
1 files changed, 5 insertions, 5 deletions
diff --git a/sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c b/sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c
index 692ae24930..6a72d6a853 100644
--- a/sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c
+++ b/sysdeps/ieee754/ldbl-128ibm/e_rem_pio2l.c
@@ -185,13 +185,13 @@ static const int32_t two_over_pi[] = {
 };
 
 static const long double c[] = {
-/* 93 bits of pi/2 */
+/* 106 bits of pi/2 */
 #define PI_2_1 c[0]
- 1.57079632679489661923132169155131424e+00L, /* 3fff921fb54442d18469898cc5100000 */
+ 0x1.921fb54442d18469898cc517018p+0L,
 
 /* pi/2 - PI_2_1 */
 #define PI_2_1t c[1]
- 8.84372056613570112025531863263659260e-29L, /* 3fa1c06e0e68948127044533e63a0106 */
+ 0x3.839a252049c1114cf98e804178p-108L,
 };
 
 int32_t __ieee754_rem_pio2l(long double x, long double *y)
@@ -216,7 +216,7 @@ int32_t __ieee754_rem_pio2l(long double x, long double *y)
     {
       if (hx > 0)
 	{
-	  /* 113 + 93 bit PI is ok */
+	  /* 106 + 106 bit PI is ok */
 	  z = x - PI_2_1;
 	  y[0] = z - PI_2_1t;
 	  y[1] = (z - y[0]) - PI_2_1t;
@@ -224,7 +224,7 @@ int32_t __ieee754_rem_pio2l(long double x, long double *y)
 	}
       else
         {
-	  /* 113 + 93 bit PI is ok */
+	  /* 106 + 106 bit PI is ok */
 	  z = x + PI_2_1;
 	  y[0] = z + PI_2_1t;
 	  y[1] = (z - y[0]) + PI_2_1t;