91 DOUBLE PRECISION FUNCTION ztzt02( M, N, AF, LDA, TAU, WORK,
100 INTEGER LDA, LWORK, M, N
103 COMPLEX*16 AF( lda, * ), TAU( * ), WORK( lwork )
109 DOUBLE PRECISION ZERO, ONE
110 parameter( zero = 0.0d0, one = 1.0d0 )
116 DOUBLE PRECISION RWORK( 1 )
119 DOUBLE PRECISION DLAMCH, ZLANGE
120 EXTERNAL dlamch, zlange
126 INTRINSIC dble, dcmplx, dconjg, max
132 IF( lwork.LT.n*n+n )
THEN
133 CALL xerbla(
'ZTZT02', 7 )
139 IF( m.LE.0 .OR. n.LE.0 )
144 CALL zlaset(
'Full', n, n, dcmplx( zero ), dcmplx( one ), work,
150 CALL zlatzm(
'Left', n-m+1, n, af( i, m+1 ), lda, tau( i ),
151 $ work( i ), work( m+1 ), n, work( n*n+1 ) )
157 CALL zlatzm(
'Left', n-m+1, n, af( i, m+1 ), lda,
158 $ dconjg( tau( i ) ), work( i ), work( m+1 ), n,
165 work( ( i-1 )*n+i ) = work( ( i-1 )*n+i ) - one
168 ztzt02 = zlange(
'One-norm', n, n, work, n, rwork ) /
169 $ ( dlamch(
'Epsilon' )*dble( max( m, n ) ) )
subroutine zlatzm(SIDE, M, N, V, INCV, TAU, C1, C2, LDC, WORK)
ZLATZM
subroutine xerbla(SRNAME, INFO)
XERBLA
double precision function ztzt02(M, N, AF, LDA, TAU, WORK, LWORK)
ZTZT02
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...