How scanf() handles floating point precision formats ??

Question

Below is an excerpt from a program -

float val=214.20;
double val1=214.20;

printf("\n\n\nfloat :: %f,  %4.6f, %4.2f \n ",val,val,val);
printf("double:: %f,  %4.6f, %4.2f \n ",val1,val1,val1);

And the output is -

float :: 214.199997,  214.199997, 214.20<--- is the correct value as we wanted 
double:: 214.200000,  214.200000, 214.20

I understand that 214.20 has infinite binary representation. The first two elements of first line has approximation of the intended value, but the last value has absolutely same data. This caused a question in my mind -

How the scanf/fscanf/sprintf functions treat the precision formats ??

With no precision provided it printed an approximate value, but with %4.2f it gave correct result. Can some friend may explain me the algorithm used by scanf when handling precisions.

Thanks, Ripunjay

Answer 1

The thing is, 214.20 cannot be expressed exactly with binary representation. Few decimal numbers can. So an approximation is stored. Now when you use printf, the binary representation is turned into a decimal representation, but it again cannot be expressed exactly and is only approximated.

As you noticed, you can give a precision to printf to tell it how to round the decimal approximation. And if you don't give it a precision then a precision of 6 is assumed (see the man page for details).

If you use %.40f for the float and %.40lf for the double in your example above, you will get these results:

214.1999969482421875000000000000000000000000
214.1999999999999886313162278383970260620117

They are different because with double, there are more bits to better approximate 214.20. But as you can see, they are still very odd when represented in decimal.

I recommend to read the Wikipedia article on floating point numbers for more insights about how floating point numbers work. An excellent read is also What Every Computer Scientist Should Know About Floating-Point Arithmetic

Answer 2

scanf will round the input value to the nearest exactly representable floating point value. As DarkDust's answer illustrates, in single precision the closest exactly representable value is below the exact value, and in double precision the closest exactly representable value is above the exact value.

Answer 3

Since you asked about scanf, one thing you should note is that POSIX requires printf and a subsequent scanf (or strtod) to reconstruct the original value exactly as long as sufficiently significant digits (at least DECIMAL_DIG, I believe) were printed. Plain C of course makes no such requirement; the C standard pretty much allows floating point operations to give whatever result the implementor likes as long as they document it. However, if your intent is to store floating point numbers in text files to be read back later, you might be better off using the C99 %a specifier to print them in hex. This way they'll be exact and there's no confusion about whether the serialize/deserialize process loses precision.

Keep in mind thatr when I said "reconstruct the original value", I mean the actual value which was held in the variable/expression passed to printf, not the original decimal you wrote in the source file which got rounded to the best binary representation by the compiler.