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

*** 16,28 **** --- 16,30 ---- * fields enclosed by brackets "[]" replaced with your own identifying * information: Portions Copyright [yyyy] [name of copyright owner] * * CDDL HEADER END */ + /* * Copyright 2011 Nexenta Systems, Inc. All rights reserved. */ + /* * Copyright 2006 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */
*** 35,73 **** #include "libm.h" #if defined(__i386) && !defined(__amd64) && (!defined(__FLT_EVAL_METHOD__) || \ __FLT_EVAL_METHOD__ != 0) extern enum fp_precision_type __swapRP(enum fp_precision_type); #define DECLRP(x) enum fp_precision_type x; #define SWAPRP(new, x) x = __swapRP(new); #define RESTRP(x) (void) __swapRP(x); #else #define DECLRP(x) #define SWAPRP(new, x) #define RESTRP(x) #endif ! static const double ! two52 = 4503599627370496.0, ! zero = 0.0, ! one = 1.0; double ! rint(double x) { DECLRP(rp) double t, w; int ix, hx; ix = ((int *)&x)[HIWORD]; hx = ix & ~0x80000000; if (hx >= 0x43300000) return (x * one); ! t = (ix < 0)? -two52 : two52; SWAPRP(fp_double, rp) /* set precision mode to double */ w = x + t; /* x+sign(x)*2**52 rounded */ RESTRP(rp) /* restore precision mode */ if (w == t) ! return ((ix < 0)? -zero : zero); return (w - t); } --- 37,78 ---- #include "libm.h" #if defined(__i386) && !defined(__amd64) && (!defined(__FLT_EVAL_METHOD__) || \ __FLT_EVAL_METHOD__ != 0) + extern enum fp_precision_type __swapRP(enum fp_precision_type); + #define DECLRP(x) enum fp_precision_type x; #define SWAPRP(new, x) x = __swapRP(new); #define RESTRP(x) (void) __swapRP(x); #else #define DECLRP(x) #define SWAPRP(new, x) #define RESTRP(x) #endif ! static const double two52 = 4503599627370496.0, zero = 0.0, one = 1.0; double ! rint(double x) ! { DECLRP(rp) double t, w; int ix, hx; ix = ((int *)&x)[HIWORD]; hx = ix & ~0x80000000; if (hx >= 0x43300000) return (x * one); ! ! t = (ix < 0) ? -two52 : two52; SWAPRP(fp_double, rp) /* set precision mode to double */ w = x + t; /* x+sign(x)*2**52 rounded */ RESTRP(rp) /* restore precision mode */ + if (w == t) ! return ((ix < 0) ? -zero : zero); ! return (w - t); }