https://gcc.gnu.org/bugzilla/show_bug.cgi?id=67419
--- Comment #6 from Francois-Xavier Coudert <fxcoudert at gcc dot gnu.org> --- (In reply to r...@cebitec.uni-bielefeld.de from comment #4) >> int main (void) >> { >> long double x; >> x = 789.1356; >> __builtin_printf ("%.30Lg\n", x); >> __builtin_printf ("%.30Lg\n", __builtin_sqrtl(x)); >> } > > I get > > 789.13559999999995397956809029 > -NaN OK, then it's definitely not Fortran-specific. The question is: is it just because of GCC's built-in handling (e.g., because the math.h header does something weird), or is it simply that the libm's sqrtl() is unusable?