… | |
… | |
99 | |
99 | |
100 | // sign0 returns -1, 0 or +1 |
100 | // sign0 returns -1, 0 or +1 |
101 | template<typename T> |
101 | template<typename T> |
102 | static inline T sign0 (T v) { return v ? sign (v) : 0; } |
102 | static inline T sign0 (T v) { return v ? sign (v) : 0; } |
103 | |
103 | |
|
|
104 | template<typename T, typename U> |
|
|
105 | static inline T copysign (T a, U b) { return a > 0 ? b : -b; } |
|
|
106 | |
104 | // div* only work correctly for div > 0 |
107 | // div* only work correctly for div > 0 |
105 | // div, with correct rounding (< 0.5 downwards, >=0.5 upwards) |
108 | // div, with correct rounding (< 0.5 downwards, >=0.5 upwards) |
106 | template<typename T> static inline T div (T val, T div) |
109 | template<typename T> static inline T div (T val, T div) |
107 | { |
110 | { |
108 | return expect_false (val < 0) ? - ((-val + (div - 1) / 2) / div) : (val + div / 2) / div; |
111 | return expect_false (val < 0) ? - ((-val + (div - 1) / 2) / div) : (val + div / 2) / div; |