projects / musl / blob
? search:
summary | shortlog | log | commit | commitdiff | tree
history | raw | HEAD
acos.s fix: use the formula acos(x) = atan2(sqrt(1-x),sqrt(1+x))
[musl] / src / math / i386 / acos.s
1 .global acosf
2 .type acosf,@function
3 acosf:
4         flds 4(%esp)
5         jmp 1f
6
7 .global acosl
8 .type acosl,@function
9 acosl:
10         fldt 4(%esp)
11         jmp 1f
12
13 .global acos
14 .type acos,@function
15 acos:
16         fldl 4(%esp)
17 1:      fld1
18         fld %st(1)
19         fld1
20         fsubp
21         fsqrt
22         fxch %st(2)
23         faddp
24         fsqrt
25         fpatan
26         fld1
27         fld1
28         faddp
29         fmulp
30         ret
musl