#include <sys/cdefs.h>
#if defined(LIBM_SCCS) && !defined(lint)
__RCSID("$NetBSD: s_scalbnf.c,v 1.7 1999/07/02 15:37:43 simonb Exp $");
#endif
#include "math.h"
#include "math_private.h"
static const float
two25 = 3.355443200e+07,
twom25 = 2.9802322388e-08,
huge = 1.0e+30,
tiny = 1.0e-30;
float scalbnf (float x, int n)
{
int32_t k,ix;
GET_FLOAT_WORD(ix,x);
k = (ix&0x7f800000)>>23;
if (k==0) {
if ((ix&0x7fffffff)==0) return x;
x *= two25;
GET_FLOAT_WORD(ix,x);
k = ((ix&0x7f800000)>>23) - 25;
if (n< -50000) return tiny*x;
}
if (k==0xff) return x+x;
k = k+n;
if (k > 0xfe) return huge*copysignf(huge,x);
if (k > 0)
{SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); return x;}
if (k <= -25) {
if (n > 50000)
return huge*copysignf(huge,x);
else return tiny*copysignf(tiny,x);
}
k += 25;
SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23));
return x*twom25;
}