FloatComparison.hpp

Go to the documentation of this file.

#pragma once
 
#include <algorithm>
#include <cmath>
#include <limits>
#include <type_traits>
 
namespace volcart
{
 
static constexpr double DEFAULT_MAX_DIFFERENCE = 1e-7;
 
template <
    typename T,
    std::enable_if_t<std::is_floating_point<T>::value, bool> = true>
inline auto AlmostEqual(
    const T lhs,
    const T rhs,
    T epsAbs = static_cast<T>(DEFAULT_MAX_DIFFERENCE),
    T epsRel = std::numeric_limits<T>::epsilon()) -> bool
{
    // Note: For some reason, std::abs resolves to integral overload? - SP
    T d = std::fabs(lhs - rhs);
    if (d <= epsAbs) {
        return true;
    }
    T l = std::fabs(lhs);
    T r = std::fabs(rhs);
 
    T largest = std::max(r, l);
    return d <= largest * epsRel;
}
 
}  // namespace volcart