utils.h

#ifndef LIB2GEOM_UTILS_HEADER
#define LIB2GEOM_UTILS_HEADER

#include <cmath>

namespace Geom {

// proper logical xor
inline bool logical_xor (bool a, bool b) { return (a || b) && !(a && b); }

template <class T> inline int sgn(const T& x) {return (x < 0 ? -1 : (x > 0 ? 1 : 0) );}

template <class T> inline T sqr(const T& x) {return x * x;}
template <class T> inline T cube(const T& x) {return x * x * x;}

template <class T> inline const T& between (const T& min, const T& max, const T& x)
    { return min < x && max > x; }

inline double round(double const x) { return std::floor(x + .5); }

inline double decimal_round(double const x, int const places) {
    //TODO: possibly implement with modulus instead?
    double const multiplier = std::pow(10.0, places);
    return round( x * multiplier ) / multiplier;
}

}

#endif