/illumos-gate/usr/src/lib/libm/common/complex/ |
H A D | ctanhl.c | 43 int hx, ix, hy, iy; local 51 iy = hy & 0x7fffffff; 58 } else if (iy >= 0x7fff0000) { /* y is inf or NaN */ 76 if (iy < 0x7ffe0000) /* t = sin(2y) */
|
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 | 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 | 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 | catanl.c | 129 int hx, hy, ix, iy; local 138 iy = hy & 0x7fffffff; 152 } else if (iy >= 0x7fff0000) { 311 if (iy >= 0x3ffe0000 && iy < 0x40000000) {
|
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 | 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 */
|
H A D | ctanh.c | 100 int hx, ix, lx, hy, iy, ly; local 110 iy = hy & 0x7fffffff; 114 if ((iy | ly) == 0) { /* ctanh(x,0) = (x,0) for x = 0 or NaN */ 117 } else if (iy >= 0x7ff00000) { /* y is inf or NaN */ 133 if (iy < 0x7fe00000) /* t = sin(2y) */
|
H A D | cacos.c | 42 * where z = x+iy, can be defined by 221 int ix, iy, hx, hy; local 232 iy = hy & 0x7fffffff; 236 if (((iy | ly) == 0) || (iy >= 0x7ff00000)) { 244 if (iy >= 0x7ff00000) { 245 if (ISINF(iy, ly)) { /* cacos(x + i inf) = pi/2 - i inf */ 274 if (iy >= 0x7ff00000) { 275 if (ISINF(iy, ly)) { 291 /* cacos(inf + iy ) [all...] |
H A D | cacosl.c | 77 int ix, iy, hx, hy; local 85 iy = hy & 0x7fffffff; 89 if (y == zero || (iy >= 0x7fff0000)) { 97 if (iy >= 0x7fff0000) { 125 if (iy >= 0x7fff0000) { 142 /* cacos(inf + iy ) = 0 - i inf */ 143 /* cacos(-inf+ iy ) = pi - i inf */ 159 if (iy >= 0x7fff0000) { 222 } else if (ix >= 0x5ffb0000 || iy >= 0x5ffb0000) { 231 if (iy > [all...] |
H A D | casin.c | 42 * where z = x+iy, can be defined by 78 * casin(inf + iy ) = pi/2 + i inf finite y 217 int ix, iy, hx, hy; local 228 iy = hy & 0x7fffffff; 238 if (iy >= 0x7ff00000) { 239 if (ISINF(iy, ly)) 244 } else /* casin(inf + iy) = pi/2 + i inf */ 247 if (iy >= 0x7ff00000) { 269 if ((ix | lx | iy | ly) == 0) 272 if (iy > [all...] |
H A D | casinl.c | 71 int ix, iy, hx, hy; local 79 iy = hy & 0x7fffffff; 89 if (iy >= 0x7fff0000) { 95 } else /* casin(inf + iy) = pi/2 + i inf */ 98 if (iy >= 0x7fff0000) { 125 if (iy >= 0x7fff0000) { /* y is inf or NaN */ 182 } else if (ix >= 0x5ffb0000 || iy >= 0x5ffb0000) { 191 if (iy >= 0x3fff8000) /* if y > Acrossover */
|
H A D | cpowf.c | 58 int ix, iy, hx, hy, hv, hu, iu, iv, j; local 69 iy = hy & 0x7fffffff; 81 } else if (iy == 0) { /* (real)**(real) */ 92 } else if (ix == 0 || ix >= hiinf || iy >= hiinf) { 93 if (ix > hiinf || iy > hiinf || iu > hiinf) { 104 } else if (ix == iy) { /* if |x| == |y| */ 127 if (iu >= hiinf || iv >= hiinf || ix >= hiinf || iy >= hiinf) { 153 else if (iy > 0x3f000000 && iy < 0x40000000) /* .5 < |y| < 2 */
|
H A D | cpow.c | 40 * Let z = x+iy, w = u+iv. 44 * log(x+iy) = log(\/ x + y ) + i tan (---) 51 * (u+iv)* log(x+iy) = --- log(x + y ) - v tan (---) + (1) 156 int ix, iy, hx, lx, hy, ly, hv, hu, iu, iv, lu, lv; local 172 iy = hy & 0x7fffffff; 184 } else if ((iy | ly) == 0) { /* (real)**(real) */ 194 } else if (((ix | lx) == 0) || ix >= hiinf || iy >= hiinf) { 207 } else if (((ix - iy) | (lx - ly)) == 0) { /* |x| = |y| */ 240 if (iu >= hiinf || iv >= hiinf || ix >= hiinf || iy >= hiinf) {
|
/illumos-gate/usr/src/lib/libm/common/C/ |
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...] |
H A D | _SVID_error.c | 127 int iy; local 928 iy = ((int *)&y)[HIWORD]; 932 else if ((x == 0.0 && iy > 0) || (!finite(x) && iy < 0)) 935 exc.retval = ieee_retval = ((iy < 0)? x / -y : x * y);
|
/illumos-gate/usr/src/lib/libm/common/R/ |
H A D | atan2f.c | 254 int i, k, hx, hy, ix, iy, sign; local 259 iy = *(int *)&fy; 261 hy = iy & ~0x80000000; 271 if (iy < 0) { 285 if (iy < 0) {
|
H A D | atanf.c | 72 * If iy is the high word of y, then 73 * single : j = (iy - 0x3e000000) >> 19 74 * (single is modified to (iy-0x3f000000)>>19) 75 * double : j = (iy - 0x3fc00000) >> 16 76 * quad : j = (iy - 0x3ffc0000) >> 12 108 int ix, iy, sign, j; local 174 iy = (ix + 0x00040000) & 0x7ff80000; 175 *(int *) &y = iy; 176 j = (iy - 0x3f000000) >> 19; 178 if (ix == iy) [all...] |
H A D | powf.c | 130 int ix, iy, jx, jy, k, iw, yisint; local 133 iy = *(int *)&y; 135 jy = iy & ~0x80000000; 176 if (iy > 0) 181 if (iy < 0) 188 if (iy < 0) 191 } else if (iy == 0x40000000) { /* y is 2 */ 193 } else if (iy == 0x3f000000) { /* y is 0.5 */ 204 if (iy < 0)
|
/illumos-gate/usr/src/lib/libm/common/m9x/ |
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...] |
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 | remquol.c | 60 int n, ix, iy, k, sx, sq, m; local 112 /* determine iy = ilogbl(y) */ 114 iy = -16382; 116 iy -= 16; 123 iy -= 1; 130 iy = (y0 >> 16) - 16383; 136 n = ix - iy; 243 iy -= 16; 253 iy -= 1; 258 if (iy > [all...] |
/illumos-gate/usr/src/lib/libm/common/Q/ |
H A D | atanl.c | 67 * If iy is the high word of y, then 68 * single : j = (iy - 0x3e000000) >> 19 69 * double : j = (iy - 0x3fc00000) >> 16 70 * quad : j = (iy - 0x3ffc0000) >> 12 115 int ix, iy, sign, j; local 187 iy = (ix + 0x00000800) & 0x7ffff000; 188 py[i0] = iy; 190 j = (iy - 0x3ffc0000) >> 12; 197 if (ix == iy)
|
H A D | fmodl.c | 58 int n, ix, iy, k, sx; local 111 /* determine iy = ilogbl(y) */ 113 iy = -16382; 115 iy -= 16; 122 iy -= 1; 129 iy = (y0 >> 16) - 16383; 134 n = ix - iy; 225 iy -= 16; 235 iy -= 1; 240 if (iy > [all...] |
H A D | log1pl.c | 144 int i, j, ix, iy, n, hx, m; local 171 iy = H0(y); 172 n = ((iy + 0x200) >> 16) - 0x3fff; 173 iy = (iy & 0x0000ffff) | 0x3fff0000; /* scale 1+x to [1,2] */ 174 H0(y) = iy; 203 i = (iy + 0x200) >> 10;
|