1 /* @(#)s_ceil.c 5.1 93/09/24 */
3 * ====================================================
4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
6 * Developed at SunPro, a Sun Microsystems, Inc. business.
7 * Permission to use, copy, modify, and distribute this
8 * software is freely granted, provided that this notice
10 * ====================================================
15 * Return x rounded toward -inf to integral value
21 #include <math_private.h>
22 #include <libm-alias-double.h>
29 EXTRACT_WORDS (i0
, i1
, x
);
30 j0
= ((i0
>> 20) & 0x7ff) - 0x3ff;
35 /* return 0*sign(x) if |x|<1 */
38 i0
= 0x80000000; i1
= 0;
40 else if ((i0
| i1
) != 0)
42 i0
= 0x3ff00000; i1
= 0;
47 i
= (0x000fffff) >> j0
;
48 if (((i0
& i
) | i1
) == 0)
49 return x
; /* x is integral */
51 i0
+= (0x00100000) >> j0
;
58 return x
+ x
; /* inf or NaN */
60 return x
; /* x is integral */
64 i
= ((uint32_t) (0xffffffff)) >> (j0
- 20);
66 return x
; /* x is integral */
73 j
= i1
+ (1 << (52 - j0
));
75 i0
+= 1; /* got a carry */
81 INSERT_WORDS (x
, i0
, i1
);
85 libm_alias_double (__ceil
, ceil
)