138 SUBROUTINE dlaqge( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX,
148 DOUBLE PRECISION AMAX, COLCND, ROWCND
151 DOUBLE PRECISION A( LDA, * ), C( * ), R( * )
157 DOUBLE PRECISION ONE, THRESH
158 parameter( one = 1.0d+0, thresh = 0.1d+0 )
162 DOUBLE PRECISION CJ, LARGE, SMALL
165 DOUBLE PRECISION DLAMCH
172 IF( m.LE.0 .OR. n.LE.0 )
THEN
179 small = dlamch(
'Safe minimum' ) / dlamch(
'Precision' )
182 IF( rowcnd.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
187 IF( colcnd.GE.thresh )
THEN
199 a( i, j ) = cj*a( i, j )
204 ELSE IF( colcnd.GE.thresh )
THEN
210 a( i, j ) = r( i )*a( i, j )
221 a( i, j ) = cj*r( i )*a( i, j )
subroutine dlaqge(m, n, a, lda, r, c, rowcnd, colcnd, amax, equed)
DLAQGE scales a general rectangular matrix, using row and column scaling factors computed by sgeequ.