3 PARAMETER(MAXDIM=4,MAXVEC=4**3*8,MAXT=20)
4 REAL(KIND=8) B,W1,W2,BNORM,BINV,WT,W0,C1,C2,R1,R2
5 DIMENSION B(MAXVEC,0:3),W1(MAXVEC,0:3),W2(MAXVEC,0:3)
6 DIMENSION BNORM(MAXVEC),BINV(MAXVEC),WT(MAXVEC),W0(MAXVEC)
7 DIMENSION C1(MAXVEC),C2(MAXVEC),R1(MAXVEC),R2(MAXVEC)
9 DO 200 I1=1,IDIM ! { dg-warning "Obsolescent feature: Shared DO termination" }
11 CALL INTACT(ILAT,I1,I1,W1)
15 WT(I)=ONE + C1(I)*LOG(EPS+R1(I))
16 IF( R2(I)**2 .LE. (ONE-WT(I)**2) )THEN