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/libc/machine/i386/f_powf.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/libc/machine/i386/f_powf.c')
-rw-r--r-- | newlib/libc/machine/i386/f_powf.c | 44 |
1 files changed, 0 insertions, 44 deletions
diff --git a/newlib/libc/machine/i386/f_powf.c b/newlib/libc/machine/i386/f_powf.c deleted file mode 100644 index 501b6b1..0000000 --- a/newlib/libc/machine/i386/f_powf.c +++ /dev/null @@ -1,44 +0,0 @@ -/* - * ==================================================== - * Copyright (C) 1998 by Cygnus Solutions. All rights reserved. - * - * Permission to use, copy, modify, and distribute this - * software is freely granted, provided that this notice - * is preserved. - * ==================================================== - */ - -/* -Fast version of pow using Intel float instructions. - - float _f_powf (float x, float y); - -Function calculates x to power of y. -The function optimizes the case where x is >0.0 and y is finite. -In such a case, there is no error checking or setting of errno. -All other cases defer to normal powf() function which will -set errno as normal. -*/ - -#include <math.h> -#include <ieeefp.h> -#include "f_math.h" - -float _f_powf (float x, float y) -{ - /* following sequence handles the majority of cases for pow() */ - if (x > 0.0 && check_finitef(y)) - { - float result; - /* calculate x ** y as 2 ** (y log2(x)). On Intel, can only - raise 2 to an integer or a small fraction, thus, we have - to perform two steps 2**integer portion * 2**fraction. */ - asm ("flds 8(%%ebp); fyl2x; fld %%st; frndint; fsub %%st,%%st(1);" \ - "fxch; fchs; f2xm1; fld1; faddp; fxch; fld1; fscale; fstp %%st(1);"\ - "fmulp" : "=t" (result) : "0" (y)); - return result; - } - else /* all other strange cases, defer to normal pow() */ - return powf (x,y); -} - |