/*! ========================================================================
** Extended Template and Library
** Misc
-** $Id: _misc.h,v 1.1.1.1 2005/01/04 01:31:47 darco Exp $
+** $Id$
**
** Copyright (c) 2002 Robert B. Quattlebaum Jr.
+** Copyright (c) 2007 Chris Moore
**
** This package is free software; you can redistribute it and/or
** modify it under the terms of the GNU General Public License as
{
#if 1
I iter(begin+(end-begin)/2);
-
+
while(end-begin>1 && !(*iter==value))
{
((*iter<value)?begin:end) = iter;
-
+
iter = begin+(end-begin)/2;
}
return iter;
#else
size_t len_(end-begin);
size_t half_(len_/2);
-
+
I iter(begin);
iter+=half_;
-
+
while(len_>1 && !(*iter==value))
{
((*iter<value)?begin:end) = iter;
-
+
len_=half_;
half_/=2;
-
+
iter=begin;
iter+=half_;
}
#endif
}
-/*! \fixme Isn't there some x86 FPU instruction for quickly
-** converting a float to a rounded integer? It's worth
-** looking into at some point... */
-inline int round_to_int(const float x) { return static_cast<int>(x+0.5f); }
-inline int round_to_int(const double x) { return static_cast<int>(x+0.5); }
+inline int round_to_int(const float x) {
+ /*! \todo Isn't there some x86 FPU instruction for quickly
+ ** converting a float to a rounded integer? It's worth
+ ** looking into at some point... */
+ // 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);
+}
inline int ceil_to_int(const float x) { return static_cast<int>(ceil(x)); }
inline int ceil_to_int(const double x) { return static_cast<int>(ceil(x)); }