Back to index

glibc  2.9
Functions | Variables
s_modf.c File Reference
#include "math.h"
#include "math_private.h"
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions

double __modf (double x, double *iptr)

Variables

static double one = 1.0

Function Documentation

double __modf ( double  x,
double*  iptr 
)

Definition at line 39 of file s_modf.c.

{
       int32_t i0,i1,j0;
       u_int32_t i;
       EXTRACT_WORDS(i0,i1,x);
       j0 = ((i0>>20)&0x7ff)-0x3ff;       /* exponent of x */
       if(j0<20) {                 /* integer part in high x */
           if(j0<0) {                     /* |x|<1 */
               INSERT_WORDS(*iptr,i0&0x80000000,0);     /* *iptr = +-0 */
              return x;
           } else {
              i = (0x000fffff)>>j0;
              if(((i0&i)|i1)==0) {        /* x is integral */
                  *iptr = x;
                  INSERT_WORDS(x,i0&0x80000000,0);      /* return +-0 */
                  return x;
              } else {
                  INSERT_WORDS(*iptr,i0&(~i),0);
                  return x - *iptr;
              }
           }
       } else if (j0>51) {         /* no fraction part */
           *iptr = x*one;
           /* We must handle NaNs separately.  */
           if (j0 == 0x400 && ((i0 & 0xfffff) | i1))
             return x*one;
           INSERT_WORDS(x,i0&0x80000000,0);      /* return +-0 */
           return x;
       } else {                    /* fraction part in low x */
           i = ((u_int32_t)(0xffffffff))>>(j0-20);
           if((i1&i)==0) {         /* x is integral */
              *iptr = x;
              INSERT_WORDS(x,i0&0x80000000,0);   /* return +-0 */
              return x;
           } else {
               INSERT_WORDS(*iptr,i0,i1&(~i));
              return x - *iptr;
           }
       }
}

Here is the call graph for this function:


Variable Documentation

double one = 1.0 [static]

Definition at line 33 of file s_modf.c.