#ifndef __MATH_AUX_H
#define __MATH_AUX_H

#include <limits>

/*
 * This function is not only used to compare floating point values with a tolerance,
 * it also leads to ambiguity problems if the called functions do not have the
 * right prototype.
 */
template <class _Tp>
bool are_equals(_Tp val, _Tp ref) {
  if (val < ref) {
    return (ref - val) <= std::numeric_limits<_Tp>::epsilon();
  }
  else {
    return (val - ref) <= std::numeric_limits<_Tp>::epsilon();
  }
}

#endif // __MATH_AUX_H