diff options
author | PikalaxALT <pikalaxalt@gmail.com> | 2018-01-03 17:39:24 -0700 |
---|---|---|
committer | PikalaxALT <pikalaxalt@gmail.com> | 2018-01-03 17:39:24 -0700 |
commit | a6c1ed4716cf02626ea035beb6dd4a921642ba80 (patch) | |
tree | ef582c1b52819e27bdd16097ec03b69799d04ede /newlib/libm/mathfp/s_ldexp.c | |
parent | f6c9a624fa8a6878a7fb2b02f55e4990a20feb59 (diff) |
Use libc from agbcc instead of standalone newlib\nYou must have AGBCC commit 80d029caec189587f8b9294b6c8a5a489b8f5f88 in order to compile pmd_red.gbalibc
Diffstat (limited to 'newlib/libm/mathfp/s_ldexp.c')
-rw-r--r-- | newlib/libm/mathfp/s_ldexp.c | 123 |
1 files changed, 0 insertions, 123 deletions
diff --git a/newlib/libm/mathfp/s_ldexp.c b/newlib/libm/mathfp/s_ldexp.c deleted file mode 100644 index f68e8e4..0000000 --- a/newlib/libm/mathfp/s_ldexp.c +++ /dev/null @@ -1,123 +0,0 @@ - -/* @(#)z_ldexp.c 1.0 98/08/13 */ - -/* -FUNCTION - <<ldexp>>, <<ldexpf>>---load exponent - -INDEX - ldexp -INDEX - ldexpf - -ANSI_SYNOPSIS - #include <math.h> - double ldexp(double <[val]>, int <[exp]>); - float ldexpf(float <[val]>, int <[exp]>); - -TRAD_SYNOPSIS - #include <math.h> - - double ldexp(<[val]>, <[exp]>) - double <[val]>; - int <[exp]>; - - float ldexpf(<[val]>, <[exp]>) - float <[val]>; - int <[exp]>; - -DESCRIPTION -<<ldexp>> calculates the value -@ifinfo -<[val]> times 2 to the power <[exp]>. -@end ifinfo -@tex -$val\times 2^{exp}$. -@end tex -<<ldexpf>> is identical, save that it takes and returns <<float>> -rather than <<double>> values. - -RETURNS -<<ldexp>> returns the calculated value. - -Underflow and overflow both set <<errno>> to <<ERANGE>>. -On underflow, <<ldexp>> and <<ldexpf>> return 0.0. -On overflow, <<ldexp>> returns plus or minus <<HUGE_VAL>>. - -PORTABILITY -<<ldexp>> is ANSI, <<ldexpf>> is an extension. - -*/ - -/****************************************************************** - * 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" - -#ifndef _DOUBLE_IS_32BITS - -double -_DEFUN (ldexp, (double, int), - double d _AND - int e) -{ - int exp; - __uint32_t hd; - - GET_HIGH_WORD (hd, d); - - /* Check for special values and then scale d by e. */ - switch (numtest (d)) - { - case NAN: - errno = EDOM; - break; - - case INF: - errno = ERANGE; - break; - - case 0: - break; - - default: - exp = (hd & 0x7ff00000) >> 20; - exp += e; - - if (exp > (DBL_MAX_EXP + 1023)) - { - errno = ERANGE; - d = z_infinity.d; - } - else if (exp < DBL_MIN_EXP) - { - errno = ERANGE; - d = -z_infinity.d; - } - else - { - hd &= 0x800fffff; - hd |= exp << 20; - SET_HIGH_WORD (d, hd); - } - } - - return (d); -} - -#endif /* _DOUBLE_IS_32BITS */ |