106 SUBROUTINE zlakf2( M, N, A, LDA, B, D, E, Z, LDZ )
114 INTEGER LDA, LDZ, M, N
117 COMPLEX*16 A( lda, * ), B( lda, * ), D( lda, * ),
118 $ e( lda, * ), z( ldz, * )
125 parameter ( zero = ( 0.0d+0, 0.0d+0 ) )
128 INTEGER I, IK, J, JK, L, MN, MN2
139 CALL zlaset(
'Full', mn2, mn2, zero, zero, z, ldz )
148 z( ik+i-1, ik+j-1 ) = a( i, j )
156 z( ik+mn+i-1, ik+j-1 ) = d( i, j )
172 z( ik+i-1, jk+i-1 ) = -b( j, l )
178 z( ik+mn+i-1, jk+i-1 ) = -e( j, l )
subroutine zlakf2(M, N, A, LDA, B, D, E, Z, LDZ)
ZLAKF2
subroutine zlaset(UPLO, M, N, ALPHA, BETA, A, LDA)
ZLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...