Print this page
5262 libm needs to be carefully unifdef'd
5268 libm doesn't need to hide symbols which are already local
@@ -25,19 +25,17 @@
/*
* Copyright 2006 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
-#if defined(ELFOBJ)
#pragma weak scalbln = __scalbln
-#endif
#include "libm.h"
#include <float.h> /* DBL_MAX, DBL_MIN */
static const double twom54 = 5.5511151231257827021181583404541015625e-17;
-#if defined(USE_FPSCALE) || defined(__x86)
+#if defined(__x86)
static const double two52 = 4503599627370496.0;
#else
/*
* Normalize non-zero subnormal x and return biased exponent of x in [-51,0]
*/
@@ -67,11 +65,11 @@
}
px[HIWORD] = (px[HIWORD] & 0x80000000) | v;
px[LOWORD] = w;
return (1 - s);
}
-#endif /* defined(USE_FPSCALE) */
+#endif /* defined(__x86) */
double
scalbln(double x, long n) {
int *px = (int *) &x, ix, k;
@@ -85,11 +83,11 @@
return (x + x);
#endif
if ((px[LOWORD] | ix) == 0 || n == 0)
return (x);
if (k == 0) {
-#if defined(USE_FPSCALE) || defined(__x86)
+#if defined(__x86)
x *= two52;
k = ((px[HIWORD] & ~0x80000000) >> 20) - 52;
#else
k = ilogb_biased((unsigned *) px);
#endif