#if defined(LIBM_SCCS) && !defined(lint)
__RCSID("$NetBSD: s_scalbnf.c,v 1.7 1999/07/02 15:37:43 simonb Exp $");
#endif
#define __DARWIN_UNIX03 1
#include "sys/cdefs.h"
#include "math.h"
#include "math_private.h"
static const float
two25 = 0x1.0p+25,
twom25 = 0x1.0p-25,
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;
}