Skip to content

Commit ad3d9b8

Browse files
xORBDB6/xUNBDB6: use safe termination criterion
With a modest amount of random testing on a work station, it is possible to detect the loss of orthogonality caused by the old early termination criterion.
1 parent bc2b374 commit ad3d9b8

File tree

4 files changed

+4
-4
lines changed

4 files changed

+4
-4
lines changed

SRC/cunbdb6.f

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -173,7 +173,7 @@ SUBROUTINE CUNBDB6( M1, M2, N, X1, INCX1, X2, INCX2, Q1, LDQ1, Q2,
173173
*
174174
* .. Parameters ..
175175
REAL ALPHA, REALONE, REALZERO
176-
PARAMETER ( ALPHA = 0.1E0, REALONE = 1.0E0,
176+
PARAMETER ( ALPHA = 0.83E0, REALONE = 1.0E0,
177177
$ REALZERO = 0.0E0 )
178178
COMPLEX NEGONE, ONE, ZERO
179179
PARAMETER ( NEGONE = (-1.0E0,0.0E0), ONE = (1.0E0,0.0E0),

SRC/dorbdb6.f

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -173,7 +173,7 @@ SUBROUTINE DORBDB6( M1, M2, N, X1, INCX1, X2, INCX2, Q1, LDQ1, Q2,
173173
*
174174
* .. Parameters ..
175175
DOUBLE PRECISION ALPHA, REALONE, REALZERO
176-
PARAMETER ( ALPHA = 0.1D0, REALONE = 1.0D0,
176+
PARAMETER ( ALPHA = 0.83D0, REALONE = 1.0D0,
177177
$ REALZERO = 0.0D0 )
178178
DOUBLE PRECISION NEGONE, ONE, ZERO
179179
PARAMETER ( NEGONE = -1.0D0, ONE = 1.0D0, ZERO = 0.0D0 )

SRC/sorbdb6.f

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -173,7 +173,7 @@ SUBROUTINE SORBDB6( M1, M2, N, X1, INCX1, X2, INCX2, Q1, LDQ1, Q2,
173173
*
174174
* .. Parameters ..
175175
REAL ALPHA, REALONE, REALZERO
176-
PARAMETER ( ALPHA = 0.1E0, REALONE = 1.0E0,
176+
PARAMETER ( ALPHA = 0.83E0, REALONE = 1.0E0,
177177
$ REALZERO = 0.0E0 )
178178
REAL NEGONE, ONE, ZERO
179179
PARAMETER ( NEGONE = -1.0E0, ONE = 1.0E0, ZERO = 0.0E0 )

SRC/zunbdb6.f

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -173,7 +173,7 @@ SUBROUTINE ZUNBDB6( M1, M2, N, X1, INCX1, X2, INCX2, Q1, LDQ1, Q2,
173173
*
174174
* .. Parameters ..
175175
DOUBLE PRECISION ALPHA, REALONE, REALZERO
176-
PARAMETER ( ALPHA = 0.1D0, REALONE = 1.0D0,
176+
PARAMETER ( ALPHA = 0.83D0, REALONE = 1.0D0,
177177
$ REALZERO = 0.0D0 )
178178
COMPLEX*16 NEGONE, ONE, ZERO
179179
PARAMETER ( NEGONE = (-1.0D0,0.0D0), ONE = (1.0D0,0.0D0),

0 commit comments

Comments
 (0)