123 $ LDAF, X, INFO, WORK,
132 INTEGER n, lda, ldaf, info
135 COMPLEX*16 a( lda, * ), af( ldaf, * ), work( * ), x( * )
136 DOUBLE PRECISION rwork( * )
143 DOUBLE PRECISION ainvnm, anorm, tmp
158 INTRINSIC abs, max, real, dimag
161 DOUBLE PRECISION cabs1
164 cabs1( zdum ) = abs( dble( zdum ) ) + abs( dimag( zdum ) )
171 upper =
lsame( uplo,
'U' )
172 IF( .NOT.upper .AND. .NOT.
lsame( uplo,
'L' ) )
THEN
174 ELSE IF ( n.LT.0 )
THEN
176 ELSE IF( lda.LT.max( 1, n ) )
THEN
178 ELSE IF( ldaf.LT.max( 1, n ) )
THEN
182 CALL xerbla(
'ZLA_PORCOND_X', -info )
186 IF (
lsame( uplo,
'U' ) ) up = .true.
195 tmp = tmp + cabs1( a( j, i ) * x( j ) )
198 tmp = tmp + cabs1( a( i, j ) * x( j ) )
201 anorm = max( anorm, tmp )
207 tmp = tmp + cabs1( a( i, j ) * x( j ) )
210 tmp = tmp + cabs1( a( j, i ) * x( j ) )
213 anorm = max( anorm, tmp )
222 ELSE IF( anorm .EQ. 0.0d+0 )
THEN
232 CALL zlacn2( n, work( n+1 ), work, ainvnm, kase, isave )
239 work( i ) = work( i ) * rwork( i )
243 CALL zpotrs(
'U', n, 1, af, ldaf,
246 CALL zpotrs(
'L', n, 1, af, ldaf,
253 work( i ) = work( i ) / x( i )
260 work( i ) = work( i ) / x( i )
264 CALL zpotrs(
'U', n, 1, af, ldaf,
267 CALL zpotrs(
'L', n, 1, af, ldaf,
274 work( i ) = work( i ) * rwork( i )
282 IF( ainvnm .NE. 0.0d+0 )
subroutine xerbla(SRNAME, INFO)
XERBLA
logical function lsame(CA, CB)
LSAME
subroutine zlacn2(N, V, X, EST, KASE, ISAVE)
ZLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...
subroutine zpotrs(UPLO, N, NRHS, A, LDA, B, LDB, INFO)
ZPOTRS
double precision function zla_porcond_x(UPLO, N, A, LDA, AF, LDAF, X, INFO, WORK, RWORK)
ZLA_PORCOND_X computes the infinity norm condition number of op(A)*diag(x) for Hermitian positive-def...