gcc/testsuite/g77.f-torture/compile/20020307-1.f

   1       SUBROUTINE SWEEP
   2       PARAMETER(MAXDIM=4,MAXVEC=4**3*8,MAXT=20)
   3       REAL*8 B,W1,W2,BNORM,BINV,WT,W0,C1,C2,R1,R2
   4       DIMENSION B(MAXVEC,0:3),W1(MAXVEC,0:3),W2(MAXVEC,0:3)
   5       DIMENSION BNORM(MAXVEC),BINV(MAXVEC),WT(MAXVEC),W0(MAXVEC)
   6       DIMENSION C1(MAXVEC),C2(MAXVEC),R1(MAXVEC),R2(MAXVEC)
   7       DO 200 ILAT=1,2**IDIM
   8       DO 200 I1=1,IDIM
   9       DO 220 I2=1,IDIM
  10       CALL INTACT(ILAT,I1,I1,W1)
  11 220   CONTINUE
  12       DO 310 IATT=1,IDIM
  13       DO 311 I=1,100
  14       WT(I)=ONE + C1(I)*LOG(EPS+R1(I))
  15       IF( R2(I)**2 .LE. (ONE-WT(I)**2) )THEN
  16       W0(I)=WT(I)
  17       ENDIF
  18 311   CONTINUE
  19 310   CONTINUE
  20 200   CONTINUE
  21       END