160 REAL FUNCTION cla_gbrcond_c( TRANS, N, KL, KU, AB, LDAB, AFB,
161 $ ldafb, ipiv, c, capply, info, work,
172 INTEGER N, KL, KU, KD, KE, LDAB, LDAFB, INFO
176 COMPLEX AB( ldab, * ), AFB( ldafb, * ), WORK( * )
177 REAL C( * ), RWORK( * )
185 REAL AINVNM, ANORM, TMP
205 cabs1( zdum ) = abs(
REAL( ZDUM ) ) + abs( AIMAG( zdum ) )
211 notrans = lsame( trans,
'N' )
212 IF ( .NOT. notrans .AND. .NOT. lsame( trans,
'T' ) .AND. .NOT.
213 $ lsame( trans,
'C' ) )
THEN
215 ELSE IF( n.LT.0 )
THEN
217 ELSE IF( kl.LT.0 .OR. kl.GT.n-1 )
THEN
219 ELSE IF( ku.LT.0 .OR. ku.GT.n-1 )
THEN
221 ELSE IF( ldab.LT.kl+ku+1 )
THEN
223 ELSE IF( ldafb.LT.2*kl+ku+1 )
THEN
227 CALL xerbla(
'CLA_GBRCOND_C', -info )
240 DO j = max( i-kl, 1 ), min( i+ku, n )
241 tmp = tmp + cabs1( ab( kd+i-j, j ) ) / c( j )
244 DO j = max( i-kl, 1 ), min( i+ku, n )
245 tmp = tmp + cabs1( ab( kd+i-j, j ) )
249 anorm = max( anorm, tmp )
255 DO j = max( i-kl, 1 ), min( i+ku, n )
256 tmp = tmp + cabs1( ab( ke-i+j, i ) ) / c( j )
259 DO j = max( i-kl, 1 ), min( i+ku, n )
260 tmp = tmp + cabs1( ab( ke-i+j, i ) )
264 anorm = max( anorm, tmp )
273 ELSE IF( anorm .EQ. 0.0e+0 )
THEN
283 CALL clacn2( n, work( n+1 ), work, ainvnm, kase, isave )
290 work( i ) = work( i ) * rwork( i )
294 CALL cgbtrs(
'No transpose', n, kl, ku, 1, afb, ldafb,
295 $ ipiv, work, n, info )
297 CALL cgbtrs(
'Conjugate transpose', n, kl, ku, 1, afb,
298 $ ldafb, ipiv, work, n, info )
305 work( i ) = work( i ) * c( i )
314 work( i ) = work( i ) * c( i )
319 CALL cgbtrs(
'Conjugate transpose', n, kl, ku, 1, afb,
320 $ ldafb, ipiv, work, n, info )
322 CALL cgbtrs(
'No transpose', n, kl, ku, 1, afb, ldafb,
323 $ ipiv, work, n, info )
329 work( i ) = work( i ) * rwork( i )
337 IF( ainvnm .NE. 0.0e+0 )
subroutine xerbla(SRNAME, INFO)
XERBLA
real function cla_gbrcond_c(TRANS, N, KL, KU, AB, LDAB, AFB, LDAFB, IPIV, C, CAPPLY, INFO, WORK, RWORK)
CLA_GBRCOND_C computes the infinity norm condition number of op(A)*inv(diag(c)) for general banded ma...
subroutine clacn2(N, V, X, EST, KASE, ISAVE)
CLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...
subroutine cgbtrs(TRANS, N, KL, KU, NRHS, AB, LDAB, IPIV, B, LDB, INFO)
CGBTRS