* sysdeps/alpha/fpu/s_ceil.c: Rewrite without branches.
[kopensolaris-gnu/glibc.git] / sysdeps / alpha / fpu / s_rintf.c
index d5d019d..9e4cbd1 100644 (file)
@@ -1,4 +1,4 @@
-/* Copyright (C) 2000 Free Software Foundation, Inc.
+/* Copyright (C) 2000, 2007 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Richard Henderson.
 
 float
 __rintf (float x)
 {
-  if (isless (fabsf (x), 16777216.0f)) /* 1 << FLT_MANT_DIG */
-    {
-      /* Note that Alpha S_Floating is stored in registers in a
-        restricted T_Floating format, so we don't even need to
-        convert back to S_Floating in the end.  The initial
-        conversion to T_Floating is needed to handle denormals.  */
-
-      float tmp1, tmp2, new_x;
-
-      __asm ("cvtst/s %3,%2\n\t"
-#ifdef _IEEE_FP_INEXACT
-            "cvttq/svid %2,%1\n\t"
-#else
-            "cvttq/svd %2,%1\n\t"
-#endif
-            "cvtqt/d %1,%0\n\t"
-            : "=f"(new_x), "=&f"(tmp1), "=&f"(tmp2)
-            : "f"(x));
-
-      /* rint(-0.1) == -0, and in general we'll always have the same
-        sign as our input.  */
-      x = copysignf(new_x, x);
-    }
-  return x;
+  float two23 = copysignf (0x1.0p23, x);
+  float r;
+
+  r = x + two23;
+  r = r - two23;
+
+  /* rint(-0.1) == -0, and in general we'll always have the same sign
+     as our input.  */
+  return copysign (r, x);
 }
 
 weak_alias (__rintf, rintf)