139 $ LDAF, IPIV, C, CAPPLY,
140 $ INFO, WORK, RWORK )
149 INTEGER n, lda, ldaf, info
153 COMPLEX*16 a( lda, * ), af( ldaf, * ), work( * )
154 DOUBLE PRECISION c( * ), rwork( * )
162 DOUBLE PRECISION ainvnm, anorm, tmp
176 INTRINSIC abs, max, real, dimag
179 DOUBLE PRECISION cabs1
182 cabs1( zdum ) = abs( dble( zdum ) ) + abs( dimag( zdum ) )
188 notrans =
lsame( trans,
'N' )
189 IF ( .NOT. notrans .AND. .NOT.
lsame( trans,
'T' ) .AND. .NOT.
190 $
lsame( trans,
'C' ) )
THEN
192 ELSE IF( n.LT.0 )
THEN
194 ELSE IF( lda.LT.max( 1, n ) )
THEN
196 ELSE IF( ldaf.LT.max( 1, n ) )
THEN
200 CALL xerbla(
'ZLA_GERCOND_C', -info )
212 tmp = tmp + cabs1( a( i, j ) ) / c( j )
216 tmp = tmp + cabs1( a( i, j ) )
220 anorm = max( anorm, tmp )
227 tmp = tmp + cabs1( a( j, i ) ) / c( j )
231 tmp = tmp + cabs1( a( j, i ) )
235 anorm = max( anorm, tmp )
244 ELSE IF( anorm .EQ. 0.0d+0 )
THEN
254 CALL zlacn2( n, work( n+1 ), work, ainvnm, kase, isave )
261 work( i ) = work( i ) * rwork( i )
265 CALL zgetrs(
'No transpose', n, 1, af, ldaf, ipiv,
268 CALL zgetrs(
'Conjugate transpose', n, 1, af, ldaf,
277 work( i ) = work( i ) * c( i )
286 work( i ) = work( i ) * c( i )
291 CALL zgetrs(
'Conjugate transpose', n, 1, af, ldaf,
295 CALL zgetrs(
'No transpose', n, 1, af, ldaf, ipiv,
302 work( i ) = work( i ) * rwork( i )
310 IF( ainvnm .NE. 0.0d+0 )
subroutine zgetrs(trans, n, nrhs, a, lda, ipiv, b, ldb, info)
ZGETRS
double precision function zla_gercond_c(trans, n, a, lda, af, ldaf, ipiv, c, capply, info, work, rwork)
ZLA_GERCOND_C computes the infinity norm condition number of op(A)*inv(diag(c)) for general matrices.
subroutine zlacn2(n, v, x, est, kase, isave)
ZLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...