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