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.
*/
*** 42,65 ****
#include <sys/isa_defs.h> /* _ILP32 */
#include "libm.h"
#if defined(_ILP32)
long
! lrint(double x) {
/*
* Note: The following code works on x86 (in the default rounding
* precision mode), but one should just use the fistpl instruction
* instead.
*/
union {
unsigned i[2];
double d;
} xx, yy;
unsigned hx;
xx.d = x;
hx = xx.i[HIWORD] & ~0x80000000;
if (hx < 0x43300000) { /* |x| < 2^52 */
/* add and subtract a power of two to round x to an integer */
#if defined(__sparc)
yy.i[HIWORD] = (xx.i[HIWORD] & 0x80000000) | 0x43300000;
#elif defined(__x86)
--- 43,69 ----
#include <sys/isa_defs.h> /* _ILP32 */
#include "libm.h"
#if defined(_ILP32)
long
! lrint(double x)
! {
/*
* Note: The following code works on x86 (in the default rounding
* precision mode), but one should just use the fistpl instruction
* instead.
*/
union {
unsigned i[2];
double d;
} xx, yy;
+
unsigned hx;
xx.d = x;
hx = xx.i[HIWORD] & ~0x80000000;
+
if (hx < 0x43300000) { /* |x| < 2^52 */
/* add and subtract a power of two to round x to an integer */
#if defined(__sparc)
yy.i[HIWORD] = (xx.i[HIWORD] & 0x80000000) | 0x43300000;
#elif defined(__x86)
*** 70,79 ****
yy.i[LOWORD] = 0;
x = (x + yy.d) - yy.d;
}
/* now x is nan, inf, or integral */
! return ((long) x);
}
#else
#error Unsupported architecture
#endif /* defined(_ILP32) */
--- 74,83 ----
yy.i[LOWORD] = 0;
x = (x + yy.d) - yy.d;
}
/* now x is nan, inf, or integral */
! return ((long)x);
}
#else
#error Unsupported architecture
#endif /* defined(_ILP32) */