125 SUBROUTINE sspcon( UPLO, N, AP, IPIV, ANORM, RCOND, WORK, IWORK,
139 INTEGER IPIV( * ), IWORK( * )
140 REAL AP( * ), WORK( * )
147 parameter ( one = 1.0e+0, zero = 0.0e+0 )
169 upper = lsame( uplo,
'U' )
170 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
172 ELSE IF( n.LT.0 )
THEN
174 ELSE IF( anorm.LT.zero )
THEN
178 CALL xerbla(
'SSPCON', -info )
188 ELSE IF( anorm.LE.zero )
THEN
200 IF( ipiv( i ).GT.0 .AND. ap( ip ).EQ.zero )
210 IF( ipiv( i ).GT.0 .AND. ap( ip ).EQ.zero )
220 CALL slacn2( n, work( n+1 ), work, iwork, ainvnm, kase, isave )
225 CALL ssptrs( uplo, n, 1, ap, ipiv, work, n, info )
232 $ rcond = ( one / ainvnm ) / anorm
subroutine ssptrs(UPLO, N, NRHS, AP, IPIV, B, LDB, INFO)
SSPTRS
subroutine sspcon(UPLO, N, AP, IPIV, ANORM, RCOND, WORK, IWORK, INFO)
SSPCON
subroutine xerbla(SRNAME, INFO)
XERBLA
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...