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_pow.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_pow.c')
-rw-r--r-- | newlib/libm/mathfp/s_pow.c | 146 |
1 files changed, 0 insertions, 146 deletions
diff --git a/newlib/libm/mathfp/s_pow.c b/newlib/libm/mathfp/s_pow.c deleted file mode 100644 index 7c0a38a..0000000 --- a/newlib/libm/mathfp/s_pow.c +++ /dev/null @@ -1,146 +0,0 @@ - -/* @(#)z_pow.c 1.0 98/08/13 */ - -/* -FUNCTION - <<pow>>, <<powf>>---x to the power y -INDEX - pow -INDEX - powf - - -ANSI_SYNOPSIS - #include <math.h> - double pow(double <[x]>, double <[y]>); - float pow(float <[x]>, float <[y]>); - -TRAD_SYNOPSIS - #include <math.h> - double pow(<[x]>, <[y]>); - double <[x]>, <[y]>; - - float pow(<[x]>, <[y]>); - float <[x]>, <[y]>; - -DESCRIPTION - <<pow>> and <<powf>> calculate <[x]> raised to the exp1.0nt <[y]>. - @tex - (That is, $x^y$.) - @end tex - -RETURNS - On success, <<pow>> and <<powf>> return the value calculated. - - When the argument values would produce overflow, <<pow>> - returns <<HUGE_VAL>> and set <<errno>> to <<ERANGE>>. If the - argument <[x]> passed to <<pow>> or <<powf>> is a negative - noninteger, and <[y]> is also not an integer, then <<errno>> - is set to <<EDOM>>. If <[x]> and <[y]> are both 0, then - <<pow>> and <<powf>> return <<1>>. - - You can modify error handling for these functions using <<matherr>>. - -PORTABILITY - <<pow>> is ANSI C. <<powf>> is an extension. */ - -#include <float.h> -#include "fdlibm.h" -#include "zmath.h" - -#ifndef _DOUBLE_IS_32BITS - -double pow (double x, double y) -{ - double d, t, r = 1.0; - int n, k, sign = 0; - __uint32_t px; - - GET_HIGH_WORD (px, x); - - k = modf (y, &d); - if (k == 0.0) - { - if (modf (ldexp (y, -1), &t)) - sign = 0; - else - sign = 1; - } - - if (x == 0.0 && y <= 0.0) - errno = EDOM; - - else if ((t = y * log (fabs (x))) >= BIGX) - { - errno = ERANGE; - if (px & 0x80000000) - { - if (!k) - { - errno = EDOM; - x = 0.0; - } - else if (sign) - x = -z_infinity.d; - else - x = z_infinity.d; - } - - else - x = z_infinity.d; - } - - else if (t < SMALLX) - { - errno = ERANGE; - x = 0.0; - } - - else - { - if ( k && fabs(d) <= 32767 ) - { - n = (int) d; - - if (sign = (n < 0)) - n = -n; - - while ( n > 0 ) - { - if ((unsigned int) n % 2) - r *= x; - x *= x; - n = (unsigned int) n / 2; - } - - if (sign) - r = 1.0 / r; - - return r; - } - - else - { - if ( px & 0x80000000 ) - { - if ( !k ) - { - errno = EDOM; - return 0.0; - } - } - - x = exp (t); - - if ( sign ) - { - px ^= 0x80000000; - SET_HIGH_WORD (x, px); - } - } - } - - return x; -} - -#endif _DOUBLE_IS_32BITS |