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) 11220 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 18311 CONTINUE 19310 CONTINUE 20200 CONTINUE 21 END 22