| ! { 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 |