141 SUBROUTINE csyt01_3( UPLO, N, A, LDA, AFAC, LDAFAC, E, IPIV, C,
142 $ ldc, rwork, resid )
151 INTEGER LDA, LDAFAC, LDC, N
157 COMPLEX A( lda, * ), AFAC( ldafac, * ), C( ldc, * ),
165 parameter ( zero = 0.0e+0, one = 1.0e+0 )
167 parameter ( czero = ( 0.0e+0, 0.0e+0 ),
168 $ cone = ( 1.0e+0, 0.0e+0 ) )
177 EXTERNAL lsame, slamch, clansy
196 CALL csyconvf_rook( uplo,
'R', n, afac, ldafac, e, ipiv, info )
200 eps = slamch(
'Epsilon' )
201 anorm = clansy(
'1', uplo, n, a, lda, rwork )
205 CALL claset(
'Full', n, n, czero, cone, c, ldc )
209 CALL clavsy_rook( uplo,
'Transpose',
'Non-unit', n, n, afac,
210 $ ldafac, ipiv, c, ldc, info )
214 CALL clavsy_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 = clansy(
'1', uplo, n, c, ldc, rwork )
237 IF( anorm.LE.zero )
THEN
241 resid = ( ( resid /
REAL( N ) ) / anorm ) / eps
247 CALL csyconvf_rook( uplo,
'C', n, afac, ldafac, e, ipiv, info )
subroutine csyconvf_rook(UPLO, WAY, N, A, LDA, E, IPIV, INFO)
CSYCONVF_ROOK
subroutine claset(UPLO, M, N, ALPHA, BETA, A, LDA)
CLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...
subroutine clavsy_rook(UPLO, TRANS, DIAG, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
CLAVSY_ROOK
subroutine csyt01_3(UPLO, N, A, LDA, AFAC, LDAFAC, E, IPIV, C, LDC, RWORK, RESID)
CSYT01_3