# use acos(x) = atan2(fabs(sqrt((1-x)*(1+x))), x).global acosf.type acosf,@functionacosf:flds 4(%esp)jmp 1f.global acosl.type acosl,@functionacosl:fldt 4(%esp)jmp 1f.global acos.type acos,@functionacos:fldl 4(%esp)1: fld %st(0)fld1fsub %st(0),%st(1)fadd %st(2)fmulpfsqrtfabs # fix sign of zero (matters in downward rounding mode)fxch %st(1)fpatanret