src/lib/libm/src/s_modf.c - diff

Return to s_modf.c CVS log

Up to [cvs.NetBSD.org] / src / lib / libm / src

Please note that diffs are not public domain; they are subject to the copyright notices on the relevant files.

version 1.11, 2002/05/26 22:01:57

version 1.11.30.1, 2008/05/18 12:30:39

Line 33 static const double one = 1.0;

double

modf(double x, double *iptr)

{

int32_t i0,i1,j0;

int32_t i0,i1,jj0;

u_int32_t i;

EXTRACT_WORDS(i0,i1,x);

j0 = ((i0>>20)&0x7ff)-0x3ff; /* exponent of x */

jj0 = ((i0>>20)&0x7ff)-0x3ff; /* exponent of x */

if(j0<20) { /* integer part in high x */

if(jj0<20) { /* integer part in high x */

if(j0<0) { /* |x|<1 */

if(jj0<0) { /* |x|<1 */

INSERT_WORDS(*iptr,i0&0x80000000,0); /* *iptr = +-0 */

return x;

} else {

i = (0x000fffff)>>j0;

i = (0x000fffff)>>jj0;

if(((i0&i)|i1)==0) { /* x is integral */

u_int32_t high;

*iptr = x;

Line 54 modf(double x, double *iptr)

return x - *iptr;

}

} else if (j0>51) { /* no fraction part */

} else if (jj0>51) { /* no fraction part */

u_int32_t high;

*iptr = x*one;

GET_HIGH_WORD(high,x);

INSERT_WORDS(x,high&0x80000000,0); /* return +-0 */

return x;

} else { /* fraction part in low x */

i = ((u_int32_t)(0xffffffff))>>(j0-20);

i = ((u_int32_t)(0xffffffff))>>(jj0-20);

if((i1&i)==0) { /* x is integral */

u_int32_t high;

*iptr = x;