107 DOUBLE PRECISION FUNCTION zla_porpvgrw( UPLO, NCOLS, A, LDA, AF,
117 INTEGER NCOLS, LDA, LDAF
120 COMPLEX*16 A( lda, * ), AF( ldaf, * )
121 DOUBLE PRECISION WORK( * )
128 DOUBLE PRECISION AMAX, UMAX, RPVGRW
137 INTRINSIC abs, max, min,
REAL, DIMAG
140 DOUBLE PRECISION CABS1
143 cabs1( zdum ) = abs( dble( zdum ) ) + abs( dimag( zdum ) )
146 upper = lsame(
'Upper', uplo )
163 $ max( cabs1( a( i, j ) ), work( ncols+j ) )
170 $ max( cabs1( a( i, j ) ), work( ncols+j ) )
178 IF ( lsame(
'Upper', uplo ) )
THEN
181 work( j ) = max( cabs1( af( i, j ) ), work( j ) )
187 work( j ) = max( cabs1( af( i, j ) ), work( j ) )
199 IF ( lsame(
'Upper', uplo ) )
THEN
202 amax = work( ncols+i )
203 IF ( umax /= 0.0d+0 )
THEN
204 rpvgrw = min( amax / umax, rpvgrw )
210 amax = work( ncols+i )
211 IF ( umax /= 0.0d+0 )
THEN
212 rpvgrw = min( amax / umax, rpvgrw )
logical function lsame(CA, CB)
LSAME
double precision function zla_porpvgrw(UPLO, NCOLS, A, LDA, AF, LDAF, WORK)
ZLA_PORPVGRW computes the reciprocal pivot growth factor norm(A)/norm(U) for a symmetric or Hermitian...