blob: 3fa5d99cbf2363089059b95daaa8b77850ec04df (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
! { dg-do run }
! { dg-require-effective-target fortran_integer_16 }
! { dg-require-effective-target fortran_real_16 }
! { dg-additional-options "-fdump-tree-original" }
! { dg-final { scan-tree-dump-times "_gfortran_stop_numeric" 2 "original" } }
!
! PR fortran/96711 - ICE on NINT() Function
program p
implicit none
real(8) :: x
real(16) :: y
! Assume radix(x) == 2
! 2/epsilon(x) = 2/(radix(x)**(1-digits(x)) = 2**digits(x) with that assumption
integer(16), parameter :: k1 = nint (2 / epsilon (x), kind(k1))
integer(16), parameter :: k2 = nint (2 / epsilon (y), kind(k2))
integer(16), parameter :: m1 = 2_16**digits(x) ! IEEE: 2**53
integer(16), parameter :: m2 = 2_16**digits(y) ! IEEE: 2**113
integer(16), volatile :: m
x = 2 / epsilon (x)
y = 2 / epsilon (y)
m = nint (x, kind(m))
! print *, m
if (k1 /= m1) stop 1
if (m /= m1) stop 2
m = nint (y, kind(m))
! print *, m
if (k2 /= m2) stop 3
if (m /= m2) stop 4
end program
|
