summaryrefslogtreecommitdiff
path: root/newlib/libm/mathfp/sf_pow.c
diff options
context:
space:
mode:
authorPikalaxALT <pikalaxalt@gmail.com>2018-01-03 17:39:24 -0700
committerPikalaxALT <pikalaxalt@gmail.com>2018-01-03 17:39:24 -0700
commita6c1ed4716cf02626ea035beb6dd4a921642ba80 (patch)
treeef582c1b52819e27bdd16097ec03b69799d04ede /newlib/libm/mathfp/sf_pow.c
parentf6c9a624fa8a6878a7fb2b02f55e4990a20feb59 (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/sf_pow.c')
-rw-r--r--newlib/libm/mathfp/sf_pow.c107
1 files changed, 0 insertions, 107 deletions
diff --git a/newlib/libm/mathfp/sf_pow.c b/newlib/libm/mathfp/sf_pow.c
deleted file mode 100644
index 2b3bed3..0000000
--- a/newlib/libm/mathfp/sf_pow.c
+++ /dev/null
@@ -1,107 +0,0 @@
-
-/* @(#)z_powf.c 1.0 98/08/13 */
-#include <float.h>
-#include "fdlibm.h"
-#include "zmath.h"
-
-float powf (float x, float y)
-{
- float d, t, r = 1.0;
- int n, k, sign = 0;
- __int32_t px;
-
- GET_FLOAT_WORD (px, x);
-
- k = modff (y, &d);
- if (k == 0.0)
- {
- if (modff (ldexpf (y, -1), &t))
- sign = 0;
- else
- sign = 1;
- }
-
- if (x == 0.0 && y <= 0.0)
- errno = EDOM;
-
- else if ((t = y * log (fabsf (x))) >= BIGX)
- {
- errno = ERANGE;
- if (px & 0x80000000)
- {
- if (!k)
- {
- errno = EDOM;
- x = 0.0;
- }
- else if (sign)
- x = -z_infinity_f.f;
- else
- x = z_infinity_f.f;
- }
-
- else
- x = z_infinity_f.f;
- }
-
- else if (t < SMALLX)
- {
- errno = ERANGE;
- x = 0.0;
- }
-
- else
- {
- if ( k && fabsf (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_FLOAT_WORD (x, px);
- }
- }
- }
-
- return x;
-}
-
-#ifdef _DOUBLE_IS_32BITS
-
-double pow (double x, double y)
-{
- return (double) powf ((float) x, (float) y);
-}
-
-#endif /* defined(_DOUBLE_IS_32BITS) */