138 SUBROUTINE dsyt01_3( UPLO, N, A, LDA, AFAC, LDAFAC, E, IPIV, C,
139 $ LDC, RWORK, RESID )
147 INTEGER LDA, LDAFAC, LDC, N
148 DOUBLE PRECISION RESID
152 DOUBLE PRECISION A( LDA, * ), AFAC( LDAFAC, * ), C( LDC, * ),
159 DOUBLE PRECISION ZERO, ONE
160 parameter( zero = 0.0d+0, one = 1.0d+0 )
164 DOUBLE PRECISION ANORM, EPS
168 DOUBLE PRECISION DLAMCH, DLANSY
169 EXTERNAL lsame, dlamch, dlansy
188 CALL dsyconvf_rook( uplo,
'R', n, afac, ldafac, e, ipiv, info )
192 eps = dlamch(
'Epsilon' )
193 anorm = dlansy(
'1', uplo, n, a, lda, rwork )
197 CALL dlaset(
'Full', n, n, zero, one, c, ldc )
201 CALL dlavsy_rook( uplo,
'Transpose',
'Non-unit', n, n, afac,
202 $ ldafac, ipiv, c, ldc, info )
206 CALL dlavsy_rook( uplo,
'No transpose',
'Unit', n, n, afac,
207 $ ldafac, ipiv, c, ldc, info )
211 IF( lsame( uplo,
'U' ) )
THEN
214 c( i, j ) = c( i, j ) - a( i, j )
220 c( i, j ) = c( i, j ) - a( i, j )
227 resid = dlansy(
'1', uplo, n, c, ldc, rwork )
229 IF( anorm.LE.zero )
THEN
233 resid = ( ( resid / dble( n ) ) / anorm ) / eps
239 CALL dsyconvf_rook( uplo,
'C', n, afac, ldafac, e, ipiv, info )
subroutine dlaset(UPLO, M, N, ALPHA, BETA, A, LDA)
DLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values.
subroutine dlavsy_rook(UPLO, TRANS, DIAG, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
DLAVSY_ROOK
subroutine dsyt01_3(UPLO, N, A, LDA, AFAC, LDAFAC, E, IPIV, C, LDC, RWORK, RESID)
DSYT01_3
subroutine dsyconvf_rook(UPLO, WAY, N, A, LDA, E, IPIV, INFO)
DSYCONVF_ROOK