1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
|
/* See LICENSE file for copyright and license details. */
#include "internals"
#define ta libzahl_tmp_divmod_a
#define tb libzahl_tmp_divmod_b
#define td libzahl_tmp_divmod_d
void
zdivmod(z_t a, z_t b, z_t c, z_t d)
{
size_t c_bits, d_bits, bit;
int sign, cmpmag;
sign = zsignum(c) * zsignum(d);
if (!sign) {
if (zzero(c)) {
if (zzero(d)) {
errno = EDOM; /* Indeterminate form: 0 divided by 0 */
FAILURE_JUMP();
} else {
SET_SIGNUM(a, 0);
SET_SIGNUM(b, 0);
}
} else {
errno = EDOM; /* Undefined form: Division by 0 */
FAILURE_JUMP();
}
return;
} else if ((cmpmag = zcmpmag(c, d)) <= 0) {
if (cmpmag == 0) {
zseti(a, sign);
SET_SIGNUM(b, 0);
return;
} else if (sign < 0) {
zsub_unsigned(b, d, c);
} else if (b != c) {
zset(b, c);
}
SET_SIGNUM(b, 1);
SET_SIGNUM(a, 0);
return;
}
c_bits = zbits(c);
d_bits = zbits(d);
bit = c_bits - d_bits;
zlsh(td, d, bit);
SET_SIGNUM(td, 1);
if (zcmpmag(td, c) > 0) {
zrsh(td, td, 1);
bit -= 1;
}
SET_SIGNUM(ta, 0);
zabs(tb, c);
for (;;) {
if (zcmpmag(td, tb) <= 0) {
zsub(tb, tb, td);
zbset(ta, ta, bit, 1);
}
if (!bit--)
break;
zrsh(td, td, 1);
}
zset(a, ta);
zset(b, tb);
SET_SIGNUM(a, sign);
}
|
