sysdeps/libm-i387/s_log1p.S

   1 /*
   2  * Written by J.T. Conklin <jtc@netbsd.org>.
   3  * Public domain.
   4  */
   5
   6 #include <machine/asm.h>
   7
   8 RCSID("$NetBSD: s_log1p.S,v 1.7 1995/05/09 00:10:58 jtc Exp $")
   9
  10 /*
  11  * Since the fyl2xp1 instruction has such a limited range:
  12  *      -(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1
  13  * it's not worth trying to use it.
  14  */
  15
  16 ENTRY(__log1p)
  17         fldln2
  18         fldl 4(%esp)
  19         fld1
  20         faddp
  21         fyl2x
  22         ret
  23 weak_alias (__log1p, log1p)