Retro68/gcc/newlib/libm/mathfp/sf_ldexp.c
2012-03-27 01:51:53 +02:00

82 lines
1.5 KiB
C

/* @(#)z_ldexpf.c 1.0 98/08/13 */
/******************************************************************
* ldexp
*
* Input:
* d - a floating point value
* e - an exponent value
*
* Output:
* A floating point value f such that f = d * 2 ^ e.
*
* Description:
* This function creates a floating point number f such that
* f = d * 2 ^ e.
*
*****************************************************************/
#include <float.h>
#include "fdlibm.h"
#include "zmath.h"
#define FLOAT_EXP_OFFS 127
float
_DEFUN (ldexpf, (float, int),
float d _AND
int e)
{
int exp;
__int32_t wd;
GET_FLOAT_WORD (wd, d);
/* Check for special values and then scale d by e. */
switch (numtestf (wd))
{
case NAN:
errno = EDOM;
break;
case INF:
errno = ERANGE;
break;
case 0:
break;
default:
exp = (wd & 0x7f800000) >> 23;
exp += e;
if (exp > FLT_MAX_EXP + FLOAT_EXP_OFFS)
{
errno = ERANGE;
d = z_infinity_f.f;
}
else if (exp < FLT_MIN_EXP + FLOAT_EXP_OFFS)
{
errno = ERANGE;
d = -z_infinity_f.f;
}
else
{
wd &= 0x807fffff;
wd |= exp << 23;
SET_FLOAT_WORD (d, wd);
}
}
return (d);
}
#ifdef _DOUBLE_IS_32BITS
double ldexp (double x, int e)
{
return (double) ldexpf ((float) x, e);
}
#endif /* defined(_DOUBLE_IS_32BITS) */