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