|
SMAX = ABS( D( M ) ) |
|
SMIN = SMAX |
|
DO 70 LLL = 1, M - 1 |
|
LL = M - LLL |
|
ABSS = ABS( D( LL ) ) |
|
ABSE = ABS( E( LL ) ) |
|
IF( TOL.LT.ZERO .AND. ABSS.LE.THRESH ) |
|
$ D( LL ) = ZERO |
|
IF( ABSE.LE.THRESH ) |
|
$ GO TO 80 |
|
SMIN = MIN( SMIN, ABSS ) |
|
SMAX = MAX( SMAX, ABSS, ABSE ) |
|
70 CONTINUE |
At the end of the loop,
SMIN is the absolute of the digonal minimum value, but its value is never used, I wonder why.
lapack/SRC/dbdsqr.f
Lines 456 to 468 in 4f7f2aa
At the end of the loop,
SMINis the absolute of the digonal minimum value, but its value is never used, I wonder why.