IEEE-correctly-rounded by mandate, so the intrinsic is already bit-identical across compliant compilers.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | x |
elemental function safe_sqrt(x) result(y) !! IEEE-correctly-rounded by mandate, so the intrinsic is !! already bit-identical across compliant compilers. real(wp), intent(in) :: x real(wp) :: y y = sqrt(x) end function safe_sqrt