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…

Reply via email to