[Bug libfortran/59836] [4.7/4.8/4.9 Regression] Wrong outputs with rounding formats for some values.

[dominiq at lps dot ens.fr]

http://gcc.gnu.org/bugzilla/show_bug.cgi?id=59836

--- Comment #3 from Dominique d'Humieres <dominiq at lps dot ens.fr> ---
> I will also test your patches and review and commit for you if you need. 
> Please let me know. 

Could you check in particular that I did not miss something when removing the
line

-      nzero--;

I have not followed in detail the use of nzero.

I don't have write access so I cannot do the commit, but I can do the packaging
and the submit if it helps (I have the FSF form signed). Feel free to use the
posted patches as you like.

> I knew you were close to the solution on this one. Great job!

Thanks! Actually I have tried

-  if (w > 0 && d == 0 && p == 0)
+  if (ft != FMT_F && nbefore == 0 && w > 0 && d == 0 && p == 0)

but gfortran.dg/round_3.f08 regressed and I did not catch the need for the
second part of the patch, hence the kludge in pr59774 comment 9.

I think the formats and values in comment 0, can be added to
gfortran.dg/round_3.f08
(as well as some tests for pr59774). I'ld like also to add some tests for RC
around 0.5.