Back to index

glibc  2.9
Functions | Variables
s_modfl.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

long double __modfl (long double x, long double *iptr)

Variables

static long double one = 1.0

Function Documentation

long double __modfl ( long double  x,
long double*  iptr 
)

Definition at line 42 of file s_modfl.c.

{
       int64_t i0,i1,j0;
       u_int64_t i;
       GET_LDOUBLE_WORDS64(i0,i1,x);
       j0 = ((i0>>48)&0x7fff)-0x3fff;     /* exponent of x */
       if(j0<48) {                 /* integer part in high x */
           if(j0<0) {                     /* |x|<1 */
              /* *iptr = +-0 */
               SET_LDOUBLE_WORDS64(*iptr,i0&0x8000000000000000ULL,0);
              return x;
           } else {
              i = (0x0000ffffffffffffLL)>>j0;
              if(((i0&i)|i1)==0) {        /* x is integral */
                  *iptr = x;
                  /* return +-0 */
                  SET_LDOUBLE_WORDS64(x,i0&0x8000000000000000ULL,0);
                  return x;
              } else {
                  SET_LDOUBLE_WORDS64(*iptr,i0&(~i),0);
                  return x - *iptr;
              }
           }
       } else if (j0>111) {        /* no fraction part */
           *iptr = x*one;
           /* We must handle NaNs separately.  */
           if (j0 == 0x4000 && ((i0 & 0x0000ffffffffffffLL) | i1))
             return x*one;
           /* return +-0 */
           SET_LDOUBLE_WORDS64(x,i0&0x8000000000000000ULL,0);
           return x;
       } else {                    /* fraction part in low x */
           i = -1ULL>>(j0-48);
           if((i1&i)==0) {         /* x is integral */
              *iptr = x;
              /* return +-0 */
              SET_LDOUBLE_WORDS64(x,i0&0x8000000000000000ULL,0);
              return x;
           } else {
              SET_LDOUBLE_WORDS64(*iptr,i0,i1&(~i));
              return x - *iptr;
           }
       }
}

Here is the call graph for this function:


Variable Documentation

long double one = 1.0 [static]

Definition at line 36 of file s_modfl.c.