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

*** 20,29 **** --- 20,30 ---- */ /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ + /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */
*** 39,83 **** #else /* * v: high part of a non-zero subnormal |x|; w: low part of |x| */ static int ! ilogb_subnormal(unsigned v, unsigned w) { int r = -1022 - 52; if (v) r += 32; else v = w; if (v & 0xffff0000) r += 16, v >>= 16; if (v & 0xff00) r += 8, v >>= 8; if (v & 0xf0) r += 4, v >>= 4; v <<= 1; return (r + ((0xffffaa50 >> v) & 0x3)); } #endif /* defined(__x86) */ double ! logb(double x) { ! int *px = (int *) &x, k = px[HIWORD] & ~0x80000000; if (k < 0x00100000) { ! if ((px[LOWORD] | k) == 0) return (_SVID_libm_err(x, x, 45)); ! else if ((__xpg6 & _C99SUSv3_logb) != 0) { #if defined(__x86) x *= two52; ! return ((double) (((px[HIWORD] & 0x7ff00000) >> 20) ! - 1075)); #else ! return ((double) ilogb_subnormal(k, px[LOWORD])); #endif ! } else return (-1022.0); ! } else if (k < 0x7ff00000) ! return ((double) ((k >> 20) - 1023)); ! else return (x * x); } --- 40,92 ---- #else /* * v: high part of a non-zero subnormal |x|; w: low part of |x| */ static int ! ilogb_subnormal(unsigned v, unsigned w) ! { int r = -1022 - 52; if (v) r += 32; else v = w; + if (v & 0xffff0000) r += 16, v >>= 16; + if (v & 0xff00) r += 8, v >>= 8; + if (v & 0xf0) r += 4, v >>= 4; + v <<= 1; return (r + ((0xffffaa50 >> v) & 0x3)); } #endif /* defined(__x86) */ double ! logb(double x) ! { ! int *px = (int *)&x, k = px[HIWORD] & ~0x80000000; if (k < 0x00100000) { ! if ((px[LOWORD] | k) == 0) { return (_SVID_libm_err(x, x, 45)); ! } else if ((__xpg6 & _C99SUSv3_logb) != 0) { #if defined(__x86) x *= two52; ! return ((double)(((px[HIWORD] & 0x7ff00000) >> 20) - ! 1075)); #else ! return ((double)ilogb_subnormal(k, px[LOWORD])); #endif ! } else { return (-1022.0); ! } ! } else if (k < 0x7ff00000) { ! return ((double)((k >> 20) - 1023)); ! } else { return (x * x); + } }