[Bug libfortran/48511] Implement Steele-White algorithm for numeric output

[burnus at gcc dot gnu.org]

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

--- Comment #8 from Tobias Burnus <burnus at gcc dot gnu.org> 2011-04-27 14:59:39 UTC ---
(In reply to comment #7)
> Gay's routines (http://www.netlib.org/fp/) can handle double, float, extended,
> quad

libquadmath initially used dtoa for __float128 [REAL(16)] I/O, but there were
several issues especially with NaN; thus, one now uses functions based on
GLIBC. (At least the code I saw in libquadmath/gdtoa was also a mess, maybe
using "dtoa.c" itself is better.)

In any case, one could think of using Fortran tailored I/O functions instead of
using the C ones; and Gay's routines might be a starting point. (Interestingly,
the netlib files are still maintained, the last changelog entry is from
2011-03-21!)