132 $ ldaf, c, capply, info,
143 INTEGER n, lda, ldaf, info
146 COMPLEX*16 a( lda, * ), af( ldaf, * ), work( * )
147 DOUBLE PRECISION c( * ), rwork( * )
154 DOUBLE PRECISION ainvnm, anorm, tmp
170 INTRINSIC abs, max,
REAL, dimag
173 DOUBLE PRECISION cabs1
176 cabs1( zdum ) = abs( dble( zdum ) ) + abs( dimag( zdum ) )
183 upper =
lsame( uplo,
'U' )
184 IF( .NOT.upper .AND. .NOT.
lsame( uplo,
'L' ) )
THEN
186 ELSE IF( n.LT.0 )
THEN
188 ELSE IF( lda.LT.max( 1, n ) )
THEN
190 ELSE IF( ldaf.LT.max( 1, n ) )
THEN
194 CALL
xerbla(
'ZLA_PORCOND_C', -info )
198 IF (
lsame( uplo,
'U' ) ) up = .true.
208 tmp = tmp + cabs1( a( j, i ) ) / c( j )
211 tmp = tmp + cabs1( a( i, j ) ) / c( j )
215 tmp = tmp + cabs1( a( j, i ) )
218 tmp = tmp + cabs1( a( i, j ) )
222 anorm = max( anorm, tmp )
229 tmp = tmp + cabs1( a( i, j ) ) / c( j )
232 tmp = tmp + cabs1( a( j, i ) ) / c( j )
236 tmp = tmp + cabs1( a( i, j ) )
239 tmp = tmp + cabs1( a( j, i ) )
243 anorm = max( anorm, tmp )
252 ELSE IF( anorm .EQ. 0.0d+0 )
THEN
262 CALL
zlacn2( n, work( n+1 ), work, ainvnm, kase, isave )
269 work( i ) = work( i ) * rwork( i )
273 CALL
zpotrs(
'U', n, 1, af, ldaf,
276 CALL
zpotrs(
'L', n, 1, af, ldaf,
284 work( i ) = work( i ) * c( i )
293 work( i ) = work( i ) * c( i )
298 CALL
zpotrs(
'U', n, 1, af, ldaf,
301 CALL
zpotrs(
'L', n, 1, af, ldaf,
308 work( i ) = work( i ) * rwork( i )
316 IF( ainvnm .NE. 0.0d+0 )