Comparing IEEE floats and doubles for equality

by

The best approach I think is to compare ULPs.

bool is_nan(float f)
{
    return (*reinterpret_cast<unsigned __int32*>(&f) & 0x7f800000) == 0x7f800000 && (*reinterpret_cast<unsigned __int32*>(&f) & 0x007fffff) != 0;
}

bool is_finite(float f)
{
    return (*reinterpret_cast<unsigned __int32*>(&f) & 0x7f800000) != 0x7f800000;
}

// if this symbol is defined, NaNs are never equal to anything (as is normal in IEEE floating point)
// if this symbol is not defined, NaNs are hugely different from regular numbers, but might be equal to each other
#define UNEQUAL_NANS 1
// if this symbol is defined, infinites are never equal to finite numbers (as they're unimaginably greater)
// if this symbol is not defined, infinities are 1 ULP away from +/- FLT_MAX
#define INFINITE_INFINITIES 1

// test whether two IEEE floats are within a specified number of representable values of each other
// This depends on the fact that IEEE floats are properly ordered when treated as signed magnitude integers
bool equal_float(float lhs, float rhs, unsigned __int32 max_ulp_difference)
{
#ifdef UNEQUAL_NANS
    if(is_nan(lhs) || is_nan(rhs))
    {
        return false;
    }
#endif
#ifdef INFINITE_INFINITIES
    if((is_finite(lhs) && !is_finite(rhs)) || (!is_finite(lhs) && is_finite(rhs)))
    {
        return false;
    }
#endif
    signed __int32 left(*reinterpret_cast<signed __int32*>(&lhs));
    // transform signed magnitude ints into 2s complement signed ints
    if(left < 0)
    {
        left = 0x80000000 - left;
    }
    signed __int32 right(*reinterpret_cast<signed __int32*>(&rhs));
    // transform signed magnitude ints into 2s complement signed ints
    if(right < 0)
    {
        right = 0x80000000 - right;
    }
    if(static_cast<unsigned __int32>(std::abs(left - right)) <= max_ulp_difference)
    {
        return true;
    }
    return false;
}

A similar technique can be used for doubles. The trick is to convert the floats so that they’re ordered (as if integers) and then just see how different they are.

I have no idea why this damn thing is screwing up my underscores. Edit: Oh, perhaps that is just an artefact of the preview. That’s OK then.

More Related Contents:

  1. In MATLAB, are variables REALLY double-precision by default?
  2. How to get variable name?
  3. Creating or referencing variables dynamically in Sass
  4. What’s the difference between placing “mut” before a variable name and after the “:”?
  5. When should an Excel VBA variable be killed or set to Nothing?
  6. Formatting a float to 2 decimal places
  7. Why are variables “i” and “j” used for counters?
  8. Which variable size to use (db, dw, dd) with x86 assembly?
  9. When should I quote CMake variables?
  10. HiveQL: Using query results as variables
  11. In Rust, what’s the difference between “shadowing” and “mutability”?
  12. setq and defvar in Lisp
  13. What is the difference between a member variable and a local variable?
  14. How to make points one color when a third column equals zero, and another color otherwise, in Gnuplot?
  15. Batch File Set Variable not working [duplicate]
  16. Determine variable names dynamically according to a string in Fortran
  17. Why there are two ways of declaring variables in Go, what’s the difference and which to use?
  18. Multiple variable assignment in Swift
  19. Why are variables “i” and “j” used for counters?
  20. Dynamic Variable Naming and Reference (ColdFusion)
  21. Windows Batch files: what is variable expansion, and what does EnableDelayedExpansion mean?
  22. What is the difference between a constant and a static variable and which should I choose?
  23. gnuplot store one number from data file into variable
  24. What is the order of local variables on the stack?
  25. What is “_,” (underscore comma) in a Go declaration?
  26. Declaring global variable with php.ini
  27. Why does Go allow compilation of unused function parameters?
  28. Create a variable in swift with dynamic name
  29. Finding out reason of Pyomo model infeasibility
  30. When should I quote CMake variable references?

Leave a Comment