1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
! { 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
! { dg-final { cleanup-tree-dump "vect" } }
|
