/titanic_44/usr/src/lib/libm/common/R/ |
H A D | fmodf.c | 46 int hx, ix, iy, iz, k, ny, nd; in fmodf() local 49 ix = hx & 0x7fffffff; in fmodf() 53 if (ix >= ii || iy > ii || iy == 0) { in fmodf() 56 } else if (ix <= iy) { in fmodf() 57 if (ix < iy) in fmodf() 70 k = ix >> 23; in fmodf() 82 while (ix < iu) { in fmodf() 84 ix += ix; in fmodf() 87 ix = iu | (ix & im); in fmodf() 91 ix = iu | (ix & im); in fmodf() [all …]
|
H A D | expf.c | 347 int hx, ix, n; in expf() local 350 ix = hx & ~0x80000000; in expf() 352 if (ix < 0x3c800000) { /* |x| < 2**-6 */ in expf() 353 if (ix < 0x38800000) /* |x| < 2**-14 */ in expf() 358 n = ix >> 23; /* biased exponent */ in expf() 372 ix -= n << 23; in expf() 374 ix += 0x800000; in expf() 376 ix = 0x800000 - ix; in expf() 378 ix <<= n - 0x7f; in expf() 379 w = ET[(ix & 0x3f) + 64] * ET[((ix >> 6) & 0x3f) + 128]; in expf() [all …]
|
H A D | nextafterf.c | 40 int ix, iy, iz; in nextafterf() local 42 ix = px[0]; in nextafterf() 44 if ((ix & ~0x80000000) > 0x7f800000) in nextafterf() 48 if (ix == iy || (ix | iy) == 0x80000000) in nextafterf() 50 if ((ix & ~0x80000000) == 0) in nextafterf() 52 else if (ix > 0) { in nextafterf() 53 if (ix > iy) in nextafterf() 54 iz = ix - 1; in nextafterf() 56 iz = ix + 1; in nextafterf() 58 if (iy < 0 && ix < iy) in nextafterf() [all …]
|
H A D | atanf.c | 108 int ix, iy, sign, j; in atanf() local 111 ix = *(int *) &x; in atanf() 112 sign = ix & 0x80000000; in atanf() 113 ix ^= sign; in atanf() 116 if (ix < 0x3e000000) { in atanf() 117 if (ix < 0x38800000) { /* if |x| < 2**(-prec/2-2) */ in atanf() 125 if (ix < 0x3c000000) { /* if |x| < 2**(-prec/4-1) */ in atanf() 135 if (ix >= 0x41000000) { in atanf() 136 *(int *) &x = ix; in atanf() 137 if (ix < 0x42820000) { /* x < 65 */ in atanf() [all …]
|
H A D | besself.c | 272 int ix, i; in __k_j0f() local 274 ix = *(int *)&fx & ~0x80000000; in __k_j0f() 276 if (ix > 0x41000000) { in __k_j0f() 287 if (ix > 0x501502f9) { in __k_j0f() 305 if (ix <= 0x3727c5ac) { in __k_j0f() 307 if (ix <= 0x219392ef) /* x <= 1.0e-18 */ in __k_j0f() 312 if (ix <= 0x3fa3d70a) { in __k_j0f() 331 int ix; in j0f() local 336 ix = *(int *)&fx & ~0x80000000; in j0f() 337 if (ix >= 0x7f800000) { /* nan or inf */ in j0f() [all …]
|
H A D | logf.c | 118 int hx, ix, i, exp, iy; in logf() local 121 ix = hx & ~0x80000000; in logf() 123 if (ix >= 0x7f800000) /* nan or inf */ in logf() 130 return ((ix == 0)? -1.0f / f : f / f); in logf() 134 f = (float)ix; in logf() 135 ix = *(int *)&f; in logf() 139 exp += (ix - 0x3f320000) >> 23; in logf() 140 ix &= 0x007fffff; in logf() 141 iy = (ix + 0x20000) & 0xfffc0000; in logf() 144 v = (double)(ix - iy) * TBL[i + 1]; in logf()
|
H A D | hypotf.c | 38 int ix, iy; in hypotf() local 40 ix = (*(int *) &x) & 0x7fffffff; in hypotf() 42 if (ix >= 0x7f800000) { in hypotf() 43 if (ix == 0x7f800000) in hypotf() 44 *(int *) &w = x == y ? iy : ix; /* w = |x| = inf */ in hypotf() 46 *(int *) &w = x == y ? ix : iy; /* w = |y| = inf */ in hypotf() 51 *(int *) &w = x == y ? ix : iy; /* w = |y| = inf */ in hypotf() 54 } else if (ix == 0) in hypotf() 57 *(int *) &w = ix; /* w = |x| */ in hypotf()
|
H A D | log1pf.c | 37 int ix; in log1pf() local 39 ix = *(int *)&x; in log1pf() 40 if (ix >= 0x7f800000) { in log1pf() 44 if (ix < 0) { in log1pf() 45 ix &= ~0x80000000; in log1pf() 46 if (ix == 0x3f800000) /* x is -1 */ in log1pf() 48 if (ix > 0x3f800000) /* x is < -1 or nan */ in log1pf()
|
/titanic_44/usr/src/lib/libm/common/m9x/ |
H A D | remquof.c | 58 int hx, ix, iy, iz, k, ny, nd, m, sq; in fmodquof() local 61 ix = hx & 0x7fffffff; in fmodquof() 68 if (ix >= ii || iy > ii || iy == 0) { in fmodquof() 71 } else if (ix <= iy) { in fmodquof() 72 if (ix < iy) in fmodquof() 87 k = ix >> 23; in fmodquof() 99 while (ix < iu) { in fmodquof() 101 ix += ix; in fmodquof() 104 ix = iu | (ix & im); in fmodquof() 107 ix = iu | (ix & im); in fmodquof() [all …]
|
H A D | tgammaf.c | 311 int k, ix, j, m; in large_gam() local 313 ix = __HI(x); in large_gam() 315 m = (ix >> 20) - 0x3ff; /* exponent of x, range:3-5 */ in large_gam() 316 ix = (ix & 0x000fffff) | 0x3ff00000; /* y = scale x to [1,2] */ in large_gam() 317 __HI(y) = ix; in large_gam() 319 __HI(z) = (ix & 0xffffc000) | 0x2000; /* z[j]=1+j/64+1/128 */ in large_gam() 321 j = (ix >> 14) & 0x3f; in large_gam() 456 int i, j, k, ix, hx, xk; in tgammaf() local 459 ix = hx & 0x7fffffff; in tgammaf() 462 if (ix < 0x33800000) in tgammaf() [all …]
|
H A D | ldexp.c | 37 int *px = (int *) &x, ix = px[HIWORD] & ~0x80000000; in ldexp() local 39 if (ix >= 0x7ff00000 || (px[LOWORD] | ix) == 0) in ldexp() 41 return (ix >= 0x7ff80000 ? x : x + x); in ldexp() 47 ix = px[HIWORD] & ~0x80000000; in ldexp() 52 if (ix >= 0x7ff00000 || (px[LOWORD] | ix) == 0) in ldexp()
|
/titanic_44/usr/src/lib/libm/common/C/ |
H A D | fmod.c | 45 long long hx, ix, iy, iz; in fmod() local 49 ix = hx & ~0x8000000000000000ull; in fmod() 56 if (ix >= 0x7ff0000000000000ll || iy > 0x7ff0000000000000ll) in fmod() 59 if (ix <= iy) in fmod() 60 return ((ix < iy)? x : x * zero); in fmod() 70 k = ix >> 52; in fmod() 81 while (ix < 0x0010000000000000ll) { in fmod() 83 ix += ix; in fmod() 86 ix = 0x0010000000000000ll | (ix & 0x000fffffffffffffll); in fmod() 91 ix = 0x0010000000000000ll | (ix & 0x000fffffffffffffll); in fmod() [all …]
|
H A D | log2.c | 153 int i, hx, ix, n, lx; in log2() local 156 hx = ((int *) &x)[HIWORD]; ix = hx & 0x7fffffff; in log2() 162 if (ix >= 0x7ff80000) /* assumes sparc-like QNaN */ in log2() 168 if (ix >= 0x7ff00000) in log2() 175 if (ix >= 0x7ff00000) /* log(NaN) = NaN */ in log2() 179 hx = ((int *) &x)[HIWORD]; ix = hx & 0x7fffffff; in log2() 184 i = ix >> 19; in log2() 187 if (ix >= 0x3fec0000 && ix < 0x3ff20000) { in log2() 191 if (((ix << 12) | lx) == 0) in log2() 200 i = (ix - 0x3fb80000) >> 15; in log2() [all …]
|
H A D | hypot.c | 81 int i, j, nx, ny, ix, iy, iscale = 0; in hypot() local 84 ix = ((int *) &x)[HIWORD] & ~0x80000000; in hypot() 91 if (iy > ix) { in hypot() 94 i = ix; in hypot() 95 ix = iy; in hypot() 104 nx = ix >> 20; in hypot() 112 if (((ix - 0x7ff00000) | lx) == 0) in hypot() 129 ix -= 768 << 20; in hypot() 136 if ((ix | lx) == 0) in hypot() 148 ix = ((int *) &ax)[HIWORD]; in hypot() [all …]
|
H A D | scalbn.c | 77 int *px, ix, hx, k; in scalbn() local 80 ix = px[HIWORD]; in scalbn() 81 hx = ix & ~0x80000000; in scalbn() 92 ix = px[HIWORD]; in scalbn() 93 k = ((ix & ~0x80000000) >> 20) - 52; in scalbn() 96 ix = px[HIWORD]; in scalbn() 110 return (huge * ((ix < 0)? -huge : huge)); in scalbn() 113 return (tiny * ((ix < 0)? -tiny : tiny)); in scalbn() 115 px[HIWORD] = (ix & ~0x7ff00000) | (k << 20); in scalbn() 118 px[HIWORD] = (ix & ~0x7ff00000) | (k << 20); in scalbn()
|
H A D | cos.c | 113 int i, j, n, hx, ix, lx; in cos() local 117 ix = hx & ~0x80000000; in cos() 119 if (ix <= 0x3fc50000) { /* |x| < 10.5/64 = 0.164062500 */ in cos() 120 if (ix < 0x3e400000) { /* |x| < 2**-27 */ in cos() 125 if (ix < 0x3f800000) /* |x| < 0.008 */ in cos() 133 n = ix >> 20; in cos() 135 i = (((ix >> 12) & 0xff) | 0x100) >> (0x401 - n); in cos() 142 i = ix - 0x3ff921fb; in cos() 150 if (((ix - 0x3ff92000) >> 12) == 0) { in cos() 163 i = ix - 0x4012D97C; in cos() [all …]
|
H A D | log.c | 159 int i, hx, ix, n, lx; in log() local 163 ix = hx & 0x7fffffff; in log() 168 if (ix > 0x7ff00000 || (ix == 0x7ff00000 && lx != 0)) /* nan */ in log() 180 ix = ((int *)&x)[HIWORD]; in log() 184 i = ix >> 19; in log() 187 if (ix >= 0x3fec0000 && ix < 0x3ff22000) { in log() 191 if (((ix - 0x3ff00000) | lx) == 0) /* x = 1 */ in log() 200 i = (ix - 0x3fb80000) >> 15; in log() 208 dn = (double)(n + ((ix >> 20) - 0x3ff)); in log() 210 i = (ix & 0x000fffff) | 0x3ff00000; /* scale x to [1,2] */ in log()
|
/titanic_44/usr/src/lib/libm/common/Q/ |
H A D | logl.c | 111 int i, j, ix, i0, i1, n; in logl() local 123 ix = px[i0]; in logl() 124 if (ix > 0x3ffee0f8) { /* if x > 31/33 */ in logl() 125 if (ix < 0x3fff1084) { /* if x < 33/31 */ in logl() 128 if (((ix - 0x3fff0000) | px[i1] | px[2] | px[1]) == 0) { in logl() 137 if (ix >= 0x7fff0000) in logl() 141 if (ix >= 0x00010000) in logl() 143 i = ix & 0x7fffffff; in logl() 148 if (ix < 0) { in logl() 149 if ((unsigned) ix >= 0xffff0000) in logl() [all …]
|
H A D | atanl.c | 115 int ix, iy, sign, j; in atanl() local 117 ix = px[i0]; in atanl() 118 sign = ix & 0x80000000; in atanl() 119 ix ^= sign; in atanl() 122 if (ix < 0x3ffc0000) { in atanl() 123 if (ix < 0x3feb0000) { /* when |x| < 2**(-prec/6-2) */ in atanl() 124 if (ix < 0x3fc50000) { /* if |x| < 2**(-prec/2-2) */ in atanl() 134 if (ix < 0x3fe20000) { /* if |x| < 2**(-prec/4-1) */ in atanl() 147 if (ix >= 0x40020000) { in atanl() 148 px[i0] = ix; in atanl() [all …]
|
/titanic_44/usr/src/lib/libm/common/complex/ |
H A D | clog.c | 74 int n, ix, iy, hx, hy; in __clog() local 83 ix = hx & 0x7fffffff; in __clog() 88 if (ix < iy || (ix == iy && lx < ly)) { in __clog() 93 n = ix, ix = iy; in __clog() 98 n = (ix - iy) >> 20; in __clog() 99 if (ix >= 0x7ff00000) { /* x or y is Inf or NaN */ in __clog() 100 if (ISINF(ix, lx)) in __clog() 107 D_RE(ans) = ((ix | lx) == 0)? -one / ax : log(ax); in __clog() 108 } else if (((0x3fffffff - ix) ^ (ix - 0x3fe00000)) >= 0) { in __clog() 110 if (ix >= 0x3ff00000) { in __clog() [all …]
|
H A D | clogl.c | 52 int n, ix, iy, hx, hy; in clogl() local 58 ix = hx & 0x7fffffff; in clogl() 63 if (ix < iy || (ix == iy && ix < 0x7fff0000 && ax < ay)) { in clogl() 68 n = ix, ix = iy; in clogl() 71 n = (ix - iy) >> 16; in clogl() 72 if (ix >= 0x7fff0000) { /* x or y is Inf or NaN */ in clogl() 81 else if (((0x3fffffff - ix) ^ (ix - 0x3ffe0000)) >= 0) { in clogl() 83 if (ix >= 0x3fff0000) { in clogl() 97 else if (ix < 0x5f3f0000 && iy >= 0x20bf0000) in clogl()
|
H A D | cabs.c | 54 int i, j, nx, ny, ix, iy, iscale = 0; in cabs() local 60 ix = ((int *)&x)[HIWORD] & ~0x80000000; in cabs() 66 if (iy > ix) { in cabs() 69 i = ix; in cabs() 70 ix = iy; in cabs() 79 nx = ix >> 20; in cabs() 87 if (((ix - 0x7ff00000) | lx) == 0) in cabs() 103 ix -= 768 << 20; in cabs() 107 if ((ix | lx) == 0) in cabs() 119 ix = ((int *)&ax)[HIWORD]; in cabs() [all …]
|
H A D | catan.c | 95 int hx, hy, ix, iy; in catan() local 106 ix = hx & 0x7fffffff; in catan() 110 if (ix >= 0x7ff00000) { in catan() 111 if (ISINF(ix, lx)) { in catan() 130 } else if ((ix | lx) == 0) { in catan() 155 } else if (iy < 0x3e200000 || ((ix - iy) >> 20) >= 30) { in catan() 172 if (ix < 0x3ff00000) in catan() 180 if (ix < 0x3e200000) in catan() 182 else if (ix > 0x41c00000) in catan() 188 } else if (iy >= 0x41c00000 && ((iy - ix) >> 20) >= 30) { in catan() [all …]
|
/titanic_44/usr/src/common/mpi/ |
H A D | mplogic.c | 128 mp_size ix; in mpl_set_bit() local 134 ix = bitNum / MP_DIGIT_BIT; in mpl_set_bit() 135 if (ix + 1 > MP_USED(a)) { in mpl_set_bit() 136 rv = s_mp_pad(a, ix + 1); in mpl_set_bit() 144 MP_DIGIT(a,ix) |= mask; in mpl_set_bit() 146 MP_DIGIT(a,ix) &= ~mask; in mpl_set_bit() 158 mp_size bit, ix; in mpl_get_bit() local 163 ix = bitNum / MP_DIGIT_BIT; in mpl_get_bit() 164 ARGCHK(ix <= MP_USED(a) - 1, MP_RANGE); in mpl_get_bit() 167 rv = (mp_err)(MP_DIGIT(a, ix) >> bit) & 1; in mpl_get_bit() [all …]
|
H A D | mpi.c | 356 int ix; in mp_set_int() local 369 for (ix = sizeof(long) - 1; ix >= 0; ix--) { in mp_set_int() 373 res = s_mp_add_d(mp, (mp_digit)((v >> (ix * CHAR_BIT)) & UCHAR_MAX)); in mp_set_int() 391 int ix; in mp_set_ulong() local 403 for (ix = sizeof(long) - 1; ix >= 0; ix--) { in mp_set_ulong() 407 res = s_mp_add_d(mp, (mp_digit)((z >> (ix * CHAR_BIT)) & UCHAR_MAX)); in mp_set_ulong() 925 mp_size ix; in mp_sqr() local 940 ix = 2 * MP_USED(a); in mp_sqr() 941 if (ix > MP_ALLOC(sqr)) { in mp_sqr() 943 MP_CHECKOK( s_mp_grow(sqr, ix) ); in mp_sqr() [all …]
|