Blame sysdeps/i386/fpu/e_acosl.c
|
Packit |
6c4009 |
/*
|
|
Packit |
6c4009 |
* Written by J.T. Conklin <jtc@netbsd.org>.
|
|
Packit |
6c4009 |
* Public domain.
|
|
Packit |
6c4009 |
*
|
|
Packit |
6c4009 |
* Adapted for `long double' by Ulrich Drepper <drepper@cygnus.com>.
|
|
Packit |
6c4009 |
*/
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
#include <math_private.h>
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
long double
|
|
Packit |
6c4009 |
__ieee754_acosl (long double x)
|
|
Packit |
6c4009 |
{
|
|
Packit |
6c4009 |
long double res;
|
|
Packit |
6c4009 |
|
|
Packit |
6c4009 |
/* acosl = atanl (sqrtl((1-x) (1+x)) / x) */
|
|
Packit |
6c4009 |
asm ( "fld %%st\n"
|
|
Packit |
6c4009 |
"fld1\n"
|
|
Packit |
6c4009 |
"fsubp\n"
|
|
Packit |
6c4009 |
"fld1\n"
|
|
Packit |
6c4009 |
"fadd %%st(2)\n"
|
|
Packit |
6c4009 |
"fmulp\n" /* 1 - x^2 */
|
|
Packit |
6c4009 |
"fsqrt\n" /* sqrtl (1 - x^2) */
|
|
Packit |
6c4009 |
"fabs\n"
|
|
Packit |
6c4009 |
"fxch %%st(1)\n"
|
|
Packit |
6c4009 |
"fpatan"
|
|
Packit |
6c4009 |
: "=t" (res) : "0" (x) : "st(1)");
|
|
Packit |
6c4009 |
return res;
|
|
Packit |
6c4009 |
}
|
|
Packit |
6c4009 |
strong_alias (__ieee754_acosl, __acosl_finite)
|