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
LOGICAL function lsame(CA, CB)
LSAME
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...
subroutine sspcon(UPLO, N, AP, IPIV, ANORM, RCOND, WORK, IWORK, INFO)
SSPCON
subroutine ssptrs(UPLO, N, NRHS, AP, IPIV, B, LDB, INFO)
SSPTRS