diff options
| author | Szabolcs Nagy <nsz@port70.net> | 2016-10-12 00:49:59 +0200 |
|---|---|---|
| committer | Rich Felker <dalias@aerifal.cx> | 2016-10-20 01:54:15 -0400 |
| commit | 51ab6db4ed115d079d7131975e4adb074ba9ef9d (patch) | |
| tree | 06791593d873bbe64a5c6c04ddf0420686e9a653 | |
| parent | 8fb28b0b3e7a5e958fb844722a4b2ef9bc244af1 (diff) | |
| download | musl-51ab6db4ed115d079d7131975e4adb074ba9ef9d.tar.gz musl-51ab6db4ed115d079d7131975e4adb074ba9ef9d.tar.xz musl-51ab6db4ed115d079d7131975e4adb074ba9ef9d.zip | |
fix float formatting of some exact halfway cases
in nearest rounding mode exact halfway cases were not following the
round to even rule if the rounding happened at a base 1000000000 digit
boundary of the internal representation and the previous digit was odd.
e.g. printf("%.0f", 1.5) printed 1 instead of 2.
| -rw-r--r-- | src/stdio/vfprintf.c | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/src/stdio/vfprintf.c b/src/stdio/vfprintf.c index e2ab2dc6..60da14d8 100644 --- a/src/stdio/vfprintf.c +++ b/src/stdio/vfprintf.c @@ -345,7 +345,8 @@ static int fmt_fp(FILE *f, long double y, int w, int p, int fl, int t) if (x || d+1!=z) { long double round = 2/LDBL_EPSILON; long double small; - if (*d/i & 1) round += 2; + if ((*d/i & 1) || (i==1000000000 && d>a && (d[-1]&1))) + round += 2; if (x<i/2) small=0x0.8p0; else if (x==i/2 && d+1==z) small=0x1.0p0; else small=0x1.8p0; |