diff options
author | Alex Auvolat <alex@adnab.me> | 2017-05-03 20:37:59 +0200 |
---|---|---|
committer | Alex Auvolat <alex@adnab.me> | 2017-05-03 20:37:59 +0200 |
commit | 1161e1d8be014945266017cb0ce735537a287677 (patch) | |
tree | 118f2201a5e12f79aefb404295794eed0d52cd6d /src/lib/libc/fdlibm/s_ceil.c | |
parent | 0b583122fb6cfcff991c54836d37cb3958c343b1 (diff) | |
download | kogata-1161e1d8be014945266017cb0ce735537a287677.tar.gz kogata-1161e1d8be014945266017cb0ce735537a287677.zip |
Truetype fonts
Diffstat (limited to 'src/lib/libc/fdlibm/s_ceil.c')
-rw-r--r-- | src/lib/libc/fdlibm/s_ceil.c | 78 |
1 files changed, 78 insertions, 0 deletions
diff --git a/src/lib/libc/fdlibm/s_ceil.c b/src/lib/libc/fdlibm/s_ceil.c new file mode 100644 index 0000000..5a2802c --- /dev/null +++ b/src/lib/libc/fdlibm/s_ceil.c @@ -0,0 +1,78 @@ + +/* @(#)s_ceil.c 1.3 95/01/18 */ +/* + * ==================================================== + * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. + * + * Developed at SunSoft, a Sun Microsystems, Inc. business. + * Permission to use, copy, modify, and distribute this + * software is freely granted, provided that this notice + * is preserved. + * ==================================================== + */ + +/* + * ceil(x) + * Return x rounded toward -inf to integral value + * Method: + * Bit twiddling. + * Exception: + * Inexact flag raised if x not equal to ceil(x). + */ + +#include "fdlibm.h" + +#ifdef __STDC__ +static const double huge = 1.0e300; +#else +static double huge = 1.0e300; +#endif + +#ifdef __STDC__ + double ceil(double x) +#else + double ceil(x) + double x; +#endif +{ + int i0,i1,j0; + unsigned i,j; + i0 = __HI(x); + i1 = __LO(x); + j0 = ((i0>>20)&0x7ff)-0x3ff; + if(j0<20) { + if(j0<0) { /* raise inexact if x != 0 */ + if(huge+x>0.0) {/* return 0*sign(x) if |x|<1 */ + if(i0<0) {i0=0x80000000;i1=0;} + else if((i0|i1)!=0) { i0=0x3ff00000;i1=0;} + } + } else { + i = (0x000fffff)>>j0; + if(((i0&i)|i1)==0) return x; /* x is integral */ + if(huge+x>0.0) { /* raise inexact flag */ + if(i0>0) i0 += (0x00100000)>>j0; + i0 &= (~i); i1=0; + } + } + } else if (j0>51) { + if(j0==0x400) return x+x; /* inf or NaN */ + else return x; /* x is integral */ + } else { + i = ((unsigned)(0xffffffff))>>(j0-20); + if((i1&i)==0) return x; /* x is integral */ + if(huge+x>0.0) { /* raise inexact flag */ + if(i0>0) { + if(j0==20) i0+=1; + else { + j = i1 + (1<<(52-j0)); + if((int)j<i1) i0+=1; /* got a carry */ + i1 = j; + } + } + i1 &= (~i); + } + } + __HI(x) = i0; + __LO(x) = i1; + return x; +} |