132 DOUBLE PRECISION FUNCTION zla_syrcond_x( UPLO, N, A, LDA, AF,
133 $ ldaf, ipiv, x, info,
143 INTEGER N, LDA, LDAF, INFO
147 COMPLEX*16 A( lda, * ), AF( ldaf, * ), WORK( * ), X( * )
148 DOUBLE PRECISION RWORK( * )
155 DOUBLE PRECISION AINVNM, ANORM, TMP
174 DOUBLE PRECISION CABS1
177 cabs1( zdum ) = abs( dble( zdum ) ) + abs( dimag( zdum ) )
184 upper = lsame( uplo,
'U' )
185 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
187 ELSE IF( n.LT.0 )
THEN
189 ELSE IF( lda.LT.max( 1, n ) )
THEN
191 ELSE IF( ldaf.LT.max( 1, n ) )
THEN
195 CALL xerbla(
'ZLA_SYRCOND_X', -info )
199 IF ( lsame( uplo,
'U' ) ) up = .true.
208 tmp = tmp + cabs1( a( j, i ) * x( j ) )
211 tmp = tmp + cabs1( a( i, j ) * x( j ) )
214 anorm = max( anorm, tmp )
220 tmp = tmp + cabs1( a( i, j ) * x( j ) )
223 tmp = tmp + cabs1( a( j, i ) * x( j ) )
226 anorm = max( anorm, tmp )
235 ELSE IF( anorm .EQ. 0.0d+0 )
THEN
245 CALL zlacn2( n, work( n+1 ), work, ainvnm, kase, isave )
252 work( i ) = work( i ) * rwork( i )
256 CALL zsytrs(
'U', n, 1, af, ldaf, ipiv,
259 CALL zsytrs(
'L', n, 1, af, ldaf, ipiv,
266 work( i ) = work( i ) / x( i )
273 work( i ) = work( i ) / x( i )
277 CALL zsytrs(
'U', n, 1, af, ldaf, ipiv,
280 CALL zsytrs(
'L', n, 1, af, ldaf, ipiv,
287 work( i ) = work( i ) * rwork( i )
295 IF( ainvnm .NE. 0.0d+0 )
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 xerbla(SRNAME, INFO)
XERBLA
subroutine zsytrs(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
ZSYTRS
double precision function zla_syrcond_x(UPLO, N, A, LDA, AF, LDAF, IPIV, X, INFO, WORK, RWORK)
ZLA_SYRCOND_X computes the infinity norm condition number of op(A)*diag(x) for symmetric indefinite m...