1! PR rtl-optimization/82973
2! { dg-do compile }
3! { dg-options "-Ofast -frounding-math" }
4
5program pr82973
6  integer, parameter :: n=16
7  real, dimension(n) :: ar, br, modulo_result, floor_result
8  integer, dimension(n) :: ai, bi , imodulo_result, ifloor_result
9  ai(1:4) = 5
10  ai(5:8) = -5
11  ai(9:12) = 1
12  ai(13:16) = -1
13  bi(1:4) = (/ 3,-3, 1, -1/)
14  bi(5:8) = bi(1:4)
15  bi(9:12) = bi(1:4)
16  bi(13:16) = bi(1:4)
17  ar = ai
18  br = bi
19  modulo_result = modulo(ar,br)
20  imodulo_result = modulo(ai,bi)
21  floor_result = ar-floor(ar/br)*br
22  ifloor_result = nint(real(ai-floor(real(ai)/real(bi))*bi))
23  do i=1,n
24    if (modulo_result(i) /= floor_result(i)) then
25      STOP 1
26    end if
27    if (imodulo_result(i) /= ifloor_result(i)) then
28      STOP 2
29    end if
30  end do
31end program pr82973
32