Handle +inf/+inf correctly.
authordrepper <drepper>
Sun, 26 Sep 2004 10:10:28 +0000 (10:10 +0000)
committerdrepper <drepper>
Sun, 26 Sep 2004 10:10:28 +0000 (10:10 +0000)
sysdeps/generic/s_fdim.c
sysdeps/generic/s_fdimf.c
sysdeps/generic/s_fdiml.c
sysdeps/i386/i686/fpu/s_fdim.S
sysdeps/i386/i686/fpu/s_fdimf.S
sysdeps/i386/i686/fpu/s_fdiml.S
sysdeps/powerpc/fpu/s_fdim.c
sysdeps/powerpc/fpu/s_fdimf.c
sysdeps/x86_64/fpu/s_fdiml.S

index 201f936..5804e63 100644 (file)
@@ -1,5 +1,5 @@
 /* Return positive difference between arguments.
-   Copyright (C) 1997 Free Software Foundation, Inc.
+   Copyright (C) 1997, 2004 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
@@ -31,7 +31,7 @@ __fdim (double x, double y)
     /* Raise invalid flag.  */
     return x - y;
 
-  return x < y ? 0 : x - y;
+  return x <= y ? 0 : x - y;
 }
 weak_alias (__fdim, fdim)
 #ifdef NO_LONG_DOUBLE
index 64d54b7..2f3ce30 100644 (file)
@@ -1,5 +1,5 @@
 /* Return positive difference between arguments.
-   Copyright (C) 1997 Free Software Foundation, Inc.
+   Copyright (C) 1997, 2004 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
@@ -31,6 +31,6 @@ __fdimf (float x, float y)
     /* Raise invalid flag.  */
     return x - y;
 
-  return x < y ? 0 : x - y;
+  return x <= y ? 0 : x - y;
 }
 weak_alias (__fdimf, fdimf)
index 83049ae..70246ba 100644 (file)
@@ -1,5 +1,5 @@
 /* Return positive difference between arguments.
-   Copyright (C) 1997 Free Software Foundation, Inc.
+   Copyright (C) 1997, 2004 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
@@ -31,6 +31,6 @@ __fdiml (long double x, long double y)
     /* Raise invalid flag.  */
     return x - y;
 
-  return x < y ? 0 : x - y;
+  return x <= y ? 0 : x - y;
 }
 weak_alias (__fdiml, fdiml)
index e610973..30ecff4 100644 (file)
@@ -1,5 +1,5 @@
 /* Compute positive difference.
-   Copyright (C) 1997, 1998 Free Software Foundation, Inc.
+   Copyright (C) 1997, 1998, 2004 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
@@ -28,12 +28,14 @@ ENTRY(__fdim)
        fucomi  %st(1), %st
        jp      1f
 
-       fsubrp  %st, %st(1)
+       jc      3f
+       fstp    %st(1)
        fldz
-       fcomi   %st(1), %st
-       fcmovb  %st(1), %st
        jmp     2f
 
+3:     fsubrp  %st, %st(1)
+       ret
+
 1:     fucomi  %st(0), %st
        fcmovnu %st(1), %st
 2:     fstp    %st(1)
index a22cbe9..888df14 100644 (file)
@@ -1,5 +1,5 @@
 /* Compute positive difference.
-   Copyright (C) 1997, 1998 Free Software Foundation, Inc.
+   Copyright (C) 1997, 1998, 2004 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
@@ -28,12 +28,14 @@ ENTRY(__fdimf)
        fucomi  %st(1), %st
        jp      1f
 
-       fsubrp  %st, %st(1)
+       jc      3f
+       fstp    %st(1)
        fldz
-       fcomi   %st(1), %st
-       fcmovb  %st(1), %st
        jmp     2f
 
+3:     fsubrp  %st, %st(1)
+       ret
+
 1:     fucomi  %st(0), %st
        fcmovnu %st(1), %st
 2:     fstp    %st(1)
index fa3c51e..cb0e26e 100644 (file)
@@ -1,5 +1,5 @@
 /* Compute positive difference.
-   Copyright (C) 1997, 1998 Free Software Foundation, Inc.
+   Copyright (C) 1997, 1998, 2004 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
@@ -28,12 +28,14 @@ ENTRY(__fdiml)
        fucomi  %st(1), %st
        jp      1f
 
-       fsubrp  %st, %st(1)
+       jc      3f
+       fstp    %st(1)
        fldz
-       fcomi   %st(1), %st
-       fcmovb  %st(1), %st
        jmp     2f
 
+3:     fsubrp  %st, %st(1)
+       ret
+
 1:     fucomi  %st(0), %st
        fcmovnu %st(1), %st
 2:     fstp    %st(1)
index 165e2ff..2b767ad 100644 (file)
@@ -1,5 +1,5 @@
 /* Return positive difference between arguments.
-   Copyright (C) 1997 Free Software Foundation, Inc.
+   Copyright (C) 1997, 2004 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -22,7 +22,7 @@
 double
 __fdim (double x, double y)
 {
-  return x < y ? 0 : x - y;
+  return x <= y ? 0 : x - y;
 }
 weak_alias (__fdim, fdim)
 #ifdef NO_LONG_DOUBLE
index 997ec89..a27c1e4 100644 (file)
@@ -1,5 +1,5 @@
 /* Return positive difference between arguments.
-   Copyright (C) 1997 Free Software Foundation, Inc.
+   Copyright (C) 1997, 2004 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -22,6 +22,6 @@
 float
 __fdimf (float x, float y)
 {
-  return x < y ? 0 : x - y;
+  return x <= y ? 0 : x - y;
 }
 weak_alias (__fdimf, fdimf)
index 3460b0f..d63ca00 100644 (file)
@@ -1,5 +1,5 @@
 /* Compute positive difference.
-   Copyright (C) 1997, 1998, 2002 Free Software Foundation, Inc.
+   Copyright (C) 1997, 1998, 2002, 2004 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Ulrich Drepper <drepper@cygnus.com>, 1997.
 
@@ -28,12 +28,14 @@ ENTRY(__fdiml)
        fucomi  %st(1), %st
        jp      1f
 
-       fsubrp  %st, %st(1)
+       jc      3f
+       fstp    %st(1)
        fldz
-       fcomi   %st(1), %st
-       fcmovb  %st(1), %st
        jmp     2f
 
+3:     fsubrp  %st, %st(1)
+       ret
+
 1:     fucomi  %st(0), %st
        fcmovnu %st(1), %st
 2:     fstp    %st(1)