gcc/testsuite/gfortran.dg/vect/pr51058.f90

   1 ! { dg-do compile }
   2
   3       SUBROUTINE MLIST(MOLsp,PBCx,PBCy,PBCz, X0)
   4       IMPLICIT NONE
   5       INTEGER, PARAMETER :: NM=16384
   6       INTEGER :: MOLsp, i
   7       REAL :: PBCx, PBCy, PBCz, boxjmp, HALf=1./2.
   8       REAL :: X0(2,-2:NM)
   9
  10          DO i = 1 , MOLsp
  11             boxjmp = PBCx*INT(X0(1,i)+SIGN(HALf,X0(1,i)))
  12             X0(1,i) = X0(1,i) - boxjmp
  13             boxjmp = PBCy*INT(X0(2,i)+SIGN(HALf,X0(2,i)))
  14             X0(2,i) = X0(2,i) - boxjmp
  15          ENDDO
  16       END
  17
  18