Double precision – decimal places

by

An IEEE double has 53 significant bits (that’s the value of DBL_MANT_DIG in <cfloat>). That’s approximately 15.95 decimal digits (log10(253)); the implementation sets DBL_DIG to 15, not 16, because it has to round down. So you have nearly an extra decimal digit of precision (beyond what’s implied by DBL_DIG==15) because of that.

The nextafter() function computes the nearest representable number to a given number; it can be used to show just how precise a given number is.

This program:

#include <cstdio>
#include <cfloat>
#include <cmath>

int main() {
    double x = 1.0/7.0;
    printf("FLT_RADIX = %d\n", FLT_RADIX);
    printf("DBL_DIG = %d\n", DBL_DIG);
    printf("DBL_MANT_DIG = %d\n", DBL_MANT_DIG);
    printf("%.17g\n%.17g\n%.17g\n", nextafter(x, 0.0), x, nextafter(x, 1.0));
}

gives me this output on my system:

FLT_RADIX = 2
DBL_DIG = 15
DBL_MANT_DIG = 53
0.14285714285714282
0.14285714285714285
0.14285714285714288

(You can replace %.17g by, say, %.64g to see more digits, none of which are significant.)

As you can see, the last displayed decimal digit changes by 3 with each consecutive value. The fact that the last displayed digit of 1.0/7.0 (5) happens to match the mathematical value is largely coincidental; it was a lucky guess. And the correct rounded digit is 6, not 5. Replacing 1.0/7.0 by 1.0/3.0 gives this output:

FLT_RADIX = 2
DBL_DIG = 15
DBL_MANT_DIG = 53
0.33333333333333326
0.33333333333333331
0.33333333333333337

which shows about 16 decimal digits of precision, as you’d expect.

More Related Contents:

  1. What is the difference between float and double?
  2. Are all integer values perfectly represented as doubles? [duplicate]
  3. Matlab vs C++ Double Precision
  4. How do I print a double value with full precision using cout?
  5. Is floating-point math consistent in C#? Can it be?
  6. C++ floating point precision [duplicate]
  7. What is a subnormal floating point number?
  8. 32-bit to 16-bit Floating Point Conversion
  9. Dividing two integers to produce a float result [duplicate]
  10. How to check if C++ compiler uses IEEE 754 floating point standard
  11. In which order should floats be added to get the most precise result?
  12. Emulate “double” using 2 “float”s
  13. How to ‘cout’ the correct number of decimal places of a double value?
  14. Portability of binary serialization of double/float type in C++
  15. What is the precision of long double in C++?
  16. Representing integers in doubles
  17. Issue when scheduling tasks using clock() function
  18. Floating point comparison revisited
  19. Set back default floating point print precision in C++
  20. How computer does floating point arithmetic?
  21. What is the purpose of max_digits10 and how is it different from digits10?
  22. Uses for negative zero floating point value?
  23. Is there a functional difference between “2.00” and “2.00f”?
  24. Are IEEE floats valid key types for std::map and std::set?
  25. Exact decimal datatype for C++?
  26. Does the C++ standard specify anything on the representation of floating point numbers?
  27. Show two digits after decimal point in c++ [duplicate]
  28. Why is unsigned integer overflow defined behavior but signed integer overflow isn’t?
  29. C++ ABI issues list
  30. How to ensure constexpr function never called at runtime?

Leave a Comment