summary refs log tree commit diff
path: root/math/s_fdim.c
diff options
context:
space:
mode:
authorJoseph Myers <joseph@codesourcery.com>2013-08-21 19:56:48 +0000
committerJoseph Myers <joseph@codesourcery.com>2013-08-21 19:56:48 +0000
commitacd06bb11f6d6436e15d0c7608fc7ea6008c224f (patch)
treec789341c0b86f19efa34ba8dbb09bef9705c2817 /math/s_fdim.c
parentc0c3f78afb6070721848574e2e5dff5cfa20e28d (diff)
downloadglibc-acd06bb11f6d6436e15d0c7608fc7ea6008c224f.tar.gz
glibc-acd06bb11f6d6436e15d0c7608fc7ea6008c224f.tar.xz
glibc-acd06bb11f6d6436e15d0c7608fc7ea6008c224f.zip
Fix fdim handling of infinities (bug 15797).
Diffstat (limited to 'math/s_fdim.c')
-rw-r--r--math/s_fdim.c8
1 files changed, 4 insertions, 4 deletions
diff --git a/math/s_fdim.c b/math/s_fdim.c
index 2f97948b2e..f8fd80490d 100644
--- a/math/s_fdim.c
+++ b/math/s_fdim.c
@@ -26,16 +26,16 @@ __fdim (double x, double y)
   int clsx = fpclassify (x);
   int clsy = fpclassify (y);
 
-  if (clsx == FP_NAN || clsy == FP_NAN
-      || (y < 0 && clsx == FP_INFINITE && clsy == FP_INFINITE))
-    /* Raise invalid flag.  */
+  if (clsx == FP_NAN || clsy == FP_NAN)
+    /* Raise invalid flag for signaling but not quiet NaN.  */
     return x - y;
 
   if (x <= y)
     return 0.0;
 
   double r = x - y;
-  if (fpclassify (r) == FP_INFINITE)
+  if (fpclassify (r) == FP_INFINITE
+      && clsx != FP_INFINITE && clsy != FP_INFINITE)
     __set_errno (ERANGE);
 
   return r;