130 SUBROUTINE ssycon( UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK,
144 INTEGER ipiv( * ), iwork( * )
145 REAL a( lda, * ), work( * )
152 parameter( one = 1.0e+0, zero = 0.0e+0 )
177 upper =
lsame( uplo,
'U' )
178 IF( .NOT.upper .AND. .NOT.
lsame( uplo,
'L' ) )
THEN
180 ELSE IF( n.LT.0 )
THEN
182 ELSE IF( lda.LT.max( 1, n ) )
THEN
184 ELSE IF( anorm.LT.zero )
THEN
188 CALL
xerbla(
'SSYCON', -info )
198 ELSE IF( anorm.LE.zero )
THEN
209 IF( ipiv( i ).GT.0 .AND. a( i, i ).EQ.zero )
217 IF( ipiv( i ).GT.0 .AND. a( i, i ).EQ.zero )
226 CALL
slacn2( n, work( n+1 ), work, iwork, ainvnm, kase, isave )
231 CALL
ssytrs( uplo, n, 1, a, lda, ipiv, work, n, info )
238 $ rcond = ( one / ainvnm ) / anorm
LOGICAL function lsame(CA, CB)
LSAME
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine ssytrs(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
SSYTRS
subroutine slacn2(N, V, X, ISGN, EST, KASE, ISAVE)
SLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...
subroutine ssycon(UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK, IWORK, INFO)
SSYCON