From a2a2059fd77eb8c30c0534cbc766c21072a9829a Mon Sep 17 00:00:00 2001
From: Bud Davis <bdavis9659@comcast.net>
Date: Sat, 12 Jun 2004 17:59:30 +0000
Subject: re PR libfortran/15292 (libgfortran depends on c99 functionality
 round and roundf)

	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.

From-SVN: r83036
---
 libgfortran/intrinsics/c99_functions.c | 90 ++++++++++++++++++++++++++++++++++
 1 file changed, 90 insertions(+)
 create mode 100644 libgfortran/intrinsics/c99_functions.c

(limited to 'libgfortran/intrinsics')

diff --git a/libgfortran/intrinsics/c99_functions.c b/libgfortran/intrinsics/c99_functions.c
new file mode 100644
index 0000000..9328387
--- /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
+
-- 
cgit v1.1