NetBSD/lib/libm/arch/i387/s_log1p.S

/*
 * 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
More i387 math support, too bad we can't enable it 1994-02-12 04:31:52 +03:00			`/*`
Winning Strategies has placed this code into the Public Domain 1995-04-29 02:43:31 +04:00			`* Written by J.T. Conklin <jtc@netbsd.org>.`
			`* Public domain.`
More i387 math support, too bad we can't enable it 1994-02-12 04:31:52 +03:00			`*/`

			`#include <machine/asm.h>`

Id -> NetBSD 1995-05-09 03:55:16 +04:00			`RCSID("$NetBSD: s_log1p.S,v 1.7 1995/05/09 00:10:58 jtc Exp $")`
Added RCS ID's. ID's use new RCSID macro from <machine/asm.h>. 1994-03-12 04:30:22 +03:00
The fyl2xp1 instruction has a limited range: -(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1 so we can't use it. Also, I'm not sure fyl2xp1's extra precision will matter once the result is converted from extended real (80 bits) back to double real (64 bits). 1994-08-19 07:55:57 +04:00			`/*`
Remove lie (in a comment) about conversion between 80 bit ot 64 bit reals, since we set up the fpu to do all calculations in 64 bit mode at this time. 1994-08-26 04:20:34 +04:00			`* Since the fyl2xp1 instruction has such a limited range:`
The fyl2xp1 instruction has a limited range: -(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1 so we can't use it. Also, I'm not sure fyl2xp1's extra precision will matter once the result is converted from extended real (80 bits) back to double real (64 bits). 1994-08-19 07:55:57 +04:00			`* -(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1`
			`* it's not worth trying to use it.`
			`*/`
Remove lie (in a comment) about conversion between 80 bit ot 64 bit reals, since we set up the fpu to do all calculations in 64 bit mode at this time. 1994-08-26 04:20:34 +04:00
More i387 math support, too bad we can't enable it 1994-02-12 04:31:52 +03:00			`ENTRY(log1p)`
			`fldln2`
Major Bogon: change %ebp to %esp, as I don't set up a stack frame. 1994-08-19 00:42:35 +04:00			`fldl 4(%esp)`
The fyl2xp1 instruction has a limited range: -(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1 so we can't use it. Also, I'm not sure fyl2xp1's extra precision will matter once the result is converted from extended real (80 bits) back to double real (64 bits). 1994-08-19 07:55:57 +04:00			`fld1`
			`faddp`
			`fyl2x`
More i387 math support, too bad we can't enable it 1994-02-12 04:31:52 +03:00			`ret`