https://gcc.gnu.org/bugzilla/show_bug.cgi?id=67419
Francois-Xavier Coudert <fxcoudert at gcc dot gnu.org> changed: What |Removed |Added ---------------------------------------------------------------------------- Status|UNCONFIRMED |NEW Last reconfirmed| |2015-09-01 Ever confirmed|0 |1 --- Comment #1 from Francois-Xavier Coudert <fxcoudert at gcc dot gnu.org> --- Can you run the following test and paste the output here? integer,parameter :: k = selected_real_kind (precision (0.0_8) + 1) real(kind=k) :: x real(kind=8) :: y print *, k x = 789.1356 y = 789.1356 print *, x print *, y print *, sqrt(x) print *, sqrt(y) print *, abs((sqrt(y)-sqrt(x))/sqrt(x)) end This should call sqrtl() and sqrt() on the same value, and check that the results are close enough to one another. But the threshold of the test case is really generous…