109 INTEGER I1, I2, LDA, N
112 COMPLEX*16 A( LDA, N )
131 upper = lsame( uplo,
'U' )
137 CALL zswap( i1-1, a(1,i1), 1, a(1,i2), 1 )
150 a(i1,i1+i)=dconjg(a(i1+i,i2))
151 a(i1+i,i2)=dconjg(tmp)
154 a(i1,i2)=dconjg(a(i1,i2))
170 CALL zswap ( i1-1, a(i1,1), lda, a(i2,1), lda )
183 a(i1+i,i1)=dconjg(a(i2,i1+i))
184 a(i2,i1+i)=dconjg(tmp)
187 a(i2,i1)=dconjg(a(i2,i1))
subroutine zswap(N, ZX, INCX, ZY, INCY)
ZSWAP
subroutine zheswapr(UPLO, N, A, LDA, I1, I2)
ZHESWAPR applies an elementary permutation on the rows and columns of a Hermitian matrix.