141 SUBROUTINE zsyt01_3( UPLO, N, A, LDA, AFAC, LDAFAC, E, IPIV, C,
142 $ ldc, rwork, resid )
151 INTEGER LDA, LDAFAC, LDC, N
152 DOUBLE PRECISION RESID
156 DOUBLE PRECISION RWORK( * )
157 COMPLEX*16 A( lda, * ), AFAC( ldafac, * ), C( ldc, * ),
164 DOUBLE PRECISION ZERO, ONE
165 parameter ( zero = 0.0d+0, one = 1.0d+0 )
166 COMPLEX*16 CZERO, CONE
167 parameter ( czero = ( 0.0d+0, 0.0d+0 ),
168 $ cone = ( 1.0d+0, 0.0d+0 ) )
172 DOUBLE PRECISION ANORM, EPS
176 DOUBLE PRECISION DLAMCH, ZLANSY
177 EXTERNAL lsame, dlamch, zlansy
196 CALL zsyconvf_rook( uplo,
'R', n, afac, ldafac, e, ipiv, info )
200 eps = dlamch(
'Epsilon' )
201 anorm = zlansy(
'1', uplo, n, a, lda, rwork )
205 CALL zlaset(
'Full', n, n, czero, cone, c, ldc )
209 CALL zlavsy_rook( uplo,
'Transpose',
'Non-unit', n, n, afac,
210 $ ldafac, ipiv, c, ldc, info )
214 CALL zlavsy_rook( uplo,
'No transpose',
'Unit', n, n, afac,
215 $ ldafac, ipiv, c, ldc, info )
219 IF( lsame( uplo,
'U' ) )
THEN
222 c( i, j ) = c( i, j ) - a( i, j )
228 c( i, j ) = c( i, j ) - a( i, j )
235 resid = zlansy(
'1', uplo, n, c, ldc, rwork )
237 IF( anorm.LE.zero )
THEN
241 resid = ( ( resid / dble( n ) ) / anorm ) / eps
247 CALL zsyconvf_rook( uplo,
'C', n, afac, ldafac, e, ipiv, info )
subroutine zlavsy_rook(UPLO, TRANS, DIAG, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
ZLAVSY_ROOK
subroutine zsyconvf_rook(UPLO, WAY, N, A, LDA, E, IPIV, INFO)
ZSYCONVF_ROOK
subroutine zsyt01_3(UPLO, N, A, LDA, AFAC, LDAFAC, E, IPIV, C, LDC, RWORK, RESID)
ZSYT01_3
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...