1 | /* |
1 | /* |
2 | * This file is part of Deliantra, the Roguelike Realtime MMORPG. |
2 | * This file is part of Deliantra, the Roguelike Realtime MMORPG. |
3 | * |
3 | * |
4 | * Copyright (©) 2005,2006,2007 Marc Alexander Lehmann / Robin Redeker / the Deliantra team |
4 | * Copyright (©) 2005,2006,2007,2008 Marc Alexander Lehmann / Robin Redeker / the Deliantra team |
5 | * |
5 | * |
6 | * Deliantra is free software: you can redistribute it and/or modify |
6 | * Deliantra is free software: you can redistribute it and/or modify |
7 | * it under the terms of the GNU General Public License as published by |
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 |
8 | * the Free Software Foundation, either version 3 of the License, or |
9 | * (at your option) any later version. |
9 | * (at your option) any later version. |
… | |
… | |
72 | #endif |
72 | #endif |
73 | |
73 | |
74 | // use C0X decltype for auto declarations until ISO C++ sanctifies them (if ever) |
74 | // use C0X decltype for auto declarations until ISO C++ sanctifies them (if ever) |
75 | #define auto(var,expr) decltype(expr) var = (expr) |
75 | #define auto(var,expr) decltype(expr) var = (expr) |
76 | |
76 | |
77 | // very ugly macro that basicaly declares and initialises a variable |
77 | // very ugly macro that basically declares and initialises a variable |
78 | // that is in scope for the next statement only |
78 | // that is in scope for the next statement only |
79 | // works only for stuff that can be assigned 0 and converts to false |
79 | // works only for stuff that can be assigned 0 and converts to false |
80 | // (note: works great for pointers) |
80 | // (note: works great for pointers) |
81 | // most ugly macro I ever wrote |
81 | // most ugly macro I ever wrote |
82 | #define statementvar(type, name, value) if (type name = 0) { } else if (((name) = (value)), 1) |
82 | #define statementvar(type, name, value) if (type name = 0) { } else if (((name) = (value)), 1) |
… | |
… | |
96 | // as a is often a constant while b is the variable. it is still a bug, though. |
96 | // as a is often a constant while b is the variable. it is still a bug, though. |
97 | template<typename T, typename U> static inline T min (T a, U b) { return (U)a < b ? (U)a : b; } |
97 | template<typename T, typename U> static inline T min (T a, U b) { return (U)a < b ? (U)a : b; } |
98 | template<typename T, typename U> static inline T max (T a, U b) { return (U)a > b ? (U)a : b; } |
98 | template<typename T, typename U> static inline T max (T a, U b) { return (U)a > b ? (U)a : b; } |
99 | 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; } |
99 | 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; } |
100 | |
100 | |
|
|
101 | template<typename T, typename U> static inline void min_it (T &v, U m) { v = min (v, (T)m); } |
|
|
102 | template<typename T, typename U> static inline void max_it (T &v, U m) { v = max (v, (T)m); } |
|
|
103 | 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); } |
|
|
104 | |
101 | template<typename T, typename U> static inline void swap (T& a, U& b) { T t=a; a=(T)b; b=(U)t; } |
105 | template<typename T, typename U> static inline void swap (T& a, U& b) { T t=a; a=(T)b; b=(U)t; } |
102 | |
106 | |
103 | template<typename T, typename U, typename V> static inline T min (T a, U b, V c) { return min (a, min (b, c)); } |
107 | template<typename T, typename U, typename V> static inline T min (T a, U b, V c) { return min (a, min (b, c)); } |
104 | template<typename T, typename U, typename V> static inline T max (T a, U b, V c) { return max (a, max (b, c)); } |
108 | template<typename T, typename U, typename V> static inline T max (T a, U b, V c) { return max (a, max (b, c)); } |
105 | |
109 | |
|
|
110 | // sign returns -1 or +1 |
|
|
111 | template<typename T> |
|
|
112 | static inline T sign (T v) { return v < 0 ? -1 : +1; } |
|
|
113 | // relies on 2c representation |
|
|
114 | template<> |
|
|
115 | inline sint8 sign (sint8 v) { return 1 - (sint8 (uint8 (v) >> 7) * 2); } |
|
|
116 | |
|
|
117 | // sign0 returns -1, 0 or +1 |
|
|
118 | template<typename T> |
|
|
119 | static inline T sign0 (T v) { return v ? sign (v) : 0; } |
|
|
120 | |
|
|
121 | // div, with correct rounding (< 0.5 downwards, >=0.5 upwards) |
|
|
122 | template<typename T> static inline T div (T val, T div) { return (val + div / 2) / div; } |
|
|
123 | // div, round-up |
|
|
124 | template<typename T> static inline T div_ru (T val, T div) { return (val + div - 1) / div; } |
|
|
125 | // div, round-down |
|
|
126 | template<typename T> static inline T div_rd (T val, T div) { return (val ) / div; } |
|
|
127 | |
106 | template<typename T> |
128 | template<typename T> |
107 | static inline T |
129 | static inline T |
108 | lerp (T val, T min_in, T max_in, T min_out, T max_out) |
130 | lerp (T val, T min_in, T max_in, T min_out, T max_out) |
109 | { |
131 | { |
110 | return (val - min_in) * (max_out - min_out) / (max_in - min_in) + min_out; |
132 | return min_out + div <T> ((val - min_in) * (max_out - min_out), max_in - min_in); |
|
|
133 | } |
|
|
134 | |
|
|
135 | // lerp, round-down |
|
|
136 | template<typename T> |
|
|
137 | static inline T |
|
|
138 | lerp_rd (T val, T min_in, T max_in, T min_out, T max_out) |
|
|
139 | { |
|
|
140 | return min_out + div_rd<T> ((val - min_in) * (max_out - min_out), max_in - min_in); |
|
|
141 | } |
|
|
142 | |
|
|
143 | // lerp, round-up |
|
|
144 | template<typename T> |
|
|
145 | static inline T |
|
|
146 | lerp_ru (T val, T min_in, T max_in, T min_out, T max_out) |
|
|
147 | { |
|
|
148 | return min_out + div_ru<T> ((val - min_in) * (max_out - min_out), max_in - min_in); |
111 | } |
149 | } |
112 | |
150 | |
113 | // lots of stuff taken from FXT |
151 | // lots of stuff taken from FXT |
114 | |
152 | |
115 | /* Rotate right. This is used in various places for checksumming */ |
153 | /* Rotate right. This is used in various places for checksumming */ |
… | |
… | |
224 | g_slice_free1 (n * sizeof (T), (void *)ptr); |
262 | g_slice_free1 (n * sizeof (T), (void *)ptr); |
225 | assert (slice_alloc >= 0);//D |
263 | assert (slice_alloc >= 0);//D |
226 | } |
264 | } |
227 | } |
265 | } |
228 | |
266 | |
|
|
267 | // nulls the pointer |
|
|
268 | template<typename T> |
|
|
269 | inline void sfree0 (T *&ptr, int n = 1) throw () |
|
|
270 | { |
|
|
271 | sfree<T> (ptr, n); |
|
|
272 | ptr = 0; |
|
|
273 | } |
|
|
274 | |
229 | // makes dynamically allocated objects zero-initialised |
275 | // makes dynamically allocated objects zero-initialised |
230 | struct zero_initialised |
276 | struct zero_initialised |
231 | { |
277 | { |
232 | void *operator new (size_t s, void *p) |
278 | void *operator new (size_t s, void *p) |
233 | { |
279 | { |
… | |
… | |
241 | } |
287 | } |
242 | |
288 | |
243 | void *operator new[] (size_t s) |
289 | void *operator new[] (size_t s) |
244 | { |
290 | { |
245 | return salloc0<char> (s); |
291 | return salloc0<char> (s); |
|
|
292 | } |
|
|
293 | |
|
|
294 | void operator delete (void *p, size_t s) |
|
|
295 | { |
|
|
296 | sfree ((char *)p, s); |
|
|
297 | } |
|
|
298 | |
|
|
299 | void operator delete[] (void *p, size_t s) |
|
|
300 | { |
|
|
301 | sfree ((char *)p, s); |
|
|
302 | } |
|
|
303 | }; |
|
|
304 | |
|
|
305 | // makes dynamically allocated objects zero-initialised |
|
|
306 | struct slice_allocated |
|
|
307 | { |
|
|
308 | void *operator new (size_t s, void *p) |
|
|
309 | { |
|
|
310 | return p; |
|
|
311 | } |
|
|
312 | |
|
|
313 | void *operator new (size_t s) |
|
|
314 | { |
|
|
315 | return salloc<char> (s); |
|
|
316 | } |
|
|
317 | |
|
|
318 | void *operator new[] (size_t s) |
|
|
319 | { |
|
|
320 | return salloc<char> (s); |
246 | } |
321 | } |
247 | |
322 | |
248 | void operator delete (void *p, size_t s) |
323 | void operator delete (void *p, size_t s) |
249 | { |
324 | { |
250 | sfree ((char *)p, s); |
325 | sfree ((char *)p, s); |
… | |
… | |
314 | // P. L'Ecuyer, “Maximally Equidistributed Combined Tausworthe Generators”, Mathematics of Computation, 65, 213 (1996), 203–213. |
389 | // P. L'Ecuyer, “Maximally Equidistributed Combined Tausworthe Generators”, Mathematics of Computation, 65, 213 (1996), 203–213. |
315 | // http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme.ps |
390 | // http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme.ps |
316 | // http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme2.ps |
391 | // http://www.iro.umontreal.ca/~lecuyer/myftp/papers/tausme2.ps |
317 | struct tausworthe_random_generator |
392 | struct tausworthe_random_generator |
318 | { |
393 | { |
319 | // generator |
|
|
320 | uint32_t state [4]; |
394 | uint32_t state [4]; |
321 | |
395 | |
322 | void operator =(const tausworthe_random_generator &src) |
396 | void operator =(const tausworthe_random_generator &src) |
323 | { |
397 | { |
324 | state [0] = src.state [0]; |
398 | state [0] = src.state [0]; |
… | |
… | |
327 | state [3] = src.state [3]; |
401 | state [3] = src.state [3]; |
328 | } |
402 | } |
329 | |
403 | |
330 | void seed (uint32_t seed); |
404 | void seed (uint32_t seed); |
331 | uint32_t next (); |
405 | uint32_t next (); |
|
|
406 | }; |
332 | |
407 | |
333 | // uniform distribution |
408 | // Xorshift RNGs, George Marsaglia |
|
|
409 | // http://www.jstatsoft.org/v08/i14/paper |
|
|
410 | // this one is about 40% faster than the tausworthe one above (i.e. not much), |
|
|
411 | // despite the inlining, and has the issue of only creating 2**32-1 numbers. |
|
|
412 | struct xorshift_random_generator |
|
|
413 | { |
|
|
414 | uint32_t x, y; |
|
|
415 | |
|
|
416 | void operator =(const xorshift_random_generator &src) |
|
|
417 | { |
|
|
418 | x = src.x; |
|
|
419 | y = src.y; |
|
|
420 | } |
|
|
421 | |
|
|
422 | void seed (uint32_t seed) |
|
|
423 | { |
|
|
424 | x = seed; |
|
|
425 | y = seed * 69069U; |
|
|
426 | } |
|
|
427 | |
|
|
428 | uint32_t next () |
|
|
429 | { |
|
|
430 | uint32_t t = x ^ (x << 10); |
|
|
431 | x = y; |
|
|
432 | y = y ^ (y >> 13) ^ t ^ (t >> 10); |
|
|
433 | return y; |
|
|
434 | } |
|
|
435 | }; |
|
|
436 | |
|
|
437 | template<class generator> |
|
|
438 | struct random_number_generator : generator |
|
|
439 | { |
|
|
440 | // uniform distribution, 0 .. max (0, num - 1) |
334 | uint32_t operator ()(uint32_t num) |
441 | uint32_t operator ()(uint32_t num) |
335 | { |
442 | { |
336 | return is_constant (num) |
443 | return !is_constant (num) ? get_range (num) // non-constant |
337 | ? (next () * (uint64_t)num) >> 32U |
444 | : num & (num - 1) ? (this->next () * (uint64_t)num) >> 32U // constant, non-power-of-two |
338 | : get_range (num); |
445 | : this->next () & (num - 1); // constant, power-of-two |
339 | } |
446 | } |
340 | |
447 | |
341 | // return a number within (min .. max) |
448 | // return a number within (min .. max) |
342 | int operator () (int r_min, int r_max) |
449 | int operator () (int r_min, int r_max) |
343 | { |
450 | { |
… | |
… | |
354 | protected: |
461 | protected: |
355 | uint32_t get_range (uint32_t r_max); |
462 | uint32_t get_range (uint32_t r_max); |
356 | int get_range (int r_min, int r_max); |
463 | int get_range (int r_min, int r_max); |
357 | }; |
464 | }; |
358 | |
465 | |
359 | typedef tausworthe_random_generator rand_gen; |
466 | typedef random_number_generator<tausworthe_random_generator> rand_gen; |
360 | |
467 | |
361 | extern rand_gen rndm; |
468 | extern rand_gen rndm, rmg_rndm; |
362 | |
469 | |
363 | INTERFACE_CLASS (attachable) |
470 | INTERFACE_CLASS (attachable) |
364 | struct refcnt_base |
471 | struct refcnt_base |
365 | { |
472 | { |
366 | typedef int refcnt_t; |
473 | typedef int refcnt_t; |