c { dg-do compile }
|
c { dg-do compile }
|
SUBROUTINE SWEEP
|
SUBROUTINE SWEEP
|
PARAMETER(MAXDIM=4,MAXVEC=4**3*8,MAXT=20)
|
PARAMETER(MAXDIM=4,MAXVEC=4**3*8,MAXT=20)
|
REAL(KIND=8) B,W1,W2,BNORM,BINV,WT,W0,C1,C2,R1,R2
|
REAL(KIND=8) B,W1,W2,BNORM,BINV,WT,W0,C1,C2,R1,R2
|
DIMENSION B(MAXVEC,0:3),W1(MAXVEC,0:3),W2(MAXVEC,0:3)
|
DIMENSION B(MAXVEC,0:3),W1(MAXVEC,0:3),W2(MAXVEC,0:3)
|
DIMENSION BNORM(MAXVEC),BINV(MAXVEC),WT(MAXVEC),W0(MAXVEC)
|
DIMENSION BNORM(MAXVEC),BINV(MAXVEC),WT(MAXVEC),W0(MAXVEC)
|
DIMENSION C1(MAXVEC),C2(MAXVEC),R1(MAXVEC),R2(MAXVEC)
|
DIMENSION C1(MAXVEC),C2(MAXVEC),R1(MAXVEC),R2(MAXVEC)
|
DO 200 ILAT=1,2**IDIM
|
DO 200 ILAT=1,2**IDIM
|
DO 200 I1=1,IDIM
|
DO 200 I1=1,IDIM
|
DO 220 I2=1,IDIM
|
DO 220 I2=1,IDIM
|
CALL INTACT(ILAT,I1,I1,W1)
|
CALL INTACT(ILAT,I1,I1,W1)
|
220 CONTINUE
|
220 CONTINUE
|
DO 310 IATT=1,IDIM
|
DO 310 IATT=1,IDIM
|
DO 311 I=1,100
|
DO 311 I=1,100
|
WT(I)=ONE + C1(I)*LOG(EPS+R1(I))
|
WT(I)=ONE + C1(I)*LOG(EPS+R1(I))
|
IF( R2(I)**2 .LE. (ONE-WT(I)**2) )THEN
|
IF( R2(I)**2 .LE. (ONE-WT(I)**2) )THEN
|
W0(I)=WT(I)
|
W0(I)=WT(I)
|
ENDIF
|
ENDIF
|
311 CONTINUE
|
311 CONTINUE
|
310 CONTINUE
|
310 CONTINUE
|
200 CONTINUE
|
200 CONTINUE
|
END
|
END
|
|
|