| /illumos-gate/usr/src/lib/libm/common/m9x/ |
| H A D | fmaf.c | 46 } xy, zz; local
56 zz.d = (double) z;
63 ez = (zz.i[0] >> 20) & 0x7ff;
66 return ((float) (xy.d + zz.d));
94 u = zz.i[1];
98 zz.i[1] = u & ~(s - 1);
100 u = zz.i[0];
102 if ((u & (s - 1)) | zz.i[1])
104 zz.i[0] = u & ~(s - 1);
105 zz
132 } xy, zz; local
[all...] |
| H A D | llroundl.c | 48 } zz; local
78 zz.i[0] = 0;
79 zz.i[1] = xx.i[0] >> (j - 96);
84 zz.i[0] = xx.i[0] >> (j - 64);
85 zz.i[1] = ((xx.i[0] << 1) << (95 - j)) | (xx.i[1] >> (j - 64));
90 zz.i[0] = ((xx.i[0] << 1) << (63 - j)) | (xx.i[1] >> (j - 32));
91 zz.i[1] = ((xx.i[1] << 1) << (63 - j)) | (xx.i[2] >> (j - 32));
99 if (++zz.i[1] == 0)
100 zz.i[0]++;
104 if (zz
[all...] |
| H A D | llrintl.c | 52 } zz; local
95 zz.i[0] = 0;
96 zz.i[1] = xx.i[0] >> (j - 96);
101 zz.i[0] = xx.i[0] >> (j - 64);
102 zz.i[1] = ((xx.i[0] << 1) << (95 - j)) | (xx.i[1] >> (j - 64));
107 zz.i[0] = ((xx.i[0] << 1) << (63 - j)) | (xx.i[1] >> (j - 32));
108 zz.i[1] = ((xx.i[1] << 1) << (63 - j)) | (xx.i[2] >> (j - 32));
116 (frac == 0x80000000 && (zz.i[1] & 1)))))) {
117 if (++zz.i[1] == 0)
118 zz
[all...] |
| H A D | fma.c | 75 } xx, yy, zz; local
84 zz.d = z;
87 hz = zz.i[0] & ~0x80000000;
97 if ((hz | zz.i[1]) == 0) /* z is zero */
114 sz = zz.i[0] & 0x80000000;
127 zz.d = z * two52;
128 ez = ((zz.i[0] & ~0x80000000) >> 20) - 52;
157 z0 = (zz.i[0] & 0xfffff) | 0x100000;
158 z1 = zz.i[1];
310 zz
406 } xx, yy, zz; local
[all...] |
| H A D | fmal.c | 87 } xx, yy, zz; local
104 zz.q = z;
107 hz = zz.i[0] & ~0x80000000;
148 if ((hz & 0xffff) | zz.i[1] | zz.i[2] | zz.i[3]) {
152 zz.i[0] |= 0x8000;
153 return (zz.q);
159 cz = (zz.i[1] | zz
878 } xx, yy, zz; local
[all...] |
| H A D | tgammaf.c | 309 double ss, zz, z, t1, t2, w, y, u; local
332 zz = z * z;
333 w = ((x - half) * ss + hln2pi) + z * (GP0 + zz * GP1 + (zz * zz) * GP2);
340 zz = S[j] * (one + z + (z * z) * (Et1 + z * Et2));
341 __HI(zz) += m << 20;
342 return (zz);
|
| H A D | tgamma.c | 1215 struct Double zz; local
1310 zz.l = S_trail[j] * (one + t3) + S[j] * t3;
1311 zz.h = S[j];
1312 return (zz);
|
| H A D | tgammal.c | 630 struct LDouble zz; local
724 zz.l = S_trail[j] * (one + t3) + S[j] * t3;
725 zz.h = S[j];
726 return (zz);
|
| /illumos-gate/usr/src/lib/libc/sparc/fp/ |
| H A D | _Q_div.c | 107 double c, d, ry, xx[4], yy[5], zz[5]; local
341 zz[0] = (ry * (xx[0] + xx[1]) + three2p27) - three2p27;
342 xx[0] = ((xx[0] - zz[0] * yy[0]) - zz[0] * yy[1]) + xx[1];
343 d = zz[0] * yy[2];
347 d = zz[0] * yy[3];
351 d = zz[0] * yy[4];
356 zz[1] = (ry * (xx[0] + xx[1]) + three2p3) - three2p3;
357 xx[0] = ((xx[0] - zz[1] * yy[0]) - zz[
[all...] |
| H A D | _Q_mul.c | 89 double xx[5], yy[5], c, d, zz[9]; local
306 zz[0] = xx[0] * yy[0];
307 zz[1] = xx[0] * yy[1] + xx[1] * yy[0];
308 zz[2] = xx[0] * yy[2] + xx[1] * yy[1] + xx[2] * yy[0];
309 zz[3] = xx[0] * yy[3] + xx[1] * yy[2] + xx[2] * yy[1] +
311 zz[4] = xx[0] * yy[4] + xx[1] * yy[3] + xx[2] * yy[2] +
313 zz[5] = xx[1] * yy[4] + xx[2] * yy[3] + xx[3] * yy[2] +
315 zz[6] = xx[2] * yy[4] + xx[3] * yy[3] + xx[4] * yy[2];
316 zz[7] = xx[3] * yy[4] + xx[4] * yy[3];
317 zz[
[all...] |
| H A D | _Q_sqrt.c | 98 double c, d, rr, r[2], tt[3], xx[4], zz[5]; local
237 zz[0] = (c + two29) - two29;
238 tt[0] = zz[0] + zz[0];
239 r[0] = (xx[0] - zz[0] * zz[0]) + xx[1];
241 zz[1] = (rr * (r[0] + xx[2]) + three2p4) - three2p4;
242 tt[1] = zz[1] + zz[1];
243 r[0] -= tt[0] * zz[
[all...] |
| /illumos-gate/usr/src/lib/libm/common/Q/ |
| H A D | sqrtl.c | 401 double xx[5], zz[3]; local
463 __q_tp_sqrt(xx, zz);
468 __q_pack(zz, ex >> 1, rm, &x, &inexact);
|
| /illumos-gate/usr/src/lib/libm/common/complex/ |
| H A D | k_atan2l.c | 574 long double y, z, r, s, t, w, s_h, s_l, x_h, x_l, zz[3], ee[2], z_h, local
608 zz[0] = z;
609 zz[1] = z_h;
610 zz[2] = z_l + x_l * (x + x_h);
614 mx_polyl(zz, pe, ee, 3);
648 zz[0] = z;
649 zz[1] = z_h;
650 zz[2] = z_l + r_l * (r + r_h);
657 mx_polyl(zz, pe, ee, 3);
664 mx_polyl(zz, q
[all...] |
| H A D | k_atan2.c | 283 double y, z, r, s, t, w, s_h, s_l, x_h, x_l, zz[3], ee[2], z_h, local
315 zz[0] = z;
316 zz[1] = z_h;
317 zz[2] = z_l + x_l * (x + x_h);
324 mx_poly(zz, p, ee, 3);
365 zz[0] = z;
366 zz[1] = z_h;
367 zz[2] = z_l + r_l * (r + r_h);
372 mx_poly(zz, p, ee, 3);
381 zz[
[all...] |
| /illumos-gate/usr/src/common/mpi/ |
| H A D | mp_gf2m.c | 357 mp_digit zz, *z, tmp; local
375 zz = z[j];
376 if (zz == 0) {
387 z[j-n] ^= (zz>>d0);
389 z[j-n-1] ^= (zz<<d1);
396 z[j-n] ^= (zz >> d0);
398 z[j-n-1] ^= (zz << d1);
406 zz = z[dN] >> d0;
407 if (zz == 0) break;
412 *z ^= zz; /* reductio
[all...] |