aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/libm-i387/s_log1pf.S
diff options
context:
space:
mode:
Diffstat (limited to 'sysdeps/libm-i387/s_log1pf.S')
-rw-r--r--sysdeps/libm-i387/s_log1pf.S54
1 files changed, 46 insertions, 8 deletions
diff --git a/sysdeps/libm-i387/s_log1pf.S b/sysdeps/libm-i387/s_log1pf.S
index 7a161f26f9..df9fdcbcfc 100644
--- a/sysdeps/libm-i387/s_log1pf.S
+++ b/sysdeps/libm-i387/s_log1pf.S
@@ -7,18 +7,56 @@
RCSID("$NetBSD: s_log1pf.S,v 1.4 1995/05/09 00:13:05 jtc Exp $")
+#ifdef __ELF__
+ .section .rodata
+#else
+ .text
+#endif
+ .align ALIGNARG(4)
+ /* The fyl2xp1 can only be used for values in
+ -1 + sqrt(2) / 2 <= x <= 1 - sqrt(2) / 2
+ 0.29 is a safe value.
+ */
+limit: .float 0.29
+one: .float 1.0
+
/*
- * 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.
+ * Use the fyl2xp1 function when the argument is in the range -0.29 to 0.29,
+ * otherwise fyl2x with the needed extra computation.
*/
-
+ .text
ENTRY(__log1pf)
fldln2
- flds 4(%esp)
- fld1
- faddp
+
+ flds 4(%esp)
+
+#ifdef PIC
+ call 1f
+1: popl %edx
+ addl $_GLOBAL_OFFSET_TABLE_+[.-1b], %edx
+#endif
+
+ fld %st
+ fabs
+#ifdef PIC
+ fcomps limit@GOTOFF(%edx)
+#else
+ fcomps limit
+#endif
+ fnstsw
+ sahf
+ jc 2f
+
+#ifdef PIC
+ fadds one@GOTOFF(%edx)
+#else
+ fadds one
+#endif
fyl2x
ret
-PSEUDO_END (__log1pf)
+
+2: fyl2xp1
+ ret
+
+END (__log1pf)
weak_alias (__log1pf, log1pf)