PR libfortran/19302

	* intrinsics/c99_functions.c (nextafterf): Special-case infinite
	numbers.

git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/trunk@95353 138bc75d-0d04-0410-961f-82ee72b054a4

diff --git a/libgfortran/ChangeLog b/libgfortran/ChangeLog
index de81b00..9c0561b 100644 (file)
--- a/libgfortran/ChangeLog
+++ b/libgfortran/ChangeLog
@@ -1,3 +1,9 @@
+2005-02-21  Eric Botcazou  <ebotcazou@libertysurf.fr>
+
+       PR libfortran/19302
+       * intrinsics/c99_functions.c (nextafterf): Special-case infinite
+       numbers.
+
 2005-02-21  Steven G. Kargl  <kargls@comcast.net>
 
        * io/write.c (output_float):  Typo in error meesage.

diff --git a/libgfortran/intrinsics/c99_functions.c b/libgfortran/intrinsics/c99_functions.c
index 7dd9591..f7dc977 100644 (file)
--- a/libgfortran/intrinsics/c99_functions.c
+++ b/libgfortran/intrinsics/c99_functions.c
@@ -218,6 +218,8 @@ nextafterf(float x, float y)
     return x + y;
   if (x == y)
     return x;
+  if (!isfinite (x))
+    return x > 0 ? __FLT_MAX__ : - __FLT_MAX__;
 
   /* absx = fabsf (x);  */
   absx = (x < 0.0) ? -x : x;