fix incorrect rounding in printf floating point corner cases
authorRich Felker <dalias@aerifal.cx>
Sun, 9 Mar 2014 07:09:49 +0000 (03:09 -0400)
committerRich Felker <dalias@aerifal.cx>
Sun, 9 Mar 2014 07:09:49 +0000 (03:09 -0400)
commit9743a399bf4d6da9a1dbdf7e8df07284c97df16f
treec9477850f2abba3b76123e1a2299b379199be9c4
parentba231cf9e5923b6216081e9a626465c6643ce4d3
fix incorrect rounding in printf floating point corner cases

the printf floating point formatting code contains an optimization to
avoid computing digits that will be thrown away by rounding at the
specified (or default) precision. while it was correctly retaining all
places up to the last decimal place to be printed, it was not
retaining enough precision to see the next nonzero decimal place in
all cases. this could cause incorrect rounding down in round-to-even
(default) rounding mode, for example, when printing 0.5+DBL_EPSILON
with "%.0f".

in the fix, LDBL_MANT_DIG/3 is a lazy (non-sharp) upper bound on the
number of zeros between any two nonzero decimal digits.
src/stdio/vfprintf.c