projects
/
musl
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
add arch_prctl syscall (amd64/x32 only)
[musl]
/
src
/
math
/
acosl.c
diff --git
a/src/math/acosl.c
b/src/math/acosl.c
index
170520f
..
83857d4
100644
(file)
--- a/
src/math/acosl.c
+++ b/
src/math/acosl.c
@@
-23,10
+23,7
@@
long double acosl(long double x)
}
#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
#include "__invtrigl.h"
}
#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
#include "__invtrigl.h"
-
-static const long double
-one = 1.00000000000000000000e+00,
-pi = 3.14159265358979323846264338327950280e+00L;
+#define ACOS_CONST (BIAS - 65) /* 2**-65 */
long double acosl(long double x)
{
long double acosl(long double x)
{
@@
-42,7
+39,8
@@
long double acosl(long double x)
if (expsign > 0)
return 0.0; /* acos(1) = 0 */
else
if (expsign > 0)
return 0.0; /* acos(1) = 0 */
else
- return pi + 2.0 * pio2_lo; /* acos(-1)= pi */
+ // FIXME
+ return pi_hi + 2.0 * pio2_lo; /* acos(-1)= pi */
}
return (x - x) / (x - x); /* acos(|x|>1) is NaN */
}
}
return (x - x) / (x - x); /* acos(|x|>1) is NaN */
}
@@
-55,15
+53,15
@@
long double acosl(long double x)
r = p / q;
return pio2_hi - (x - (pio2_lo - x * r));
} else if (expsign < 0) { /* x < -0.5 */
r = p / q;
return pio2_hi - (x - (pio2_lo - x * r));
} else if (expsign < 0) { /* x < -0.5 */
- z = (
one
+ x) * 0.5;
+ z = (
1.0
+ x) * 0.5;
p = P(z);
q = Q(z);
s = sqrtl(z);
r = p / q;
w = r * s - pio2_lo;
p = P(z);
q = Q(z);
s = sqrtl(z);
r = p / q;
w = r * s - pio2_lo;
- return pi - 2.0 * (s + w);
+ return pi
_hi
- 2.0 * (s + w);
} else { /* x > 0.5 */
} else { /* x > 0.5 */
- z = (
one
- x) * 0.5;
+ z = (
1.0
- x) * 0.5;
s = sqrtl(z);
u.e = s;
u.bits.manl = 0;
s = sqrtl(z);
u.e = s;
u.bits.manl = 0;