sysdeps/libm-i387/e_asinf.S

   1 /*
   2  * Written by J.T. Conklin <jtc@netbsd.org>.
   3  * Public domain.
   4  * Adapted for float type by Ulrich Drepper <drepper@cygnus.com>.
   5  */
   6
   7 #include <machine/asm.h>
   8
   9 RCSID("$NetBSD: $")
  10
  11 /* asin = atan (x / sqrt(1 - x^2)) */
  12 ENTRY(__ieee754_asinf)
  13         flds    4(%esp)                 /* x */
  14         fld     %st
  15         fmul    %st(0)                  /* x^2 */
  16         fld1
  17         fsubp                           /* 1 - x^2 */
  18         fsqrt                           /* sqrt (1 - x^2) */
  19         fpatan
  20         ret
  21 END (__ieee754_asinf)