Home | History | Annotate | Line # | Download | only in i387
s_log1p.S revision 1.9 
      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.9    fvdl #include "abi.h"
      9  1.9    fvdl 
     10  1.9    fvdl RCSID("$NetBSD: s_log1p.S,v 1.9 2001/06/19 00:26:30 fvdl Exp $")
     11  1.2     jtc 
     12  1.4     jtc /*
     13  1.5     jtc  * Since the fyl2xp1 instruction has such a limited range:
     14  1.4     jtc  *	-(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1
     15  1.8  simonb  * it's not worth trying to use it.
     16  1.4     jtc  */
     17  1.5     jtc 
     18  1.1     jtc ENTRY(log1p)
     19  1.9    fvdl 	XMM_ONE_ARG_DOUBLE_PROLOGUE
     20  1.1     jtc 	fldln2
     21  1.9    fvdl 	fldl	ARG_DOUBLE_ONE
     22  1.4     jtc 	fld1
     23  1.4     jtc 	faddp
     24  1.4     jtc 	fyl2x
     25  1.9    fvdl 	XMM_DOUBLE_EPILOGUE
     26  1.1     jtc 	ret
     27