Back to index

lightning-sunbird  0.9+nobinonly
Functions | Variables
s_ceil.c File Reference
#include "fdlibm.h"

Go to the source code of this file.

Functions

double fd_ceil (double x)

Variables

static double really_big = 1.0e300

Function Documentation

double fd_ceil ( double  x)

Definition at line 72 of file s_ceil.c.

{
        fd_twoints u;
       int i0,i1,j0;
       unsigned i,j;
        u.d = x;
       i0 =  __HI(u);
       i1 =  __LO(u);
       j0 = ((i0>>20)&0x7ff)-0x3ff;
       if(j0<20) {
           if(j0<0) {       /* raise inexact if x != 0 */
              if(really_big+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(really_big+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(really_big+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);
           }
       }
        u.d = x;
       __HI(u) = i0;
       __LO(u) = i1;
        x = u.d;
       return x;
}

Here is the caller graph for this function:


Variable Documentation

double really_big = 1.0e300 [static]

Definition at line 66 of file s_ceil.c.