140 SUBROUTINE ssyt01_3( UPLO, N, A, LDA, AFAC, LDAFAC, E, IPIV, C,
141 $ ldc, rwork, resid )
150 INTEGER LDA, LDAFAC, LDC, N
155 REAL A( lda, * ), AFAC( ldafac, * ), C( ldc, * ),
163 parameter ( zero = 0.0e+0, one = 1.0e+0 )
172 EXTERNAL lsame, slamch, slansy
191 CALL ssyconvf_rook( uplo,
'R', n, afac, ldafac, e, ipiv, info )
195 eps = slamch(
'Epsilon' )
196 anorm = slansy(
'1', uplo, n, a, lda, rwork )
200 CALL slaset(
'Full', n, n, zero, one, c, ldc )
204 CALL slavsy_rook( uplo,
'Transpose',
'Non-unit', n, n, afac,
205 $ ldafac, ipiv, c, ldc, info )
209 CALL slavsy_rook( uplo,
'No transpose',
'Unit', n, n, afac,
210 $ ldafac, ipiv, c, ldc, info )
214 IF( lsame( uplo,
'U' ) )
THEN
217 c( i, j ) = c( i, j ) - a( i, j )
223 c( i, j ) = c( i, j ) - a( i, j )
230 resid = slansy(
'1', uplo, n, c, ldc, rwork )
232 IF( anorm.LE.zero )
THEN
236 resid = ( ( resid /
REAL( N ) ) / anorm ) / eps
242 CALL ssyconvf_rook( uplo,
'C', n, afac, ldafac, e, ipiv, info )
subroutine ssyconvf_rook(UPLO, WAY, N, A, LDA, E, IPIV, INFO)
SSYCONVF_ROOK
subroutine slaset(UPLO, M, N, ALPHA, BETA, A, LDA)
SLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...
subroutine ssyt01_3(UPLO, N, A, LDA, AFAC, LDAFAC, E, IPIV, C, LDC, RWORK, RESID)
SSYT01_3
subroutine slavsy_rook(UPLO, TRANS, DIAG, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
SLAVSY_ROOK