Print this page
11210 libm should be cstyle(1ONBLD) clean
*** 20,45 ****
*/
/*
* Copyright 2011 Nexenta Systems, Inc. All rights reserved.
*/
/*
* Copyright 2006 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
#pragma weak __remquo = remquo
! /* INDENT OFF */
/*
* double remquo(double x, double y, int *quo) return remainder(x,y) and an
* integer pointer quo such that *quo = N mod {2**31}, where N is the
* exact integral part of x/y rounded to nearest even.
*
* remquo call internal fmodquo
*/
- /* INDENT ON */
#include "libm.h"
#include "libm_protos.h"
#include <math.h> /* fabs() */
#include <sys/isa_defs.h>
--- 20,45 ----
*/
/*
* Copyright 2011 Nexenta Systems, Inc. All rights reserved.
*/
+
/*
* Copyright 2006 Sun Microsystems, Inc. All rights reserved.
* Use is subject to license terms.
*/
#pragma weak __remquo = remquo
!
/*
* double remquo(double x, double y, int *quo) return remainder(x,y) and an
* integer pointer quo such that *quo = N mod {2**31}, where N is the
* exact integral part of x/y rounded to nearest even.
*
* remquo call internal fmodquo
*/
#include "libm.h"
#include "libm_protos.h"
#include <math.h> /* fabs() */
#include <sys/isa_defs.h>
*** 49,65 ****
#define LOWORD 1
#else
#define HIWORD 1
#define LOWORD 0
#endif
! #define __HI(x) ((int *) &x)[HIWORD]
! #define __LO(x) ((int *) &x)[LOWORD]
!
! static const double one = 1.0, Zero[] = {0.0, -0.0};
static double
! fmodquo(double x, double y, int *quo) {
int n, hx, hy, hz, ix, iy, sx, sq, i, m;
unsigned lx, ly, lz;
hx = __HI(x); /* high word of x */
lx = __LO(x); /* low word of x */
--- 49,65 ----
#define LOWORD 1
#else
#define HIWORD 1
#define LOWORD 0
#endif
! #define __HI(x) ((int *)&x)[HIWORD]
! #define __LO(x) ((int *)&x)[LOWORD]
+ static const double one = 1.0, Zero[] = { 0.0, -0.0 };
static double
! fmodquo(double x, double y, int *quo)
! {
int n, hx, hy, hz, ix, iy, sx, sq, i, m;
unsigned lx, ly, lz;
hx = __HI(x); /* high word of x */
lx = __LO(x); /* low word of x */
*** 70,89 ****
hx ^= sx; /* |x| */
hy &= 0x7fffffff; /* |y| */
/* purge off exception values */
*quo = 0;
if ((hy | ly) == 0 || hx >= 0x7ff00000 || /* y=0, or x !finite */
(hy | ((ly | -ly) >> 31)) > 0x7ff00000) /* or y is NaN */
return ((x * y) / (x * y));
if (hx <= hy) {
if (hx < hy || lx < ly)
return (x); /* |x|<|y| return x */
if (lx == ly) {
*quo = 1 + (sq >> 30);
/* |x|=|y| return x*0 */
! return (Zero[(unsigned) sx >> 31]);
}
}
/* determine ix = ilogb(x) */
if (hx < 0x00100000) { /* subnormal x */
--- 70,92 ----
hx ^= sx; /* |x| */
hy &= 0x7fffffff; /* |y| */
/* purge off exception values */
*quo = 0;
+
if ((hy | ly) == 0 || hx >= 0x7ff00000 || /* y=0, or x !finite */
(hy | ((ly | -ly) >> 31)) > 0x7ff00000) /* or y is NaN */
return ((x * y) / (x * y));
+
if (hx <= hy) {
if (hx < hy || lx < ly)
return (x); /* |x|<|y| return x */
+
if (lx == ly) {
*quo = 1 + (sq >> 30);
/* |x|=|y| return x*0 */
! return (Zero[(unsigned)sx >> 31]);
}
}
/* determine ix = ilogb(x) */
if (hx < 0x00100000) { /* subnormal x */
*** 92,133 ****
ix -= 1;
} else {
for (ix = -1022, i = (hx << 11); i > 0; i <<= 1)
ix -= 1;
}
! } else
ix = (hx >> 20) - 1023;
/* determine iy = ilogb(y) */
if (hy < 0x00100000) { /* subnormal y */
if (hy == 0) {
for (iy = -1043, i = ly; i > 0; i <<= 1)
iy -= 1;
} else {
for (iy = -1022, i = (hy << 11); i > 0; i <<= 1)
iy -= 1;
}
! } else
iy = (hy >> 20) - 1023;
/* set up {hx,lx}, {hy,ly} and align y to x */
! if (ix >= -1022)
hx = 0x00100000 | (0x000fffff & hx);
! else { /* subnormal x, shift x to normal */
n = -1022 - ix;
if (n <= 31) {
hx = (hx << n) | (lx >> (32 - n));
lx <<= n;
} else {
hx = lx << (n - 32);
lx = 0;
}
}
! if (iy >= -1022)
hy = 0x00100000 | (0x000fffff & hy);
! else { /* subnormal y, shift y to normal */
n = -1022 - iy;
if (n <= 31) {
hy = (hy << n) | (ly >> (32 - n));
ly <<= n;
} else {
hy = ly << (n - 32);
--- 95,141 ----
ix -= 1;
} else {
for (ix = -1022, i = (hx << 11); i > 0; i <<= 1)
ix -= 1;
}
! } else {
ix = (hx >> 20) - 1023;
+ }
/* determine iy = ilogb(y) */
if (hy < 0x00100000) { /* subnormal y */
if (hy == 0) {
for (iy = -1043, i = ly; i > 0; i <<= 1)
iy -= 1;
} else {
for (iy = -1022, i = (hy << 11); i > 0; i <<= 1)
iy -= 1;
}
! } else {
iy = (hy >> 20) - 1023;
+ }
/* set up {hx,lx}, {hy,ly} and align y to x */
! if (ix >= -1022) {
hx = 0x00100000 | (0x000fffff & hx);
! } else { /* subnormal x, shift x to normal */
n = -1022 - ix;
+
if (n <= 31) {
hx = (hx << n) | (lx >> (32 - n));
lx <<= n;
} else {
hx = lx << (n - 32);
lx = 0;
}
}
!
! if (iy >= -1022) {
hy = 0x00100000 | (0x000fffff & hy);
! } else { /* subnormal y, shift y to normal */
n = -1022 - iy;
+
if (n <= 31) {
hy = (hy << n) | (ly >> (32 - n));
ly <<= n;
} else {
hy = ly << (n - 32);
*** 136,217 ****
}
/* fix point fmod */
n = ix - iy;
m = 0;
while (n--) {
hz = hx - hy;
lz = lx - ly;
if (lx < ly)
hz -= 1;
if (hz < 0) {
hx = hx + hx + (lx >> 31);
lx = lx + lx;
} else {
m += 1;
if ((hz | lz) == 0) { /* return sign(x)*0 */
if (n < 31)
m <<= 1 + n;
else
m = 0;
m &= 0x7fffffff;
*quo = sq >= 0 ? m : -m;
! return (Zero[(unsigned) sx >> 31]);
}
hx = hz + hz + (lz >> 31);
lx = lz + lz;
}
m += m;
}
hz = hx - hy;
lz = lx - ly;
if (lx < ly)
hz -= 1;
if (hz >= 0) {
hx = hz;
lx = lz;
m += 1;
}
m &= 0x7fffffff;
*quo = sq >= 0 ? m : -m;
/* convert back to floating value and restore the sign */
! if ((hx | lx) == 0) { /* return sign(x)*0 */
! return (Zero[(unsigned) sx >> 31]);
! }
while (hx < 0x00100000) { /* normalize x */
hx = hx + hx + (lx >> 31);
lx = lx + lx;
iy -= 1;
}
if (iy >= -1022) { /* normalize output */
hx = (hx - 0x00100000) | ((iy + 1023) << 20);
__HI(x) = hx | sx;
__LO(x) = lx;
} else { /* subnormal output */
n = -1022 - iy;
if (n <= 20) {
! lx = (lx >> n) | ((unsigned) hx << (32 - n));
hx >>= n;
} else if (n <= 31) {
lx = (hx << (32 - n)) | (lx >> n);
hx = sx;
} else {
lx = hx >> (n - 32);
hx = sx;
}
__HI(x) = hx | sx;
__LO(x) = lx;
x *= one; /* create necessary signal */
}
return (x); /* exact output */
}
#define zero Zero[0]
double
! remquo(double x, double y, int *quo) {
int hx, hy, sx, sq;
double v;
unsigned ly;
hx = __HI(x); /* high word of x */
--- 144,241 ----
}
/* fix point fmod */
n = ix - iy;
m = 0;
+
while (n--) {
hz = hx - hy;
lz = lx - ly;
+
if (lx < ly)
hz -= 1;
+
if (hz < 0) {
hx = hx + hx + (lx >> 31);
lx = lx + lx;
} else {
m += 1;
+
if ((hz | lz) == 0) { /* return sign(x)*0 */
if (n < 31)
m <<= 1 + n;
else
m = 0;
+
m &= 0x7fffffff;
*quo = sq >= 0 ? m : -m;
! return (Zero[(unsigned)sx >> 31]);
}
+
hx = hz + hz + (lz >> 31);
lx = lz + lz;
}
+
m += m;
}
+
hz = hx - hy;
lz = lx - ly;
+
if (lx < ly)
hz -= 1;
+
if (hz >= 0) {
hx = hz;
lx = lz;
m += 1;
}
+
m &= 0x7fffffff;
*quo = sq >= 0 ? m : -m;
/* convert back to floating value and restore the sign */
! if ((hx | lx) == 0) /* return sign(x)*0 */
! return (Zero[(unsigned)sx >> 31]);
!
while (hx < 0x00100000) { /* normalize x */
hx = hx + hx + (lx >> 31);
lx = lx + lx;
iy -= 1;
}
+
if (iy >= -1022) { /* normalize output */
hx = (hx - 0x00100000) | ((iy + 1023) << 20);
__HI(x) = hx | sx;
__LO(x) = lx;
} else { /* subnormal output */
n = -1022 - iy;
+
if (n <= 20) {
! lx = (lx >> n) | ((unsigned)hx << (32 - n));
hx >>= n;
} else if (n <= 31) {
lx = (hx << (32 - n)) | (lx >> n);
hx = sx;
} else {
lx = hx >> (n - 32);
hx = sx;
}
+
__HI(x) = hx | sx;
__LO(x) = lx;
x *= one; /* create necessary signal */
}
+
return (x); /* exact output */
}
#define zero Zero[0]
double
! remquo(double x, double y, int *quo)
! {
int hx, hy, sx, sq;
double v;
unsigned ly;
hx = __HI(x); /* high word of x */
*** 222,266 ****
--- 246,300 ----
hx ^= sx; /* |x| */
hy &= 0x7fffffff; /* |y| */
/* purge off exception values */
*quo = 0;
+
if ((hy | ly) == 0 || hx >= 0x7ff00000 || /* y=0, or x !finite */
(hy | ((ly | -ly) >> 31)) > 0x7ff00000) /* or y is NaN */
return ((x * y) / (x * y));
y = fabs(y);
x = fabs(x);
+
if (hy <= 0x7fdfffff) {
x = fmodquo(x, y + y, quo);
*quo = ((*quo) & 0x3fffffff) << 1;
}
+
if (hy < 0x00200000) {
if (x + x > y) {
*quo += 1;
+
if (x == y)
x = zero;
else
x -= y;
+
if (x + x >= y) {
x -= y;
*quo += 1;
}
}
} else {
v = 0.5 * y;
+
if (x > v) {
*quo += 1;
+
if (x == y)
x = zero;
else
x -= y;
+
if (x >= v) {
x -= y;
*quo += 1;
}
}
}
+
if (sq != 0)
*quo = -(*quo);
+
return (sx == 0 ? x : -x);
}