138 REAL FUNCTION cla_syrcond_c( UPLO, N, A, LDA, AF, LDAF, IPIV, C,
139 $ capply, info, work, rwork )
149 INTEGER N, LDA, LDAF, INFO
153 COMPLEX A( lda, * ), AF( ldaf, * ), WORK( * )
154 REAL C( * ), RWORK( * )
161 REAL AINVNM, ANORM, TMP
183 cabs1( zdum ) = abs(
REAL( ZDUM ) ) + abs( AIMAG( zdum ) )
190 upper = lsame( uplo,
'U' )
191 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
193 ELSE IF( n.LT.0 )
THEN
195 ELSE IF( lda.LT.max( 1, n ) )
THEN
197 ELSE IF( ldaf.LT.max( 1, n ) )
THEN
201 CALL xerbla(
'CLA_SYRCOND_C', -info )
205 IF ( lsame( uplo,
'U' ) ) up = .true.
215 tmp = tmp + cabs1( a( j, i ) ) / c( j )
218 tmp = tmp + cabs1( a( i, j ) ) / c( j )
222 tmp = tmp + cabs1( a( j, i ) )
225 tmp = tmp + cabs1( a( i, j ) )
229 anorm = max( anorm, tmp )
236 tmp = tmp + cabs1( a( i, j ) ) / c( j )
239 tmp = tmp + cabs1( a( j, i ) ) / c( j )
243 tmp = tmp + cabs1( a( i, j ) )
246 tmp = tmp + cabs1( a( j, i ) )
250 anorm = max( anorm, tmp )
259 ELSE IF( anorm .EQ. 0.0e+0 )
THEN
269 CALL clacn2( n, work( n+1 ), work, ainvnm, kase, isave )
276 work( i ) = work( i ) * rwork( i )
280 CALL csytrs(
'U', n, 1, af, ldaf, ipiv,
283 CALL csytrs(
'L', n, 1, af, ldaf, ipiv,
291 work( i ) = work( i ) * c( i )
300 work( i ) = work( i ) * c( i )
305 CALL csytrs(
'U', n, 1, af, ldaf, ipiv,
308 CALL csytrs(
'L', n, 1, af, ldaf, ipiv,
315 work( i ) = work( i ) * rwork( i )
323 IF( ainvnm .NE. 0.0e+0 )
subroutine csytrs(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
CSYTRS
real function cla_syrcond_c(UPLO, N, A, LDA, AF, LDAF, IPIV, C, CAPPLY, INFO, WORK, RWORK)
CLA_SYRCOND_C computes the infinity norm condition number of op(A)*inv(diag(c)) for symmetric indefin...
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine clacn2(N, V, X, EST, KASE, ISAVE)
CLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...