119 SUBROUTINE cspcon( UPLO, N, AP, IPIV, ANORM, RCOND, WORK, INFO )
133 COMPLEX AP( * ), WORK( * )
140 parameter ( one = 1.0e+0, zero = 0.0e+0 )
162 upper = lsame( uplo,
'U' )
163 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
165 ELSE IF( n.LT.0 )
THEN
167 ELSE IF( anorm.LT.zero )
THEN
171 CALL xerbla(
'CSPCON', -info )
181 ELSE IF( anorm.LE.zero )
THEN
193 IF( ipiv( i ).GT.0 .AND. ap( ip ).EQ.zero )
203 IF( ipiv( i ).GT.0 .AND. ap( ip ).EQ.zero )
213 CALL clacn2( n, work( n+1 ), work, ainvnm, kase, isave )
218 CALL csptrs( uplo, n, 1, ap, ipiv, work, n, info )
225 $ rcond = ( one / ainvnm ) / anorm
subroutine cspcon(UPLO, N, AP, IPIV, ANORM, RCOND, WORK, INFO)
CSPCON
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine csptrs(UPLO, N, NRHS, AP, IPIV, B, LDB, INFO)
CSPTRS
subroutine clacn2(N, V, X, EST, KASE, ISAVE)
CLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...