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. */
*** 32,83 **** #include "libm.h" #include <float.h> /* LDBL_MAX, LDBL_MIN */ #include <stdlib.h> /* abs */ #if defined(__sparc) ! #define XSET_EXP(k, x) ((int *) &x)[0] = (((int *) &x)[0] & ~0x7fff0000) | \ ! (k << 16) #define ISINFNANL(k, x) (k == 0x7fff) #define XTWOT_OFFSET 113 ! static const long double xtwot = 10384593717069655257060992658440192.0L, ! /* 2^113 */ ! twomtm1 = 4.814824860968089632639944856462318296E-35L; /* 2^-114 */ #elif defined(__x86) ! #define XSET_EXP(k, x) ((int *) &x)[2] = (((int *) &x)[2] & ~0x7fff) | k #if defined(HANDLE_UNSUPPORTED) ! #define ISINFNANL(k, x) (k == 0x7fff || k != 0 && \ ! (((int *) &x)[1] & 0x80000000) == 0) #else #define ISINFNANL(k, x) (k == 0x7fff) #endif #define XTWOT_OFFSET 64 static const long double xtwot = 18446744073709551616.0L, /* 2^64 */ twomtm1 = 2.7105054312137610850186E-20L; /* 2^-65 */ #endif long double ! scalbnl(long double x, int n) { int k = XBIASED_EXP(x); if (ISINFNANL(k, x)) return (x + x); if (ISZEROL(x) || n == 0) return (x); if (k == 0) { x *= xtwot; k = XBIASED_EXP(x) - XTWOT_OFFSET; } ! if ((unsigned) abs(n) >= 131072) /* cast to unsigned for -2^31 */ n >>= 1; /* avoid subsequent integer overflow */ k += n; if (k > 0x7ffe) return (LDBL_MAX * copysignl(LDBL_MAX, x)); if (k <= -XTWOT_OFFSET - 1) return (LDBL_MIN * copysignl(LDBL_MIN, x)); if (k > 0) { XSET_EXP(k, x); return (x); } k += XTWOT_OFFSET + 1; XSET_EXP(k, x); return (x * twomtm1); } --- 33,96 ---- #include "libm.h" #include <float.h> /* LDBL_MAX, LDBL_MIN */ #include <stdlib.h> /* abs */ #if defined(__sparc) ! #define XSET_EXP(k, x) ((int *)&x)[0] = (((int *)&x)[0] & \ ! ~0x7fff0000) | (k << 16) #define ISINFNANL(k, x) (k == 0x7fff) #define XTWOT_OFFSET 113 ! ! /* 2^113 */ ! static const long double xtwot = 10384593717069655257060992658440192.0L; ! /* 2^-114 */ ! static const long double twomtm1 = 4.814824860968089632639944856462318296E-35L; #elif defined(__x86) ! #define XSET_EXP(k, x) ((int *)&x)[2] = (((int *)&x)[2] & ~0x7fff) | k #if defined(HANDLE_UNSUPPORTED) ! #define ISINFNANL(k, x) (k == 0x7fff || k != 0 && (((int *)&x)[1] & \ ! 0x80000000) == 0) #else #define ISINFNANL(k, x) (k == 0x7fff) #endif #define XTWOT_OFFSET 64 + static const long double xtwot = 18446744073709551616.0L, /* 2^64 */ twomtm1 = 2.7105054312137610850186E-20L; /* 2^-65 */ #endif long double ! scalbnl(long double x, int n) ! { int k = XBIASED_EXP(x); if (ISINFNANL(k, x)) return (x + x); + if (ISZEROL(x) || n == 0) return (x); + if (k == 0) { x *= xtwot; k = XBIASED_EXP(x) - XTWOT_OFFSET; } ! ! if ((unsigned)abs(n) >= 131072) /* cast to unsigned for -2^31 */ n >>= 1; /* avoid subsequent integer overflow */ + k += n; + if (k > 0x7ffe) return (LDBL_MAX * copysignl(LDBL_MAX, x)); + if (k <= -XTWOT_OFFSET - 1) return (LDBL_MIN * copysignl(LDBL_MIN, x)); + if (k > 0) { XSET_EXP(k, x); return (x); } + k += XTWOT_OFFSET + 1; XSET_EXP(k, x); return (x * twomtm1); }