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!)
