| 1 | /* |
1 | /* |
| 2 | * This file is part of Crossfire TRT, the Multiplayer Online Role Playing Game. |
2 | * This file is part of Deliantra, the Roguelike Realtime MMORPG. |
| 3 | * |
3 | * |
| 4 | * Copyright (©) 2005,2006,2007 Marc Alexander Lehmann / Robin Redeker / the Crossfire TRT team |
4 | * Copyright (©) 2005,2006,2007,2008 Marc Alexander Lehmann / Robin Redeker / the Deliantra team |
| 5 | * |
5 | * |
| 6 | * Crossfire TRT is free software; you can redistribute it and/or modify it |
6 | * Deliantra is free software: you can redistribute it and/or modify it under |
| 7 | * under the terms of the GNU General Public License as published by the Free |
7 | * the terms of the Affero GNU General Public License as published by the |
| 8 | * Software Foundation; either version 2 of the License, or (at your option) |
8 | * Free Software Foundation, either version 3 of the License, or (at your |
| 9 | * any later version. |
9 | * option) any later version. |
| 10 | * |
10 | * |
| 11 | * This program is distributed in the hope that it will be useful, but |
11 | * This program is distributed in the hope that it will be useful, |
| 12 | * WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY |
12 | * but WITHOUT ANY WARRANTY; without even the implied warranty of |
| 13 | * or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
13 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
| 14 | * for more details. |
14 | * GNU General Public License for more details. |
| 15 | * |
15 | * |
| 16 | * You should have received a copy of the GNU General Public License along |
16 | * You should have received a copy of the Affero GNU General Public License |
| 17 | * with Crossfire TRT; if not, write to the Free Software Foundation, Inc. 51 |
17 | * and the GNU General Public License along with this program. If not, see |
| 18 | * Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA |
18 | * <http://www.gnu.org/licenses/>. |
| 19 | * |
19 | * |
| 20 | * The authors can be reached via e-mail to <crossfire@schmorp.de> |
20 | * The authors can be reached via e-mail to <support@deliantra.net> |
| 21 | */ |
21 | */ |
| 22 | |
22 | |
| 23 | #ifndef UTIL_H__ |
23 | #ifndef UTIL_H__ |
| 24 | #define UTIL_H__ |
24 | #define UTIL_H__ |
| 25 | |
25 | |
| 26 | //#define PREFER_MALLOC |
26 | #include <compiler.h> |
| 27 | |
27 | |
| 28 | #if __GNUC__ >= 3 |
28 | #define DEBUG_POISON 0x00 // poison memory before freeing it if != 0 |
| 29 | # define is_constant(c) __builtin_constant_p (c) |
29 | #define DEBUG_SALLOC 0 // add a debug wrapper around all sallocs |
| 30 | # define expect(expr,value) __builtin_expect ((expr),(value)) |
30 | #define PREFER_MALLOC 0 // use malloc and not the slice allocator |
| 31 | # define prefetch(addr,rw,locality) __builtin_prefetch (addr, rw, locality) |
|
|
| 32 | #else |
|
|
| 33 | # define is_constant(c) 0 |
|
|
| 34 | # define expect(expr,value) (expr) |
|
|
| 35 | # define prefetch(addr,rw,locality) |
|
|
| 36 | #endif |
|
|
| 37 | |
31 | |
| 38 | // put into ifs if you are very sure that the expression |
32 | #include <pthread.h> |
| 39 | // is mostly true or mosty false. note that these return |
|
|
| 40 | // booleans, not the expression. |
|
|
| 41 | #define expect_false(expr) expect ((expr) != 0, 0) |
|
|
| 42 | #define expect_true(expr) expect ((expr) != 0, 1) |
|
|
| 43 | |
33 | |
| 44 | #include <cstddef> |
34 | #include <cstddef> |
| 45 | #include <cmath> |
35 | #include <cmath> |
| 46 | #include <new> |
36 | #include <new> |
| 47 | #include <vector> |
37 | #include <vector> |
| … | |
… | |
| 49 | #include <glib.h> |
39 | #include <glib.h> |
| 50 | |
40 | |
| 51 | #include <shstr.h> |
41 | #include <shstr.h> |
| 52 | #include <traits.h> |
42 | #include <traits.h> |
| 53 | |
43 | |
|
|
44 | #if DEBUG_SALLOC |
|
|
45 | # define g_slice_alloc0(s) debug_slice_alloc0(s) |
|
|
46 | # define g_slice_alloc(s) debug_slice_alloc(s) |
|
|
47 | # define g_slice_free1(s,p) debug_slice_free1(s,p) |
|
|
48 | void *g_slice_alloc (unsigned long size); |
|
|
49 | void *g_slice_alloc0 (unsigned long size); |
|
|
50 | void g_slice_free1 (unsigned long size, void *ptr); |
|
|
51 | #elif PREFER_MALLOC |
|
|
52 | # define g_slice_alloc0(s) calloc (1, (s)) |
|
|
53 | # define g_slice_alloc(s) malloc ((s)) |
|
|
54 | # define g_slice_free1(s,p) free ((p)) |
|
|
55 | #endif |
|
|
56 | |
| 54 | // use a gcc extension for auto declarations until ISO C++ sanctifies them |
57 | // use C0X decltype for auto declarations until ISO C++ sanctifies them (if ever) |
| 55 | #define auto(var,expr) typeof(expr) var = (expr) |
58 | #define auto(var,expr) decltype(expr) var = (expr) |
| 56 | |
59 | |
| 57 | // very ugly macro that basicaly declares and initialises a variable |
60 | // very ugly macro that basically declares and initialises a variable |
| 58 | // that is in scope for the next statement only |
61 | // that is in scope for the next statement only |
| 59 | // works only for stuff that can be assigned 0 and converts to false |
62 | // works only for stuff that can be assigned 0 and converts to false |
| 60 | // (note: works great for pointers) |
63 | // (note: works great for pointers) |
| 61 | // most ugly macro I ever wrote |
64 | // most ugly macro I ever wrote |
| 62 | #define declvar(type, name, value) if (type name = 0) { } else if (((name) = (value)), 1) |
65 | #define statementvar(type, name, value) if (type name = 0) { } else if (((name) = (value)), 1) |
| 63 | |
66 | |
| 64 | // in range including end |
67 | // in range including end |
| 65 | #define IN_RANGE_INC(val,beg,end) \ |
68 | #define IN_RANGE_INC(val,beg,end) \ |
| 66 | ((unsigned int)(val) - (unsigned int)(beg) <= (unsigned int)(end) - (unsigned int)(beg)) |
69 | ((unsigned int)(val) - (unsigned int)(beg) <= (unsigned int)(end) - (unsigned int)(beg)) |
| 67 | |
70 | |
| 68 | // in range excluding end |
71 | // in range excluding end |
| 69 | #define IN_RANGE_EXC(val,beg,end) \ |
72 | #define IN_RANGE_EXC(val,beg,end) \ |
| 70 | ((unsigned int)(val) - (unsigned int)(beg) < (unsigned int)(end) - (unsigned int)(beg)) |
73 | ((unsigned int)(val) - (unsigned int)(beg) < (unsigned int)(end) - (unsigned int)(beg)) |
| 71 | |
74 | |
|
|
75 | void cleanup (const char *cause, bool make_core = false); |
| 72 | void fork_abort (const char *msg); |
76 | void fork_abort (const char *msg); |
| 73 | |
77 | |
| 74 | // rationale for using (U) not (T) is to reduce signed/unsigned issues, |
78 | // rationale for using (U) not (T) is to reduce signed/unsigned issues, |
| 75 | // as a is often a constant while b is the variable. it is still a bug, though. |
79 | // as a is often a constant while b is the variable. it is still a bug, though. |
| 76 | template<typename T, typename U> static inline T min (T a, U b) { return (U)a < b ? (U)a : b; } |
80 | template<typename T, typename U> static inline T min (T a, U b) { return (U)a < b ? (U)a : b; } |
| 77 | template<typename T, typename U> static inline T max (T a, U b) { return (U)a > b ? (U)a : b; } |
81 | template<typename T, typename U> static inline T max (T a, U b) { return (U)a > b ? (U)a : b; } |
| 78 | template<typename T, typename U, typename V> static inline T clamp (T v, U a, V b) { return v < (T)a ? (T)a : v >(T)b ? (T)b : v; } |
82 | template<typename T, typename U, typename V> static inline T clamp (T v, U a, V b) { return v < (T)a ? (T)a : v >(T)b ? (T)b : v; } |
| 79 | |
83 | |
|
|
84 | template<typename T, typename U> static inline void min_it (T &v, U m) { v = min (v, (T)m); } |
|
|
85 | template<typename T, typename U> static inline void max_it (T &v, U m) { v = max (v, (T)m); } |
|
|
86 | template<typename T, typename U, typename V> static inline void clamp_it (T &v, U a, V b) { v = clamp (v, (T)a, (T)b); } |
|
|
87 | |
| 80 | template<typename T, typename U> static inline void swap (T& a, U& b) { T t=a; a=(T)b; b=(U)t; } |
88 | template<typename T, typename U> static inline void swap (T& a, U& b) { T t=a; a=(T)b; b=(U)t; } |
| 81 | |
89 | |
|
|
90 | template<typename T, typename U, typename V> static inline T min (T a, U b, V c) { return min (a, min (b, c)); } |
|
|
91 | template<typename T, typename U, typename V> static inline T max (T a, U b, V c) { return max (a, max (b, c)); } |
|
|
92 | |
|
|
93 | // sign returns -1 or +1 |
|
|
94 | template<typename T> |
|
|
95 | static inline T sign (T v) { return v < 0 ? -1 : +1; } |
|
|
96 | // relies on 2c representation |
|
|
97 | template<> |
|
|
98 | inline sint8 sign (sint8 v) { return 1 - (sint8 (uint8 (v) >> 7) * 2); } |
|
|
99 | |
|
|
100 | // sign0 returns -1, 0 or +1 |
|
|
101 | template<typename T> |
|
|
102 | static inline T sign0 (T v) { return v ? sign (v) : 0; } |
|
|
103 | |
|
|
104 | // div* only work correctly for div > 0 |
|
|
105 | // div, with correct rounding (< 0.5 downwards, >=0.5 upwards) |
|
|
106 | template<typename T> static inline T div (T val, T div) |
|
|
107 | { |
|
|
108 | return expect_false (val < 0) ? - ((-val + (div - 1) / 2) / div) : (val + div / 2) / div; |
|
|
109 | } |
|
|
110 | // div, round-up |
|
|
111 | template<typename T> static inline T div_ru (T val, T div) |
|
|
112 | { |
|
|
113 | return expect_false (val < 0) ? - ((-val ) / div) : (val + div - 1) / div; |
|
|
114 | } |
|
|
115 | // div, round-down |
|
|
116 | template<typename T> static inline T div_rd (T val, T div) |
|
|
117 | { |
|
|
118 | return expect_false (val < 0) ? - ((-val + (div - 1) ) / div) : (val ) / div; |
|
|
119 | } |
|
|
120 | |
|
|
121 | // lerp* only work correctly for min_in < max_in |
|
|
122 | // Linear intERPolate, scales val from min_in..max_in to min_out..max_out |
| 82 | template<typename T> |
123 | template<typename T> |
| 83 | static inline T |
124 | static inline T |
| 84 | lerp (T val, T min_in, T max_in, T min_out, T max_out) |
125 | lerp (T val, T min_in, T max_in, T min_out, T max_out) |
| 85 | { |
126 | { |
| 86 | return (val - min_in) * (max_out - min_out) / (max_in - min_in) + min_out; |
127 | return min_out + div <T> ((val - min_in) * (max_out - min_out), max_in - min_in); |
|
|
128 | } |
|
|
129 | |
|
|
130 | // lerp, round-down |
|
|
131 | template<typename T> |
|
|
132 | static inline T |
|
|
133 | lerp_rd (T val, T min_in, T max_in, T min_out, T max_out) |
|
|
134 | { |
|
|
135 | return min_out + div_rd<T> ((val - min_in) * (max_out - min_out), max_in - min_in); |
|
|
136 | } |
|
|
137 | |
|
|
138 | // lerp, round-up |
|
|
139 | template<typename T> |
|
|
140 | static inline T |
|
|
141 | lerp_ru (T val, T min_in, T max_in, T min_out, T max_out) |
|
|
142 | { |
|
|
143 | return min_out + div_ru<T> ((val - min_in) * (max_out - min_out), max_in - min_in); |
| 87 | } |
144 | } |
| 88 | |
145 | |
| 89 | // lots of stuff taken from FXT |
146 | // lots of stuff taken from FXT |
| 90 | |
147 | |
| 91 | /* Rotate right. This is used in various places for checksumming */ |
148 | /* Rotate right. This is used in various places for checksumming */ |
| … | |
… | |
| 129 | int32_t d = b - a; |
186 | int32_t d = b - a; |
| 130 | d &= d >> 31; |
187 | d &= d >> 31; |
| 131 | return b - d; |
188 | return b - d; |
| 132 | } |
189 | } |
| 133 | |
190 | |
| 134 | // this is much faster than crossfires original algorithm |
191 | // this is much faster than crossfire's original algorithm |
| 135 | // on modern cpus |
192 | // on modern cpus |
| 136 | inline int |
193 | inline int |
| 137 | isqrt (int n) |
194 | isqrt (int n) |
| 138 | { |
195 | { |
| 139 | return (int)sqrtf ((float)n); |
196 | return (int)sqrtf ((float)n); |
|
|
197 | } |
|
|
198 | |
|
|
199 | // this is kind of like the ^^ operator, if it would exist, without sequence point. |
|
|
200 | // more handy than it looks like, due to the implicit !! done on its arguments |
|
|
201 | inline bool |
|
|
202 | logical_xor (bool a, bool b) |
|
|
203 | { |
|
|
204 | return a != b; |
|
|
205 | } |
|
|
206 | |
|
|
207 | inline bool |
|
|
208 | logical_implies (bool a, bool b) |
|
|
209 | { |
|
|
210 | return a <= b; |
| 140 | } |
211 | } |
| 141 | |
212 | |
| 142 | // this is only twice as fast as naive sqrtf (dx*dy+dy*dy) |
213 | // this is only twice as fast as naive sqrtf (dx*dy+dy*dy) |
| 143 | #if 0 |
214 | #if 0 |
| 144 | // and has a max. error of 6 in the range -100..+100. |
215 | // and has a max. error of 6 in the range -100..+100. |
| … | |
… | |
| 169 | absdir (int d) |
240 | absdir (int d) |
| 170 | { |
241 | { |
| 171 | return ((d - 1) & 7) + 1; |
242 | return ((d - 1) & 7) + 1; |
| 172 | } |
243 | } |
| 173 | |
244 | |
|
|
245 | extern ssize_t slice_alloc; // statistics |
|
|
246 | |
|
|
247 | void *salloc_ (int n) throw (std::bad_alloc); |
|
|
248 | void *salloc_ (int n, void *src) throw (std::bad_alloc); |
|
|
249 | |
|
|
250 | // strictly the same as g_slice_alloc, but never returns 0 |
|
|
251 | template<typename T> |
|
|
252 | inline T *salloc (int n = 1) throw (std::bad_alloc) { return (T *)salloc_ (n * sizeof (T)); } |
|
|
253 | |
|
|
254 | // also copies src into the new area, like "memdup" |
|
|
255 | // if src is 0, clears the memory |
|
|
256 | template<typename T> |
|
|
257 | inline T *salloc (int n, T *src) throw (std::bad_alloc) { return (T *)salloc_ (n * sizeof (T), (void *)src); } |
|
|
258 | |
|
|
259 | // clears the memory |
|
|
260 | template<typename T> |
|
|
261 | inline T *salloc0(int n = 1) throw (std::bad_alloc) { return (T *)salloc_ (n * sizeof (T), 0); } |
|
|
262 | |
|
|
263 | // for symmetry |
|
|
264 | template<typename T> |
|
|
265 | inline void sfree (T *ptr, int n = 1) throw () |
|
|
266 | { |
|
|
267 | if (expect_true (ptr)) |
|
|
268 | { |
|
|
269 | slice_alloc -= n * sizeof (T); |
|
|
270 | if (DEBUG_POISON) memset (ptr, DEBUG_POISON, n * sizeof (T)); |
|
|
271 | g_slice_free1 (n * sizeof (T), (void *)ptr); |
|
|
272 | assert (slice_alloc >= 0);//D |
|
|
273 | } |
|
|
274 | } |
|
|
275 | |
|
|
276 | // nulls the pointer |
|
|
277 | template<typename T> |
|
|
278 | inline void sfree0 (T *&ptr, int n = 1) throw () |
|
|
279 | { |
|
|
280 | sfree<T> (ptr, n); |
|
|
281 | ptr = 0; |
|
|
282 | } |
|
|
283 | |
| 174 | // makes dynamically allocated objects zero-initialised |
284 | // makes dynamically allocated objects zero-initialised |
| 175 | struct zero_initialised |
285 | struct zero_initialised |
| 176 | { |
286 | { |
| 177 | void *operator new (size_t s, void *p) |
287 | void *operator new (size_t s, void *p) |
| 178 | { |
288 | { |
| … | |
… | |
| 180 | return p; |
290 | return p; |
| 181 | } |
291 | } |
| 182 | |
292 | |
| 183 | void *operator new (size_t s) |
293 | void *operator new (size_t s) |
| 184 | { |
294 | { |
| 185 | return g_slice_alloc0 (s); |
295 | return salloc0<char> (s); |
| 186 | } |
296 | } |
| 187 | |
297 | |
| 188 | void *operator new[] (size_t s) |
298 | void *operator new[] (size_t s) |
| 189 | { |
299 | { |
| 190 | return g_slice_alloc0 (s); |
300 | return salloc0<char> (s); |
| 191 | } |
301 | } |
| 192 | |
302 | |
| 193 | void operator delete (void *p, size_t s) |
303 | void operator delete (void *p, size_t s) |
| 194 | { |
304 | { |
| 195 | g_slice_free1 (s, p); |
305 | sfree ((char *)p, s); |
| 196 | } |
306 | } |
| 197 | |
307 | |
| 198 | void operator delete[] (void *p, size_t s) |
308 | void operator delete[] (void *p, size_t s) |
| 199 | { |
309 | { |
| 200 | g_slice_free1 (s, p); |
310 | sfree ((char *)p, s); |
| 201 | } |
311 | } |
| 202 | }; |
312 | }; |
| 203 | |
313 | |
| 204 | void *salloc_ (int n) throw (std::bad_alloc); |
314 | // makes dynamically allocated objects zero-initialised |
| 205 | void *salloc_ (int n, void *src) throw (std::bad_alloc); |
315 | struct slice_allocated |
| 206 | |
|
|
| 207 | // strictly the same as g_slice_alloc, but never returns 0 |
|
|
| 208 | template<typename T> |
|
|
| 209 | inline T *salloc (int n = 1) throw (std::bad_alloc) { return (T *)salloc_ (n * sizeof (T)); } |
|
|
| 210 | |
|
|
| 211 | // also copies src into the new area, like "memdup" |
|
|
| 212 | // if src is 0, clears the memory |
|
|
| 213 | template<typename T> |
|
|
| 214 | inline T *salloc (int n, T *src) throw (std::bad_alloc) { return (T *)salloc_ (n * sizeof (T), (void *)src); } |
|
|
| 215 | |
|
|
| 216 | // clears the memory |
|
|
| 217 | template<typename T> |
|
|
| 218 | inline T *salloc0(int n = 1) throw (std::bad_alloc) { return (T *)salloc_ (n * sizeof (T), 0); } |
|
|
| 219 | |
|
|
| 220 | // for symmetry |
|
|
| 221 | template<typename T> |
|
|
| 222 | inline void sfree (T *ptr, int n = 1) throw () |
|
|
| 223 | { |
316 | { |
| 224 | #ifdef PREFER_MALLOC |
317 | void *operator new (size_t s, void *p) |
| 225 | free (ptr); |
318 | { |
| 226 | #else |
319 | return p; |
| 227 | g_slice_free1 (n * sizeof (T), (void *)ptr); |
320 | } |
| 228 | #endif |
321 | |
| 229 | } |
322 | void *operator new (size_t s) |
|
|
323 | { |
|
|
324 | return salloc<char> (s); |
|
|
325 | } |
|
|
326 | |
|
|
327 | void *operator new[] (size_t s) |
|
|
328 | { |
|
|
329 | return salloc<char> (s); |
|
|
330 | } |
|
|
331 | |
|
|
332 | void operator delete (void *p, size_t s) |
|
|
333 | { |
|
|
334 | sfree ((char *)p, s); |
|
|
335 | } |
|
|
336 | |
|
|
337 | void operator delete[] (void *p, size_t s) |
|
|
338 | { |
|
|
339 | sfree ((char *)p, s); |
|
|
340 | } |
|
|
341 | }; |
| 230 | |
342 | |
| 231 | // a STL-compatible allocator that uses g_slice |
343 | // a STL-compatible allocator that uses g_slice |
| 232 | // boy, this is verbose |
344 | // boy, this is verbose |
| 233 | template<typename Tp> |
345 | template<typename Tp> |
| 234 | struct slice_allocator |
346 | struct slice_allocator |
| … | |
… | |
| 246 | { |
358 | { |
| 247 | typedef slice_allocator<U> other; |
359 | typedef slice_allocator<U> other; |
| 248 | }; |
360 | }; |
| 249 | |
361 | |
| 250 | slice_allocator () throw () { } |
362 | slice_allocator () throw () { } |
| 251 | slice_allocator (const slice_allocator &o) throw () { } |
363 | slice_allocator (const slice_allocator &) throw () { } |
| 252 | template<typename Tp2> |
364 | template<typename Tp2> |
| 253 | slice_allocator (const slice_allocator<Tp2> &) throw () { } |
365 | slice_allocator (const slice_allocator<Tp2> &) throw () { } |
| 254 | |
366 | |
| 255 | ~slice_allocator () { } |
367 | ~slice_allocator () { } |
| 256 | |
368 | |
| … | |
… | |
| 265 | void deallocate (pointer p, size_type n) |
377 | void deallocate (pointer p, size_type n) |
| 266 | { |
378 | { |
| 267 | sfree<Tp> (p, n); |
379 | sfree<Tp> (p, n); |
| 268 | } |
380 | } |
| 269 | |
381 | |
| 270 | size_type max_size ()const throw () |
382 | size_type max_size () const throw () |
| 271 | { |
383 | { |
| 272 | return size_t (-1) / sizeof (Tp); |
384 | return size_t (-1) / sizeof (Tp); |
| 273 | } |
385 | } |
| 274 | |
386 | |
| 275 | void construct (pointer p, const Tp &val) |
387 | void construct (pointer p, const Tp &val) |
| … | |
… | |
| 286 | // P. L'Ecuyer, “Maximally Equidistributed Combined Tausworthe Generators”, Mathematics of Computation, 65, 213 (1996), 203–213. |
398 | // P. L'Ecuyer, “Maximally Equidistributed Combined Tausworthe Generators”, Mathematics of Computation, 65, 213 (1996), 203–213. |
| 287 | // http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme.ps |
399 | // http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme.ps |
| 288 | // http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme2.ps |
400 | // http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme2.ps |
| 289 | struct tausworthe_random_generator |
401 | struct tausworthe_random_generator |
| 290 | { |
402 | { |
| 291 | // generator |
|
|
| 292 | uint32_t state [4]; |
403 | uint32_t state [4]; |
| 293 | |
404 | |
| 294 | void operator =(const tausworthe_random_generator &src) |
405 | void operator =(const tausworthe_random_generator &src) |
| 295 | { |
406 | { |
| 296 | state [0] = src.state [0]; |
407 | state [0] = src.state [0]; |
| … | |
… | |
| 299 | state [3] = src.state [3]; |
410 | state [3] = src.state [3]; |
| 300 | } |
411 | } |
| 301 | |
412 | |
| 302 | void seed (uint32_t seed); |
413 | void seed (uint32_t seed); |
| 303 | uint32_t next (); |
414 | uint32_t next (); |
|
|
415 | }; |
| 304 | |
416 | |
| 305 | // uniform distribution |
417 | // Xorshift RNGs, George Marsaglia |
|
|
418 | // http://www.jstatsoft.org/v08/i14/paper |
|
|
419 | // this one is about 40% faster than the tausworthe one above (i.e. not much), |
|
|
420 | // despite the inlining, and has the issue of only creating 2**32-1 numbers. |
|
|
421 | // see also http://www.iro.umontreal.ca/~lecuyer/myftp/papers/xorshift.pdf |
|
|
422 | struct xorshift_random_generator |
|
|
423 | { |
|
|
424 | uint32_t x, y; |
|
|
425 | |
|
|
426 | void operator =(const xorshift_random_generator &src) |
|
|
427 | { |
|
|
428 | x = src.x; |
|
|
429 | y = src.y; |
|
|
430 | } |
|
|
431 | |
|
|
432 | void seed (uint32_t seed) |
|
|
433 | { |
|
|
434 | x = seed; |
|
|
435 | y = seed * 69069U; |
|
|
436 | } |
|
|
437 | |
|
|
438 | uint32_t next () |
|
|
439 | { |
|
|
440 | uint32_t t = x ^ (x << 10); |
|
|
441 | x = y; |
|
|
442 | y = y ^ (y >> 13) ^ t ^ (t >> 10); |
|
|
443 | return y; |
|
|
444 | } |
|
|
445 | }; |
|
|
446 | |
|
|
447 | template<class generator> |
|
|
448 | struct random_number_generator : generator |
|
|
449 | { |
|
|
450 | // uniform distribution, 0 .. max (0, num - 1) |
| 306 | uint32_t operator ()(uint32_t num) |
451 | uint32_t operator ()(uint32_t num) |
| 307 | { |
452 | { |
| 308 | return is_constant (num) |
453 | return !is_constant (num) ? get_range (num) // non-constant |
| 309 | ? (next () * (uint64_t)num) >> 32U |
454 | : num & (num - 1) ? (this->next () * (uint64_t)num) >> 32U // constant, non-power-of-two |
| 310 | : get_range (num); |
455 | : this->next () & (num - 1); // constant, power-of-two |
| 311 | } |
456 | } |
| 312 | |
457 | |
| 313 | // return a number within (min .. max) |
458 | // return a number within (min .. max) |
| 314 | int operator () (int r_min, int r_max) |
459 | int operator () (int r_min, int r_max) |
| 315 | { |
460 | { |
| … | |
… | |
| 326 | protected: |
471 | protected: |
| 327 | uint32_t get_range (uint32_t r_max); |
472 | uint32_t get_range (uint32_t r_max); |
| 328 | int get_range (int r_min, int r_max); |
473 | int get_range (int r_min, int r_max); |
| 329 | }; |
474 | }; |
| 330 | |
475 | |
| 331 | typedef tausworthe_random_generator rand_gen; |
476 | typedef random_number_generator<tausworthe_random_generator> rand_gen; |
| 332 | |
477 | |
| 333 | extern rand_gen rndm; |
478 | extern rand_gen rndm, rmg_rndm; |
|
|
479 | |
|
|
480 | INTERFACE_CLASS (attachable) |
|
|
481 | struct refcnt_base |
|
|
482 | { |
|
|
483 | typedef int refcnt_t; |
|
|
484 | mutable refcnt_t ACC (RW, refcnt); |
|
|
485 | |
|
|
486 | MTH void refcnt_inc () const { ++refcnt; } |
|
|
487 | MTH void refcnt_dec () const { --refcnt; } |
|
|
488 | |
|
|
489 | refcnt_base () : refcnt (0) { } |
|
|
490 | }; |
|
|
491 | |
|
|
492 | // to avoid branches with more advanced compilers |
|
|
493 | extern refcnt_base::refcnt_t refcnt_dummy; |
| 334 | |
494 | |
| 335 | template<class T> |
495 | template<class T> |
| 336 | struct refptr |
496 | struct refptr |
| 337 | { |
497 | { |
|
|
498 | // p if not null |
|
|
499 | refcnt_base::refcnt_t *refcnt_ref () { return p ? &p->refcnt : &refcnt_dummy; } |
|
|
500 | |
|
|
501 | void refcnt_dec () |
|
|
502 | { |
|
|
503 | if (!is_constant (p)) |
|
|
504 | --*refcnt_ref (); |
|
|
505 | else if (p) |
|
|
506 | --p->refcnt; |
|
|
507 | } |
|
|
508 | |
|
|
509 | void refcnt_inc () |
|
|
510 | { |
|
|
511 | if (!is_constant (p)) |
|
|
512 | ++*refcnt_ref (); |
|
|
513 | else if (p) |
|
|
514 | ++p->refcnt; |
|
|
515 | } |
|
|
516 | |
| 338 | T *p; |
517 | T *p; |
| 339 | |
518 | |
| 340 | refptr () : p(0) { } |
519 | refptr () : p(0) { } |
| 341 | refptr (const refptr<T> &p) : p(p.p) { if (p) p->refcnt_inc (); } |
520 | refptr (const refptr<T> &p) : p(p.p) { refcnt_inc (); } |
| 342 | refptr (T *p) : p(p) { if (p) p->refcnt_inc (); } |
521 | refptr (T *p) : p(p) { refcnt_inc (); } |
| 343 | ~refptr () { if (p) p->refcnt_dec (); } |
522 | ~refptr () { refcnt_dec (); } |
| 344 | |
523 | |
| 345 | const refptr<T> &operator =(T *o) |
524 | const refptr<T> &operator =(T *o) |
| 346 | { |
525 | { |
|
|
526 | // if decrementing ever destroys we need to reverse the order here |
| 347 | if (p) p->refcnt_dec (); |
527 | refcnt_dec (); |
| 348 | p = o; |
528 | p = o; |
| 349 | if (p) p->refcnt_inc (); |
529 | refcnt_inc (); |
| 350 | |
|
|
| 351 | return *this; |
530 | return *this; |
| 352 | } |
531 | } |
| 353 | |
532 | |
| 354 | const refptr<T> &operator =(const refptr<T> o) |
533 | const refptr<T> &operator =(const refptr<T> &o) |
| 355 | { |
534 | { |
| 356 | *this = o.p; |
535 | *this = o.p; |
| 357 | return *this; |
536 | return *this; |
| 358 | } |
537 | } |
| 359 | |
538 | |
| 360 | T &operator * () const { return *p; } |
539 | T &operator * () const { return *p; } |
| 361 | T *operator ->() const { return p; } |
540 | T *operator ->() const { return p; } |
| 362 | |
541 | |
| 363 | operator T *() const { return p; } |
542 | operator T *() const { return p; } |
| 364 | }; |
543 | }; |
| 365 | |
544 | |
| 366 | typedef refptr<maptile> maptile_ptr; |
545 | typedef refptr<maptile> maptile_ptr; |
| … | |
… | |
| 371 | |
550 | |
| 372 | struct str_hash |
551 | struct str_hash |
| 373 | { |
552 | { |
| 374 | std::size_t operator ()(const char *s) const |
553 | std::size_t operator ()(const char *s) const |
| 375 | { |
554 | { |
| 376 | unsigned long hash = 0; |
555 | #if 0 |
|
|
556 | uint32_t hash = 0; |
| 377 | |
557 | |
| 378 | /* use the one-at-a-time hash function, which supposedly is |
558 | /* use the one-at-a-time hash function, which supposedly is |
| 379 | * better than the djb2-like one used by perl5.005, but |
559 | * better than the djb2-like one used by perl5.005, but |
| 380 | * certainly is better then the bug used here before. |
560 | * certainly is better then the bug used here before. |
| 381 | * see http://burtleburtle.net/bob/hash/doobs.html |
561 | * see http://burtleburtle.net/bob/hash/doobs.html |
| … | |
… | |
| 388 | } |
568 | } |
| 389 | |
569 | |
| 390 | hash += hash << 3; |
570 | hash += hash << 3; |
| 391 | hash ^= hash >> 11; |
571 | hash ^= hash >> 11; |
| 392 | hash += hash << 15; |
572 | hash += hash << 15; |
|
|
573 | #else |
|
|
574 | // use FNV-1a hash (http://isthe.com/chongo/tech/comp/fnv/) |
|
|
575 | // it is about twice as fast as the one-at-a-time one, |
|
|
576 | // with good distribution. |
|
|
577 | // FNV-1a is faster on many cpus because the multiplication |
|
|
578 | // runs concurrent with the looping logic. |
|
|
579 | uint32_t hash = 2166136261; |
|
|
580 | |
|
|
581 | while (*s) |
|
|
582 | hash = (hash ^ *s++) * 16777619; |
|
|
583 | #endif |
| 393 | |
584 | |
| 394 | return hash; |
585 | return hash; |
| 395 | } |
586 | } |
| 396 | }; |
587 | }; |
| 397 | |
588 | |
| … | |
… | |
| 401 | { |
592 | { |
| 402 | return !strcmp (a, b); |
593 | return !strcmp (a, b); |
| 403 | } |
594 | } |
| 404 | }; |
595 | }; |
| 405 | |
596 | |
|
|
597 | // Mostly the same as std::vector, but insert/erase can reorder |
|
|
598 | // the elements, making append(=insert)/remove O(1) instead of O(n). |
|
|
599 | // |
|
|
600 | // NOTE: only some forms of erase are available |
| 406 | template<class T> |
601 | template<class T> |
| 407 | struct unordered_vector : std::vector<T, slice_allocator<T> > |
602 | struct unordered_vector : std::vector<T, slice_allocator<T> > |
| 408 | { |
603 | { |
| 409 | typedef typename unordered_vector::iterator iterator; |
604 | typedef typename unordered_vector::iterator iterator; |
| 410 | |
605 | |
| … | |
… | |
| 420 | { |
615 | { |
| 421 | erase ((unsigned int )(i - this->begin ())); |
616 | erase ((unsigned int )(i - this->begin ())); |
| 422 | } |
617 | } |
| 423 | }; |
618 | }; |
| 424 | |
619 | |
| 425 | template<class T, int T::* index> |
620 | // This container blends advantages of linked lists |
|
|
621 | // (efficiency) with vectors (random access) by |
|
|
622 | // by using an unordered vector and storing the vector |
|
|
623 | // index inside the object. |
|
|
624 | // |
|
|
625 | // + memory-efficient on most 64 bit archs |
|
|
626 | // + O(1) insert/remove |
|
|
627 | // + free unique (but varying) id for inserted objects |
|
|
628 | // + cache-friendly iteration |
|
|
629 | // - only works for pointers to structs |
|
|
630 | // |
|
|
631 | // NOTE: only some forms of erase/insert are available |
|
|
632 | typedef int object_vector_index; |
|
|
633 | |
|
|
634 | template<class T, object_vector_index T::*indexmember> |
| 426 | struct object_vector : std::vector<T *, slice_allocator<T *> > |
635 | struct object_vector : std::vector<T *, slice_allocator<T *> > |
| 427 | { |
636 | { |
|
|
637 | typedef typename object_vector::iterator iterator; |
|
|
638 | |
|
|
639 | bool contains (const T *obj) const |
|
|
640 | { |
|
|
641 | return obj->*indexmember; |
|
|
642 | } |
|
|
643 | |
|
|
644 | iterator find (const T *obj) |
|
|
645 | { |
|
|
646 | return obj->*indexmember |
|
|
647 | ? this->begin () + obj->*indexmember - 1 |
|
|
648 | : this->end (); |
|
|
649 | } |
|
|
650 | |
|
|
651 | void push_back (T *obj) |
|
|
652 | { |
|
|
653 | std::vector<T *, slice_allocator<T *> >::push_back (obj); |
|
|
654 | obj->*indexmember = this->size (); |
|
|
655 | } |
|
|
656 | |
| 428 | void insert (T *obj) |
657 | void insert (T *obj) |
| 429 | { |
658 | { |
| 430 | assert (!(obj->*index)); |
|
|
| 431 | push_back (obj); |
659 | push_back (obj); |
| 432 | obj->*index = this->size (); |
|
|
| 433 | } |
660 | } |
| 434 | |
661 | |
| 435 | void insert (T &obj) |
662 | void insert (T &obj) |
| 436 | { |
663 | { |
| 437 | insert (&obj); |
664 | insert (&obj); |
| 438 | } |
665 | } |
| 439 | |
666 | |
| 440 | void erase (T *obj) |
667 | void erase (T *obj) |
| 441 | { |
668 | { |
| 442 | assert (obj->*index); |
|
|
| 443 | unsigned int pos = obj->*index; |
669 | unsigned int pos = obj->*indexmember; |
| 444 | obj->*index = 0; |
670 | obj->*indexmember = 0; |
| 445 | |
671 | |
| 446 | if (pos < this->size ()) |
672 | if (pos < this->size ()) |
| 447 | { |
673 | { |
| 448 | (*this)[pos - 1] = (*this)[this->size () - 1]; |
674 | (*this)[pos - 1] = (*this)[this->size () - 1]; |
| 449 | (*this)[pos - 1]->*index = pos; |
675 | (*this)[pos - 1]->*indexmember = pos; |
| 450 | } |
676 | } |
| 451 | |
677 | |
| 452 | this->pop_back (); |
678 | this->pop_back (); |
| 453 | } |
679 | } |
| 454 | |
680 | |
| 455 | void erase (T &obj) |
681 | void erase (T &obj) |
| 456 | { |
682 | { |
| 457 | errase (&obj); |
683 | erase (&obj); |
| 458 | } |
684 | } |
| 459 | }; |
685 | }; |
| 460 | |
686 | |
| 461 | // basically does what strncpy should do, but appends "..." to strings exceeding length |
687 | // basically does what strncpy should do, but appends "..." to strings exceeding length |
|
|
688 | // returns the number of bytes actually used (including \0) |
| 462 | void assign (char *dst, const char *src, int maxlen); |
689 | int assign (char *dst, const char *src, int maxsize); |
| 463 | |
690 | |
| 464 | // type-safe version of assign |
691 | // type-safe version of assign |
| 465 | template<int N> |
692 | template<int N> |
| 466 | inline void assign (char (&dst)[N], const char *src) |
693 | inline int assign (char (&dst)[N], const char *src) |
| 467 | { |
694 | { |
| 468 | assign ((char *)&dst, src, N); |
695 | return assign ((char *)&dst, src, N); |
| 469 | } |
696 | } |
| 470 | |
697 | |
| 471 | typedef double tstamp; |
698 | typedef double tstamp; |
| 472 | |
699 | |
| 473 | // return current time as timestampe |
700 | // return current time as timestamp |
| 474 | tstamp now (); |
701 | tstamp now (); |
| 475 | |
702 | |
| 476 | int similar_direction (int a, int b); |
703 | int similar_direction (int a, int b); |
| 477 | |
704 | |
| 478 | // like printf, but returns a std::string |
705 | // like v?sprintf, but returns a "static" buffer |
| 479 | const std::string format (const char *format, ...); |
706 | char *vformat (const char *format, va_list ap); |
|
|
707 | char *format (const char *format, ...) attribute ((format (printf, 1, 2))); |
| 480 | |
708 | |
|
|
709 | // safety-check player input which will become object->msg |
|
|
710 | bool msg_is_safe (const char *msg); |
|
|
711 | |
|
|
712 | ///////////////////////////////////////////////////////////////////////////// |
|
|
713 | // threads, very very thin wrappers around pthreads |
|
|
714 | |
|
|
715 | struct thread |
|
|
716 | { |
|
|
717 | pthread_t id; |
|
|
718 | |
|
|
719 | void start (void *(*start_routine)(void *), void *arg = 0); |
|
|
720 | |
|
|
721 | void cancel () |
|
|
722 | { |
|
|
723 | pthread_cancel (id); |
|
|
724 | } |
|
|
725 | |
|
|
726 | void *join () |
|
|
727 | { |
|
|
728 | void *ret; |
|
|
729 | |
|
|
730 | if (pthread_join (id, &ret)) |
|
|
731 | cleanup ("pthread_join failed", 1); |
|
|
732 | |
|
|
733 | return ret; |
|
|
734 | } |
|
|
735 | }; |
|
|
736 | |
|
|
737 | // note that mutexes are not classes |
|
|
738 | typedef pthread_mutex_t smutex; |
|
|
739 | |
|
|
740 | #if __linux && defined (PTHREAD_ADAPTIVE_MUTEX_INITIALIZER_NP) |
|
|
741 | #define SMUTEX_INITIALISER PTHREAD_ADAPTIVE_MUTEX_INITIALIZER_NP |
|
|
742 | #else |
|
|
743 | #define SMUTEX_INITIALISER PTHREAD_MUTEX_INITIALIZER |
| 481 | #endif |
744 | #endif |
| 482 | |
745 | |
|
|
746 | #define SMUTEX(name) smutex name = SMUTEX_INITIALISER |
|
|
747 | #define SMUTEX_LOCK(name) pthread_mutex_lock (&(name)) |
|
|
748 | #define SMUTEX_UNLOCK(name) pthread_mutex_unlock (&(name)) |
|
|
749 | |
|
|
750 | typedef pthread_cond_t scond; |
|
|
751 | |
|
|
752 | #define SCOND(name) scond name = PTHREAD_COND_INITIALIZER |
|
|
753 | #define SCOND_SIGNAL(name) pthread_cond_signal (&(name)) |
|
|
754 | #define SCOND_BROADCAST(name) pthread_cond_broadcast (&(name)) |
|
|
755 | #define SCOND_WAIT(name,mutex) pthread_cond_wait (&(name), &(mutex)) |
|
|
756 | |
|
|
757 | #endif |
|
|
758 | |
