123 REAL FUNCTION cla_porcond_x( UPLO, N, A, LDA, AF, LDAF, X, INFO,
133 INTEGER N, LDA, LDAF, INFO
136 COMPLEX A( lda, * ), AF( ldaf, * ), WORK( * ), X( * )
144 REAL AINVNM, ANORM, TMP
159 INTRINSIC abs, max,
REAL, AIMAG
165 cabs1( zdum ) = abs(
REAL( ZDUM ) ) + abs( AIMAG( zdum ) )
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( ldaf.LT.max( 1, n ) )
THEN
183 CALL xerbla(
'CLA_PORCOND_X', -info )
187 IF ( lsame( uplo,
'U' ) ) up = .true.
196 tmp = tmp + cabs1( a( j, i ) * x( j ) )
199 tmp = tmp + cabs1( a( i, j ) * x( j ) )
202 anorm = max( anorm, tmp )
208 tmp = tmp + cabs1( a( i, j ) * x( j ) )
211 tmp = tmp + cabs1( a( j, i ) * x( j ) )
214 anorm = max( anorm, tmp )
223 ELSE IF( anorm .EQ. 0.0e+0 )
THEN
233 CALL clacn2( n, work( n+1 ), work, ainvnm, kase, isave )
240 work( i ) = work( i ) * rwork( i )
244 CALL cpotrs(
'U', n, 1, af, ldaf,
247 CALL cpotrs(
'L', n, 1, af, ldaf,
254 work( i ) = work( i ) / x( i )
261 work( i ) = work( i ) / x( i )
265 CALL cpotrs(
'U', n, 1, af, ldaf,
268 CALL cpotrs(
'L', n, 1, af, ldaf,
275 work( i ) = work( i ) * rwork( i )
283 IF( ainvnm .NE. 0.0e+0 )
subroutine cpotrs(UPLO, N, NRHS, A, LDA, B, LDB, INFO)
CPOTRS
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...
real function cla_porcond_x(UPLO, N, A, LDA, AF, LDAF, X, INFO, WORK, RWORK)
CLA_PORCOND_X computes the infinity norm condition number of op(A)*diag(x) for Hermitian positive-def...