![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to asm0.h CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / pari-2.2 / src / kernel / alpha|
Return to asm0.h CVS log | Up to [local] / OpenXM_contrib / pari-2.2 / src / kernel / alpha |
1.1 noro 1: ulong overflow;
2: ulong hiremainder;
3:
4: #define LOCAL_OVERFLOW
5: #define LOCAL_HIREMAINDER
6: #define SAVE_OVERFLOW
7: #define SAVE_HIREMAINDER
8:
9:
10: /* From the PARI source, using gcc __asm__ format. */
11:
12: #define addll(a, b)\
13: ({ register ulong __value, __arg1 = (a), __arg2 = (b); \
14: __asm__ volatile ("addq %2,%3,%0\n\tcmpult %4,%2,%1" \
15: : "=r" (__value), "=r" (overflow) \
16: : "r" (__arg1), "r" (__arg2), "0" ((ulong) 0)); \
17: __value; \
18: })
19:
20: #define addllx(a, b)\
21: ({ register ulong __value, __arg1 = (a), __arg2 = (b), __temp; \
22: __asm__ volatile ("addq %3,%4,%0\n\tcmpult %5,%3,%2\n\taddq %5,%6,%0\n\tcmpult %5,%6,%1\n\taddq %6,%7,%1\n\t" \
23: : "=r" (__value), "=r" (overflow), "=r" (__temp) \
24: : "r" (__arg1), "r" (__arg2), "0" ((ulong) 0), "1" (overflow), "2" ((ulong) 0)); \
25: __value; \
26: })
27:
28: #define subll(a, b)\
29: ({ register ulong __value, __arg1 = (a), __arg2 = (b); \
30: __asm__ volatile ("subq %2,%3,%0\n\tcmpult %2,%4,%1" \
31: : "=r" (__value), "=r" (overflow) \
32: : "r" (__arg1), "r" (__arg2), "0" ((ulong)0)); \
33: __value; \
34: })
35:
36: #define subllx(a, b)\
37: ({ register ulong __value, __arg1 = (a), __arg2 = (b), __temp1, __temp2; \
38: __asm__ volatile ("subq %4,%5,%2\n\tcmpult %4,%8,%3\n\tsubq %8,%7,%0\n\tcmpult %8,%6,%1\n\taddq %7,%9,%1\n\t" \
39: : "=r" (__value), "=r" (overflow), "=r" (__temp1), "=r" (__temp2) \
40: : "r" (__arg1), "r" (__arg2), "0" ((ulong)0), "1" (overflow), "2" ((ulong)0), "3" ((ulong)0)); \
41: __value; \
42: })
43:
44: #define shiftl(a, b) \
45: ({ register ulong __value, __arg1 = (a), __arg2 = (b), __temp; \
46: __asm__ volatile ("subq %5,%4,%2\n\tsll %3,%4,%0\n\tsrl %3,%6,%1\n\t" \
47: : "=r" (__value), "=r" (hiremainder), "=r" (__temp) \
48: : "r" (__arg1), "r" (__arg2), "n" ((ulong) 64), "2" ((ulong)0)); \
49: __value; \
50: })
51:
52: #define shiftlr(a, b) \
53: ({ register ulong __value, __arg1 = (a), __arg2 = (b), __temp; \
54: __asm__ volatile ("subq %5,%4,%2\n\tsrl %3,%4,%0\n\tsll %3,%6,%1\n\t" \
55: : "=r" (__value), "=r" (hiremainder), "=r" (__temp) \
56: : "r" (__arg1), "r" (__arg2), "n" ((ulong) 64), "2" ((ulong)0)); \
57: __value; \
58: })
59:
60: #define mulll(a, b) \
61: ({ register ulong __value, __arg1 = (a), __arg2 = (b); \
62: __asm__ volatile ("umulh %2,%3,%1\n\tmulq %2,%3,%0\n\t" \
63: : "=r" (__value), "=r" (hiremainder) \
64: : "r" (__arg1), "r" (__arg2)); \
65: __value; \
66: })
67:
68: #define addmul(a, b) \
69: ({ register ulong __value, __arg1 = (a), __arg2 = (b), __temp; \
70: __asm__ volatile ("mulq %3,%4,%0\n\tumulh %3,%4,%2\n\taddq %5,%6,%0\n\tcmpult %5,%6,%1\n\taddq %7,%6,%1\n\t" \
71: : "=r" (__value), "=r" (hiremainder), "=r" (__temp) \
72: : "r" (__arg1), "r" (__arg2), "0" ((ulong) 0), "1" (hiremainder), "2" ((ulong) 0)); \
73: __value; \
74: })
75:
76: /*
77: The end of the present file is a slight adaptation of source code
78: extracted from gmp-3.1.1 (from T. Granlund), files longlong.h and
79: gmp-impl.h
80:
81: Copyright (C) 1991, 1992, 1993, 1994, 1996, 1997, 1999, 2000 Free Software
82: Foundation, Inc.
83: */
84:
85: extern const unsigned char __clz_tab[];
86: extern ulong invert_word(ulong);
87:
88: #define bfffo(x) \
89: ({ \
90: ulong __xr = (x); \
91: ulong __a; \
92: \
93: for (__a = 56; __a > 0; __a -= 8) \
94: if (((__xr >> __a) & 0xff) != 0) \
95: break; \
96: 64 - (__clz_tab[__xr >> __a] + __a); \
97: })
98:
99: #define sub_ddmmss(sh, sl, ah, al, bh, bl) \
100: do { \
101: ulong __x; \
102: __x = (al) - (bl); \
103: (sh) = (ah) - (bh) - (__x > (al)); \
104: (sl) = __x; \
105: } while (0)
106:
107: #define divll(x, y) \
108: ({ \
109: register ulong _di, _x = (x), _y = (y), _q, _ql, _r; \
110: register ulong _xh, _xl, _k, __hire; \
111: \
112: if (_y & 0x8000000000000000UL) \
113: { _k = 0; __hire = hiremainder; } \
114: else \
115: { \
116: _k = bfffo(_y); \
117: __hire = (hiremainder << _k) | (_x >> (64 - _k)); \
118: _x <<= _k; _y <<= _k; \
119: } \
120: _di = invert_word(_y); \
121: _ql = mulll (__hire, _di); \
122: _q = __hire + hiremainder; \
123: _xl = mulll(_q, _y); _xh = hiremainder; \
124: sub_ddmmss (_xh, _r, __hire, _x, _xh, _xl); \
125: if (_xh != 0) \
126: { \
127: sub_ddmmss (_xh, _r, _xh, _r, 0, _y); _q += 1; \
128: if (_xh != 0) \
129: { sub_ddmmss (_xh, _r, _xh, _r, 0, _y); _q += 1; } \
130: } \
131: if (_r >= _y) \
132: { _r -= _y; _q += 1; } \
133: hiremainder = _r >> _k; \
134: _q; \
135: })
136:
137: