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 9RCSID("$NetBSD: $") 10 11/* acos = atan (sqrt(1 - x^2) / x) */ 12ENTRY(__ieee754_acosf) 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 fxch %st(1) 20 fpatan 21 ret 22END (__ieee754_acosf) 23