Print this page
11210 libm should be cstyle(1ONBLD) clean

*** 20,37 **** */ /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #pragma weak __asin = asin ! /* INDENT OFF */ /* * asin(x) * Method : * Since asin(x) = x + x^3/6 + x^5*3/40 + x^7*15/336 + ... * we approximate asin(x) on [0,0.5] by --- 20,38 ---- */ /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ + /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #pragma weak __asin = asin ! /* * asin(x) * Method : * Since asin(x) = x + x^3/6 + x^5*3/40 + x^7*15/336 + ... * we approximate asin(x) on [0,0.5] by
*** 58,92 **** * Special cases: * if x is NaN, return x itself; * if |x|>1, return NaN with invalid signal. * */ - /* INDENT ON */ #include "libm_protos.h" /* _SVID_libm_error */ #include "libm_macros.h" #include <math.h> - /* INDENT OFF */ static const double xxx[] = { ! /* one */ 1.00000000000000000000e+00, /* 3FF00000, 00000000 */ /* huge */ 1.000e+300, /* pio2_hi */ 1.57079632679489655800e+00, /* 3FF921FB, 54442D18 */ /* pio2_lo */ 6.12323399573676603587e-17, /* 3C91A626, 33145C07 */ /* pio4_hi */ 7.85398163397448278999e-01, /* 3FE921FB, 54442D18 */ ! /* coefficient for R(x^2) */ ! /* pS0 */ 1.66666666666666657415e-01, /* 3FC55555, 55555555 */ /* pS1 */ -3.25565818622400915405e-01, /* BFD4D612, 03EB6F7D */ /* pS2 */ 2.01212532134862925881e-01, /* 3FC9C155, 0E884455 */ /* pS3 */ -4.00555345006794114027e-02, /* BFA48228, B5688F3B */ /* pS4 */ 7.91534994289814532176e-04, /* 3F49EFE0, 7501B288 */ /* pS5 */ 3.47933107596021167570e-05, /* 3F023DE1, 0DFDF709 */ /* qS1 */ -2.40339491173441421878e+00, /* C0033A27, 1C8A2D4B */ /* qS2 */ 2.02094576023350569471e+00, /* 40002AE5, 9C598AC8 */ /* qS3 */ -6.88283971605453293030e-01, /* BFE6066C, 1B8D0159 */ /* qS4 */ 7.70381505559019352791e-02 /* 3FB3B8C5, B12E9282 */ }; #define one xxx[0] #define huge xxx[1] #define pio2_hi xxx[2] #define pio2_lo xxx[3] #define pio4_hi xxx[4] --- 59,97 ---- * Special cases: * if x is NaN, return x itself; * if |x|>1, return NaN with invalid signal. * */ #include "libm_protos.h" /* _SVID_libm_error */ #include "libm_macros.h" #include <math.h> static const double xxx[] = { ! /* one */ ! 1.00000000000000000000e+00, /* 3FF00000, 00000000 */ /* huge */ 1.000e+300, /* pio2_hi */ 1.57079632679489655800e+00, /* 3FF921FB, 54442D18 */ /* pio2_lo */ 6.12323399573676603587e-17, /* 3C91A626, 33145C07 */ /* pio4_hi */ 7.85398163397448278999e-01, /* 3FE921FB, 54442D18 */ ! ! /* ! * coefficient for R(x^2) ! * pS0 ! */ ! 1.66666666666666657415e-01, /* 3FC55555, 55555555 */ /* pS1 */ -3.25565818622400915405e-01, /* BFD4D612, 03EB6F7D */ /* pS2 */ 2.01212532134862925881e-01, /* 3FC9C155, 0E884455 */ /* pS3 */ -4.00555345006794114027e-02, /* BFA48228, B5688F3B */ /* pS4 */ 7.91534994289814532176e-04, /* 3F49EFE0, 7501B288 */ /* pS5 */ 3.47933107596021167570e-05, /* 3F023DE1, 0DFDF709 */ /* qS1 */ -2.40339491173441421878e+00, /* C0033A27, 1C8A2D4B */ /* qS2 */ 2.02094576023350569471e+00, /* 40002AE5, 9C598AC8 */ /* qS3 */ -6.88283971605453293030e-01, /* BFE6066C, 1B8D0159 */ /* qS4 */ 7.70381505559019352791e-02 /* 3FB3B8C5, B12E9282 */ }; + #define one xxx[0] #define huge xxx[1] #define pio2_hi xxx[2] #define pio2_lo xxx[3] #define pio4_hi xxx[4]
*** 98,157 **** #define pS5 xxx[10] #define qS1 xxx[11] #define qS2 xxx[12] #define qS3 xxx[13] #define qS4 xxx[14] - /* INDENT ON */ double ! asin(double x) { double t, w, p, q, c, r, s; int hx, ix, i; ! hx = ((int *) &x)[HIWORD]; ix = hx & 0x7fffffff; if (ix >= 0x3ff00000) { /* |x| >= 1 */ ! if (((ix - 0x3ff00000) | ((int *) &x)[LOWORD]) == 0) /* asin(1)=+-pi/2 with inexact */ return (x * pio2_hi + x * pio2_lo); else if (isnan(x)) #if defined(FPADD_TRAPS_INCOMPLETE_ON_NAN) return (ix >= 0x7ff80000 ? x : (x - x) / (x - x)); /* assumes sparc-like QNaN */ #else ! return (x - x) / (x - x); /* asin(|x|>1) is NaN */ #endif else return (_SVID_libm_err(x, x, 2)); } else if (ix < 0x3fe00000) { /* |x| < 0.5 */ if (ix < 0x3e400000) { /* if |x| < 2**-27 */ ! if ((i = (int) x) == 0) /* return x with inexact if x != 0 */ return (x); } t = x * x; ! p = t * (pS0 + t * (pS1 + t * (pS2 + t * (pS3 + ! t * (pS4 + t * pS5))))); q = one + t * (qS1 + t * (qS2 + t * (qS3 + t * qS4))); w = p / q; return (x + x * w); } /* 1 > |x| >= 0.5 */ w = one - fabs(x); t = w * 0.5; p = t * (pS0 + t * (pS1 + t * (pS2 + t * (pS3 + t * (pS4 + t * pS5))))); q = one + t * (qS1 + t * (qS2 + t * (qS3 + t * qS4))); s = sqrt(t); if (ix >= 0x3FEF3333) { /* if |x| > 0.975 */ w = p / q; t = pio2_hi - (2.0 * (s + s * w) - pio2_lo); } else { w = s; ! ((int *) &w)[LOWORD] = 0; c = (t - w * w) / (s + w); r = p / q; p = 2.0 * s * r - (pio2_lo - 2.0 * c); q = pio4_hi - 2.0 * w; t = pio4_hi - (p - q); } return (hx > 0 ? t : -t); } --- 103,168 ---- #define pS5 xxx[10] #define qS1 xxx[11] #define qS2 xxx[12] #define qS3 xxx[13] #define qS4 xxx[14] double ! asin(double x) ! { double t, w, p, q, c, r, s; int hx, ix, i; ! hx = ((int *)&x)[HIWORD]; ix = hx & 0x7fffffff; + if (ix >= 0x3ff00000) { /* |x| >= 1 */ ! if (((ix - 0x3ff00000) | ((int *)&x)[LOWORD]) == 0) /* asin(1)=+-pi/2 with inexact */ return (x * pio2_hi + x * pio2_lo); else if (isnan(x)) #if defined(FPADD_TRAPS_INCOMPLETE_ON_NAN) return (ix >= 0x7ff80000 ? x : (x - x) / (x - x)); + /* assumes sparc-like QNaN */ #else ! return ((x - x) / (x - x)); /* asin(|x|>1) is NaN */ #endif else return (_SVID_libm_err(x, x, 2)); } else if (ix < 0x3fe00000) { /* |x| < 0.5 */ if (ix < 0x3e400000) { /* if |x| < 2**-27 */ ! if ((i = (int)x) == 0) /* return x with inexact if x != 0 */ return (x); } + t = x * x; ! p = t * (pS0 + t * (pS1 + t * (pS2 + t * (pS3 + t * (pS4 + t * ! pS5))))); q = one + t * (qS1 + t * (qS2 + t * (qS3 + t * qS4))); w = p / q; return (x + x * w); } + /* 1 > |x| >= 0.5 */ w = one - fabs(x); t = w * 0.5; p = t * (pS0 + t * (pS1 + t * (pS2 + t * (pS3 + t * (pS4 + t * pS5))))); q = one + t * (qS1 + t * (qS2 + t * (qS3 + t * qS4))); s = sqrt(t); + if (ix >= 0x3FEF3333) { /* if |x| > 0.975 */ w = p / q; t = pio2_hi - (2.0 * (s + s * w) - pio2_lo); } else { w = s; ! ((int *)&w)[LOWORD] = 0; c = (t - w * w) / (s + w); r = p / q; p = 2.0 * s * r - (pio2_lo - 2.0 * c); q = pio4_hi - 2.0 * w; t = pio4_hi - (p - q); } + return (hx > 0 ? t : -t); }