From 8e27e3cc453849158e069dbed0eea7b8113ab8ff Mon Sep 17 00:00:00 2001 From: Joseph Myers Date: Thu, 22 Nov 2012 14:57:22 +0000 Subject: Fix ldbl-128ibm hypotl inaccuracy for arguments with large ratio (bug 14868). --- sysdeps/ieee754/ldbl-128ibm/e_hypotl.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'sysdeps/ieee754/ldbl-128ibm/e_hypotl.c') diff --git a/sysdeps/ieee754/ldbl-128ibm/e_hypotl.c b/sysdeps/ieee754/ldbl-128ibm/e_hypotl.c index 00bfb15..ce21194 100644 --- a/sysdeps/ieee754/ldbl-128ibm/e_hypotl.c +++ b/sysdeps/ieee754/ldbl-128ibm/e_hypotl.c @@ -61,7 +61,7 @@ __ieee754_hypotl(long double x, long double y) if(hb > ha) {a=y;b=x;j=ha; ha=hb;hb=j;} else {a=x;b=y;} a = fabsl(a); /* a <- |a| */ b = fabsl(b); /* b <- |b| */ - if((ha-hb)>0x3c0000000000000LL) {return a+b;} /* x/y > 2**60 */ + if((ha-hb)>0x780000000000000LL) {return a+b;} /* x/y > 2**120 */ k=0; kld = 1.0L; if(ha > 0x5f30000000000000LL) { /* a>2**500 */ -- cgit v1.1