Back to index

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

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 43 of file s_modfl.c.

{
       int32_t i0,i1,j0;
       u_int32_t i,se;
       GET_LDOUBLE_WORDS(se,i0,i1,x);
       j0 = (se&0x7fff)-0x3fff;    /* exponent of x */
       if(j0<32) {                 /* integer part in high x */
           if(j0<0) {                     /* |x|<1 */
               SET_LDOUBLE_WORDS(*iptr,se&0x8000,0,0);  /* *iptr = +-0 */
              return x;
           } else {
              i = (0x7fffffff)>>j0;
              if(((i0&i)|i1)==0) {        /* x is integral */
                  *iptr = x;
                  SET_LDOUBLE_WORDS(x,se&0x8000,0,0);   /* return +-0 */
                  return x;
              } else {
                  SET_LDOUBLE_WORDS(*iptr,se,i0&(~i),0);
                  return x - *iptr;
              }
           }
       } else if (j0>63) {         /* no fraction part */
           *iptr = x*one;
           /* We must handle NaNs separately.  */
           if (j0 == 0x4000 && ((i0 & 0x7fffffff) | i1))
             return x*one;
           SET_LDOUBLE_WORDS(x,se&0x8000,0,0);   /* return +-0 */
           return x;
       } else {                    /* fraction part in low x */
           i = ((u_int32_t)(0x7fffffff))>>(j0-32);
           if((i1&i)==0) {         /* x is integral */
              *iptr = x;
              SET_LDOUBLE_WORDS(x,se&0x8000,0,0);       /* return +-0 */
              return x;
           } else {
              SET_LDOUBLE_WORDS(*iptr,se,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 37 of file s_modfl.c.