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 )
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...
logical function lsame(CA, CB)
LSAME
subroutine zlaset(UPLO, M, N, ALPHA, BETA, A, LDA)
ZLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j