125 SUBROUTINE csycon( UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK,
140 COMPLEX A( lda, * ), WORK( * )
147 parameter( one = 1.0e+0, zero = 0.0e+0 )
172 upper = lsame( uplo,
'U' )
173 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
175 ELSE IF( n.LT.0 )
THEN
177 ELSE IF( lda.LT.max( 1, n ) )
THEN
179 ELSE IF( anorm.LT.zero )
THEN
183 CALL
xerbla(
'CSYCON', -info )
193 ELSE IF( anorm.LE.zero )
THEN
204 IF( ipiv( i ).GT.0 .AND. a( i, i ).EQ.zero )
212 IF( ipiv( i ).GT.0 .AND. a( i, i ).EQ.zero )
221 CALL
clacn2( n, work( n+1 ), work, ainvnm, kase, isave )
226 CALL
csytrs( uplo, n, 1, a, lda, ipiv, work, n, info )
233 $ rcond = ( one / ainvnm ) / anorm
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine csytrs(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
CSYTRS
subroutine csycon(UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK, INFO)
CSYCON
subroutine clacn2(N, V, X, EST, KASE, ISAVE)
CLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...