1 SUBROUTINE dlaqsy( UPLO, N, A, LDA, S, SCOND, AMAX, EQUED )
10 DOUBLE PRECISION amax, scond
13 DOUBLE PRECISION a( lda, * ), s( * )
78 DOUBLE PRECISION one, thresh
79 parameter( one = 1.0d+0, thresh = 0.1d+0 )
83 DOUBLE PRECISION cj, large, small
87 DOUBLE PRECISION dlamch
88 EXTERNAL lsame, dlamch
101 small = dlamch(
'Safe minimum' ) / dlamch(
'Precision' )
104 IF( scond.GE.thresh .AND. amax.GE.small .AND. amax.LE.large )
THEN
113 IF( lsame(
uplo,
'U' ) )
THEN
120 a( i, j ) = cj*s( i )*a( i, j )
130 a( i, j ) = cj*s( i )*a( i, j )