106 DOUBLE PRECISION FUNCTION dla_porpvgrw( UPLO, NCOLS, A, LDA, AF,
116 INTEGER NCOLS, LDA, LDAF
119 DOUBLE PRECISION A( lda, * ), AF( ldaf, * ), WORK( * )
126 DOUBLE PRECISION AMAX, UMAX, RPVGRW
130 INTRINSIC abs, max, min
138 upper = lsame(
'Upper', uplo )
155 $ max( abs( a( i, j ) ), work( ncols+j ) )
162 $ max( abs( a( i, j ) ), work( ncols+j ) )
170 IF ( lsame(
'Upper', uplo ) )
THEN
173 work( j ) = max( abs( af( i, j ) ), work( j ) )
179 work( j ) = max( abs( af( i, j ) ), work( j ) )
191 IF ( lsame(
'Upper', uplo ) )
THEN
194 amax = work( ncols+i )
195 IF ( umax /= 0.0d+0 )
THEN
196 rpvgrw = min( amax / umax, rpvgrw )
202 amax = work( ncols+i )
203 IF ( umax /= 0.0d+0 )
THEN
204 rpvgrw = min( amax / umax, rpvgrw )
logical function lsame(CA, CB)
LSAME
double precision function dla_porpvgrw(UPLO, NCOLS, A, LDA, AF, LDAF, WORK)
DLA_PORPVGRW computes the reciprocal pivot growth factor norm(A)/norm(U) for a symmetric or Hermitian...