1 | /* s_ceilf.c -- float version of s_ceil.c. |
2 | * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. |
3 | */ |
4 | |
5 | /* |
6 | * ==================================================== |
7 | * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. |
8 | * |
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 |
12 | * is preserved. |
13 | * ==================================================== |
14 | */ |
15 | |
16 | #define NO_MATH_REDIRECT |
17 | #include <math.h> |
18 | #include <math_private.h> |
19 | #include <libm-alias-float.h> |
20 | #include <math-use-builtins.h> |
21 | |
22 | float |
23 | __ceilf (float x) |
24 | { |
25 | #if USE_CEILF_BUILTIN |
26 | return __builtin_ceilf (x); |
27 | #else |
28 | /* Use generic implementation. */ |
29 | int32_t i0, j0; |
30 | uint32_t i; |
31 | |
32 | GET_FLOAT_WORD (i0, x); |
33 | j0 = ((i0 >> 23) & 0xff) - 0x7f; |
34 | if (j0 < 23) |
35 | { |
36 | if (j0 < 0) |
37 | { |
38 | /* return 0 * sign (x) if |x| < 1 */ |
39 | if (i0 < 0) |
40 | i0 = 0x80000000; |
41 | else if (i0 != 0) |
42 | i0 = 0x3f800000; |
43 | } |
44 | else |
45 | { |
46 | i = (0x007fffff) >> j0; |
47 | if ((i0 & i) == 0) |
48 | return x; /* x is integral */ |
49 | if (i0 > 0) |
50 | i0 += (0x00800000) >> j0; |
51 | i0 &= (~i); |
52 | } |
53 | } |
54 | else |
55 | { |
56 | if (__glibc_unlikely (j0 == 0x80)) |
57 | return x + x; /* inf or NaN */ |
58 | else |
59 | return x; /* x is integral */ |
60 | } |
61 | SET_FLOAT_WORD (x, i0); |
62 | return x; |
63 | #endif /* ! USE_CEILF_BUILTIN */ |
64 | } |
65 | #ifndef __ceilf |
66 | libm_alias_float (__ceil, ceil) |
67 | #endif |
68 | |