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