12 template <
unsigned int BITS>
15 static_assert(BITS/32 > 0 && BITS%32 == 0,
"Template parameter BITS must be a positive multiple of 32.");
20 template <
unsigned int BITS>
24 for (
int i = 0; i < WIDTH; i++)
28 for (
int i = 0; i < WIDTH; i++) {
29 if (i + k + 1 < WIDTH && shift != 0)
30 pn[i + k + 1] |= (a.
pn[i] >> (32 - shift));
32 pn[i + k] |= (a.
pn[i] << shift);
37 template <
unsigned int BITS>
41 for (
int i = 0; i < WIDTH; i++)
45 for (
int i = 0; i < WIDTH; i++) {
46 if (i - k - 1 >= 0 && shift != 0)
47 pn[i - k - 1] |= (a.
pn[i] << (32 - shift));
49 pn[i - k] |= (a.
pn[i] >> shift);
54 template <
unsigned int BITS>
58 for (
int i = 0; i < WIDTH; i++) {
59 uint64_t n = carry + (uint64_t)b32 * pn[i];
60 pn[i] = n & 0xffffffff;
66 template <
unsigned int BITS>
70 for (
int j = 0; j < WIDTH; j++) {
72 for (
int i = 0; i + j < WIDTH; i++) {
73 uint64_t n = carry + a.
pn[i + j] + (uint64_t)pn[j] * b.
pn[i];
74 a.
pn[i + j] = n & 0xffffffff;
82 template <
unsigned int BITS>
88 int num_bits = num.
bits();
89 int div_bits = div.
bits();
92 if (div_bits > num_bits)
94 int shift = num_bits - div_bits;
99 pn[shift / 32] |= (1 << (shift & 31));
108 template <
unsigned int BITS>
111 for (
int i = WIDTH - 1; i >= 0; i--) {
120 template <
unsigned int BITS>
123 for (
int i = WIDTH - 1; i >= 2; i--) {
127 if (pn[1] != (b >> 32))
129 if (pn[0] != (b & 0xfffffffful))
134 template <
unsigned int BITS>
139 for (
int i = 0; i < WIDTH; i++) {
141 fact *= 4294967296.0;
146 template <
unsigned int BITS>
152 template <
unsigned int BITS>
158 template <
unsigned int BITS>
164 template <
unsigned int BITS>
170 template <
unsigned int BITS>
173 for (
int pos = WIDTH - 1; pos >= 0; pos--) {
175 for (
int nbits = 31; nbits > 0; nbits--) {
176 if (pn[pos] & 1U << nbits)
177 return 32 * pos + nbits + 1;
205 int nSize = nCompact >> 24;
206 uint32_t nWord = nCompact & 0x007fffff;
208 nWord >>= 8 * (3 - nSize);
212 *
this <<= 8 * (nSize - 3);
215 *pfNegative = nWord != 0 && (nCompact & 0x00800000) != 0;
217 *pfOverflow = nWord != 0 && ((nSize > 34) ||
218 (nWord > 0xff && nSize > 33) ||
219 (nWord > 0xffff && nSize > 32));
225 int nSize = (
bits() + 7) / 8;
226 uint32_t nCompact = 0;
228 nCompact =
GetLow64() << 8 * (3 - nSize);
235 if (nCompact & 0x00800000) {
239 assert((nCompact & ~0x007fffff) == 0);
241 nCompact |= nSize << 24;
242 nCompact |= (fNegative && (nCompact & 0x007fffff) ? 0x00800000 : 0);
249 for(
int x=0; x<a.
WIDTH; ++x)
256 for(
int x=0; x<b.
WIDTH; ++x)