blob: 2daa5b03815d76154c2660f10919f1f61468c3ff (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
/*
* ====================================================
* 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.
* ====================================================
*/
#if !defined(_SOFT_FLOAT)
/*
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);
}
#endif
|
