musl - musl - an implementation of the standard library for Linux-based systems

acos.s fix: use the formula acos(x) = atan2(sqrt(1-x),sqrt(1+x))

the old formula atan2(1,sqrt((1+x)/(1-x))) was faster but could give nan result at x=1 when the rounding mode is FE_DOWNWARD (so 1-1 == -0 and 2/-0 == -inf), the new formula gives -0 at x=+-1 with downward rounding.

Diffstat (limited to 'include/complex.h')

0 files changed, 0 insertions, 0 deletions