/illumos-gate/usr/src/lib/libm/common/R/ |
H A D | hypotf.c | 38 int ix, iy; local 41 iy = (*(int *) &y) & 0x7fffffff; 44 *(int *) &w = x == y ? iy : ix; /* w = |x| = inf */ 45 else if (iy == 0x7f800000) 46 *(int *) &w = x == y ? ix : iy; /* w = |y| = inf */ 49 } else if (iy >= 0x7f800000) { 50 if (iy == 0x7f800000) 51 *(int *) &w = x == y ? ix : iy; /* w = |y| = inf */ 55 *(int *) &w = iy; /* w = |y| */ 56 else if (iy [all...] |
H A D | nextafterf.c | 40 int ix, iy, iz; local 43 iy = py[0]; 46 if ((iy & ~0x80000000) > 0x7f800000) 48 if (ix == iy || (ix | iy) == 0x80000000) 51 iz = 1 | (iy & 0x80000000); 53 if (ix > iy) 58 if (iy < 0 && ix < iy)
|
H A D | atan2pif.c | 35 int ix, iy, hx, hy; local 38 iy = *(int *)&y; 40 hy = iy & ~0x80000000; 47 return ((iy == 0)? 1.0f : -1.0f);
|
H A D | scalbf.c | 35 int ix, iy, hx, hy, n; local 38 iy = *(int *)&y; 40 hy = iy & ~0x80000000; 44 return ((iy < 0)? x / -y : x * y); 50 n = (iy < 0)? -65000 : 65000;
|
H A D | fmodf.c | 46 int hx, ix, iy, iz, k, ny, nd; local 50 iy = *(int *)&y & 0x7fffffff; 53 if (ix >= ii || iy > ii || iy == 0) { 56 } else if (ix <= iy) { 57 if (ix < iy) 69 ny = iy >> 23; 75 while (iy < iu) { 77 iy += iy; [all...] |
/illumos-gate/usr/src/lib/libm/common/C/ |
H A D | fmod.c | 45 long long hx, ix, iy, iz; local 50 iy = *(long long *)&y & ~0x8000000000000000ull; 53 if (iy == 0ll) 56 if (ix >= 0x7ff0000000000000ll || iy > 0x7ff0000000000000ll) 59 if (ix <= iy) 60 return ((ix < iy)? x : x * zero); 67 * iy = normalized significand of y 69 ny = iy >> 52; 74 while (iy < 0x0010000000000000ll) { 76 iy [all...] |
H A D | atan2pi.c | 39 int ix, iy; local 43 iy = ((int *)&y)[HIWORD]; 46 return ((iy >= 0)? 1.0 : -1.0);
|
H A D | hypot.c | 81 int i, j, nx, ny, ix, iy, iscale = 0; local 86 iy = ((int *) &y)[HIWORD] & ~0x80000000; 91 if (iy > ix) { 95 ix = iy; 96 iy = i; 105 ny = iy >> 20; 114 else if (((iy - 0x7ff00000) | ly) == 0) 130 iy -= 768 << 20; 138 if ((iy | ly) == 0) 154 iy [all...] |
/illumos-gate/usr/src/lib/libm/common/complex/ |
H A D | carg.c | 40 int ix, iy; local 43 iy = ((int *)&(D_IM(z)))[HIWORD]; 44 if ((((ix | iy) & ~0x80000000) | ((int *)&(D_RE(z)))[LOWORD] | 49 return ((iy == 0)? pi + pi_lo : -pi - pi_lo);
|
H A D | cprojf.c | 42 int ix, iy, hx, hy; local 49 iy = hy & 0x7fffffff; 50 if (iy == 0x7f800000) {
|
H A D | clogf.c | 43 int ix, iy, hx, hy; local 50 iy = hy & 0x7fffffff; 54 if (ix >= 0x7f800000 || iy >= 0x7f800000) { 56 if (iy == 0x7f800000) 70 else if (iy == 0x3f800000) 72 else if ((ix | iy) == 0)
|
H A D | clog.c | 37 * log(x+iy) = log(\/ x + y ) + i tan (---) 52 * clog(-inf + iy )= +inf + i pi, for finite positive-signed y 53 * clog(+inf + iy )= +inf + i 0 , for finite positive-signed y 74 int n, ix, iy, hx, hy; local 84 iy = hy & 0x7fffffff; 88 if (ix < iy || (ix == iy && lx < ly)) { 93 n = ix, ix = iy; 94 iy = n; 98 n = (ix - iy) >> 2 [all...] |
H A D | cabs.c | 54 int i, j, nx, ny, ix, iy, iscale = 0; local 62 iy = ((int *)&y)[HIWORD] & ~0x80000000; 66 if (iy > ix) { 70 ix = iy; 71 iy = i; 80 ny = iy >> 20; 89 else if (((iy - 0x7ff00000) | ly) == 0) 104 iy -= 768 << 20; 109 if ((iy | ly) == 0) 126 iy [all...] |
H A D | cproj.c | 51 int ix, iy, hx, hy, lx, ly; local 60 iy = hy & 0x7fffffff; 61 if (ISINF(iy, ly)) {
|
H A D | catanf.c | 50 int hx, hy, ix, iy; local 59 iy = hy & 0x7fffffff; 67 if (iy == 0 || iy == 0x7f800000) 72 } else if (iy >= 0x7f800000) { /* y is inf or NaN */ 73 if (iy == 0x7f800000) { 94 if (iy == 0x3f800000) { 98 } else if (iy > 0x3f800000) { /* y>1 */
|
H A D | clogl.c | 52 int n, ix, iy, hx, hy; local 59 iy = hy & 0x7fffffff; 63 if (ix < iy || (ix == iy && ix < 0x7fff0000 && ax < ay)) { 68 n = ix, ix = iy; 69 iy = n; 71 n = (ix - iy) >> 16; 97 else if (ix < 0x5f3f0000 && iy >= 0x20bf0000)
|
H A D | csqrt.c | 37 * Let w=r+i*s = sqrt(x+iy). Then (r + i s) = r - s + i 2sr = x + i y. 96 * csqrt(-inf+ iy ) = 0 + i inf for finite positive-signed y 97 * csqrt(+inf+ iy ) = inf + i 0 for finite positive-signed y 123 int n, ix, iy, hx, hy, lx, ly; local 132 iy = hy & 0x7fffffff; 135 if (ix >= 0x7ff00000 || iy >= 0x7ff00000) { 137 if (ISINF(iy, ly)) 149 } else if ((iy | ly) == 0) { /* y = 0 */ 157 } else if (ix >= iy) { 158 n = (ix - iy) >> 2 [all...] |
H A D | catan.c | 37 * z = x + iy, 95 int hx, hy, ix, iy; local 107 iy = hy & 0x7fffffff; 116 if ((iy | ly) == 0 || (ISINF(iy, ly))) 121 } else if (iy >= 0x7ff00000) { 123 if (ISINF(iy, ly)) { 144 if (((iy - 0x3ff00000) | ly) == 0) { 148 } else if (iy >= 0x3ff00000) { /* y>1 */ 155 } else if (iy < [all...] |
H A D | cexpl.c | 46 int n, ix, iy, hx, hy; local 53 iy = hy & 0x7fffffff; 59 if (iy >= 0x7fff0000) { 68 if (iy >= 0x7fff0000) {
|
H A D | cexp.c | 36 * x+iy x 43 * Thus if exp(x+iy) = A + Bi and t = __k_cexp(x,&n), then 70 int n, ix, iy, hx, hy, lx, ly; local 79 iy = hy & 0x7fffffff; 80 if ((iy | ly) == 0) { /* y = 0 */ 85 if (iy >= 0x7ff00000) { 94 if (iy >= 0x7ff00000) {
|
H A D | cexpf.c | 45 int n, ix, iy, hx, hy; local 52 iy = hy & 0x7fffffff; 53 if (iy == 0) { /* y = 0 */ 58 if (iy >= 0x7f800000) { 67 if (iy >= 0x7f800000) {
|
H A D | csqrtf.c | 44 int ix, iy, hx, hy; local 51 iy = hy & 0x7fffffff; 54 if (ix >= 0x7f800000 || iy >= 0x7f800000) { 56 if (iy == 0x7f800000) 68 } else if (iy == 0) {
|
H A D | csqrtl.c | 51 int n, ix, iy, hx, hy; local 58 iy = hy & 0x7fffffff; 61 if (ix >= 0x7fff0000 || iy >= 0x7fff0000) { 83 } else if (ix >= iy) { 84 n = (ix - iy) >> 16; 95 } else if (iy <= 0x20bf0000) { /* y < 2**-8000 */ 110 n = (iy - ix) >> 16; 118 else if (iy >= 0x7ffe0000) 124 } else if (iy >= 0x5f3f0000) { /* y > 2**8000 */
|
/illumos-gate/usr/src/lib/libm/common/m9x/ |
H A D | remquof.c | 58 int hx, ix, iy, iz, k, ny, nd, m, sq; local 62 iy = *(int *) &y; 63 sq = (iy ^ hx) & is; /* sign of x/y */ 64 iy &= 0x7fffffff; 68 if (ix >= ii || iy > ii || iy == 0) { 71 } else if (ix <= iy) { 72 if (ix < iy) 86 ny = iy >> 23; 92 while (iy < i [all...] |
H A D | remquo.c | 61 int n, hx, hy, hz, ix, iy, sx, sq, i, m; local 100 /* determine iy = ilogb(y) */ 103 for (iy = -1043, i = ly; i > 0; i <<= 1) 104 iy -= 1; 106 for (iy = -1022, i = (hy << 11); i > 0; i <<= 1) 107 iy -= 1; 110 iy = (hy >> 20) - 1023; 125 if (iy >= -1022) 128 n = -1022 - iy; 139 n = ix - iy; [all...] |