diff options
Diffstat (limited to 'gcc-4.8/libgcc/config/rs6000/ibm-ldouble.c')
-rw-r--r-- | gcc-4.8/libgcc/config/rs6000/ibm-ldouble.c | 11 |
1 files changed, 10 insertions, 1 deletions
diff --git a/gcc-4.8/libgcc/config/rs6000/ibm-ldouble.c b/gcc-4.8/libgcc/config/rs6000/ibm-ldouble.c index 28e02e995..574e395f8 100644 --- a/gcc-4.8/libgcc/config/rs6000/ibm-ldouble.c +++ b/gcc-4.8/libgcc/config/rs6000/ibm-ldouble.c @@ -188,7 +188,16 @@ __gcc_qdiv (double a, double b, double c, double d) || nonfinite (t)) return t; - /* Finite nonzero result requires corrections to the highest order term. */ + /* Finite nonzero result requires corrections to the highest order + term. These corrections require the low part of c * t to be + exactly represented in double. */ + if (fabs (a) <= 0x1p-969) + { + a *= 0x1p106; + b *= 0x1p106; + c *= 0x1p106; + d *= 0x1p106; + } s = c * t; /* (s,sigma) = c*t exactly. */ w = -(-b + d * t); /* Written to get fnmsub for speed, but not |