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 csytrs(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
CSYTRS
subroutine xerbla(SRNAME, INFO)
XERBLA
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...