Home | History | Annotate | Line # | Download | only in i387
s_log1p.S revision 1.6
      1  1.1  jtc /*
      2  1.6  jtc  * Written by J.T. Conklin <jtc (at) netbsd.org>.
      3  1.6  jtc  * Public domain.
      4  1.1  jtc  */
      5  1.1  jtc 
      6  1.1  jtc #include <machine/asm.h>
      7  1.1  jtc 
      8  1.6  jtc RCSID("$Id: s_log1p.S,v 1.6 1995/04/28 22:44:10 jtc Exp $")
      9  1.2  jtc 
     10  1.4  jtc /*
     11  1.5  jtc  * Since the fyl2xp1 instruction has such a limited range:
     12  1.4  jtc  *	-(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1
     13  1.4  jtc  * it's not worth trying to use it.
     14  1.4  jtc  */
     15  1.5  jtc 
     16  1.1  jtc ENTRY(log1p)
     17  1.1  jtc 	fldln2
     18  1.3  jtc 	fldl 4(%esp)
     19  1.4  jtc 	fld1
     20  1.4  jtc 	faddp
     21  1.4  jtc 	fyl2x
     22  1.1  jtc 	ret
     23