… | |
… | |
55 | #endif |
55 | #endif |
56 | |
56 | |
57 | // use C0X decltype for auto declarations until ISO C++ sanctifies them (if ever) |
57 | // use C0X decltype for auto declarations until ISO C++ sanctifies them (if ever) |
58 | #define auto(var,expr) decltype(expr) var = (expr) |
58 | #define auto(var,expr) decltype(expr) var = (expr) |
59 | |
59 | |
60 | // could use the sizeof (arr) /( sizeof (arr [0]) here, but C++ is |
60 | #if cplusplus_does_not_suck |
61 | // much more obfuscated... :) |
61 | // does not work for local types (http://www.open-std.org/jtc1/sc22/wg21/docs/papers/2008/n2657.htm) |
62 | |
|
|
63 | template<typename T, int N> |
62 | template<typename T, int N> |
64 | inline int array_length (const T (&arr)[N]) |
63 | static inline int array_length (const T (&arr)[N]) |
65 | { |
64 | { |
66 | return N; |
65 | return N; |
67 | } |
66 | } |
|
|
67 | #else |
|
|
68 | #define array_length(name) (sizeof (name) / sizeof (name [0])) |
|
|
69 | #endif |
68 | |
70 | |
69 | // very ugly macro that basically declares and initialises a variable |
71 | // very ugly macro that basically declares and initialises a variable |
70 | // that is in scope for the next statement only |
72 | // that is in scope for the next statement only |
71 | // works only for stuff that can be assigned 0 and converts to false |
73 | // works only for stuff that can be assigned 0 and converts to false |
72 | // (note: works great for pointers) |
74 | // (note: works great for pointers) |
… | |
… | |
121 | // div, with correct rounding (< 0.5 downwards, >=0.5 upwards) |
123 | // div, with correct rounding (< 0.5 downwards, >=0.5 upwards) |
122 | template<typename T> static inline T div (T val, T div) |
124 | template<typename T> static inline T div (T val, T div) |
123 | { |
125 | { |
124 | return expect_false (val < 0) ? - ((-val + (div - 1) / 2) / div) : (val + div / 2) / div; |
126 | return expect_false (val < 0) ? - ((-val + (div - 1) / 2) / div) : (val + div / 2) / div; |
125 | } |
127 | } |
|
|
128 | |
|
|
129 | template<> inline float div (float val, float div) { return val / div; } |
|
|
130 | template<> inline double div (double val, double div) { return val / div; } |
|
|
131 | |
126 | // div, round-up |
132 | // div, round-up |
127 | template<typename T> static inline T div_ru (T val, T div) |
133 | template<typename T> static inline T div_ru (T val, T div) |
128 | { |
134 | { |
129 | return expect_false (val < 0) ? - ((-val ) / div) : (val + div - 1) / div; |
135 | return expect_false (val < 0) ? - ((-val ) / div) : (val + div - 1) / div; |
130 | } |
136 | } |
… | |
… | |
475 | }; |
481 | }; |
476 | |
482 | |
477 | template<class generator> |
483 | template<class generator> |
478 | struct random_number_generator : generator |
484 | struct random_number_generator : generator |
479 | { |
485 | { |
480 | // uniform distribution, 0 .. max (0, num - 1) |
486 | // uniform distribution, [0 .. num - 1] |
481 | uint32_t operator ()(uint32_t num) |
487 | uint32_t operator ()(uint32_t num) |
482 | { |
488 | { |
483 | return !is_constant (num) ? get_range (num) // non-constant |
489 | return !is_constant (num) ? get_range (num) // non-constant |
484 | : num & (num - 1) ? (this->next () * (uint64_t)num) >> 32U // constant, non-power-of-two |
490 | : num & (num - 1) ? (this->next () * (uint64_t)num) >> 32U // constant, non-power-of-two |
485 | : this->next () & (num - 1); // constant, power-of-two |
491 | : this->next () & (num - 1); // constant, power-of-two |
486 | } |
492 | } |
487 | |
493 | |
488 | // return a number within the closed interval [min .. max] |
494 | // return a number within the closed interval [min .. max], max can be >, < or == min. |
489 | int operator () (int r_min, int r_max) |
495 | int operator () (int r_min, int r_max) |
490 | { |
496 | { |
491 | return is_constant (r_min <= r_max) && r_min <= r_max |
497 | return is_constant (r_min <= r_max) && r_min <= r_max |
492 | ? r_min + operator ()(r_max - r_min + 1) |
498 | ? r_min + operator ()(r_max - r_min + 1) |
493 | : get_range (r_min, r_max); |
499 | : get_range (r_min, r_max); |
494 | } |
500 | } |
495 | |
501 | |
496 | // return a number within the closed interval [0..1] |
502 | // return a number within the half-open interval [0..1[ |
497 | double operator ()() |
503 | double operator ()() |
498 | { |
504 | { |
499 | return this->next () / (double)0xFFFFFFFFU; |
505 | return this->next () / (double)0x100000000ULL; |
500 | } |
506 | } |
501 | |
507 | |
502 | protected: |
508 | protected: |
503 | uint32_t get_range (uint32_t r_max); |
509 | uint32_t get_range (uint32_t r_max); |
504 | int get_range (int r_min, int r_max); |
510 | int get_range (int r_min, int r_max); |