Why is this C program giving a 'wrong' output?

Question

Why is this C program giving a 'wrong' output?

#include<stdio.h>

void main()
{
    float f = 12345.054321;

    printf("%f", f);

    getch();
}

Output:

12345.054688

But the output should be, 12345.054321.

I am using VC++ in VS2008.

Answer 1

The decimal number 12345.054321 cannot be represented accurately as a float on your platform. The result that you are seeing is a decimal approximation to the closest number that can be represented as a float.

Answer 2

It's giving the "wrong" answer simply because not all real values are representable by floats (or doubles, for that matter). What you'll get is an approximation based on the underlying encoding.

In order to represent every real value, even between 1.0x10^-100 and 1.1x10^-100 (a truly minuscule range), you still require an infinite number of bits.

Single-precision IEEE754 values have only 32 bits available (some of which are tasked to other things such as exponent and NaN/Inf representations) and cannot therefore give you infinite precision. They actually have 23 bits available giving precision of about 2²⁴ (there's an extra implicit bit) or just over 7 decimal digits (log₁₀(2²⁴) is roughly 7.2).

I enclose the word "wrong" in quotes because it's not actually wrong. What's wrong is your understanding about how computers represent numbers (don't be offended though, you're not alone in this misapprehension).

Head on over to http://www.h-schmidt.net/FloatApplet/IEEE754.html and type your number into the "Decimal representation" box to see this in action.

If you want a more accurate number, use doubles instead of floats - these have double the number of bits available for representing values (assuming your C implementation is using IEEE754 single and double precision data types for float and double respectively).

If you want arbitrary precision, you'll need to use a "bignum" library like GMP although that's somewhat slower than native types so make sure you understand the trade-offs.

Answer 3

Single-precision floating point values can only represent about eight to nine significant (decimal) digits. Beyond that point, you're seeing quantization error.

Answer 4

It's all to do with precision. Your number cannot be stored accurately in a float.

Answer 5

floats are about convenience and speed, and use a binary representation - if you care about precision use a decimal type.

To understand the problem, read What Every Computer Scientist Should Know About Floating-Point Arithmetic: http://docs.sun.com/source/806-3568/ncg_goldberg.html

For a solution, see the Decimal Arithmetic FAQ: http://speleotrove.com/decimal/decifaq.html