1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
! { dg-do compile }
! { dg-options "-O2 -fprefetch-loop-arrays -w" }
SUBROUTINE GHDSYM(IZ,IS,LMMAX,S,LMS,Y,L2M,DRL,NLAY2,K0,DCUT)!,
!
COMPLEX Y(L2M,L2M),H(33),S(LMS)
COMPLEX RU,CI,CZ,K0,FF,Z,Z1,Z2,Z3,ST
!
DO 140 KK=1,4
DO 130 L=1,L2M
L1=L*L-L
DO 120 M=1,L
IPM=L1+M
IMM=L1-M+2
S(IPM)=S(IPM)+Z3*Y(L,M)
IF (M.NE.1) S(IMM)=S(IMM)+Z3*Y(M-1,L)*CSGN
120 CONTINUE
130 CONTINUE
140 CONTINUE
END
|
