about summary refs log tree commit diff
path: root/REORG.TODO/sysdeps/i386/fpu/e_acosl.c
diff options
context:
space:
mode:
Diffstat (limited to 'REORG.TODO/sysdeps/i386/fpu/e_acosl.c')
-rw-r--r--REORG.TODO/sysdeps/i386/fpu/e_acosl.c29
1 files changed, 29 insertions, 0 deletions
diff --git a/REORG.TODO/sysdeps/i386/fpu/e_acosl.c b/REORG.TODO/sysdeps/i386/fpu/e_acosl.c
new file mode 100644
index 0000000000..ab08931924
--- /dev/null
+++ b/REORG.TODO/sysdeps/i386/fpu/e_acosl.c
@@ -0,0 +1,29 @@
+/*
+ * Written by J.T. Conklin <jtc@netbsd.org>.
+ * Public domain.
+ *
+ * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
+ */
+
+#include <math_private.h>
+
+long double
+__ieee754_acosl (long double x)
+{
+  long double res;
+
+  /* acosl = atanl (sqrtl((1-x) (1+x)) / x) */
+  asm (	"fld	%%st\n"
+	"fld1\n"
+	"fsubp\n"
+	"fld1\n"
+	"fadd	%%st(2)\n"
+	"fmulp\n"			/* 1 - x^2 */
+	"fsqrt\n"			/* sqrtl (1 - x^2) */
+	"fabs\n"
+	"fxch	%%st(1)\n"
+	"fpatan"
+	: "=t" (res) : "0" (x) : "st(1)");
+  return res;
+}
+strong_alias (__ieee754_acosl, __acosl_finite)