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