37 static constexpr double minMagnitude = std::numeric_limits<double>::min();
42 static constexpr double ulpFactor = std::numeric_limits<double>::epsilon();
47 const double tolerance = maxMagnitude * (ulpFactor * (maxUlpsDiff + 1));
51 return a ==
b || std::abs(
b -
a) < tolerance;