111
|
1 ! { dg-do compile }
|
|
2 ! { dg-options "-Ofast -fwrapv" }
|
|
3 ! { dg-additional-options "-march=broadwell" { target x86_64-*-* i?86-*-* } }
|
|
4 SUBROUTINE ECPDRA(IC4C,FP,FQ,G)
|
|
5 IMPLICIT DOUBLE PRECISION (A-H,O-Z)
|
|
6 DIMENSION FP(*),FQ(*),G(*)
|
|
7 DIMENSION CKLU(23,12,12)
|
|
8 !
|
|
9 DO 240 I=IAMIN,IAMAX
|
|
10 DO 240 J=JAMIN,MMAX
|
|
11 DO 230 K=1,NPNP
|
|
12 DO 230 L=1,K
|
|
13 DO 230 MU=1,2*L-1
|
|
14 CKLTEM= CKLU(MU,L,K)
|
|
15 IF(IC4C.LE.0) THEN
|
|
16 IF(ABS(CKLTEM).GT.TOL) SUM= SUM+FP(N)*CKLTEM
|
|
17 ELSE
|
|
18 IF(ABS(CKLTEM).GT.TOL) SUM= SUM+FQ(N)*CKLTEM
|
|
19 END IF
|
|
20 230 N= N+1
|
|
21 G(NN)= G(NN)+DUMJ*SUM
|
|
22 240 NN= NN+1
|
|
23 END
|