BLI_math_rotation: properly name the quaternion power function.
[blender.git] / source / blender / blenlib / intern / hash_mm2a.c
1 /*
2  * ***** BEGIN GPL LICENSE BLOCK *****
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  *
18  * ***** END GPL LICENSE BLOCK *****
19  *
20  * Copyright (C) 2014 Blender Foundation.
21  *
22  */
23
24 /** \file blender/blenlib/intern/hash_mm2a.c
25  *  \ingroup bli
26  *
27  *  Functions to compute Murmur2A hash key.
28  *
29  * A very fast hash generating int32 result, with few collisions and good repartition.
30  *
31  * See also:
32  *     reference implementation: https://smhasher.googlecode.com/svn-history/r130/trunk/MurmurHash2.cpp
33  *     and http://programmers.stackexchange.com/questions/49550/which-hashing-algorithm-is-best-for-uniqueness-and-speed
34  *
35  * \warning Do not store that hash in files or such, it is not endian-agnostic, so you should only use it
36  *          for temporary data.
37  */
38
39 #include "BLI_compiler_attrs.h"
40
41 #include "BLI_hash_mm2a.h"  /* own include */
42
43 /* Helpers. */
44 #define MM2A_M 0x5bd1e995
45
46 #define MM2A_MIX(h, k)           \
47 {                                \
48         (k) *= MM2A_M;               \
49         (k) ^= (k) >> 24;            \
50         (k) *= MM2A_M;               \
51         (h) = ((h) * MM2A_M) ^ (k);  \
52 } (void)0
53
54 #define MM2A_MIX_FINALIZE(h)     \
55 {                                \
56         (h) ^= (h) >> 13;            \
57         (h) *= MM2A_M;               \
58         (h) ^= (h) >> 15;            \
59 } (void)0
60
61 static void mm2a_mix_tail(BLI_HashMurmur2A *mm2, const unsigned char **data, size_t *len)
62 {
63         while (*len && ((*len < 4) || mm2->count)) {
64                 mm2->tail |= (uint32_t)(**data) << (mm2->count * 8);
65
66                 mm2->count++;
67                 (*len)--;
68                 (*data)++;
69
70                 if (mm2->count == 4) {
71                         MM2A_MIX(mm2->hash, mm2->tail);
72                         mm2->tail = 0;
73                         mm2->count = 0;
74                 }
75         }
76 }
77
78 void BLI_hash_mm2a_init(BLI_HashMurmur2A *mm2, uint32_t seed)
79 {
80         mm2->hash  = seed;
81         mm2->tail  = 0;
82         mm2->count = 0;
83         mm2->size  = 0;
84 }
85
86 void BLI_hash_mm2a_add(BLI_HashMurmur2A *mm2, const unsigned char *data, size_t len)
87 {
88         mm2->size += (uint32_t)len;
89
90         mm2a_mix_tail(mm2, &data, &len);
91
92         for (; len >= 4; data += 4, len -= 4) {
93                 uint32_t k = *(const uint32_t *)data;
94
95                 MM2A_MIX(mm2->hash, k);
96         }
97
98         mm2a_mix_tail(mm2, &data, &len);
99 }
100
101 void BLI_hash_mm2a_add_int(BLI_HashMurmur2A *mm2, int data)
102 {
103         BLI_hash_mm2a_add(mm2, (const unsigned char *)&data, sizeof(data));
104 }
105
106 uint32_t BLI_hash_mm2a_end(BLI_HashMurmur2A *mm2)
107 {
108         MM2A_MIX(mm2->hash, mm2->tail);
109         MM2A_MIX(mm2->hash, mm2->size);
110
111         MM2A_MIX_FINALIZE(mm2->hash);
112
113         return mm2->hash;
114 }
115
116 /* Non-incremental version, quicker for small keys. */
117 uint32_t BLI_hash_mm2(const unsigned char *data, size_t len, uint32_t seed)
118 {
119         /* Initialize the hash to a 'random' value */
120         uint32_t h = seed ^ len;
121
122         /* Mix 4 bytes at a time into the hash */
123         for (; len >= 4; data += 4, len -= 4) {
124                 uint32_t k = *(uint32_t *)data;
125
126                 MM2A_MIX(h, k);
127         }
128
129         /* Handle the last few bytes of the input array */
130         switch (len) {
131                 case 3:
132                         h ^= data[2] << 16;
133                         ATTR_FALLTHROUGH;
134                 case 2:
135                         h ^= data[1] << 8;
136                         ATTR_FALLTHROUGH;
137                 case 1:
138                         h ^= data[0];
139                         h *= MM2A_M;
140         }
141
142         /* Do a few final mixes of the hash to ensure the last few bytes are well-incorporated. */
143         MM2A_MIX_FINALIZE(h);
144
145         return h;
146 }