Back to index

glibc  2.9
Functions | Variables
s_modfl.c File Reference
#include "math.h"
#include "math_private.h"
#include <math_ldbl_opt.h>

Go to the source code of this file.

Functions

long double __modfl (long double x, long double *iptr)
 long_double_symbol (libc, __modfl, modfl)

Variables

static long double one = 1.0

Function Documentation

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

Definition at line 43 of file s_modfl.c.

{
       int64_t i0,i1,j0;
       u_int64_t i;
       GET_LDOUBLE_WORDS64(i0,i1,x);
       i1 &= 0x000fffffffffffffLL;
       j0 = ((i0>>52)&0x7ff)-0x3ff;       /* exponent of x */
       if(j0<52) {                 /* integer part in high x */
           if(j0<0) {                     /* |x|<1 */
              /* *iptr = +-0 */
               SET_LDOUBLE_WORDS64(*iptr,i0&0x8000000000000000ULL,0);
              return x;
           } else {
              i = (0x000fffffffffffffLL)>>j0;
              if(((i0&i)|(i1&0x7fffffffffffffffLL))==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>103) {        /* no fraction part */
           *iptr = x*one;
           /* We must handle NaNs separately.  */
           if (j0 == 0x400 && ((i0 & 0x000fffffffffffffLL) | i1))
             return x*one;
           /* return +-0 */
           SET_LDOUBLE_WORDS64(x,i0&0x8000000000000000ULL,0);
           return x;
       } else {                    /* fraction part in low x */
           i = -1ULL>>(j0-52);
           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:

long_double_symbol ( libc  ,
__modfl  ,
modfl   
)

Variable Documentation

long double one = 1.0 [static]

Definition at line 37 of file s_modfl.c.