1! { dg-do run }
2! testcase from PR 19032 adapted for testsuite
3! Our implementation of modulo was wrong for P = 1 and P = -1,
4! both in the real and the integer case
5program main
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
10  ai(1:4) = 5
11  ai(5:8) = -5
12  ai(9:12) = 1
13  ai(13:16) = -1
14  bi(1:4) = (/ 3,-3, 1, -1/)
15  bi(5:8) = bi(1:4)
16  bi(9:12) = bi(1:4)
17  bi(13:16) = bi(1:4)
18  ar = ai
19  br = bi
20  modulo_result = modulo(ar,br)
21  imodulo_result = modulo(ai,bi)
22  floor_result = ar-floor(ar/br)*br
23  ifloor_result = nint(real(ai-floor(real(ai)/real(bi))*bi))
24
25  do i=1,n
26     if (modulo_result(i) /= floor_result(i) ) then
27!        print "(A,4F5.0)" ,"real case failed: ", &
28!             ar(i),br(i), modulo_result(i), floor_result(i)
29        call abort()
30    end if
31    if (imodulo_result(i) /= ifloor_result(i)) then
32!       print "(A,4I5)", "int case failed: ", &
33!            ai(i), bi(i), imodulo_result(i), ifloor_result(i)
34       call abort ()
35     end if
36  end do
37end program main
38