about summary refs log tree commit diff
path: root/src/math/s_rintf.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/math/s_rintf.c')
-rw-r--r--src/math/s_rintf.c45
1 files changed, 0 insertions, 45 deletions
diff --git a/src/math/s_rintf.c b/src/math/s_rintf.c
deleted file mode 100644
index c441870d..00000000
--- a/src/math/s_rintf.c
+++ /dev/null
@@ -1,45 +0,0 @@
-/* s_rintf.c -- float version of s_rint.c.
- * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
- */
-
-/*
- * ====================================================
- * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
- *
- * Developed at SunPro, a Sun Microsystems, Inc. business.
- * Permission to use, copy, modify, and distribute this
- * software is freely granted, provided that this notice
- * is preserved.
- * ====================================================
- */
-
-#include <math.h>
-#include "math_private.h"
-
-static const float
-TWO23[2]={
-  8.3886080000e+06, /* 0x4b000000 */
- -8.3886080000e+06, /* 0xcb000000 */
-};
-
-float
-rintf(float x)
-{
-        int32_t i0,j0,sx;
-        volatile float w,t;     /* volatile works around gcc bug */
-        GET_FLOAT_WORD(i0,x);
-        sx = (i0>>31)&1;
-        j0 = ((i0>>23)&0xff)-0x7f;
-        if(j0<23) {
-            if(j0<0) {
-                if((i0&0x7fffffff)==0) return x;
-                w = TWO23[sx]+x;
-                t =  w-TWO23[sx];
-                return t;
-            }
-            w = TWO23[sx]+x;
-            return w-TWO23[sx];
-        }
-        if(j0==0x80) return x+x;        /* inf or NaN */
-        else return x;                  /* x is integral */
-}