3 SUBROUTINE PROPAGATE
(ICI1
,ICI2
,I
,J
,J1
,ELEM
,NHSO
,HSO
5 IMPLICIT DOUBLE PRECISION(A
-H
,O
-Z
)
7 DIMENSION HSO
(NHSO
,NHSO
),MULST
(*),IROOTS
(*)
8 ISHIFT
=MULST
(ICI1
)*(I
-1)+1
9 JSHIFT
=MULST
(ICI2
)*(J
-1)+1
11 ISHIFT
=ISHIFT
+MULST
(ICI
)*IROOTS
(ICI
)
14 JSHIFT
=JSHIFT
+MULST
(ICI
)*IROOTS
(ICI
)
20 HSO
(IND2
,IND1
)=DCONJG
(ELEM
)