Lines Matching refs:x
314 #define word0(x) (x)->L[1] argument
315 #define word1(x) (x)->L[0] argument
317 #define word0(x) (x)->L[0] argument
318 #define word1(x) (x)->L[1] argument
320 #define dval(x) (x)->d argument
541 ULong x[1]; member
586 int x; local
598 x = 1 << k;
600 rv = (Bigint *)MALLOC(sizeof(Bigint) + (x-1)*sizeof(ULong));
606 len = (sizeof(Bigint) + (x-1)*sizeof(ULong) + sizeof(double) - 1)
620 rv->maxwds = x;
651 #define Bcopy(x,y) memcpy((char *)&x->sign, (char *)&y->sign, \ argument
664 ULong *x; local
667 ULong carry, *x, y; local
675 x = b->x;
680 y = *x * (ULLong)m + carry;
682 *x++ = y & FFFFFFFF;
685 xi = *x;
689 *x++ = (z << 16) + (y & 0xffff);
691 y = *x * m + carry;
693 *x++ = y & 0xffff;
705 b->x[wds++] = carry;
721 Long x, y; local
723 x = (nd + 8) / 9;
724 for(k = 0, y = 1; x > y; y <<= 1, k++) ;
727 b->x[0] = y9;
731 b->x[0] = y9 & 0xffff;
732 b->wds = (b->x[1] = y9 >> 16) ? 2 : 1;
752 (x) ULong x; in hi0bits() argument
754 (ULong x)
759 if (!(x & 0xffff0000)) {
761 x <<= 16;
763 if (!(x & 0xff000000)) {
765 x <<= 8;
767 if (!(x & 0xf0000000)) {
769 x <<= 4;
771 if (!(x & 0xc0000000)) {
773 x <<= 2;
775 if (!(x & 0x80000000)) {
777 if (!(x & 0x40000000))
792 ULong x = *y; local
794 if (x & 7) {
795 if (x & 1)
797 if (x & 2) {
798 *y = x >> 1;
801 *y = x >> 2;
805 if (!(x & 0xffff)) {
807 x >>= 16;
809 if (!(x & 0xff)) {
811 x >>= 8;
813 if (!(x & 0xf)) {
815 x >>= 4;
817 if (!(x & 0x3)) {
819 x >>= 2;
821 if (!(x & 1)) {
823 x >>= 1;
824 if (!x)
827 *y = x;
842 b->x[0] = i;
857 ULong *x, *xa, *xae, *xb, *xbe, *xc, *xc0; local
880 for(x = c->x, xa = x + wc; x < xa; x++)
881 *x = 0;
882 xa = a->x;
884 xb = b->x;
886 xc0 = c->x;
890 x = xa;
894 z = *x++ * (ULLong)y + *xc + carry;
898 while(x < xae);
906 x = xa;
910 z = (*x & 0xffff) * y + (*xc & 0xffff) + carry;
912 z2 = (*x++ >> 16) * y + (*xc >> 16) + carry;
916 while(x < xae);
920 x = xa;
925 z = (*x & 0xffff) * y + (*xc >> 16) + carry;
928 z2 = (*x++ >> 16) * y + (*xc & 0xffff) + carry;
931 while(x < xae);
938 x = xa;
942 z = *x++ * y + *xc + carry;
946 while(x < xae);
952 for(xc0 = c->x, xc = xc0 + wc; wc > 0 && !*--xc; --wc) ;
1026 ULong *x, *x1, *xe, z; local
1038 x1 = b1->x;
1041 x = b->x;
1042 xe = x + b->wds;
1048 *x1++ = *x << k | z;
1049 z = *x++ >> k1;
1051 while(x < xe);
1060 *x1++ = *x << k & 0xffff | z;
1061 z = *x++ >> k1;
1063 while(x < xe);
1069 *x1++ = *x++;
1070 while(x < xe);
1090 if (i > 1 && !a->x[i-1])
1092 if (j > 1 && !b->x[j-1])
1097 xa0 = a->x;
1099 xb0 = b->x;
1134 c->x[0] = 0;
1148 xa = a->x;
1151 xb = b->x;
1153 xc = c->x;
1207 (x) U *x; in ulp() argument
1209 (U *x)
1215 L = (word0(x) & Exp_mask) - (P-1)*Exp_msk1;
1264 xa0 = a->x;
1326 ULong *x, y, z; local
1344 x = b->x;
1360 x[0] = y | z << (32 - k);
1364 x[0] = y;
1368 b->wds = (x[1] = z) ? 2 : 1;
1372 x[0] = z;
1383 x[0] = y | z << 32 - k & 0xffff;
1384 x[1] = z >> k - 16 & 0xffff;
1385 x[2] = z >> k;
1389 x[0] = y & 0xffff;
1390 x[1] = y >> 16 | z << 16 - k & 0xffff;
1391 x[2] = z >> k & 0xffff;
1392 x[3] = z >> k+16;
1396 x[0] = y & 0xffff;
1397 x[1] = y >> 16;
1398 x[2] = z & 0xffff;
1399 x[3] = z >> 16;
1410 x[0] = z;
1414 x[0] = z & 0xffff;
1415 x[1] = z >> 16;
1420 while(!x[i])
1439 *bits = 32*i - hi0bits(x[i-1]);
1441 *bits = (i+2)*16 - hi0bits(x[i]);
1624 ULong c, x[2]; local
1629 x[0] = x[1] = 0;
1669 x[0] = x[1];
1670 x[1] = 0;
1673 x[0] = (x[0] << 4) | (x[1] >> 28);
1674 x[1] = (x[1] << 4) | c;
1676 if ((x[0] &= 0xfffff) || x[1]) {
1677 word0(rvp) = Exp_mask | x[0];
1678 word1(rvp) = x[1];
1702 ULong *x, *xe; local
1705 x = b->x;
1706 xe = x + b->wds;
1708 if (*x < (ULong)0xffffffffL) {
1709 ++*x;
1712 *x++ = 0;
1713 } while(x < xe);
1721 b->x[b->wds++] = 1;
1737 ULong *x, *x1, *xe, y; local
1740 x = x1 = b->x;
1743 xe = x + b->wds;
1744 x += n;
1747 y = *x++ >> k;
1748 while(x < xe) {
1749 *x1++ = (y | (*x << n)) & 0xffffffff;
1750 y = *x++ >> k;
1756 while(x < xe)
1757 *x1++ = *x++;
1759 if ((b->wds = x1 - b->x) == 0)
1760 b->x[0] = 0;
1771 ULong *x, *x0, x1, x2; local
1773 x = b->x;
1779 x1 = x2 = x[n];
1785 x0 = x;
1786 x += n;
1787 while(x > x0)
1788 if (*--x)
1811 ULong L, lostbits, *x; local
1982 x = b->x;
1999 *x++ = L;
2006 *x++ = L;
2007 b->wds = n = x - b->x;
2011 x = b->x;
2017 if (x[k>>kshift] & 1 << (k & kmask)) {
2030 x = b->x;
2077 if (x[k>>kshift] & 1 << (k & kmask))
2090 && (lostbits & 1) | (x[0] & 1))
2102 x = b->x;
2106 && x[nbits >> kshift] & 1 << (nbits & kmask))
2112 && hi0bits(x[k-1]) < 32-n)) {
2121 word0(rvp) = b->wds > 1 ? b->x[1] & ~0x100000 : 0;
2123 word0(rvp) = (b->x[1] & ~0x100000) | ((e + 0x3ff + 52) << 20);
2124 word1(rvp) = b->x[0];
2128 k = b->x[0] & ((1 << j) - 1);
2148 word0(rvp) = b->x[1] | ((e + 65 + 13) << 24);
2149 word1(rvp) = b->x[0];
2155 word0(rvp) = ((b->x[1] & ~0x800000) >> 16) | ((e + 129 + 55) << 7) | (b->x[1] << 16);
2156 word1(rvp) = (b->x[0] >> 16) | (b->x[0] << 16);
2169 int rv = hi0bits(b->x[b->wds-1]) - 4;
2201 sx = S->x;
2203 bx = b->x;
2248 bx = b->x;
2258 bx = b->x;
2259 sx = S->x;
2288 bx = b->x;
2303 (x, bc) U *x; BCinfo *bc; in sulp() argument
2305 (U *x, BCinfo *bc)
2312 rv = ulp(x);
2313 if (!bc->scale || (i = 2*P + 1 - ((word0(x) & Exp_mask) >> Exp_shift)) <= 0)
2397 b->x[0] |= 1;
2438 if (!b->x[0] && b->wds == 1) {
2449 if (!b->x[0] && b->wds == 1) {
2457 if (dig > 0 || b->x[0] || b->wds > 1)
3143 if (!delta->x[0] && delta->wds <= 1) {
3247 if (!delta->x[0] && delta->wds <= 1)
3252 if (!delta->x[0] && delta->wds <= 1) {
3759 ULong x; local
3889 x = i > 32 ? word0(&u) << (64 - i) | word1(&u) >> (i - 32)
3891 dval(&d2) = x;
4279 else if (!b->x[0] && b->wds <= 1)
4291 if (!b->x[0] && b->wds <= 1) {
4350 if (!b->x[0] && b->wds <= 1) {
