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. */
*** 31,80 **** #include "libm.h" #include <float.h> /* LDBL_MAX, LDBL_MIN */ #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 ! scalblnl(long double x, long 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; } ! k += (int) n; if (n > 50000 || k > 0x7ffe) return (LDBL_MAX * copysignl(LDBL_MAX, x)); if (n < -50000 || 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); } --- 32,92 ---- #include "libm.h" #include <float.h> /* LDBL_MAX, LDBL_MIN */ #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 ! scalblnl(long double x, long 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; } ! ! k += (int)n; ! if (n > 50000 || k > 0x7ffe) return (LDBL_MAX * copysignl(LDBL_MAX, x)); + if (n < -50000 || 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); }