Skip to content

cleanup: normalize named constants in LAPACK *equb routines #1227

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
nakatamaho wants to merge 1 commit into
base: master
Choose a base branch
from
Open

cleanup: normalize named constants in LAPACK *equb routines #1227

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions SRC/cgbequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -175,8 +175,8 @@ SUBROUTINE CGBEQUB( M, N, KL, KU, AB, LDAB, R, C, ROWCND,
* =====================================================================
*
* .. Parameters ..
REAL ONE, ZERO
PARAMETER ( ONE = 1.0E+0, ZERO = 0.0E+0 )
REAL ZERO, ONE
PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )
* ..
* .. Local Scalars ..
INTEGER I, J, KD
Expand Down
4 changes: 2 additions & 2 deletions SRC/cgeequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -160,8 +160,8 @@ SUBROUTINE CGEEQUB( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX,
* =====================================================================
*
* .. Parameters ..
REAL ONE, ZERO
PARAMETER ( ONE = 1.0E+0, ZERO = 0.0E+0 )
REAL ZERO, ONE
PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )
* ..
* .. Local Scalars ..
INTEGER I, J
Expand Down
16 changes: 8 additions & 8 deletions SRC/cheequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -147,8 +147,8 @@ SUBROUTINE CHEEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
* =====================================================================
*
* .. Parameters ..
REAL ONE, ZERO
PARAMETER ( ONE = 1.0E0, ZERO = 0.0E0 )
REAL ZERO, ONE, TWO
PARAMETER ( ZERO = 0.0E0, ONE = 1.0E0, TWO = 2.0E0 )
INTEGER MAX_ITER
PARAMETER ( MAX_ITER = 100 )
* ..
Expand Down Expand Up @@ -231,14 +231,14 @@ SUBROUTINE CHEEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END DO
END IF
DO J = 1, N
S( J ) = 1.0E0 / S( J )
S( J ) = ONE / S( J )
END DO

TOL = ONE / SQRT( 2.0E0 * REAL( N ) )
TOL = ONE / SQRT( TWO * REAL( N ) )

DO ITER = 1, MAX_ITER
SCALE = 0.0E0
SUMSQ = 0.0E0
SCALE = ZERO
SUMSQ = ZERO
* beta = |A|s
DO I = 1, N
WORK( I ) = ZERO
Expand All @@ -262,13 +262,13 @@ SUBROUTINE CHEEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END IF

* avg = s^T beta / n
AVG = 0.0E0
AVG = ZERO
DO I = 1, N
AVG = AVG + REAL( S( I )*WORK( I ) )
END DO
AVG = AVG / REAL( N )

STD = 0.0E0
STD = ZERO
DO I = N+1, 2*N
WORK( I ) = S( I-N ) * WORK( I-N ) - AVG
END DO
Expand Down
6 changes: 3 additions & 3 deletions SRC/cpoequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -132,8 +132,8 @@ SUBROUTINE CPOEQUB( N, A, LDA, S, SCOND, AMAX, INFO )
* =====================================================================
*
* .. Parameters ..
REAL ZERO, ONE
PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )
REAL ZERO, HALF, ONE
PARAMETER ( ZERO = 0.0E+0, HALF = 0.5E+0, ONE = 1.0E+0 )
* ..
* .. Local Scalars ..
INTEGER I
Expand Down Expand Up @@ -175,7 +175,7 @@ SUBROUTINE CPOEQUB( N, A, LDA, S, SCOND, AMAX, INFO )
END IF

BASE = SLAMCH( 'B' )
TMP = -0.5 / LOG ( BASE )
TMP = -HALF / LOG ( BASE )
*
* Find the minimum and maximum diagonal elements.
*
Expand Down
16 changes: 8 additions & 8 deletions SRC/csyequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -147,8 +147,8 @@ SUBROUTINE CSYEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
* =====================================================================
*
* .. Parameters ..
REAL ONE, ZERO
PARAMETER ( ONE = 1.0E0, ZERO = 0.0E0 )
REAL ZERO, ONE, TWO
PARAMETER ( ZERO = 0.0E0, ONE = 1.0E0, TWO = 2.0E0 )
INTEGER MAX_ITER
PARAMETER ( MAX_ITER = 100 )
* ..
Expand Down Expand Up @@ -231,14 +231,14 @@ SUBROUTINE CSYEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END DO
END IF
DO J = 1, N
S( J ) = 1.0 / S( J )
S( J ) = ONE / S( J )
END DO

TOL = ONE / SQRT( 2.0E0 * REAL( N ) )
TOL = ONE / SQRT( TWO * REAL( N ) )

DO ITER = 1, MAX_ITER
SCALE = 0.0E0
SUMSQ = 0.0E0
SCALE = ZERO
SUMSQ = ZERO
* beta = |A|s
DO I = 1, N
WORK( I ) = ZERO
Expand All @@ -262,13 +262,13 @@ SUBROUTINE CSYEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END IF

* avg = s^T beta / n
AVG = 0.0E0
AVG = ZERO
DO I = 1, N
AVG = AVG + REAL( S( I )*WORK( I ) )
END DO
AVG = AVG / REAL( N )

STD = 0.0E0
STD = ZERO
DO I = N+1, 2*N
WORK( I ) = S( I-N ) * WORK( I-N ) - AVG
END DO
Expand Down
4 changes: 2 additions & 2 deletions SRC/dgbequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -173,8 +173,8 @@ SUBROUTINE DGBEQUB( M, N, KL, KU, AB, LDAB, R, C, ROWCND,
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION ONE, ZERO
PARAMETER ( ONE = 1.0D+0, ZERO = 0.0D+0 )
DOUBLE PRECISION ZERO, ONE
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
* ..
* .. Local Scalars ..
INTEGER I, J, KD
Expand Down
4 changes: 2 additions & 2 deletions SRC/dgeequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -158,8 +158,8 @@ SUBROUTINE DGEEQUB( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX,
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION ONE, ZERO
PARAMETER ( ONE = 1.0D+0, ZERO = 0.0D+0 )
DOUBLE PRECISION ZERO, ONE
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
* ..
* .. Local Scalars ..
INTEGER I, J
Expand Down
6 changes: 3 additions & 3 deletions SRC/dpoequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,8 +130,8 @@ SUBROUTINE DPOEQUB( N, A, LDA, S, SCOND, AMAX, INFO )
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION ZERO, ONE
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
DOUBLE PRECISION ZERO, HALF, ONE
PARAMETER ( ZERO = 0.0D+0, HALF = 0.5D+0, ONE = 1.0D+0 )
* ..
* .. Local Scalars ..
INTEGER I
Expand Down Expand Up @@ -173,7 +173,7 @@ SUBROUTINE DPOEQUB( N, A, LDA, S, SCOND, AMAX, INFO )
END IF

BASE = DLAMCH( 'B' )
TMP = -0.5D+0 / LOG ( BASE )
TMP = -HALF / LOG ( BASE )
*
* Find the minimum and maximum diagonal elements.
*
Expand Down
16 changes: 8 additions & 8 deletions SRC/dsyequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -145,8 +145,8 @@ SUBROUTINE DSYEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION ONE, ZERO
PARAMETER ( ONE = 1.0D0, ZERO = 0.0D0 )
DOUBLE PRECISION ZERO, ONE, TWO
PARAMETER ( ZERO = 0.0D0, ONE = 1.0D0, TWO = 2.0D0 )
INTEGER MAX_ITER
PARAMETER ( MAX_ITER = 100 )
* ..
Expand Down Expand Up @@ -222,14 +222,14 @@ SUBROUTINE DSYEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END DO
END IF
DO J = 1, N
S( J ) = 1.0D0 / S( J )
S( J ) = ONE / S( J )
END DO

TOL = ONE / SQRT( 2.0D0 * N )
TOL = ONE / SQRT( TWO * N )

DO ITER = 1, MAX_ITER
SCALE = 0.0D0
SUMSQ = 0.0D0
SCALE = ZERO
SUMSQ = ZERO
* beta = |A|s
DO I = 1, N
WORK( I ) = ZERO
Expand All @@ -253,13 +253,13 @@ SUBROUTINE DSYEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END IF

* avg = s^T beta / n
AVG = 0.0D0
AVG = ZERO
DO I = 1, N
AVG = AVG + S( I )*WORK( I )
END DO
AVG = AVG / N

STD = 0.0D0
STD = ZERO
DO I = N+1, 2*N
WORK( I ) = S( I-N ) * WORK( I-N ) - AVG
END DO
Expand Down
4 changes: 2 additions & 2 deletions SRC/sgbequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -173,8 +173,8 @@ SUBROUTINE SGBEQUB( M, N, KL, KU, AB, LDAB, R, C, ROWCND,
* =====================================================================
*
* .. Parameters ..
REAL ONE, ZERO
PARAMETER ( ONE = 1.0E+0, ZERO = 0.0E+0 )
REAL ZERO, ONE
PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )
* ..
* .. Local Scalars ..
INTEGER I, J, KD
Expand Down
4 changes: 2 additions & 2 deletions SRC/sgeequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -158,8 +158,8 @@ SUBROUTINE SGEEQUB( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX,
* =====================================================================
*
* .. Parameters ..
REAL ONE, ZERO
PARAMETER ( ONE = 1.0E+0, ZERO = 0.0E+0 )
REAL ZERO, ONE
PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )
* ..
* .. Local Scalars ..
INTEGER I, J
Expand Down
6 changes: 3 additions & 3 deletions SRC/spoequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,8 +130,8 @@ SUBROUTINE SPOEQUB( N, A, LDA, S, SCOND, AMAX, INFO )
* =====================================================================
*
* .. Parameters ..
REAL ZERO, ONE
PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )
REAL ZERO, HALF, ONE
PARAMETER ( ZERO = 0.0E+0, HALF = 0.5E+0, ONE = 1.0E+0 )
* ..
* .. Local Scalars ..
INTEGER I
Expand Down Expand Up @@ -173,7 +173,7 @@ SUBROUTINE SPOEQUB( N, A, LDA, S, SCOND, AMAX, INFO )
END IF

BASE = SLAMCH( 'B' )
TMP = -0.5 / LOG ( BASE )
TMP = -HALF / LOG ( BASE )
*
* Find the minimum and maximum diagonal elements.
*
Expand Down
16 changes: 8 additions & 8 deletions SRC/ssyequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -145,8 +145,8 @@ SUBROUTINE SSYEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
* =====================================================================
*
* .. Parameters ..
REAL ONE, ZERO
PARAMETER ( ONE = 1.0E0, ZERO = 0.0E0 )
REAL ZERO, ONE, TWO
PARAMETER ( ZERO = 0.0E0, ONE = 1.0E0, TWO = 2.0E0 )
INTEGER MAX_ITER
PARAMETER ( MAX_ITER = 100 )
* ..
Expand Down Expand Up @@ -222,14 +222,14 @@ SUBROUTINE SSYEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END DO
END IF
DO J = 1, N
S( J ) = 1.0E0 / S( J )
S( J ) = ONE / S( J )
END DO

TOL = ONE / SQRT( 2.0E0 * REAL( N ) )
TOL = ONE / SQRT( TWO * REAL( N ) )

DO ITER = 1, MAX_ITER
SCALE = 0.0E0
SUMSQ = 0.0E0
SCALE = ZERO
SUMSQ = ZERO
* beta = |A|s
DO I = 1, N
WORK( I ) = ZERO
Expand All @@ -253,13 +253,13 @@ SUBROUTINE SSYEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END IF

* avg = s^T beta / n
AVG = 0.0E0
AVG = ZERO
DO I = 1, N
AVG = AVG + S( I )*WORK( I )
END DO
AVG = AVG / REAL( N )

STD = 0.0E0
STD = ZERO
DO I = N+1, 2*N
WORK( I ) = S( I-N ) * WORK( I-N ) - AVG
END DO
Expand Down
4 changes: 2 additions & 2 deletions SRC/zgbequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -175,8 +175,8 @@ SUBROUTINE ZGBEQUB( M, N, KL, KU, AB, LDAB, R, C, ROWCND,
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION ONE, ZERO
PARAMETER ( ONE = 1.0D+0, ZERO = 0.0D+0 )
DOUBLE PRECISION ZERO, ONE
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
* ..
* .. Local Scalars ..
INTEGER I, J, KD
Expand Down
4 changes: 2 additions & 2 deletions SRC/zgeequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -160,8 +160,8 @@ SUBROUTINE ZGEEQUB( M, N, A, LDA, R, C, ROWCND, COLCND, AMAX,
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION ONE, ZERO
PARAMETER ( ONE = 1.0D+0, ZERO = 0.0D+0 )
DOUBLE PRECISION ZERO, ONE
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
* ..
* .. Local Scalars ..
INTEGER I, J
Expand Down
16 changes: 8 additions & 8 deletions SRC/zheequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -147,8 +147,8 @@ SUBROUTINE ZHEEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION ONE, ZERO
PARAMETER ( ONE = 1.0D0, ZERO = 0.0D0 )
DOUBLE PRECISION ZERO, ONE, TWO
PARAMETER ( ZERO = 0.0D0, ONE = 1.0D0, TWO = 2.0D0 )
INTEGER MAX_ITER
PARAMETER ( MAX_ITER = 100 )
* ..
Expand Down Expand Up @@ -231,14 +231,14 @@ SUBROUTINE ZHEEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END DO
END IF
DO J = 1, N
S( J ) = 1.0D0 / S( J )
S( J ) = ONE / S( J )
END DO

TOL = ONE / SQRT( 2.0D0 * N )
TOL = ONE / SQRT( TWO * N )

DO ITER = 1, MAX_ITER
SCALE = 0.0D0
SUMSQ = 0.0D0
SCALE = ZERO
SUMSQ = ZERO
* beta = |A|s
DO I = 1, N
WORK( I ) = ZERO
Expand All @@ -262,13 +262,13 @@ SUBROUTINE ZHEEQUB( UPLO, N, A, LDA, S, SCOND, AMAX, WORK,
END IF

* avg = s^T beta / n
AVG = 0.0D0
AVG = ZERO
DO I = 1, N
AVG = AVG + DBLE( S( I )*WORK( I ) )
END DO
AVG = AVG / N

STD = 0.0D0
STD = ZERO
DO I = N+1, 2*N
WORK( I ) = S( I-N ) * WORK( I-N ) - AVG
END DO
Expand Down
6 changes: 3 additions & 3 deletions SRC/zpoequb.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -132,8 +132,8 @@ SUBROUTINE ZPOEQUB( N, A, LDA, S, SCOND, AMAX, INFO )
* =====================================================================
*
* .. Parameters ..
DOUBLE PRECISION ZERO, ONE
PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 )
DOUBLE PRECISION ZERO, HALF, ONE
PARAMETER ( ZERO = 0.0D+0, HALF = 0.5D+0, ONE = 1.0D+0 )
* ..
* .. Local Scalars ..
INTEGER I
Expand Down Expand Up @@ -175,7 +175,7 @@ SUBROUTINE ZPOEQUB( N, A, LDA, S, SCOND, AMAX, INFO )
END IF

BASE = DLAMCH( 'B' )
TMP = -0.5D+0 / LOG ( BASE )
TMP = -HALF / LOG ( BASE )
*
* Find the minimum and maximum diagonal elements.
*
Expand Down
Loading
Loading