fluencelabs/musl
1
0
Fork 0
mirror of https://github.com/fluencelabs/musl synced 2025-07-02 16:12:02 +00:00
Code Issues Projects Releases Wiki Activity

math: sin cos cleanup

Browse Source 
* use unsigned arithmetics
* use unsigned to store arg reduction quotient (so n&3 is understood)
* remove z=0.0 variables, use literal 0
* raise underflow and inexact exceptions properly when x is small
* fix spurious underflow in tanl
This commit is contained in:
Szabolcs Nagy
2013-05-18 14:40:22 +00:00
parent 1d5ba3bb5a
commit bfda379358
10 changed files with 128 additions and 112 deletions
Download Patch File Download Diff File Expand all files Collapse all files

33
src/math/cosf.c
View File

@ -26,34 +26,39 @@ c4pio2 = 4*M_PI_2; /* 0x401921FB, 0x54442D18 */
float cosf(float x)
{
double y;
int32_t n, hx, ix;
uint32_t ix;
unsigned n, sign;
GET_FLOAT_WORD(ix, x);
sign = ix >> 31;
ix &= 0x7fffffff;
GET_FLOAT_WORD(hx, x);
ix = hx & 0x7fffffff;
if (ix <= 0x3f490fda) { /* |x| ~<= pi/4 */
if (ix < 0x39800000) /* |x| < 2**-12 */
if ((int)x == 0) /* raise inexact if x != 0 */
return 1.0;
if (ix < 0x39800000) { /* |x| < 2**-12 */
/* raise inexact if x != 0 */
FORCE_EVAL(x + 0x1p120f);
return 1.0f;
}
return __cosdf(x);
}
if (ix <= 0x407b53d1) { /* |x| ~<= 5*pi/4 */
if (ix > 0x4016cbe3) /* |x| ~> 3*pi/4 */
return -__cosdf(hx > 0 ? x-c2pio2 : x+c2pio2);
return -__cosdf(sign ? x+c2pio2 : x-c2pio2);
else {
if (hx > 0)
return __sindf(c1pio2 - x);
else
if (sign)
return __sindf(x + c1pio2);
else
return __sindf(c1pio2 - x);
}
}
if (ix <= 0x40e231d5) { /* |x| ~<= 9*pi/4 */
if (ix > 0x40afeddf) /* |x| ~> 7*pi/4 */
return __cosdf(hx > 0 ? x-c4pio2 : x+c4pio2);
return __cosdf(sign ? x+c4pio2 : x-c4pio2);
else {
if (hx > 0)
return __sindf(x - c3pio2);
if (sign)
return __sindf(-x - c3pio2);
else
return __sindf(-c3pio2 - x);
return __sindf(x - c3pio2);
}
}