summaryrefslogtreecommitdiff
path: root/src/math/i386/acos.s
diff options
context:
space:
mode:
authorVladimir Azarov <avm@intermediate-node.net>2024-10-01 15:47:05 +0200
committerVladimir Azarov <avm@intermediate-node.net>2024-10-01 15:47:05 +0200
commit4abab5ad6c8465a7528ccdd5f49367da05f78bbd (patch)
treeebf009bf1376a5a223a915bc27cbbd791a1316bc /src/math/i386/acos.s
Initial version
Diffstat (limited to 'src/math/i386/acos.s')
-rw-r--r--src/math/i386/acos.s18
1 files changed, 18 insertions, 0 deletions
diff --git a/src/math/i386/acos.s b/src/math/i386/acos.s
new file mode 100644
index 0000000..af423a2
--- /dev/null
+++ b/src/math/i386/acos.s
@@ -0,0 +1,18 @@
+# use acos(x) = atan2(fabs(sqrt((1-x)*(1+x))), x)
+
+.global acos
+.type acos,@function
+acos:
+ fldl 4(%esp)
+ fld %st(0)
+ fld1
+ fsub %st(0),%st(1)
+ fadd %st(2)
+ fmulp
+ fsqrt
+ fabs # fix sign of zero (matters in downward rounding mode)
+ fxch %st(1)
+ fpatan
+ fstpl 4(%esp)
+ fldl 4(%esp)
+ ret