129 $ info, work, rwork )
137 INTEGER n, lda, ldaf, info
141 COMPLEX a( lda, * ), af( ldaf, * ), work( * ), x( * )
149 REAL ainvnm, anorm, tmp
170 cabs1( zdum ) = abs( real( zdum ) ) + abs( aimag( zdum ) )
177 upper =
lsame( uplo,
'U' )
178 IF( .NOT.upper .AND. .NOT.
lsame( uplo,
'L' ) )
THEN
180 ELSE IF ( n.LT.0 )
THEN
182 ELSE IF( lda.LT.max( 1, n ) )
THEN
184 ELSE IF( ldaf.LT.max( 1, n ) )
THEN
188 CALL xerbla(
'CLA_HERCOND_X', -info )
192 IF (
lsame( uplo,
'U' ) ) up = .true.
201 tmp = tmp + cabs1( a( j, i ) * x( j ) )
204 tmp = tmp + cabs1( a( i, j ) * x( j ) )
207 anorm = max( anorm, tmp )
213 tmp = tmp + cabs1( a( i, j ) * x( j ) )
216 tmp = tmp + cabs1( a( j, i ) * x( j ) )
219 anorm = max( anorm, tmp )
228 ELSE IF( anorm .EQ. 0.0e+0 )
THEN
238 CALL clacn2( n, work( n+1 ), work, ainvnm, kase, isave )
245 work( i ) = work( i ) * rwork( i )
249 CALL chetrs(
'U', n, 1, af, ldaf, ipiv,
252 CALL chetrs(
'L', n, 1, af, ldaf, ipiv,
259 work( i ) = work( i ) / x( i )
266 work( i ) = work( i ) / x( i )
270 CALL chetrs(
'U', n, 1, af, ldaf, ipiv,
273 CALL chetrs(
'L', n, 1, af, ldaf, ipiv,
280 work( i ) = work( i ) * rwork( i )
288 IF( ainvnm .NE. 0.0e+0 )
subroutine chetrs(uplo, n, nrhs, a, lda, ipiv, b, ldb, info)
CHETRS
real function cla_hercond_x(uplo, n, a, lda, af, ldaf, ipiv, x, info, work, rwork)
CLA_HERCOND_X computes the infinity norm condition number of op(A)*diag(x) for Hermitian indefinite m...
subroutine clacn2(n, v, x, est, kase, isave)
CLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...