114 INTEGER ncols, lda, ldaf
117 REAL a( lda, * ), af( ldaf, * ), work( * )
124 REAL amax, umax, rpvgrw
128 INTRINSIC abs, max, min
136 upper =
lsame(
'Upper', uplo )
153 $ max( abs( a( i,
j ) ), work( ncols+
j ) )
160 $ max( abs( a( i,
j ) ), work( ncols+
j ) )
168 IF (
lsame(
'Upper', uplo ) )
THEN
171 work(
j ) = max( abs( af( i,
j ) ), work(
j ) )
177 work(
j ) = max( abs( af( i,
j ) ), work(
j ) )
189 IF (
lsame(
'Upper', uplo ) )
THEN
192 amax = work( ncols+i )
193 IF ( umax /= 0.0 )
THEN
194 rpvgrw = min( amax / umax, rpvgrw )
200 amax = work( ncols+i )
201 IF ( umax /= 0.0 )
THEN
202 rpvgrw = min( amax / umax, rpvgrw )
subroutine slaset(UPLO, M, N, ALPHA, BETA, A, LDA)
SLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values...
logical function lsame(CA, CB)
LSAME
real function sla_porpvgrw(UPLO, NCOLS, A, LDA, AF, LDAF, WORK)
SLA_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