diff options
Diffstat (limited to 'sysdeps/libm-i387/s_log1p.S')
| -rw-r--r-- | sysdeps/libm-i387/s_log1p.S | 23 |
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) |