#include <sys/cdefs.h>
#if defined(LIBM_SCCS) && !defined(lint)
__RCSID("$NetBSD: s_modff.c,v 1.6 1999/07/02 15:37:43 simonb Exp $");
#endif
#include "math.h"
#include "math_private.h"
static const float one = 1.0;
float modff(float x, float *iptr)
{
int32_t i0,j0;
u_int32_t i;
if (x != x) {
*iptr = x;
return x;
}
GET_FLOAT_WORD(i0,x);
j0 = ((i0>>23)&0xff)-0x7f;
if(j0<23) {
if(j0<0) {
SET_FLOAT_WORD(*iptr,i0&0x80000000);
return x;
} else {
i = (0x007fffff)>>j0;
if((i0&i)==0) {
u_int32_t ix;
*iptr = x;
GET_FLOAT_WORD(ix,x);
SET_FLOAT_WORD(x,ix&0x80000000);
return x;
} else {
SET_FLOAT_WORD(*iptr,i0&(~i));
return x - *iptr;
}
}
} else {
u_int32_t ix;
*iptr = x*one;
GET_FLOAT_WORD(ix,x);
SET_FLOAT_WORD(x,ix&0x80000000);
return x;
}
}