Home | History | Annotate | Line # | Download | only in i387
s_log1p.S revision 1.10 
      1   1.1     jtc /*
      2  1.10    salo  * 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.10    salo RCSID("$NetBSD: s_log1p.S,v 1.10 2003/07/26 19:25:02 salo 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