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
|
# File 'ext/sumbur/sumbur.c', line 4
static VALUE
rb_sumbur(VALUE self, VALUE hashed_int, VALUE capacity)
{
unsigned int h = NUM2UINT(hashed_int);
unsigned int capa = NUM2UINT(capacity);
unsigned int part, n, i, c;
if (capacity == 0) {
rb_raise(rb_eArgError, "Sumbur is not applicable to empty cluster");
}
part = L / capa;
if (L - h <= part) return INT2FIX(0);
n = 1;
do {
if (h >= L / 2) h -= L / 2;
else {
n = 2;
if (L / 2 - h < part) break;
}
if (capa == 2) break;
#define curslice(i) (L / (i * (i - 1)))
#define unroll(i) \
if (curslice(i) <= h) h -= curslice(i); \
else { \
h += curslice(i) * (i - n - 1); \
n = i; \
if (L / i - h < part) break; \
} \
if (capa == i) break
unroll(3); unroll(4); unroll(5); unroll(6); unroll(7);
unroll(8); unroll(9); unroll(10); unroll(11); unroll(12);
unroll(13); unroll(14); unroll(15); unroll(16); unroll(17);
unroll(18); unroll(19); unroll(20); unroll(21); unroll(22);
unroll(23); unroll(24);
for(i = 25; i <= capa; i++) {
c = L / (i * (i - 1));
if (c <= h) {
h -= c;
}
else {
h += c * (i - n - 1);
n = i;
if (L / i - h < part) break;
}
}
} while(0);
return INT2FIX(n - 1);
}
|