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 xerbla(SRNAME, INFO)
XERBLA
logical function lsame(CA, CB)
LSAME
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...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine cpotrs(UPLO, N, NRHS, A, LDA, B, LDB, INFO)
CPOTRS
subroutine clacn2(N, V, X, EST, KASE, ISAVE)
CLACN2 estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...