! { dg-do compile } | |
SUBROUTINE PROPAGATE(ICI1,ICI2,I,J,J1,ELEM,NHSO,HSO | |
* ,MULST,IROOTS) | |
IMPLICIT DOUBLE PRECISION(A-H,O-Z) | |
COMPLEX*16 HSO,ELEM | |
DIMENSION HSO(NHSO,NHSO),MULST(*),IROOTS(*) | |
ISHIFT=MULST(ICI1)*(I-1)+1 | |
JSHIFT=MULST(ICI2)*(J-1)+1 | |
DO 200 ICI=1,ICI1-1 | |
ISHIFT=ISHIFT+MULST(ICI)*IROOTS(ICI) | |
200 CONTINUE | |
DO 220 ICI=1,ICI2-1 | |
JSHIFT=JSHIFT+MULST(ICI)*IROOTS(ICI) | |
220 CONTINUE | |
DO 150 MSS=MS,-MS,-2 | |
IND1=ISHIFT+K | |
IND2=JSHIFT+K | |
HSO(IND1,IND2)=ELEM | |
HSO(IND2,IND1)=DCONJG(ELEM) | |
150 CONTINUE | |
END | |