summaryrefslogtreecommitdiff
path: root/sysdeps/libm-i387/s_log1p.S
diff options
context:
space:
mode:
Diffstat (limited to 'sysdeps/libm-i387/s_log1p.S')
-rw-r--r--sysdeps/libm-i387/s_log1p.S23
1 files changed, 23 insertions, 0 deletions
diff --git a/sysdeps/libm-i387/s_log1p.S b/sysdeps/libm-i387/s_log1p.S
new file mode 100644
index 0000000000..c4aee4aa9c
--- /dev/null
+++ b/sysdeps/libm-i387/s_log1p.S
@@ -0,0 +1,23 @@
+/*
+ * Written by J.T. Conklin <jtc@netbsd.org>.
+ * Public domain.
+ */
+
+#include <machine/asm.h>
+
+RCSID("$NetBSD: s_log1p.S,v 1.7 1995/05/09 00:10:58 jtc Exp $")
+
+/*
+ * Since the fyl2xp1 instruction has such a limited range:
+ * -(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1
+ * it's not worth trying to use it.
+ */
+
+ENTRY(__log1p)
+ fldln2
+ fldl 4(%esp)
+ fld1
+ faddp
+ fyl2x
+ ret
+weak_alias (__log1p, log1p)