| /* origin: FreeBSD /usr/src/lib/msun/src/e_acosf.c */ |
| /* |
| * Conversion to float by Ian Lance Taylor, Cygnus Support, ian@cygnus.com. |
| */ |
| /* |
| * ==================================================== |
| * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved. |
| * |
| * Developed at SunPro, a Sun Microsystems, Inc. business. |
| * Permission to use, copy, modify, and distribute this |
| * software is freely granted, provided that this notice |
| * is preserved. |
| * ==================================================== |
| */ |
| |
| #include "libm.h" |
| |
| static const float |
| pio2_hi = 1.5707962513e+00, /* 0x3fc90fda */ |
| pio2_lo = 7.5497894159e-08, /* 0x33a22168 */ |
| pS0 = 1.6666586697e-01, |
| pS1 = -4.2743422091e-02, |
| pS2 = -8.6563630030e-03, |
| qS1 = -7.0662963390e-01; |
| |
| static float R(float z) |
| { |
| float_t p, q; |
| p = z*(pS0+z*(pS1+z*pS2)); |
| q = 1.0f+z*qS1; |
| return p/q; |
| } |
| |
| float acosf(float x) |
| { |
| float z,w,s,c,df; |
| uint32_t hx,ix; |
| |
| GET_FLOAT_WORD(hx, x); |
| ix = hx & 0x7fffffff; |
| /* |x| >= 1 or nan */ |
| if (ix >= 0x3f800000) { |
| if (ix == 0x3f800000) { |
| if (hx >> 31) |
| return 2*pio2_hi + 0x1p-120f; |
| return 0; |
| } |
| return 0/(x-x); |
| } |
| /* |x| < 0.5 */ |
| if (ix < 0x3f000000) { |
| if (ix <= 0x32800000) /* |x| < 2**-26 */ |
| return pio2_hi + 0x1p-120f; |
| return pio2_hi - (x - (pio2_lo-x*R(x*x))); |
| } |
| /* x < -0.5 */ |
| if (hx >> 31) { |
| z = (1+x)*0.5f; |
| s = sqrtf(z); |
| w = R(z)*s-pio2_lo; |
| return 2*(pio2_hi - (s+w)); |
| } |
| /* x > 0.5 */ |
| z = (1-x)*0.5f; |
| s = sqrtf(z); |
| GET_FLOAT_WORD(hx,s); |
| SET_FLOAT_WORD(df,hx&0xfffff000); |
| c = (z-df*df)/(s+df); |
| w = R(z)*s+c; |
| return 2*(df+w); |
| } |
