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

*** 16,49 **** * 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 2005 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #pragma weak __fmodf = fmodf #include "libm.h" ! /* INDENT OFF */ ! static const int ! is = (int)0x80000000, im = 0x007fffff, ii = 0x7f800000, iu = 0x00800000; - /* INDENT ON */ static const float zero = 0.0; float ! fmodf(float x, float y) { float w; int hx, ix, iy, iz, k, ny, nd; hx = *(int *)&x; ix = hx & 0x7fffffff; --- 16,49 ---- * 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 2005 Sun Microsystems, Inc. All rights reserved. * Use is subject to license terms. */ #pragma weak __fmodf = fmodf #include "libm.h" ! static const int is = (int)0x80000000, im = 0x007fffff, ii = 0x7f800000, iu = 0x00800000; static const float zero = 0.0; float ! fmodf(float x, float y) ! { float w; int hx, ix, iy, iz, k, ny, nd; hx = *(int *)&x; ix = hx & 0x7fffffff;
*** 57,86 **** if (ix < iy) w = x; /* return x if |x|<|y| */ else w = zero * x; /* return sign(x)*0.0 */ } else { ! /* INDENT OFF */ /* * scale x,y to "normal" with * ny = exponent of y * nd = exponent of x minus exponent of y */ - /* INDENT ON */ ny = iy >> 23; k = ix >> 23; /* special case for subnormal y or x */ if (ny == 0) { ny = 1; while (iy < iu) { ny -= 1; iy += iy; } nd = k - ny; if (k == 0) { nd += 1; while (ix < iu) { nd -= 1; ix += ix; } } else { --- 57,89 ---- if (ix < iy) w = x; /* return x if |x|<|y| */ else w = zero * x; /* return sign(x)*0.0 */ } else { ! /* * scale x,y to "normal" with * ny = exponent of y * nd = exponent of x minus exponent of y */ ny = iy >> 23; k = ix >> 23; /* special case for subnormal y or x */ if (ny == 0) { ny = 1; + while (iy < iu) { ny -= 1; iy += iy; } + nd = k - ny; + if (k == 0) { nd += 1; + while (ix < iu) { nd -= 1; ix += ix; } } else {
*** 90,101 **** nd = k - ny; ix = iu | (ix & im); iy = iu | (iy & im); } ! /* fix point fmod for normalized ix and iy */ ! /* INDENT OFF */ /* * while (nd--) { * iz = ix - iy; * if (iz < 0) * ix = ix + ix; --- 93,107 ---- nd = k - ny; ix = iu | (ix & im); iy = iu | (iy & im); } ! /* ! * fix point fmod for normalized ix and iy ! */ ! ! /* BEGIN CSTYLED */ /* * while (nd--) { * iz = ix - iy; * if (iz < 0) * ix = ix + ix;
*** 105,176 **** * } * else * ix = iz + iz; * } */ ! /* INDENT ON */ ! /* unroll the above loop 4 times to gain performance */ k = nd >> 2; nd -= k << 2; while (k--) { iz = ix - iy; if (iz >= 0) ix = iz + iz; else ix += ix; iz = ix - iy; if (iz >= 0) ix = iz + iz; else ix += ix; iz = ix - iy; if (iz >= 0) ix = iz + iz; else ix += ix; iz = ix - iy; if (iz >= 0) ix = iz + iz; else ix += ix; if (iz == 0) { *(int *)&w = is & hx; return (w); } } while (nd--) { iz = ix - iy; if (iz >= 0) ix = iz + iz; else ix += ix; } /* end of unrolling */ iz = ix - iy; if (iz >= 0) ix = iz; /* convert back to floating value and restore the sign */ if (ix == 0) { *(int *)&w = is & hx; return (w); } while (ix < iu) { ix += ix; ny -= 1; } while (ix > (iu + iu)) { ny += 1; ix >>= 1; } if (ny > 0) { *(int *)&w = (is & hx) | (ix & im) | (ny << 23); } else { /* subnormal output */ k = -ny + 1; ix >>= k; *(int *)&w = (is & hx) | ix; } } return (w); } --- 111,202 ---- * } * else * ix = iz + iz; * } */ ! /* END CSTYLED */ ! ! /* ! * unroll the above loop 4 times to gain performance ! */ k = nd >> 2; nd -= k << 2; + while (k--) { iz = ix - iy; + if (iz >= 0) ix = iz + iz; else ix += ix; + iz = ix - iy; + if (iz >= 0) ix = iz + iz; else ix += ix; + iz = ix - iy; + if (iz >= 0) ix = iz + iz; else ix += ix; + iz = ix - iy; + if (iz >= 0) ix = iz + iz; else ix += ix; + if (iz == 0) { *(int *)&w = is & hx; return (w); } } + while (nd--) { iz = ix - iy; + if (iz >= 0) ix = iz + iz; else ix += ix; } + /* end of unrolling */ iz = ix - iy; + if (iz >= 0) ix = iz; /* convert back to floating value and restore the sign */ if (ix == 0) { *(int *)&w = is & hx; return (w); } + while (ix < iu) { ix += ix; ny -= 1; } + while (ix > (iu + iu)) { ny += 1; ix >>= 1; } + if (ny > 0) { *(int *)&w = (is & hx) | (ix & im) | (ny << 23); } else { /* subnormal output */ k = -ny + 1; ix >>= k; *(int *)&w = (is & hx) | ix; } } + return (w); }