sysdeps/libm-i387/s_log1pl.S

   1 /*
   2  * Written by J.T. Conklin <jtc@netbsd.org>.
   3  * Public domain.
   4  *
   5  * Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
   6  */
   7
   8 #include <machine/asm.h>
   9
  10 RCSID("$NetBSD: $")
  11
  12 /*
  13  * Since the fyl2xp1 instruction has such a limited range:
  14  *      -(1 - (sqrt(2) / 2)) <= x <= sqrt(2) - 1
  15  * it's not worth trying to use it.
  16  */
  17
  18 ENTRY(__log1pl)
  19         fldln2
  20         fldt 4(%esp)
  21         fld1
  22         faddp
  23         fyl2x
  24         ret
  25 PSEUDO_END (__log1pl)
  26 weak_alias (__log1pl, log1pl)