+ // return static_cast<int>(x+0.5f); // <-- (a) fast, but rounds -1.333 to 0!
+ // return static_cast<int>(rintf(x)); // <-- (b) slow, but correct
+ if (x>=0) return static_cast<int>(x + 0.5); // <-- slower than (a), but correct, and faster than (b)
+ else return static_cast<int>(x - 0.5);
+}
+inline int round_to_int(const double x) {
+ // return static_cast<int>(x+0.5);
+ // return static_cast<int>(rint(x));
+ if (x>=0) return static_cast<int>(x + 0.5);
+ else return static_cast<int>(x - 0.5);