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