2 * Routines to provide a memory-efficient hashtable.
4 * Copyright (C) 2007-2019 Wayne Davison
6 * This program is free software; you can redistribute it and/or modify
7 * it under the terms of the GNU General Public License as published by
8 * the Free Software Foundation; either version 3 of the License, or
9 * (at your option) any later version.
11 * This program is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 * GNU General Public License for more details.
16 * You should have received a copy of the GNU General Public License along
17 * with this program; if not, visit the http://fsf.org website.
22 #define HASH_LOAD_LIMIT(size) ((size)*3/4)
24 struct hashtable
*hashtable_create(int size
, int key64
)
27 struct hashtable
*tbl
;
28 int node_size
= key64
? sizeof (struct ht_int64_node
)
29 : sizeof (struct ht_int32_node
);
31 /* Pick a power of 2 that can hold the requested size. */
32 if (size
& (size
-1) || size
< 16) {
38 if (!(tbl
= new(struct hashtable
))
39 || !(tbl
->nodes
= new_array0(char, size
* node_size
)))
40 out_of_memory("hashtable_create");
43 tbl
->node_size
= node_size
;
44 tbl
->key64
= key64
? 1 : 0;
46 if (DEBUG_GTE(HASH
, 1)) {
49 snprintf(buf
, sizeof buf
, "req: %d, ", req
);
52 rprintf(FINFO
, "[%s] created hashtable %lx (%ssize: %d, keys: %d-bit)\n",
53 who_am_i(), (long)tbl
, buf
, size
, key64
? 64 : 32);
59 void hashtable_destroy(struct hashtable
*tbl
)
61 if (DEBUG_GTE(HASH
, 1)) {
62 rprintf(FINFO
, "[%s] destroyed hashtable %lx (size: %d, keys: %d-bit)\n",
63 who_am_i(), (long)tbl
, tbl
->size
, tbl
->key64
? 64 : 32);
69 /* This returns the node for the indicated key, either newly created or
70 * already existing. Returns NULL if not allocating and not found. */
71 void *hashtable_find(struct hashtable
*tbl
, int64 key
, int allocate_if_missing
)
73 int key64
= tbl
->key64
;
74 struct ht_int32_node
*node
;
77 if (key64
? key
== 0 : (int32
)key
== 0) {
78 rprintf(FERROR
, "Internal hashtable error: illegal key supplied!\n");
79 exit_cleanup(RERR_MESSAGEIO
);
82 if (allocate_if_missing
&& tbl
->entries
> HASH_LOAD_LIMIT(tbl
->size
)) {
83 void *old_nodes
= tbl
->nodes
;
84 int size
= tbl
->size
* 2;
87 if (!(tbl
->nodes
= new_array0(char, size
* tbl
->node_size
)))
88 out_of_memory("hashtable_node");
92 if (DEBUG_GTE(HASH
, 1)) {
93 rprintf(FINFO
, "[%s] growing hashtable %lx (size: %d, keys: %d-bit)\n",
94 who_am_i(), (long)tbl
, size
, key64
? 64 : 32);
97 for (i
= size
/ 2; i
-- > 0; ) {
98 struct ht_int32_node
*move_node
= HT_NODE(tbl
, old_nodes
, i
);
99 int64 move_key
= HT_KEY(move_node
, key64
);
102 node
= hashtable_find(tbl
, move_key
, 1);
103 node
->data
= move_node
->data
;
110 /* Based on Jenkins One-at-a-time hash. */
111 uchar buf
[4], *keyp
= buf
;
115 for (ndx
= 0, i
= 0; i
< 4; i
++) {
124 /* Based on Jenkins hashword() from lookup3.c. */
127 /* Set up the internal state */
128 a
= b
= c
= 0xdeadbeef + (8 << 2);
130 #define rot(x,k) (((x)<<(k)) ^ ((x)>>(32-(k))))
131 #if SIZEOF_INT64 >= 8
132 b
+= (uint32
)(key
>> 32);
135 c
^= b
; c
-= rot(b
, 14);
136 a
^= c
; a
-= rot(c
, 11);
137 b
^= a
; b
-= rot(a
, 25);
138 c
^= b
; c
-= rot(b
, 16);
139 a
^= c
; a
-= rot(c
, 4);
140 b
^= a
; b
-= rot(a
, 14);
141 c
^= b
; c
-= rot(b
, 24);
146 /* If it already exists, return the node. If we're not
147 * allocating, return NULL if the key is not found. */
151 ndx
&= tbl
->size
- 1;
152 node
= HT_NODE(tbl
, tbl
->nodes
, ndx
);
153 nkey
= HT_KEY(node
, key64
);
158 if (!allocate_if_missing
)
165 /* Take over this empty spot and then return the node. */
167 ((struct ht_int64_node
*)node
)->key
= key
;
169 node
->key
= (int32
)key
;
174 #ifndef WORDS_BIGENDIAN
175 # define HASH_LITTLE_ENDIAN 1
176 # define HASH_BIG_ENDIAN 0
178 # define HASH_LITTLE_ENDIAN 0
179 # define HASH_BIG_ENDIAN 1
183 -------------------------------------------------------------------------------
184 lookup3.c, by Bob Jenkins, May 2006, Public Domain.
186 These are functions for producing 32-bit hashes for hash table lookup.
187 hash_word(), hashlittle(), hashlittle2(), hashbig(), mix(), and final()
188 are externally useful functions. Routines to test the hash are included
189 if SELF_TEST is defined. You can use this free for any purpose. It's in
190 the public domain. It has no warranty.
192 You probably want to use hashlittle(). hashlittle() and hashbig()
193 hash byte arrays. hashlittle() is is faster than hashbig() on
194 little-endian machines. Intel and AMD are little-endian machines.
195 On second thought, you probably want hashlittle2(), which is identical to
196 hashlittle() except it returns two 32-bit hashes for the price of one.
197 You could implement hashbig2() if you wanted but I haven't bothered here.
199 If you want to find a hash of, say, exactly 7 integers, do
200 a = i1; b = i2; c = i3;
202 a += i4; b += i5; c += i6;
206 then use c as the hash value. If you have a variable length array of
207 4-byte integers to hash, use hash_word(). If you have a byte array (like
208 a character string), use hashlittle(). If you have several byte arrays, or
209 a mix of things, see the comments above hashlittle().
211 Why is this so big? I read 12 bytes at a time into 3 4-byte integers,
212 then mix those integers. This is fast (you can do a lot more thorough
213 mixing with 12*3 instructions on 3 integers than you can with 3 instructions
214 on 1 byte), but shoehorning those bytes into integers efficiently is messy.
217 #define hashsize(n) ((uint32_t)1<<(n))
218 #define hashmask(n) (hashsize(n)-1)
219 #define rot(x,k) (((x)<<(k)) | ((x)>>(32-(k))))
222 -------------------------------------------------------------------------------
223 mix -- mix 3 32-bit values reversibly.
225 This is reversible, so any information in (a,b,c) before mix() is
226 still in (a,b,c) after mix().
228 If four pairs of (a,b,c) inputs are run through mix(), or through
229 mix() in reverse, there are at least 32 bits of the output that
230 are sometimes the same for one pair and different for another pair.
232 * pairs that differed by one bit, by two bits, in any combination
233 of top bits of (a,b,c), or in any combination of bottom bits of
235 * "differ" is defined as +, -, ^, or ~^. For + and -, I transformed
236 the output delta to a Gray code (a^(a>>1)) so a string of 1's (as
237 is commonly produced by subtraction) look like a single 1-bit
239 * the base values were pseudorandom, all zero but one bit set, or
240 all zero plus a counter that starts at zero.
242 Some k values for my "a-=c; a^=rot(c,k); c+=b;" arrangement that
247 Well, "9 15 3 18 27 15" didn't quite get 32 bits diffing
248 for "differ" defined as + with a one-bit base and a two-bit delta. I
249 used http://burtleburtle.net/bob/hash/avalanche.html to choose
250 the operations, constants, and arrangements of the variables.
252 This does not achieve avalanche. There are input bits of (a,b,c)
253 that fail to affect some output bits of (a,b,c), especially of a. The
254 most thoroughly mixed value is c, but it doesn't really even achieve
257 This allows some parallelism. Read-after-writes are good at doubling
258 the number of bits affected, so the goal of mixing pulls in the opposite
259 direction as the goal of parallelism. I did what I could. Rotates
260 seem to cost as much as shifts on every machine I could lay my hands
261 on, and rotates are much kinder to the top and bottom bits, so I used
263 -------------------------------------------------------------------------------
267 a -= c; a ^= rot(c, 4); c += b; \
268 b -= a; b ^= rot(a, 6); a += c; \
269 c -= b; c ^= rot(b, 8); b += a; \
270 a -= c; a ^= rot(c,16); c += b; \
271 b -= a; b ^= rot(a,19); a += c; \
272 c -= b; c ^= rot(b, 4); b += a; \
276 -------------------------------------------------------------------------------
277 final -- final mixing of 3 32-bit values (a,b,c) into c
279 Pairs of (a,b,c) values differing in only a few bits will usually
280 produce values of c that look totally different. This was tested for
281 * pairs that differed by one bit, by two bits, in any combination
282 of top bits of (a,b,c), or in any combination of bottom bits of
284 * "differ" is defined as +, -, ^, or ~^. For + and -, I transformed
285 the output delta to a Gray code (a^(a>>1)) so a string of 1's (as
286 is commonly produced by subtraction) look like a single 1-bit
288 * the base values were pseudorandom, all zero but one bit set, or
289 all zero plus a counter that starts at zero.
291 These constants passed:
294 and these came close:
298 -------------------------------------------------------------------------------
300 #define final(a,b,c) \
302 c ^= b; c -= rot(b,14); \
303 a ^= c; a -= rot(c,11); \
304 b ^= a; b -= rot(a,25); \
305 c ^= b; c -= rot(b,16); \
306 a ^= c; a -= rot(c,4); \
307 b ^= a; b -= rot(a,14); \
308 c ^= b; c -= rot(b,24); \
313 -------------------------------------------------------------------------------
314 hashlittle() -- hash a variable-length key into a 32-bit value
315 k : the key (the unaligned variable-length array of bytes)
316 length : the length of the key, counting by bytes
317 val2 : IN: can be any 4-byte value OUT: second 32 bit hash.
318 Returns a 32-bit value. Every bit of the key affects every bit of
319 the return value. Two keys differing by one or two bits will have
320 totally different hash values. Note that the return value is better
321 mixed than val2, so use that first.
323 The best hash table sizes are powers of 2. There is no need to do
324 mod a prime (mod is sooo slow!). If you need less than 32 bits,
325 use a bitmask. For example, if you need only 10 bits, do
326 h = (h & hashmask(10));
327 In which case, the hash table should have hashsize(10) elements.
329 If you are hashing n strings (uint8_t **)k, do it like this:
330 for (i=0, h=0; i<n; ++i) h = hashlittle( k[i], len[i], h);
332 By Bob Jenkins, 2006. bob_jenkins@burtleburtle.net. You may use this
333 code any way you wish, private, educational, or commercial. It's free.
335 Use for hash table lookup, or anything where one collision in 2^^32 is
336 acceptable. Do NOT use for cryptographic purposes.
337 -------------------------------------------------------------------------------
340 uint32_t hashlittle(const void *key
, size_t length
)
342 uint32_t a
,b
,c
; /* internal state */
343 union { const void *ptr
; size_t i
; } u
; /* needed for Mac Powerbook G4 */
345 /* Set up the internal state */
346 a
= b
= c
= 0xdeadbeef + ((uint32_t)length
);
349 if (HASH_LITTLE_ENDIAN
&& ((u
.i
& 0x3) == 0)) {
350 const uint32_t *k
= (const uint32_t *)key
; /* read 32-bit chunks */
353 /*------ all but last block: aligned reads and affect 32 bits of (a,b,c) */
364 /*----------------------------- handle the last (probably partial) block */
365 k8
= (const uint8_t *)k
;
368 case 12: c
+=k
[2]; b
+=k
[1]; a
+=k
[0]; break;
369 case 11: c
+=((uint32_t)k8
[10])<<16; /* fall through */
370 case 10: c
+=((uint32_t)k8
[9])<<8; /* fall through */
371 case 9 : c
+=k8
[8]; /* fall through */
372 case 8 : b
+=k
[1]; a
+=k
[0]; break;
373 case 7 : b
+=((uint32_t)k8
[6])<<16; /* fall through */
374 case 6 : b
+=((uint32_t)k8
[5])<<8; /* fall through */
375 case 5 : b
+=k8
[4]; /* fall through */
376 case 4 : a
+=k
[0]; break;
377 case 3 : a
+=((uint32_t)k8
[2])<<16; /* fall through */
378 case 2 : a
+=((uint32_t)k8
[1])<<8; /* fall through */
379 case 1 : a
+=k8
[0]; break;
382 } else if (HASH_LITTLE_ENDIAN
&& ((u
.i
& 0x1) == 0)) {
383 const uint16_t *k
= (const uint16_t *)key
; /* read 16-bit chunks */
386 /*--------------- all but last block: aligned reads and different mixing */
389 a
+= k
[0] + (((uint32_t)k
[1])<<16);
390 b
+= k
[2] + (((uint32_t)k
[3])<<16);
391 c
+= k
[4] + (((uint32_t)k
[5])<<16);
397 /*----------------------------- handle the last (probably partial) block */
398 k8
= (const uint8_t *)k
;
401 case 12: c
+=k
[4]+(((uint32_t)k
[5])<<16);
402 b
+=k
[2]+(((uint32_t)k
[3])<<16);
403 a
+=k
[0]+(((uint32_t)k
[1])<<16);
405 case 11: c
+=((uint32_t)k8
[10])<<16; /* fall through */
407 b
+=k
[2]+(((uint32_t)k
[3])<<16);
408 a
+=k
[0]+(((uint32_t)k
[1])<<16);
410 case 9 : c
+=k8
[8]; /* fall through */
411 case 8 : b
+=k
[2]+(((uint32_t)k
[3])<<16);
412 a
+=k
[0]+(((uint32_t)k
[1])<<16);
414 case 7 : b
+=((uint32_t)k8
[6])<<16; /* fall through */
416 a
+=k
[0]+(((uint32_t)k
[1])<<16);
418 case 5 : b
+=k8
[4]; /* fall through */
419 case 4 : a
+=k
[0]+(((uint32_t)k
[1])<<16);
421 case 3 : a
+=((uint32_t)k8
[2])<<16; /* fall through */
426 case 0 : return c
; /* zero length requires no mixing */
429 } else { /* need to read the key one byte at a time */
430 const uint8_t *k
= (const uint8_t *)key
;
432 /*--------------- all but the last block: affect some 32 bits of (a,b,c) */
436 a
+= ((uint32_t)k
[1])<<8;
437 a
+= ((uint32_t)k
[2])<<16;
438 a
+= ((uint32_t)k
[3])<<24;
440 b
+= ((uint32_t)k
[5])<<8;
441 b
+= ((uint32_t)k
[6])<<16;
442 b
+= ((uint32_t)k
[7])<<24;
444 c
+= ((uint32_t)k
[9])<<8;
445 c
+= ((uint32_t)k
[10])<<16;
446 c
+= ((uint32_t)k
[11])<<24;
452 /*-------------------------------- last block: affect all 32 bits of (c) */
453 switch(length
) /* all the case statements fall through */
455 case 12: c
+=((uint32_t)k
[11])<<24;
457 case 11: c
+=((uint32_t)k
[10])<<16;
459 case 10: c
+=((uint32_t)k
[9])<<8;
463 case 8 : b
+=((uint32_t)k
[7])<<24;
465 case 7 : b
+=((uint32_t)k
[6])<<16;
467 case 6 : b
+=((uint32_t)k
[5])<<8;
471 case 4 : a
+=((uint32_t)k
[3])<<24;
473 case 3 : a
+=((uint32_t)k
[2])<<16;
475 case 2 : a
+=((uint32_t)k
[1])<<8;