projects
/
oweals
/
musl.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
math: add fma TODO comments about the underflow issue
[oweals/musl.git]
/
src
/
math
/
acos.c
diff --git
a/src/math/acos.c
b/src/math/acos.c
index be95d25e24a14fda2d7e672dd2424585266a851c..ea9c87bf087809f9a6b1fa6a2ed2f065d8184f72 100644
(file)
--- a/
src/math/acos.c
+++ b/
src/math/acos.c
@@
-51,7
+51,7
@@
qS4 = 7.70381505559019352791e-02; /* 0x3FB3B8C5, 0xB12E9282 */
static double R(double z)
{
static double R(double z)
{
- double p, q;
+ double
_t
p, q;
p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5)))));
q = 1.0+z*(qS1+z*(qS2+z*(qS3+z*qS4)));
return p/q;
p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5)))));
q = 1.0+z*(qS1+z*(qS2+z*(qS3+z*qS4)));
return p/q;
@@
-72,7
+72,7
@@
double acos(double x)
if ((ix-0x3ff00000 | lx) == 0) {
/* acos(1)=0, acos(-1)=pi */
if (hx >> 31)
if ((ix-0x3ff00000 | lx) == 0) {
/* acos(1)=0, acos(-1)=pi */
if (hx >> 31)
- return 2*pio2_hi + 0x1p-1
000
;
+ return 2*pio2_hi + 0x1p-1
20f
;
return 0;
}
return 0/(x-x);
return 0;
}
return 0/(x-x);
@@
-80,7
+80,7
@@
double acos(double x)
/* |x| < 0.5 */
if (ix < 0x3fe00000) {
if (ix <= 0x3c600000) /* |x| < 2**-57 */
/* |x| < 0.5 */
if (ix < 0x3fe00000) {
if (ix <= 0x3c600000) /* |x| < 2**-57 */
- return pio2_hi + 0x1p-1
000
;
+ return pio2_hi + 0x1p-1
20f
;
return pio2_hi - (x - (pio2_lo-x*R(x*x)));
}
/* x < -0.5 */
return pio2_hi - (x - (pio2_lo-x*R(x*x)));
}
/* x < -0.5 */