How should I do floating point comparison?

by

Comparing for greater/smaller is not really a problem unless you’re working right at the edge of the float/double precision limit.

For a “fuzzy equals” comparison, this (Java code, should be easy to adapt) is what I came up with for The Floating-Point Guide after a lot of work and taking into account lots of criticism:

public static boolean nearlyEqual(float a, float b, float epsilon) {
    final float absA = Math.abs(a);
    final float absB = Math.abs(b);
    final float diff = Math.abs(a - b);

    if (a == b) { // shortcut, handles infinities
        return true;
    } else if (a == 0 || b == 0 || diff < Float.MIN_NORMAL) {
        // a or b is zero or both are extremely close to it
        // relative error is less meaningful here
        return diff < (epsilon * Float.MIN_NORMAL);
    } else { // use relative error
        return diff / (absA + absB) < epsilon;
    }
}

It comes with a test suite. You should immediately dismiss any solution that doesn’t, because it is virtually guaranteed to fail in some edge cases like having one value 0, two very small values opposite of zero, or infinities.

An alternative (see link above for more details) is to convert the floats’ bit patterns to integer and accept everything within a fixed integer distance.

In any case, there probably isn’t any solution that is perfect for all applications. Ideally, you’d develop/adapt your own with a test suite covering your actual use cases.

More Related Contents:

  1. Should we compare floating point numbers for equality against a *relative* error?
  2. How do I display the binary representation of a float or double?
  3. Why does division by zero in IEEE754 standard results in Infinite value?
  4. String Compare “Logic”
  5. What’s the difference between passing by reference vs. passing by value?
  6. Shortest distance between a point and a line segment
  7. What’s your most controversial programming opinion?
  8. What is a callback function?
  9. Strangest language feature
  10. Runtime vs. Compile time
  11. De Morgan’s rules explained
  12. What is boxing and unboxing and what are the trade offs?
  13. Code Golf: Number to Words
  14. Is it OK to use DYLD_LIBRARY_PATH on Mac OS X? And, what’s the dynamic library search algorithm with it?
  15. Algorithm to convert an IEEE 754 double to a string?
  16. How can I find the difference between two angles?
  17. What’s the point of OOP?
  18. Does float have a negative zero? (-0f)
  19. Acronyms in CamelCase [closed]
  20. What are the advantages of memory-mapped files?
  21. What is the best way to replace or substitute if..else if..else trees in programs?
  22. MySQL floating point comparison issues
  23. How best to sum up lots of floating point numbers?
  24. Is Sleep() evil?
  25. Recursion or iteration?
  26. While vs. Do While
  27. What are “magic numbers” in computer programming?
  28. Find all complete sub-graphs within a graph
  29. How to design extensible software (plugin architecture)? [closed]
  30. XOR of three values

Leave a Comment