142 SUBROUTINE dlaqge( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX,
153 DOUBLE PRECISION AMAX, COLCND, ROWCND
156 DOUBLE PRECISION A( lda, * ), C( * ), R( * )
162 DOUBLE PRECISION ONE, THRESH
163 parameter ( one = 1.0d+0, thresh = 0.1d+0 )
167 DOUBLE PRECISION CJ, LARGE, SMALL
170 DOUBLE PRECISION DLAMCH
177 IF( m.LE.0 .OR. n.LE.0 )
THEN
184 small = dlamch(
'Safe minimum' ) / dlamch(
'Precision' )
187 IF( rowcnd.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
192 IF( colcnd.GE.thresh )
THEN
204 a( i, j ) = cj*a( i, j )
209 ELSE IF( colcnd.GE.thresh )
THEN
215 a( i, j ) = r( i )*a( i, j )
226 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...