Skip to content

Commit a015b21

Browse files
SORBDB6: replace integer literal with constant
This change makes it easier to port the recent changes to SORBDB6 to (complex) double precision.
1 parent deed8f3 commit a015b21

File tree

1 file changed

+4
-3
lines changed

1 file changed

+4
-3
lines changed

SRC/sorbdb6.f

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -173,8 +173,9 @@ SUBROUTINE SORBDB6( M1, M2, N, X1, INCX1, X2, INCX2, Q1, LDQ1, Q2,
173173
* =====================================================================
174174
*
175175
* .. Parameters ..
176-
REAL ALPHA, REALZERO
177-
PARAMETER ( ALPHA = 0.1E0, REALZERO = 0.0E0 )
176+
REAL ALPHA, REALONE, REALZERO
177+
PARAMETER ( ALPHA = 0.01E0, REALONE = 1.0E0,
178+
$ REALZERO = 0.0E0 )
178179
REAL NEGONE, ONE, ZERO
179180
PARAMETER ( NEGONE = -1.0E0, ONE = 1.0E0, ZERO = 0.0E0 )
180181
* ..
@@ -228,7 +229,7 @@ SUBROUTINE SORBDB6( M1, M2, N, X1, INCX1, X2, INCX2, Q1, LDQ1, Q2,
228229
* and an assertion added comparing the norm with one. Alas, Fortran
229230
* never made it into 1989 when assert() was introduced into the C
230231
* programming language.
231-
NORM = 1.0E0
232+
NORM = REALONE
232233
*
233234
IF( M1 .EQ. 0 ) THEN
234235
DO I = 1, N

0 commit comments

Comments
 (0)