1 /* origin: FreeBSD /usr/src/lib/msun/src/s_rintf.c */
3 * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com.
6 * ====================================================
7 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
9 * Developed at SunPro, a Sun Microsystems, Inc. business.
10 * Permission to use, copy, modify, and distribute this
11 * software is freely granted, provided that this notice
13 * ====================================================
20 8.3886080000e+06, /* 0x4b000000 */
21 -8.3886080000e+06, /* 0xcb000000 */
29 GET_FLOAT_WORD(i0, x);
31 j0 = ((i0>>23)&0xff) - 0x7f;
34 if ((i0&0x7fffffff) == 0)
36 STRICT_ASSIGN(float, w, TWO23[sx] + x);
38 GET_FLOAT_WORD(i0, t);
39 SET_FLOAT_WORD(t, (i0&0x7fffffff)|(sx<<31));
42 STRICT_ASSIGN(float, w, TWO23[sx] + x);
46 return x+x; /* inf or NaN */
47 return x; /* x is integral */