PR gfortran/15292

	* intrinsics/c99_functions.c: New file.
	* Makefile.am: Add new file.
	* configure.ac: Added test for round/roundf.
	* Makefile.in: Regenerate.
	* configure: Regenerate.
        * configure.h.in: Regenerate.


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

1 files changed, 90 insertions, 0 deletions

diff --git a/libgfortran/intrinsics/c99_functions.c b/libgfortran/intrinsics/c99_functions.c
new file mode 100644
index 00000000000..932838780b8
--- /dev/null
+++ b/libgfortran/intrinsics/c99_functions.c
@@ -0,0 +1,90 @@
+/* Implementation of various C99 functions 
+   Copyright (C) 2004 Free Software Foundation, Inc.
+
+This file is part of the GNU Fortran 95 runtime library (libgfortran).
+
+Libgfortran is free software; you can redistribute it and/or
+modify it under the terms of the GNU Lesser General Public
+License as published by the Free Software Foundation; either
+version 2.1 of the License, or (at your option) any later version.
+
+Libgfortran is distributed in the hope that it will be useful,
+but WITHOUT ANY WARRANTY; without even the implied warranty of
+MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+GNU Lesser General Public License for more details.
+
+You should have received a copy of the GNU Lesser General Public
+License along with libgfortran; see the file COPYING.LIB.  If not,
+write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
+Boston, MA 02111-1307, USA.  */
+
+#include "config.h"
+#include <sys/types.h>
+#include <math.h>
+#include "libgfortran.h"
+
+
+/* Algorithm by Steven G. Kargl.  */
+
+#ifndef HAVE_ROUND
+/* Round to nearest integral value.  If the argument is halfway between two
+   integral values then round away from zero.  */
+
+double
+round(double x)
+{
+   double t;
+   int i;
+
+   i = fpclassify(x);
+   if (i == FP_INFINITE || i == FP_NAN)
+     return (x);
+
+   if (x >= 0.0) 
+    {
+      t = ceil(x);
+      if (t - x > 0.5)
+	t -= 1.0;
+      return (t);
+    } 
+   else 
+    {
+      t = ceil(-x);
+      if (t + x > 0.5)
+	t -= 1.0;
+      return (-t);
+    }
+}
+#endif
+
+#ifndef HAVE_ROUNDF
+/* Round to nearest integral value.  If the argument is halfway between two
+   integral values then round away from zero.  */
+
+float
+roundf(float x)
+{
+   float t;
+   int i;
+
+   i = fpclassify(x);
+   if (i == FP_INFINITE || i == FP_NAN)
+     return (x);
+
+   if (x >= 0.0) 
+    {
+      t = ceilf(x);
+      if (t - x > 0.5)
+	t -= 1.0;
+      return (t);
+    } 
+   else 
+    {
+      t = ceilf(-x);
+      if (t + x > 0.5)
+	t -= 1.0;
+      return (-t);
+    }
+}
+#endif
+