1/* @(#)s_ceil.c 5.1 93/09/24 */ 2/* 3 * ==================================================== 4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 5 * 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 9 * is preserved. 10 * ==================================================== 11 */ 12 13#ifndef lint
| 1/* @(#)s_ceil.c 5.1 93/09/24 */ 2/* 3 * ==================================================== 4 * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. 5 * 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 9 * is preserved. 10 * ==================================================== 11 */ 12 13#ifndef lint
|
14static char rcsid[] = "$FreeBSD: head/lib/msun/src/s_ceil.c 97409 2002-05-28 17:51:46Z alfred $";
| 14static char rcsid[] = "$FreeBSD: head/lib/msun/src/s_ceil.c 97413 2002-05-28 18:15:04Z alfred $";
|
15#endif 16 17/* 18 * ceil(x) 19 * Return x rounded toward -inf to integral value 20 * Method: 21 * Bit twiddling. 22 * Exception: 23 * Inexact flag raised if x not equal to ceil(x). 24 */ 25 26#include "math.h" 27#include "math_private.h" 28 29static const double huge = 1.0e300; 30
| 15#endif 16 17/* 18 * ceil(x) 19 * Return x rounded toward -inf to integral value 20 * Method: 21 * Bit twiddling. 22 * Exception: 23 * Inexact flag raised if x not equal to ceil(x). 24 */ 25 26#include "math.h" 27#include "math_private.h" 28 29static const double huge = 1.0e300; 30
|
31 double __generic_ceil(double x)
| 31double 32__generic_ceil(double x)
|
32{ 33 int32_t i0,i1,j0; 34 u_int32_t i,j; 35 EXTRACT_WORDS(i0,i1,x); 36 j0 = ((i0>>20)&0x7ff)-0x3ff; 37 if(j0<20) { 38 if(j0<0) { /* raise inexact if x != 0 */ 39 if(huge+x>0.0) {/* return 0*sign(x) if |x|<1 */
--- 32 unchanged lines hidden --- | 33{ 34 int32_t i0,i1,j0; 35 u_int32_t i,j; 36 EXTRACT_WORDS(i0,i1,x); 37 j0 = ((i0>>20)&0x7ff)-0x3ff; 38 if(j0<20) { 39 if(j0<0) { /* raise inexact if x != 0 */ 40 if(huge+x>0.0) {/* return 0*sign(x) if |x|<1 */
--- 32 unchanged lines hidden --- |