decreased EPS by a factor of 0.9 to correct problem with orthogonality of U in certain cases

james · james · commit 8b5dd5433b20 · 2011-04-07T12:30:43.000Z
diff --git a/SRC/dbdsdc.f b/SRC/dbdsdc.f
@@ -287,7 +287,7 @@ SUBROUTINE DBDSDC( UPLO, COMPQ, N, D, E, U, LDU, VT, LDVT, Q, IQ,
       CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, N, 1, D, N, IERR )
       CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, NM1, 1, E, NM1, IERR )
 *
-      EPS = DLAMCH( 'Epsilon' )
+      EPS = (0.9D+0)*DLAMCH( 'Epsilon' )
 *
       MLVL = INT( LOG( DBLE( N ) / DBLE( SMLSIZ+1 ) ) / LOG( TWO ) ) + 1
       SMLSZP = SMLSIZ + 1

Original file line number	Diff line number	Diff line change
`@@ -287,7 +287,7 @@ SUBROUTINE DBDSDC( UPLO, COMPQ, N, D, E, U, LDU, VT, LDVT, Q, IQ,`
`287`	`287`	`CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, N, 1, D, N, IERR )`
`288`	`288`	`CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, NM1, 1, E, NM1, IERR )`
`289`	`289`	`*`
`290`		`- EPS = DLAMCH( 'Epsilon' )`
	`290`	`+ EPS = (0.9D+0)*DLAMCH( 'Epsilon' )`
`291`	`291`	`*`
`292`	`292`	`MLVL = INT( LOG( DBLE( N ) / DBLE( SMLSIZ+1 ) ) / LOG( TWO ) ) + 1`
`293`	`293`	`SMLSZP = SMLSIZ + 1`