From e3950c04711f692553e9f80d80272ed899bb3c82 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Tue, 1 Jun 2021 17:55:40 +1000 Subject: [PATCH 01/16] comments fixed: epsilon was mistakenly called a machine precision in xGBT02 --- TESTING/LIN/cgbt02.f | 2 +- TESTING/LIN/dgbt02.f | 2 +- TESTING/LIN/sgbt02.f | 2 +- TESTING/LIN/zgbt02.f | 2 +- 4 files changed, 4 insertions(+), 4 deletions(-) diff --git a/TESTING/LIN/cgbt02.f b/TESTING/LIN/cgbt02.f index 8c51e8b2a2..3204a89867 100644 --- a/TESTING/LIN/cgbt02.f +++ b/TESTING/LIN/cgbt02.f @@ -31,7 +31,7 @@ *> equations op(A)*X = B: *> RESID = norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ), *> where op(A) = A, A**T, or A**H, depending on TRANS, and EPS is the -*> machine precision. +*> machine epsilon. *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/dgbt02.f b/TESTING/LIN/dgbt02.f index ef7e00c3d6..1184587a9c 100644 --- a/TESTING/LIN/dgbt02.f +++ b/TESTING/LIN/dgbt02.f @@ -31,7 +31,7 @@ *> equations op(A)*X = B: *> RESID = norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ), *> where op(A) = A or A**T, depending on TRANS, and EPS is the -*> machine precision. +*> machine epsilon. *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/sgbt02.f b/TESTING/LIN/sgbt02.f index 2c0734d150..82a5555490 100644 --- a/TESTING/LIN/sgbt02.f +++ b/TESTING/LIN/sgbt02.f @@ -31,7 +31,7 @@ *> equations op(A)*X = B: *> RESID = norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ), *> where op(A) = A or A**T, depending on TRANS, and EPS is the -*> machine precision. +*> machine epsilon. *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/zgbt02.f b/TESTING/LIN/zgbt02.f index 3547ea0b5e..0d13b571ad 100644 --- a/TESTING/LIN/zgbt02.f +++ b/TESTING/LIN/zgbt02.f @@ -31,7 +31,7 @@ *> equations op(A)*X = B: *> RESID = norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ), *> where op(A) = A, A**T, or A**H, depending on TRANS, and EPS is the -*> machine precision. +*> machine epsilon. *> \endverbatim * * Arguments: From 457afcf7dfe149fee23fc82e0b26492cb7318a66 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Wed, 2 Jun 2021 21:25:18 +1000 Subject: [PATCH 02/16] comments unified: homgenize notation for transpose (**T) and conjugate transpose (**H) This commit continues changes made in: * Reference-LAPACK/lapack@f2953573ede24d7f8c01fdb18de48f65f00a9943 2011-04-02 * Reference-LAPACK/lapack@d9d50d133bdf00870f37a31faefea58f0926a0a6 2011-04-02 * Reference-LAPACK/lapack@16973f0c15a4327f8a78b9652f178246895778c1 2011-04-08 --- TESTING/EIG/cbdt01.f | 22 +++++++++++----------- TESTING/EIG/cbdt02.f | 5 +++-- TESTING/EIG/dbdt01.f | 23 ++++++++++++----------- TESTING/EIG/dbdt02.f | 5 +++-- TESTING/EIG/dlarhs.f | 6 +++--- TESTING/EIG/sbdt01.f | 23 ++++++++++++----------- TESTING/EIG/sbdt02.f | 5 +++-- TESTING/EIG/slarhs.f | 6 +++--- TESTING/EIG/zbdt01.f | 23 ++++++++++++----------- TESTING/EIG/zbdt02.f | 5 +++-- TESTING/LIN/cpot01.f | 14 +++++++------- TESTING/LIN/cqrt17.f | 8 ++++---- TESTING/LIN/ctbt02.f | 2 +- TESTING/LIN/dgtt02.f | 4 ++-- TESTING/LIN/dlarhs.f | 6 +++--- TESTING/LIN/dpot01.f | 14 +++++++------- TESTING/LIN/dqrt17.f | 8 ++++---- TESTING/LIN/dtbt02.f | 6 +++--- TESTING/LIN/sgtt02.f | 4 ++-- TESTING/LIN/slarhs.f | 6 +++--- TESTING/LIN/spot01.f | 14 +++++++------- TESTING/LIN/sqrt17.f | 8 ++++---- TESTING/LIN/stbt02.f | 13 ++++++------- TESTING/LIN/zpot01.f | 14 +++++++------- TESTING/LIN/zqrt17.f | 8 ++++---- TESTING/LIN/ztbt02.f | 2 +- 26 files changed, 130 insertions(+), 124 deletions(-) diff --git a/TESTING/EIG/cbdt01.f b/TESTING/EIG/cbdt01.f index edcbcbb39f..c1bfbffa53 100644 --- a/TESTING/EIG/cbdt01.f +++ b/TESTING/EIG/cbdt01.f @@ -28,13 +28,13 @@ *> \verbatim *> *> CBDT01 reconstructs a general matrix A from its bidiagonal form -*> A = Q * B * P' -*> where Q (m by min(m,n)) and P' (min(m,n) by n) are unitary +*> A = Q * B * P**H +*> where Q (m by min(m,n)) and P**H (min(m,n) by n) are unitary *> matrices and B is bidiagonal. *> *> The test ratio to test the reduction is -*> RESID = norm( A - Q * B * PT ) / ( n * norm(A) * EPS ) -*> where PT = P' and EPS is the machine precision. +*> RESID = norm( A - Q * B * P**H ) / ( n * norm(A) * EPS ) +*> where EPS is the machine precision. *> \endverbatim * * Arguments: @@ -49,7 +49,7 @@ *> \param[in] N *> \verbatim *> N is INTEGER -*> The number of columns of the matrices A and P'. +*> The number of columns of the matrices A and P**H. *> \endverbatim *> *> \param[in] KD @@ -78,7 +78,7 @@ *> \verbatim *> Q is COMPLEX array, dimension (LDQ,N) *> The m by min(m,n) unitary matrix Q in the reduction -*> A = Q * B * P'. +*> A = Q * B * P**H. *> \endverbatim *> *> \param[in] LDQ @@ -103,8 +103,8 @@ *> \param[in] PT *> \verbatim *> PT is COMPLEX array, dimension (LDPT,N) -*> The min(m,n) by n unitary matrix P' in the reduction -*> A = Q * B * P'. +*> The min(m,n) by n unitary matrix P**H in the reduction +*> A = Q * B * P**H. *> \endverbatim *> *> \param[in] LDPT @@ -127,7 +127,7 @@ *> \param[out] RESID *> \verbatim *> RESID is REAL -*> The test ratio: norm(A - Q * B * P') / ( n * norm(A) * EPS ) +*> The test ratio: norm(A - Q * B * P**H) / ( n * norm(A) * EPS ) *> \endverbatim * * Authors: @@ -187,7 +187,7 @@ SUBROUTINE CBDT01( M, N, KD, A, LDA, Q, LDQ, D, E, PT, LDPT, WORK, RETURN END IF * -* Compute A - Q * B * P' one column at a time. +* Compute A - Q * B * P**H one column at a time. * RESID = ZERO IF( KD.NE.0 ) THEN @@ -265,7 +265,7 @@ SUBROUTINE CBDT01( M, N, KD, A, LDA, Q, LDQ, D, E, PT, LDPT, WORK, END IF END IF * -* Compute norm(A - Q * B * P') / ( n * norm(A) * EPS ) +* Compute norm(A - Q * B * P**H) / ( n * norm(A) * EPS ) * ANORM = CLANGE( '1', M, N, A, LDA, RWORK ) EPS = SLAMCH( 'Precision' ) diff --git a/TESTING/EIG/cbdt02.f b/TESTING/EIG/cbdt02.f index 87fda4a71d..8aafefd227 100644 --- a/TESTING/EIG/cbdt02.f +++ b/TESTING/EIG/cbdt02.f @@ -27,7 +27,8 @@ *> *> \verbatim *> -*> CBDT02 tests the change of basis C = U' * B by computing the residual +*> CBDT02 tests the change of basis C = U**H * B by computing the +*> residual *> *> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), *> @@ -66,7 +67,7 @@ *> \param[in] C *> \verbatim *> C is COMPLEX array, dimension (LDC,N) -*> The m by n matrix C, assumed to contain U' * B. +*> The m by n matrix C, assumed to contain U**H * B. *> \endverbatim *> *> \param[in] LDC diff --git a/TESTING/EIG/dbdt01.f b/TESTING/EIG/dbdt01.f index 40ac373ba8..f978493ea4 100644 --- a/TESTING/EIG/dbdt01.f +++ b/TESTING/EIG/dbdt01.f @@ -27,13 +27,13 @@ *> \verbatim *> *> DBDT01 reconstructs a general matrix A from its bidiagonal form -*> A = Q * B * P' -*> where Q (m by min(m,n)) and P' (min(m,n) by n) are orthogonal +*> A = Q * B * P**T +*> where Q (m by min(m,n)) and P**T (min(m,n) by n) are orthogonal *> matrices and B is bidiagonal. *> *> The test ratio to test the reduction is -*> RESID = norm( A - Q * B * PT ) / ( n * norm(A) * EPS ) -*> where PT = P' and EPS is the machine precision. +*> RESID = norm( A - Q * B * P**T ) / ( n * norm(A) * EPS ) +*> where EPS is the machine precision. *> \endverbatim * * Arguments: @@ -48,7 +48,7 @@ *> \param[in] N *> \verbatim *> N is INTEGER -*> The number of columns of the matrices A and P'. +*> The number of columns of the matrices A and P**T. *> \endverbatim *> *> \param[in] KD @@ -77,7 +77,7 @@ *> \verbatim *> Q is DOUBLE PRECISION array, dimension (LDQ,N) *> The m by min(m,n) orthogonal matrix Q in the reduction -*> A = Q * B * P'. +*> A = Q * B * P**T. *> \endverbatim *> *> \param[in] LDQ @@ -102,8 +102,8 @@ *> \param[in] PT *> \verbatim *> PT is DOUBLE PRECISION array, dimension (LDPT,N) -*> The min(m,n) by n orthogonal matrix P' in the reduction -*> A = Q * B * P'. +*> The min(m,n) by n orthogonal matrix P**T in the reduction +*> A = Q * B * P**T. *> \endverbatim *> *> \param[in] LDPT @@ -121,7 +121,8 @@ *> \param[out] RESID *> \verbatim *> RESID is DOUBLE PRECISION -*> The test ratio: norm(A - Q * B * P') / ( n * norm(A) * EPS ) +*> The test ratio: +*> norm(A - Q * B * P**T) / ( n * norm(A) * EPS ) *> \endverbatim * * Authors: @@ -180,7 +181,7 @@ SUBROUTINE DBDT01( M, N, KD, A, LDA, Q, LDQ, D, E, PT, LDPT, WORK, RETURN END IF * -* Compute A - Q * B * P' one column at a time. +* Compute A - Q * B * P**T one column at a time. * RESID = ZERO IF( KD.NE.0 ) THEN @@ -258,7 +259,7 @@ SUBROUTINE DBDT01( M, N, KD, A, LDA, Q, LDQ, D, E, PT, LDPT, WORK, END IF END IF * -* Compute norm(A - Q * B * P') / ( n * norm(A) * EPS ) +* Compute norm(A - Q * B * P**T) / ( n * norm(A) * EPS ) * ANORM = DLANGE( '1', M, N, A, LDA, WORK ) EPS = DLAMCH( 'Precision' ) diff --git a/TESTING/EIG/dbdt02.f b/TESTING/EIG/dbdt02.f index fd1b030e0c..76e3ff713d 100644 --- a/TESTING/EIG/dbdt02.f +++ b/TESTING/EIG/dbdt02.f @@ -25,7 +25,8 @@ *> *> \verbatim *> -*> DBDT02 tests the change of basis C = U' * B by computing the residual +*> DBDT02 tests the change of basis C = U**H * B by computing the +*> residual *> *> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), *> @@ -64,7 +65,7 @@ *> \param[in] C *> \verbatim *> C is DOUBLE PRECISION array, dimension (LDC,N) -*> The m by n matrix C, assumed to contain U' * B. +*> The m by n matrix C, assumed to contain U**H * B. *> \endverbatim *> *> \param[in] LDC diff --git a/TESTING/EIG/dlarhs.f b/TESTING/EIG/dlarhs.f index 2505ae227b..a5ba318bed 100644 --- a/TESTING/EIG/dlarhs.f +++ b/TESTING/EIG/dlarhs.f @@ -30,7 +30,7 @@ *> DLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system *> op( A ) * X = B, -*> where op( A ) may be A or A' (transpose of A). +*> where op( A ) = A or A**T, depending on TRANS. *> \endverbatim * * Arguments: @@ -81,8 +81,8 @@ *> TRANS is CHARACTER*1 *> Specifies the operation applied to the matrix A. *> = 'N': System is A * x = b -*> = 'T': System is A'* x = b -*> = 'C': System is A'* x = b +*> = 'T': B := A**T * X (Transpose) +*> = 'C': B := A**H * X (Conjugate transpose = Transpose) *> \endverbatim *> *> \param[in] M diff --git a/TESTING/EIG/sbdt01.f b/TESTING/EIG/sbdt01.f index 2e4e48df33..0861818506 100644 --- a/TESTING/EIG/sbdt01.f +++ b/TESTING/EIG/sbdt01.f @@ -27,13 +27,13 @@ *> \verbatim *> *> SBDT01 reconstructs a general matrix A from its bidiagonal form -*> A = Q * B * P' -*> where Q (m by min(m,n)) and P' (min(m,n) by n) are orthogonal +*> A = Q * B * P**T +*> where Q (m by min(m,n)) and P**T (min(m,n) by n) are orthogonal *> matrices and B is bidiagonal. *> *> The test ratio to test the reduction is -*> RESID = norm( A - Q * B * PT ) / ( n * norm(A) * EPS ) -*> where PT = P' and EPS is the machine precision. +*> RESID = norm( A - Q * B * P**T ) / ( n * norm(A) * EPS ) +*> where EPS is the machine precision. *> \endverbatim * * Arguments: @@ -48,7 +48,7 @@ *> \param[in] N *> \verbatim *> N is INTEGER -*> The number of columns of the matrices A and P'. +*> The number of columns of the matrices A and P**T. *> \endverbatim *> *> \param[in] KD @@ -77,7 +77,7 @@ *> \verbatim *> Q is REAL array, dimension (LDQ,N) *> The m by min(m,n) orthogonal matrix Q in the reduction -*> A = Q * B * P'. +*> A = Q * B * P**T. *> \endverbatim *> *> \param[in] LDQ @@ -102,8 +102,8 @@ *> \param[in] PT *> \verbatim *> PT is REAL array, dimension (LDPT,N) -*> The min(m,n) by n orthogonal matrix P' in the reduction -*> A = Q * B * P'. +*> The min(m,n) by n orthogonal matrix P**T in the reduction +*> A = Q * B * P**T. *> \endverbatim *> *> \param[in] LDPT @@ -121,7 +121,8 @@ *> \param[out] RESID *> \verbatim *> RESID is REAL -*> The test ratio: norm(A - Q * B * P') / ( n * norm(A) * EPS ) +*> The test ratio: +*> norm(A - Q * B * P**T) / ( n * norm(A) * EPS ) *> \endverbatim * * Authors: @@ -180,7 +181,7 @@ SUBROUTINE SBDT01( M, N, KD, A, LDA, Q, LDQ, D, E, PT, LDPT, WORK, RETURN END IF * -* Compute A - Q * B * P' one column at a time. +* Compute A - Q * B * P**T one column at a time. * RESID = ZERO IF( KD.NE.0 ) THEN @@ -258,7 +259,7 @@ SUBROUTINE SBDT01( M, N, KD, A, LDA, Q, LDQ, D, E, PT, LDPT, WORK, END IF END IF * -* Compute norm(A - Q * B * P') / ( n * norm(A) * EPS ) +* Compute norm(A - Q * B * P**T) / ( n * norm(A) * EPS ) * ANORM = SLANGE( '1', M, N, A, LDA, WORK ) EPS = SLAMCH( 'Precision' ) diff --git a/TESTING/EIG/sbdt02.f b/TESTING/EIG/sbdt02.f index cbce82ba7c..1cf8a22725 100644 --- a/TESTING/EIG/sbdt02.f +++ b/TESTING/EIG/sbdt02.f @@ -25,7 +25,8 @@ *> *> \verbatim *> -*> SBDT02 tests the change of basis C = U' * B by computing the residual +*> SBDT02 tests the change of basis C = U**H * B by computing the +*> residual *> *> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), *> @@ -64,7 +65,7 @@ *> \param[in] C *> \verbatim *> C is REAL array, dimension (LDC,N) -*> The m by n matrix C, assumed to contain U' * B. +*> The m by n matrix C, assumed to contain U**H * B. *> \endverbatim *> *> \param[in] LDC diff --git a/TESTING/EIG/slarhs.f b/TESTING/EIG/slarhs.f index 0d736f3505..0e9d309479 100644 --- a/TESTING/EIG/slarhs.f +++ b/TESTING/EIG/slarhs.f @@ -30,7 +30,7 @@ *> SLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system *> op( A ) * X = B, -*> where op( A ) may be A or A' (transpose of A). +*> where op( A ) = A or A**T, depending on TRANS. *> \endverbatim * * Arguments: @@ -81,8 +81,8 @@ *> TRANS is CHARACTER*1 *> Specifies the operation applied to the matrix A. *> = 'N': System is A * x = b -*> = 'T': System is A'* x = b -*> = 'C': System is A'* x = b +*> = 'T': B := A**T * X (Transpose) +*> = 'C': B := A**H * X (Conjugate transpose = Transpose) *> \endverbatim *> *> \param[in] M diff --git a/TESTING/EIG/zbdt01.f b/TESTING/EIG/zbdt01.f index 78d771f198..cb62993ade 100644 --- a/TESTING/EIG/zbdt01.f +++ b/TESTING/EIG/zbdt01.f @@ -28,13 +28,13 @@ *> \verbatim *> *> ZBDT01 reconstructs a general matrix A from its bidiagonal form -*> A = Q * B * P' -*> where Q (m by min(m,n)) and P' (min(m,n) by n) are unitary +*> A = Q * B * P**H +*> where Q (m by min(m,n)) and P**H (min(m,n) by n) are unitary *> matrices and B is bidiagonal. *> *> The test ratio to test the reduction is -*> RESID = norm( A - Q * B * PT ) / ( n * norm(A) * EPS ) -*> where PT = P' and EPS is the machine precision. +*> RESID = norm( A - Q * B * P**H ) / ( n * norm(A) * EPS ) +*> where EPS is the machine precision. *> \endverbatim * * Arguments: @@ -49,7 +49,7 @@ *> \param[in] N *> \verbatim *> N is INTEGER -*> The number of columns of the matrices A and P'. +*> The number of columns of the matrices A and P**H. *> \endverbatim *> *> \param[in] KD @@ -78,7 +78,7 @@ *> \verbatim *> Q is COMPLEX*16 array, dimension (LDQ,N) *> The m by min(m,n) unitary matrix Q in the reduction -*> A = Q * B * P'. +*> A = Q * B * P**H. *> \endverbatim *> *> \param[in] LDQ @@ -103,8 +103,8 @@ *> \param[in] PT *> \verbatim *> PT is COMPLEX*16 array, dimension (LDPT,N) -*> The min(m,n) by n unitary matrix P' in the reduction -*> A = Q * B * P'. +*> The min(m,n) by n unitary matrix P**H in the reduction +*> A = Q * B * P**H. *> \endverbatim *> *> \param[in] LDPT @@ -127,7 +127,8 @@ *> \param[out] RESID *> \verbatim *> RESID is DOUBLE PRECISION -*> The test ratio: norm(A - Q * B * P') / ( n * norm(A) * EPS ) +*> The test ratio: +*> norm(A - Q * B * P**H) / ( n * norm(A) * EPS ) *> \endverbatim * * Authors: @@ -187,7 +188,7 @@ SUBROUTINE ZBDT01( M, N, KD, A, LDA, Q, LDQ, D, E, PT, LDPT, WORK, RETURN END IF * -* Compute A - Q * B * P' one column at a time. +* Compute A - Q * B * P**H one column at a time. * RESID = ZERO IF( KD.NE.0 ) THEN @@ -265,7 +266,7 @@ SUBROUTINE ZBDT01( M, N, KD, A, LDA, Q, LDQ, D, E, PT, LDPT, WORK, END IF END IF * -* Compute norm(A - Q * B * P') / ( n * norm(A) * EPS ) +* Compute norm(A - Q * B * P**H) / ( n * norm(A) * EPS ) * ANORM = ZLANGE( '1', M, N, A, LDA, RWORK ) EPS = DLAMCH( 'Precision' ) diff --git a/TESTING/EIG/zbdt02.f b/TESTING/EIG/zbdt02.f index 72fa78a37f..27117a7d1f 100644 --- a/TESTING/EIG/zbdt02.f +++ b/TESTING/EIG/zbdt02.f @@ -27,7 +27,8 @@ *> *> \verbatim *> -*> ZBDT02 tests the change of basis C = U' * B by computing the residual +*> ZBDT02 tests the change of basis C = U**H * B by computing the +*> residual *> *> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), *> @@ -66,7 +67,7 @@ *> \param[in] C *> \verbatim *> C is COMPLEX*16 array, dimension (LDC,N) -*> The m by n matrix C, assumed to contain U' * B. +*> The m by n matrix C, assumed to contain U**H * B. *> \endverbatim *> *> \param[in] LDC diff --git a/TESTING/LIN/cpot01.f b/TESTING/LIN/cpot01.f index e9b19184f2..b80e0253b0 100644 --- a/TESTING/LIN/cpot01.f +++ b/TESTING/LIN/cpot01.f @@ -67,10 +67,10 @@ *> \param[in,out] AFAC *> \verbatim *> AFAC is COMPLEX array, dimension (LDAFAC,N) -*> On entry, the factor L or U from the L*L' or U'*U +*> On entry, the factor L or U from the L * L**H or U**H * U *> factorization of A. *> Overwritten with the reconstructed matrix, and then with the -*> difference L*L' - A (or U'*U - A). +*> difference L * L**H - A (or U**H * U - A). *> \endverbatim *> *> \param[in] LDAFAC @@ -87,8 +87,8 @@ *> \param[out] RESID *> \verbatim *> RESID is REAL -*> If UPLO = 'L', norm(L*L' - A) / ( N * norm(A) * EPS ) -*> If UPLO = 'U', norm(U'*U - A) / ( N * norm(A) * EPS ) +*> If UPLO = 'L', norm(L * L**H - A) / ( N * norm(A) * EPS ) +*> If UPLO = 'U', norm(U**H * U - A) / ( N * norm(A) * EPS ) *> \endverbatim * * Authors: @@ -169,7 +169,7 @@ SUBROUTINE CPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) END IF 10 CONTINUE * -* Compute the product U'*U, overwriting U. +* Compute the product U**H * U, overwriting U. * IF( LSAME( UPLO, 'U' ) ) THEN DO 20 K = N, 1, -1 @@ -186,7 +186,7 @@ SUBROUTINE CPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) * 20 CONTINUE * -* Compute the product L*L', overwriting L. +* Compute the product L * L**H, overwriting L. * ELSE DO 30 K = N, 1, -1 @@ -206,7 +206,7 @@ SUBROUTINE CPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 30 CONTINUE END IF * -* Compute the difference L*L' - A (or U'*U - A). +* Compute the difference L * L**H - A (or U**H * U - A). * IF( LSAME( UPLO, 'U' ) ) THEN DO 50 J = 1, N diff --git a/TESTING/LIN/cqrt17.f b/TESTING/LIN/cqrt17.f index c4e5dd6f09..d38486e3fb 100644 --- a/TESTING/LIN/cqrt17.f +++ b/TESTING/LIN/cqrt17.f @@ -28,9 +28,9 @@ *> *> CQRT17 computes the ratio *> -*> || R'*op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) +*> || R**H * op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) *> -*> where R = op(A)*X - B, op(A) is A or A', and +*> where R = op(A)*X - B, op(A) is A or A**H, and *> *> alpha = ||B|| if IRESID = 1 (zero-residual problem) *> alpha = ||R|| if IRESID = 2 (otherwise). @@ -44,7 +44,7 @@ *> TRANS is CHARACTER*1 *> Specifies whether or not the transpose of A is used. *> = 'N': No transpose, op(A) = A. -*> = 'C': Conjugate transpose, op(A) = A'. +*> = 'C': Conjugate transpose, op(A) = A**H. *> \endverbatim *> *> \param[in] IRESID @@ -224,7 +224,7 @@ REAL FUNCTION CQRT17( TRANS, IRESID, M, N, NRHS, A, $ INFO ) END IF * -* compute R'*op(A) +* compute R**H * op(A) * CALL CGEMM( 'Conjugate transpose', TRANS, NRHS, NCOLS, NROWS, $ CMPLX( ONE ), C, LDB, A, LDA, CMPLX( ZERO ), WORK, diff --git a/TESTING/LIN/ctbt02.f b/TESTING/LIN/ctbt02.f index 8f95c95af3..fc6f2dab7f 100644 --- a/TESTING/LIN/ctbt02.f +++ b/TESTING/LIN/ctbt02.f @@ -204,7 +204,7 @@ SUBROUTINE CTBT02( UPLO, TRANS, DIAG, N, KD, NRHS, AB, LDAB, X, RETURN END IF * -* Compute the 1-norm of A or A'. +* Compute the 1-norm of op(A). * IF( LSAME( TRANS, 'N' ) ) THEN ANORM = CLANTB( '1', UPLO, DIAG, N, KD, AB, LDAB, RWORK ) diff --git a/TESTING/LIN/dgtt02.f b/TESTING/LIN/dgtt02.f index 1cc17215a3..30de915b73 100644 --- a/TESTING/LIN/dgtt02.f +++ b/TESTING/LIN/dgtt02.f @@ -41,8 +41,8 @@ *> TRANS is CHARACTER *> Specifies the form of the residual. *> = 'N': B - A * X (No transpose) -*> = 'T': B - A'* X (Transpose) -*> = 'C': B - A'* X (Conjugate transpose = Transpose) +*> = 'T': B - A**T * X (Transpose) +*> = 'C': B - A**H * X (Conjugate transpose = Transpose) *> \endverbatim *> *> \param[in] N diff --git a/TESTING/LIN/dlarhs.f b/TESTING/LIN/dlarhs.f index b9ce743386..b9e6c789c1 100644 --- a/TESTING/LIN/dlarhs.f +++ b/TESTING/LIN/dlarhs.f @@ -30,7 +30,7 @@ *> DLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system *> op( A ) * X = B, -*> where op( A ) may be A or A' (transpose of A). +*> where op( A ) = A or A**T, depending on TRANS. *> \endverbatim * * Arguments: @@ -81,8 +81,8 @@ *> TRANS is CHARACTER*1 *> Specifies the operation applied to the matrix A. *> = 'N': System is A * x = b -*> = 'T': System is A'* x = b -*> = 'C': System is A'* x = b +*> = 'T': B := A**T * X (Transpose) +*> = 'C': B := A**H * X (Conjugate transpose = Transpose) *> \endverbatim *> *> \param[in] M diff --git a/TESTING/LIN/dpot01.f b/TESTING/LIN/dpot01.f index f48a1a4fc5..049cc81dd3 100644 --- a/TESTING/LIN/dpot01.f +++ b/TESTING/LIN/dpot01.f @@ -65,10 +65,10 @@ *> \param[in,out] AFAC *> \verbatim *> AFAC is DOUBLE PRECISION array, dimension (LDAFAC,N) -*> On entry, the factor L or U from the L*L' or U'*U +*> On entry, the factor L or U from the L * L**T or U**T * U *> factorization of A. *> Overwritten with the reconstructed matrix, and then with the -*> difference L*L' - A (or U'*U - A). +*> difference L * L**T - A (or U**T * U - A). *> \endverbatim *> *> \param[in] LDAFAC @@ -85,8 +85,8 @@ *> \param[out] RESID *> \verbatim *> RESID is DOUBLE PRECISION -*> If UPLO = 'L', norm(L*L' - A) / ( N * norm(A) * EPS ) -*> If UPLO = 'U', norm(U'*U - A) / ( N * norm(A) * EPS ) +*> If UPLO = 'L', norm(L * L**T - A) / ( N * norm(A) * EPS ) +*> If UPLO = 'U', norm(U**T * U - A) / ( N * norm(A) * EPS ) *> \endverbatim * * Authors: @@ -154,7 +154,7 @@ SUBROUTINE DPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) RETURN END IF * -* Compute the product U'*U, overwriting U. +* Compute the product U**T * U, overwriting U. * IF( LSAME( UPLO, 'U' ) ) THEN DO 10 K = N, 1, -1 @@ -171,7 +171,7 @@ SUBROUTINE DPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) * 10 CONTINUE * -* Compute the product L*L', overwriting L. +* Compute the product L * L**T, overwriting L. * ELSE DO 20 K = N, 1, -1 @@ -191,7 +191,7 @@ SUBROUTINE DPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 20 CONTINUE END IF * -* Compute the difference L*L' - A (or U'*U - A). +* Compute the difference L * L**T - A (or U**T * U - A). * IF( LSAME( UPLO, 'U' ) ) THEN DO 40 J = 1, N diff --git a/TESTING/LIN/dqrt17.f b/TESTING/LIN/dqrt17.f index 592904a497..e34ce3381b 100644 --- a/TESTING/LIN/dqrt17.f +++ b/TESTING/LIN/dqrt17.f @@ -28,9 +28,9 @@ *> *> DQRT17 computes the ratio *> -*> || R'*op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) +*> || R**T * op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) *> -*> where R = op(A)*X - B, op(A) is A or A', and +*> where R = op(A)*X - B, op(A) is A or A**T, and *> *> alpha = ||B|| if IRESID = 1 (zero-residual problem) *> alpha = ||R|| if IRESID = 2 (otherwise). @@ -44,7 +44,7 @@ *> TRANS is CHARACTER*1 *> Specifies whether or not the transpose of A is used. *> = 'N': No transpose, op(A) = A. -*> = 'T': Transpose, op(A) = A'. +*> = 'T': Transpose, op(A) = A**T. *> \endverbatim *> *> \param[in] IRESID @@ -225,7 +225,7 @@ DOUBLE PRECISION FUNCTION DQRT17( TRANS, IRESID, M, N, NRHS, A, $ INFO ) END IF * -* compute R'*op(A) +* compute R**T * op(A) * CALL DGEMM( 'Transpose', TRANS, NRHS, NCOLS, NROWS, ONE, C, LDB, $ A, LDA, ZERO, WORK, NRHS ) diff --git a/TESTING/LIN/dtbt02.f b/TESTING/LIN/dtbt02.f index 2870b32677..429800c868 100644 --- a/TESTING/LIN/dtbt02.f +++ b/TESTING/LIN/dtbt02.f @@ -52,8 +52,8 @@ *> TRANS is CHARACTER*1 *> Specifies the operation applied to A. *> = 'N': A *x = b (No transpose) -*> = 'T': A'*x = b (Transpose) -*> = 'C': A'*x = b (Conjugate transpose = Transpose) +*> = 'T': A**T * X = B (Transpose) +*> = 'C': A**H * X = B (Conjugate transpose = Transpose) *> \endverbatim *> *> \param[in] DIAG @@ -196,7 +196,7 @@ SUBROUTINE DTBT02( UPLO, TRANS, DIAG, N, KD, NRHS, AB, LDAB, X, RETURN END IF * -* Compute the 1-norm of A or A'. +* Compute the 1-norm of op(A). * IF( LSAME( TRANS, 'N' ) ) THEN ANORM = DLANTB( '1', UPLO, DIAG, N, KD, AB, LDAB, WORK ) diff --git a/TESTING/LIN/sgtt02.f b/TESTING/LIN/sgtt02.f index 52e2d443e7..0d19b70f24 100644 --- a/TESTING/LIN/sgtt02.f +++ b/TESTING/LIN/sgtt02.f @@ -41,8 +41,8 @@ *> TRANS is CHARACTER *> Specifies the form of the residual. *> = 'N': B - A * X (No transpose) -*> = 'T': B - A'* X (Transpose) -*> = 'C': B - A'* X (Conjugate transpose = Transpose) +*> = 'T': B - A**T * X (Transpose) +*> = 'C': B - A**H * X (Conjugate transpose = Transpose) *> \endverbatim *> *> \param[in] N diff --git a/TESTING/LIN/slarhs.f b/TESTING/LIN/slarhs.f index b126ca72a4..857640e35a 100644 --- a/TESTING/LIN/slarhs.f +++ b/TESTING/LIN/slarhs.f @@ -30,7 +30,7 @@ *> SLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system *> op( A ) * X = B, -*> where op( A ) may be A or A' (transpose of A). +*> where op(A) = A or A**T. *> \endverbatim * * Arguments: @@ -81,8 +81,8 @@ *> TRANS is CHARACTER*1 *> Specifies the operation applied to the matrix A. *> = 'N': System is A * x = b -*> = 'T': System is A'* x = b -*> = 'C': System is A'* x = b +*> = 'T': B := A**T * X (Transpose) +*> = 'C': B := A**H * X (Conjugate transpose = Transpose) *> \endverbatim *> *> \param[in] M diff --git a/TESTING/LIN/spot01.f b/TESTING/LIN/spot01.f index a0176c822e..7e2a7ca30e 100644 --- a/TESTING/LIN/spot01.f +++ b/TESTING/LIN/spot01.f @@ -65,10 +65,10 @@ *> \param[in,out] AFAC *> \verbatim *> AFAC is REAL array, dimension (LDAFAC,N) -*> On entry, the factor L or U from the L*L' or U'*U +*> On entry, the factor L or U from the L * L**T or U**T * U *> factorization of A. *> Overwritten with the reconstructed matrix, and then with the -*> difference L*L' - A (or U'*U - A). +*> difference L * L**T - A (or U**T * U - A). *> \endverbatim *> *> \param[in] LDAFAC @@ -85,8 +85,8 @@ *> \param[out] RESID *> \verbatim *> RESID is REAL -*> If UPLO = 'L', norm(L*L' - A) / ( N * norm(A) * EPS ) -*> If UPLO = 'U', norm(U'*U - A) / ( N * norm(A) * EPS ) +*> If UPLO = 'L', norm(L * L**T - A) / ( N * norm(A) * EPS ) +*> If UPLO = 'U', norm(U**T * U - A) / ( N * norm(A) * EPS ) *> \endverbatim * * Authors: @@ -154,7 +154,7 @@ SUBROUTINE SPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) RETURN END IF * -* Compute the product U'*U, overwriting U. +* Compute the product U**T * U, overwriting U. * IF( LSAME( UPLO, 'U' ) ) THEN DO 10 K = N, 1, -1 @@ -171,7 +171,7 @@ SUBROUTINE SPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) * 10 CONTINUE * -* Compute the product L*L', overwriting L. +* Compute the product L * L**T, overwriting L. * ELSE DO 20 K = N, 1, -1 @@ -191,7 +191,7 @@ SUBROUTINE SPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 20 CONTINUE END IF * -* Compute the difference L*L' - A (or U'*U - A). +* Compute the difference L * L**T - A (or U**T * U - A). * IF( LSAME( UPLO, 'U' ) ) THEN DO 40 J = 1, N diff --git a/TESTING/LIN/sqrt17.f b/TESTING/LIN/sqrt17.f index 1a80b8c784..fe2914aba3 100644 --- a/TESTING/LIN/sqrt17.f +++ b/TESTING/LIN/sqrt17.f @@ -28,9 +28,9 @@ *> *> SQRT17 computes the ratio *> -*> || R'*op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) +*> || R**T * op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) *> -*> where R = op(A)*X - B, op(A) is A or A', and +*> where R = op(A)*X - B, op(A) is A or A**T, and *> *> alpha = ||B|| if IRESID = 1 (zero-residual problem) *> alpha = ||R|| if IRESID = 2 (otherwise). @@ -44,7 +44,7 @@ *> TRANS is CHARACTER*1 *> Specifies whether or not the transpose of A is used. *> = 'N': No transpose, op(A) = A. -*> = 'T': Transpose, op(A) = A'. +*> = 'T': Transpose, op(A) = A**T. *> \endverbatim *> *> \param[in] IRESID @@ -225,7 +225,7 @@ REAL FUNCTION SQRT17( TRANS, IRESID, M, N, NRHS, A, $ INFO ) END IF * -* compute R'*op(A) +* compute R**T * op(A) * CALL SGEMM( 'Transpose', TRANS, NRHS, NCOLS, NROWS, ONE, C, LDB, $ A, LDA, ZERO, WORK, NRHS ) diff --git a/TESTING/LIN/stbt02.f b/TESTING/LIN/stbt02.f index 96ca3b6baa..3b602112d3 100644 --- a/TESTING/LIN/stbt02.f +++ b/TESTING/LIN/stbt02.f @@ -28,10 +28,9 @@ *> \verbatim *> *> STBT02 computes the residual for the computed solution to a -*> triangular system of linear equations A*x = b or A' *x = b when -*> A is a triangular band matrix. Here A' is the transpose of A and -*> x and b are N by NRHS matrices. The test ratio is the maximum over -*> the number of right hand sides of +*> triangular system of linear equations op(A)*X = B when A is a +*> triangular band matrix. The test ratio is the maximum over the +*> number of right hand sides of *> norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ), *> where op(A) denotes A or A' and EPS is the machine epsilon. *> \endverbatim @@ -52,8 +51,8 @@ *> TRANS is CHARACTER*1 *> Specifies the operation applied to A. *> = 'N': A *x = b (No transpose) -*> = 'T': A'*x = b (Transpose) -*> = 'C': A'*x = b (Conjugate transpose = Transpose) +*> = 'T': A**T * X = B (Transpose) +*> = 'C': A**H * X = B (Conjugate transpose = Transpose) *> \endverbatim *> *> \param[in] DIAG @@ -196,7 +195,7 @@ SUBROUTINE STBT02( UPLO, TRANS, DIAG, N, KD, NRHS, AB, LDAB, X, RETURN END IF * -* Compute the 1-norm of A or A'. +* Compute the 1-norm of op(A). * IF( LSAME( TRANS, 'N' ) ) THEN ANORM = SLANTB( '1', UPLO, DIAG, N, KD, AB, LDAB, WORK ) diff --git a/TESTING/LIN/zpot01.f b/TESTING/LIN/zpot01.f index c08a27bf71..5a89f33169 100644 --- a/TESTING/LIN/zpot01.f +++ b/TESTING/LIN/zpot01.f @@ -67,10 +67,10 @@ *> \param[in,out] AFAC *> \verbatim *> AFAC is COMPLEX*16 array, dimension (LDAFAC,N) -*> On entry, the factor L or U from the L*L' or U'*U +*> On entry, the factor L or U from the L * L**H or U**H * U *> factorization of A. *> Overwritten with the reconstructed matrix, and then with the -*> difference L*L' - A (or U'*U - A). +*> the difference L * L**H - A (or U**H * U - A). *> \endverbatim *> *> \param[in] LDAFAC @@ -87,8 +87,8 @@ *> \param[out] RESID *> \verbatim *> RESID is DOUBLE PRECISION -*> If UPLO = 'L', norm(L*L' - A) / ( N * norm(A) * EPS ) -*> If UPLO = 'U', norm(U'*U - A) / ( N * norm(A) * EPS ) +*> If UPLO = 'L', norm(L * L**H - A) / ( N * norm(A) * EPS ) +*> If UPLO = 'U', norm(U**H * U - A) / ( N * norm(A) * EPS ) *> \endverbatim * * Authors: @@ -169,7 +169,7 @@ SUBROUTINE ZPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) END IF 10 CONTINUE * -* Compute the product U'*U, overwriting U. +* Compute the product U**H * U, overwriting U. * IF( LSAME( UPLO, 'U' ) ) THEN DO 20 K = N, 1, -1 @@ -186,7 +186,7 @@ SUBROUTINE ZPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) * 20 CONTINUE * -* Compute the product L*L', overwriting L. +* Compute the product L * L**H, overwriting L. * ELSE DO 30 K = N, 1, -1 @@ -206,7 +206,7 @@ SUBROUTINE ZPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 30 CONTINUE END IF * -* Compute the difference L*L' - A (or U'*U - A). +* Compute the difference L * L**H - A (or U**H * U - A). * IF( LSAME( UPLO, 'U' ) ) THEN DO 50 J = 1, N diff --git a/TESTING/LIN/zqrt17.f b/TESTING/LIN/zqrt17.f index 77dd77b417..42c658b293 100644 --- a/TESTING/LIN/zqrt17.f +++ b/TESTING/LIN/zqrt17.f @@ -28,9 +28,9 @@ *> *> ZQRT17 computes the ratio *> -*> || R'*op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) +*> || R**H * op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) *> -*> where R = op(A)*X - B, op(A) is A or A', and +*> where R = op(A)*X - B, op(A) is A or A**H, and *> *> alpha = ||B|| if IRESID = 1 (zero-residual problem) *> alpha = ||R|| if IRESID = 2 (otherwise). @@ -44,7 +44,7 @@ *> TRANS is CHARACTER*1 *> Specifies whether or not the transpose of A is used. *> = 'N': No transpose, op(A) = A. -*> = 'C': Conjugate transpose, op(A) = A'. +*> = 'C': Conjugate transpose, op(A) = A**H. *> \endverbatim *> *> \param[in] IRESID @@ -225,7 +225,7 @@ DOUBLE PRECISION FUNCTION ZQRT17( TRANS, IRESID, M, N, NRHS, A, $ INFO ) END IF * -* compute R'*op(A) +* compute R**H * op(A) * CALL ZGEMM( 'Conjugate transpose', TRANS, NRHS, NCOLS, NROWS, $ DCMPLX( ONE ), C, LDB, A, LDA, DCMPLX( ZERO ), WORK, diff --git a/TESTING/LIN/ztbt02.f b/TESTING/LIN/ztbt02.f index fdfa35e6be..2553d9eaca 100644 --- a/TESTING/LIN/ztbt02.f +++ b/TESTING/LIN/ztbt02.f @@ -204,7 +204,7 @@ SUBROUTINE ZTBT02( UPLO, TRANS, DIAG, N, KD, NRHS, AB, LDAB, X, RETURN END IF * -* Compute the 1-norm of A or A'. +* Compute the 1-norm of op(A). * IF( LSAME( TRANS, 'N' ) ) THEN ANORM = ZLANTB( '1', UPLO, DIAG, N, KD, AB, LDAB, RWORK ) From 76bcd0cd4de7b2140a4787e7b926a80366d61a9f Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Wed, 2 Jun 2021 21:56:03 +1000 Subject: [PATCH 03/16] comments unified: remove spaces around expr in norm( expr ) --- TESTING/EIG/cbdt01.f | 2 +- TESTING/EIG/cbdt02.f | 6 +++--- TESTING/EIG/dbdt01.f | 2 +- TESTING/EIG/dbdt02.f | 6 +++--- TESTING/EIG/sbdt01.f | 2 +- TESTING/EIG/sbdt02.f | 6 +++--- TESTING/EIG/zbdt01.f | 2 +- TESTING/EIG/zbdt02.f | 6 +++--- TESTING/LIN/cgbt01.f | 2 +- TESTING/LIN/cpot01.f | 2 +- TESTING/LIN/dgbt01.f | 2 +- TESTING/LIN/dpot01.f | 2 +- TESTING/LIN/sgbt01.f | 2 +- TESTING/LIN/spot01.f | 2 +- TESTING/LIN/zgbt01.f | 2 +- TESTING/LIN/zpot01.f | 2 +- 16 files changed, 24 insertions(+), 24 deletions(-) diff --git a/TESTING/EIG/cbdt01.f b/TESTING/EIG/cbdt01.f index c1bfbffa53..2f678ef233 100644 --- a/TESTING/EIG/cbdt01.f +++ b/TESTING/EIG/cbdt01.f @@ -33,7 +33,7 @@ *> matrices and B is bidiagonal. *> *> The test ratio to test the reduction is -*> RESID = norm( A - Q * B * P**H ) / ( n * norm(A) * EPS ) +*> RESID = norm(A - Q * B * P**H) / ( n * norm(A) * EPS ) *> where EPS is the machine precision. *> \endverbatim * diff --git a/TESTING/EIG/cbdt02.f b/TESTING/EIG/cbdt02.f index 8aafefd227..6ef39a96b3 100644 --- a/TESTING/EIG/cbdt02.f +++ b/TESTING/EIG/cbdt02.f @@ -30,7 +30,7 @@ *> CBDT02 tests the change of basis C = U**H * B by computing the *> residual *> -*> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), +*> RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ), *> *> where B and C are M by N matrices, U is an M by M orthogonal matrix, *> and EPS is the machine precision. @@ -101,7 +101,7 @@ *> \param[out] RESID *> \verbatim *> RESID is REAL -*> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), +*> RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ), *> \endverbatim * * Authors: @@ -162,7 +162,7 @@ SUBROUTINE CBDT02( M, N, B, LDB, C, LDC, U, LDU, WORK, RWORK, REALMN = REAL( MAX( M, N ) ) EPS = SLAMCH( 'Precision' ) * -* Compute norm( B - U * C ) +* Compute norm(B - U * C) * DO 10 J = 1, N CALL CCOPY( M, B( 1, J ), 1, WORK, 1 ) diff --git a/TESTING/EIG/dbdt01.f b/TESTING/EIG/dbdt01.f index f978493ea4..69e0b5ffc2 100644 --- a/TESTING/EIG/dbdt01.f +++ b/TESTING/EIG/dbdt01.f @@ -32,7 +32,7 @@ *> matrices and B is bidiagonal. *> *> The test ratio to test the reduction is -*> RESID = norm( A - Q * B * P**T ) / ( n * norm(A) * EPS ) +*> RESID = norm(A - Q * B * P**T) / ( n * norm(A) * EPS ) *> where EPS is the machine precision. *> \endverbatim * diff --git a/TESTING/EIG/dbdt02.f b/TESTING/EIG/dbdt02.f index 76e3ff713d..ac0291219b 100644 --- a/TESTING/EIG/dbdt02.f +++ b/TESTING/EIG/dbdt02.f @@ -28,7 +28,7 @@ *> DBDT02 tests the change of basis C = U**H * B by computing the *> residual *> -*> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), +*> RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ), *> *> where B and C are M by N matrices, U is an M by M orthogonal matrix, *> and EPS is the machine precision. @@ -94,7 +94,7 @@ *> \param[out] RESID *> \verbatim *> RESID is DOUBLE PRECISION -*> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), +*> RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ), *> \endverbatim * * Authors: @@ -153,7 +153,7 @@ SUBROUTINE DBDT02( M, N, B, LDB, C, LDC, U, LDU, WORK, RESID ) REALMN = DBLE( MAX( M, N ) ) EPS = DLAMCH( 'Precision' ) * -* Compute norm( B - U * C ) +* Compute norm(B - U * C) * DO 10 J = 1, N CALL DCOPY( M, B( 1, J ), 1, WORK, 1 ) diff --git a/TESTING/EIG/sbdt01.f b/TESTING/EIG/sbdt01.f index 0861818506..fdc30a51ea 100644 --- a/TESTING/EIG/sbdt01.f +++ b/TESTING/EIG/sbdt01.f @@ -32,7 +32,7 @@ *> matrices and B is bidiagonal. *> *> The test ratio to test the reduction is -*> RESID = norm( A - Q * B * P**T ) / ( n * norm(A) * EPS ) +*> RESID = norm(A - Q * B * P**T) / ( n * norm(A) * EPS ) *> where EPS is the machine precision. *> \endverbatim * diff --git a/TESTING/EIG/sbdt02.f b/TESTING/EIG/sbdt02.f index 1cf8a22725..e86863f78b 100644 --- a/TESTING/EIG/sbdt02.f +++ b/TESTING/EIG/sbdt02.f @@ -28,7 +28,7 @@ *> SBDT02 tests the change of basis C = U**H * B by computing the *> residual *> -*> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), +*> RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ), *> *> where B and C are M by N matrices, U is an M by M orthogonal matrix, *> and EPS is the machine precision. @@ -94,7 +94,7 @@ *> \param[out] RESID *> \verbatim *> RESID is REAL -*> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), +*> RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ), *> \endverbatim * * Authors: @@ -153,7 +153,7 @@ SUBROUTINE SBDT02( M, N, B, LDB, C, LDC, U, LDU, WORK, RESID ) REALMN = REAL( MAX( M, N ) ) EPS = SLAMCH( 'Precision' ) * -* Compute norm( B - U * C ) +* Compute norm(B - U * C) * DO 10 J = 1, N CALL SCOPY( M, B( 1, J ), 1, WORK, 1 ) diff --git a/TESTING/EIG/zbdt01.f b/TESTING/EIG/zbdt01.f index cb62993ade..afe2fa89a5 100644 --- a/TESTING/EIG/zbdt01.f +++ b/TESTING/EIG/zbdt01.f @@ -33,7 +33,7 @@ *> matrices and B is bidiagonal. *> *> The test ratio to test the reduction is -*> RESID = norm( A - Q * B * P**H ) / ( n * norm(A) * EPS ) +*> RESID = norm(A - Q * B * P**H) / ( n * norm(A) * EPS ) *> where EPS is the machine precision. *> \endverbatim * diff --git a/TESTING/EIG/zbdt02.f b/TESTING/EIG/zbdt02.f index 27117a7d1f..f0da37222b 100644 --- a/TESTING/EIG/zbdt02.f +++ b/TESTING/EIG/zbdt02.f @@ -30,7 +30,7 @@ *> ZBDT02 tests the change of basis C = U**H * B by computing the *> residual *> -*> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), +*> RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ), *> *> where B and C are M by N matrices, U is an M by M orthogonal matrix, *> and EPS is the machine precision. @@ -101,7 +101,7 @@ *> \param[out] RESID *> \verbatim *> RESID is DOUBLE PRECISION -*> RESID = norm( B - U * C ) / ( max(m,n) * norm(B) * EPS ), +*> RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ), *> \endverbatim * * Authors: @@ -162,7 +162,7 @@ SUBROUTINE ZBDT02( M, N, B, LDB, C, LDC, U, LDU, WORK, RWORK, REALMN = DBLE( MAX( M, N ) ) EPS = DLAMCH( 'Precision' ) * -* Compute norm( B - U * C ) +* Compute norm(B - U * C) * DO 10 J = 1, N CALL ZCOPY( M, B( 1, J ), 1, WORK, 1 ) diff --git a/TESTING/LIN/cgbt01.f b/TESTING/LIN/cgbt01.f index 277b0b1ab8..375fbc1e3c 100644 --- a/TESTING/LIN/cgbt01.f +++ b/TESTING/LIN/cgbt01.f @@ -226,7 +226,7 @@ SUBROUTINE CGBT01( M, N, KL, KU, A, LDA, AFAC, LDAFAC, IPIV, WORK, END IF 40 CONTINUE * -* Compute norm( L*U - A ) / ( N * norm(A) * EPS ) +* Compute norm(L*U - A) / ( N * norm(A) * EPS ) * IF( ANORM.LE.ZERO ) THEN IF( RESID.NE.ZERO ) diff --git a/TESTING/LIN/cpot01.f b/TESTING/LIN/cpot01.f index b80e0253b0..30dc55d884 100644 --- a/TESTING/LIN/cpot01.f +++ b/TESTING/LIN/cpot01.f @@ -224,7 +224,7 @@ SUBROUTINE CPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 70 CONTINUE END IF * -* Compute norm( L*U - A ) / ( N * norm(A) * EPS ) +* Compute norm(L*U - A) / ( N * norm(A) * EPS ) * RESID = CLANHE( '1', UPLO, N, AFAC, LDAFAC, RWORK ) * diff --git a/TESTING/LIN/dgbt01.f b/TESTING/LIN/dgbt01.f index 28f7200ddb..bf0e77a4fa 100644 --- a/TESTING/LIN/dgbt01.f +++ b/TESTING/LIN/dgbt01.f @@ -225,7 +225,7 @@ SUBROUTINE DGBT01( M, N, KL, KU, A, LDA, AFAC, LDAFAC, IPIV, WORK, END IF 40 CONTINUE * -* Compute norm( L*U - A ) / ( N * norm(A) * EPS ) +* Compute norm(L*U - A) / ( N * norm(A) * EPS ) * IF( ANORM.LE.ZERO ) THEN IF( RESID.NE.ZERO ) diff --git a/TESTING/LIN/dpot01.f b/TESTING/LIN/dpot01.f index 049cc81dd3..5c463c12c2 100644 --- a/TESTING/LIN/dpot01.f +++ b/TESTING/LIN/dpot01.f @@ -207,7 +207,7 @@ SUBROUTINE DPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 60 CONTINUE END IF * -* Compute norm( L*U - A ) / ( N * norm(A) * EPS ) +* Compute norm(L*U - A) / ( N * norm(A) * EPS ) * RESID = DLANSY( '1', UPLO, N, AFAC, LDAFAC, RWORK ) * diff --git a/TESTING/LIN/sgbt01.f b/TESTING/LIN/sgbt01.f index 0f6129240a..9abcf56e28 100644 --- a/TESTING/LIN/sgbt01.f +++ b/TESTING/LIN/sgbt01.f @@ -225,7 +225,7 @@ SUBROUTINE SGBT01( M, N, KL, KU, A, LDA, AFAC, LDAFAC, IPIV, WORK, END IF 40 CONTINUE * -* Compute norm( L*U - A ) / ( N * norm(A) * EPS ) +* Compute norm(L*U - A) / ( N * norm(A) * EPS ) * IF( ANORM.LE.ZERO ) THEN IF( RESID.NE.ZERO ) diff --git a/TESTING/LIN/spot01.f b/TESTING/LIN/spot01.f index 7e2a7ca30e..05c85815b7 100644 --- a/TESTING/LIN/spot01.f +++ b/TESTING/LIN/spot01.f @@ -207,7 +207,7 @@ SUBROUTINE SPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 60 CONTINUE END IF * -* Compute norm( L*U - A ) / ( N * norm(A) * EPS ) +* Compute norm(L*U - A) / ( N * norm(A) * EPS ) * RESID = SLANSY( '1', UPLO, N, AFAC, LDAFAC, RWORK ) * diff --git a/TESTING/LIN/zgbt01.f b/TESTING/LIN/zgbt01.f index 9219193713..67efa58f8e 100644 --- a/TESTING/LIN/zgbt01.f +++ b/TESTING/LIN/zgbt01.f @@ -226,7 +226,7 @@ SUBROUTINE ZGBT01( M, N, KL, KU, A, LDA, AFAC, LDAFAC, IPIV, WORK, END IF 40 CONTINUE * -* Compute norm( L*U - A ) / ( N * norm(A) * EPS ) +* Compute norm(L*U - A) / ( N * norm(A) * EPS ) * IF( ANORM.LE.ZERO ) THEN IF( RESID.NE.ZERO ) diff --git a/TESTING/LIN/zpot01.f b/TESTING/LIN/zpot01.f index 5a89f33169..2c48ff001f 100644 --- a/TESTING/LIN/zpot01.f +++ b/TESTING/LIN/zpot01.f @@ -224,7 +224,7 @@ SUBROUTINE ZPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 70 CONTINUE END IF * -* Compute norm( L*U - A ) / ( N * norm(A) * EPS ) +* Compute norm(L*U - A) / ( N * norm(A) * EPS ) * RESID = ZLANHE( '1', UPLO, N, AFAC, LDAFAC, RWORK ) * From bf04fe9cd862be0cfa5a58c9af3119df66114ef1 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Wed, 2 Jun 2021 22:03:28 +1000 Subject: [PATCH 04/16] comments unified: remove spaces around expr in op( expr ) in xLARHS --- TESTING/EIG/clarhs.f | 6 +++--- TESTING/EIG/dlarhs.f | 6 +++--- TESTING/EIG/slarhs.f | 6 +++--- TESTING/EIG/zlarhs.f | 6 +++--- TESTING/LIN/clarhs.f | 6 +++--- TESTING/LIN/dlarhs.f | 6 +++--- TESTING/LIN/slarhs.f | 4 ++-- TESTING/LIN/zlarhs.f | 6 +++--- 8 files changed, 23 insertions(+), 23 deletions(-) diff --git a/TESTING/EIG/clarhs.f b/TESTING/EIG/clarhs.f index e1d7ebdddc..bba91f8913 100644 --- a/TESTING/EIG/clarhs.f +++ b/TESTING/EIG/clarhs.f @@ -29,8 +29,8 @@ *> *> CLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system -*> op( A ) * X = B, -*> where op( A ) may be A, A**T (transpose of A), or A**H (conjugate +*> op(A) * X = B, +*> where op(A) = A, A**T (transpose of A) or A**H (conjugate *> transpose of A). *> \endverbatim * @@ -313,7 +313,7 @@ SUBROUTINE CLARHS( PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, 10 CONTINUE END IF * -* Multiply X by op( A ) using an appropriate +* Multiply X by op(A) using an appropriate * matrix multiply routine. * IF( LSAMEN( 2, C2, 'GE' ) .OR. LSAMEN( 2, C2, 'QR' ) .OR. diff --git a/TESTING/EIG/dlarhs.f b/TESTING/EIG/dlarhs.f index a5ba318bed..f391fafc0b 100644 --- a/TESTING/EIG/dlarhs.f +++ b/TESTING/EIG/dlarhs.f @@ -29,8 +29,8 @@ *> *> DLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system -*> op( A ) * X = B, -*> where op( A ) = A or A**T, depending on TRANS. +*> op(A) * X = B, +*> where op(A) = A or A**T, depending on TRANS. *> \endverbatim * * Arguments: @@ -305,7 +305,7 @@ SUBROUTINE DLARHS( PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, 10 CONTINUE END IF * -* Multiply X by op( A ) using an appropriate +* Multiply X by op(A) using an appropriate * matrix multiply routine. * IF( LSAMEN( 2, C2, 'GE' ) .OR. LSAMEN( 2, C2, 'QR' ) .OR. diff --git a/TESTING/EIG/slarhs.f b/TESTING/EIG/slarhs.f index 0e9d309479..c32c1bf24e 100644 --- a/TESTING/EIG/slarhs.f +++ b/TESTING/EIG/slarhs.f @@ -29,8 +29,8 @@ *> *> SLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system -*> op( A ) * X = B, -*> where op( A ) = A or A**T, depending on TRANS. +*> op(A) * X = B, +*> where op(A) = A or A**T, depending on TRANS. *> \endverbatim * * Arguments: @@ -305,7 +305,7 @@ SUBROUTINE SLARHS( PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, 10 CONTINUE END IF * -* Multiply X by op( A ) using an appropriate +* Multiply X by op(A) using an appropriate * matrix multiply routine. * IF( LSAMEN( 2, C2, 'GE' ) .OR. LSAMEN( 2, C2, 'QR' ) .OR. diff --git a/TESTING/EIG/zlarhs.f b/TESTING/EIG/zlarhs.f index e3107f231e..d9ab7cc9fd 100644 --- a/TESTING/EIG/zlarhs.f +++ b/TESTING/EIG/zlarhs.f @@ -29,8 +29,8 @@ *> *> ZLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system -*> op( A ) * X = B, -*> where op( A ) may be A, A**T (transpose of A), or A**H (conjugate +*> op(A) * X = B, +*> where op(A) = A, A**T (transpose of A), or A**H (conjugate *> transpose of A). *> \endverbatim * @@ -313,7 +313,7 @@ SUBROUTINE ZLARHS( PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, 10 CONTINUE END IF * -* Multiply X by op( A ) using an appropriate +* Multiply X by op(A) using an appropriate * matrix multiply routine. * IF( LSAMEN( 2, C2, 'GE' ) .OR. LSAMEN( 2, C2, 'QR' ) .OR. diff --git a/TESTING/LIN/clarhs.f b/TESTING/LIN/clarhs.f index 0052d3667e..fcb196e1bf 100644 --- a/TESTING/LIN/clarhs.f +++ b/TESTING/LIN/clarhs.f @@ -29,8 +29,8 @@ *> *> CLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system -*> op( A ) * X = B, -*> where op( A ) may be A, A**T (transpose of A), or A**H (conjugate +*> op(A) * X = B, +*> where op(A) = A, A**T (transpose of A), or A**H (conjugate *> transpose of A). *> \endverbatim * @@ -313,7 +313,7 @@ SUBROUTINE CLARHS( PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, 10 CONTINUE END IF * -* Multiply X by op( A ) using an appropriate +* Multiply X by op(A) using an appropriate * matrix multiply routine. * IF( LSAMEN( 2, C2, 'GE' ) .OR. LSAMEN( 2, C2, 'QR' ) .OR. diff --git a/TESTING/LIN/dlarhs.f b/TESTING/LIN/dlarhs.f index b9e6c789c1..8e63387e8c 100644 --- a/TESTING/LIN/dlarhs.f +++ b/TESTING/LIN/dlarhs.f @@ -29,8 +29,8 @@ *> *> DLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system -*> op( A ) * X = B, -*> where op( A ) = A or A**T, depending on TRANS. +*> op(A) * X = B, +*> where op(A) = A or A**T, depending on TRANS. *> \endverbatim * * Arguments: @@ -305,7 +305,7 @@ SUBROUTINE DLARHS( PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, 10 CONTINUE END IF * -* Multiply X by op( A ) using an appropriate +* Multiply X by op(A) using an appropriate * matrix multiply routine. * IF( LSAMEN( 2, C2, 'GE' ) .OR. LSAMEN( 2, C2, 'QR' ) .OR. diff --git a/TESTING/LIN/slarhs.f b/TESTING/LIN/slarhs.f index 857640e35a..2dc3361f5e 100644 --- a/TESTING/LIN/slarhs.f +++ b/TESTING/LIN/slarhs.f @@ -29,7 +29,7 @@ *> *> SLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system -*> op( A ) * X = B, +*> op(A) * X = B, *> where op(A) = A or A**T. *> \endverbatim * @@ -305,7 +305,7 @@ SUBROUTINE SLARHS( PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, 10 CONTINUE END IF * -* Multiply X by op( A ) using an appropriate +* Multiply X by op(A) using an appropriate * matrix multiply routine. * IF( LSAMEN( 2, C2, 'GE' ) .OR. LSAMEN( 2, C2, 'QR' ) .OR. diff --git a/TESTING/LIN/zlarhs.f b/TESTING/LIN/zlarhs.f index 92d746fc45..87e463c4d6 100644 --- a/TESTING/LIN/zlarhs.f +++ b/TESTING/LIN/zlarhs.f @@ -29,8 +29,8 @@ *> *> ZLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system -*> op( A ) * X = B, -*> where op( A ) may be A, A**T (transpose of A), or A**H (conjugate +*> op(A) * X = B, +*> where op(A) = A, A**T (transpose of A), or A**H (conjugate *> transpose of A). *> \endverbatim * @@ -313,7 +313,7 @@ SUBROUTINE ZLARHS( PATH, XTYPE, UPLO, TRANS, M, N, KL, KU, NRHS, 10 CONTINUE END IF * -* Multiply X by op( A ) using an appropriate +* Multiply X by op(A) using an appropriate * matrix multiply routine. * IF( LSAMEN( 2, C2, 'GE' ) .OR. LSAMEN( 2, C2, 'QR' ) .OR. From 2f9f42bf2c703e9ed690e9dd4f19300f98cec970 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Wed, 2 Jun 2021 22:15:53 +1000 Subject: [PATCH 05/16] comments unified and enhanced: TRANS parameter --- TESTING/EIG/clarhs.f | 6 +++--- TESTING/EIG/dlarhs.f | 5 +++-- TESTING/EIG/slarhs.f | 5 +++-- TESTING/EIG/zlarhs.f | 6 +++--- TESTING/LIN/cgbt05.f | 2 +- TESTING/LIN/cgtt02.f | 2 +- TESTING/LIN/clarhs.f | 6 +++--- TESTING/LIN/ctbt02.f | 6 +++--- TESTING/LIN/dgbt05.f | 2 +- TESTING/LIN/dgtt02.f | 2 +- TESTING/LIN/dlarhs.f | 5 +++-- TESTING/LIN/dtbt02.f | 2 +- TESTING/LIN/sgbt05.f | 2 +- TESTING/LIN/sgtt02.f | 2 +- TESTING/LIN/slarhs.f | 5 +++-- TESTING/LIN/stbt02.f | 2 +- TESTING/LIN/zgbt05.f | 2 +- TESTING/LIN/zgtt02.f | 2 +- TESTING/LIN/zlarhs.f | 6 +++--- TESTING/LIN/ztbt02.f | 6 +++--- 20 files changed, 40 insertions(+), 36 deletions(-) diff --git a/TESTING/EIG/clarhs.f b/TESTING/EIG/clarhs.f index bba91f8913..4505ab34bd 100644 --- a/TESTING/EIG/clarhs.f +++ b/TESTING/EIG/clarhs.f @@ -85,9 +85,9 @@ *> TRANS is CHARACTER*1 *> Used only if A is nonsymmetric; specifies the operation *> applied to the matrix A. -*> = 'N': B := A * X -*> = 'T': B := A**T * X -*> = 'C': B := A**H * X +*> = 'N': B := A * X (No transpose) +*> = 'T': B := A**T * X (Transpose) +*> = 'C': B := A**H * X (Conjugate transpose) *> \endverbatim *> *> \param[in] M diff --git a/TESTING/EIG/dlarhs.f b/TESTING/EIG/dlarhs.f index f391fafc0b..bc048bbcfe 100644 --- a/TESTING/EIG/dlarhs.f +++ b/TESTING/EIG/dlarhs.f @@ -79,8 +79,9 @@ *> \param[in] TRANS *> \verbatim *> TRANS is CHARACTER*1 -*> Specifies the operation applied to the matrix A. -*> = 'N': System is A * x = b +*> Used only if A is nonsymmetric; specifies the operation +*> applied to the matrix A. +*> = 'N': B := A * X (No transpose) *> = 'T': B := A**T * X (Transpose) *> = 'C': B := A**H * X (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/EIG/slarhs.f b/TESTING/EIG/slarhs.f index c32c1bf24e..bdb74724e8 100644 --- a/TESTING/EIG/slarhs.f +++ b/TESTING/EIG/slarhs.f @@ -79,8 +79,9 @@ *> \param[in] TRANS *> \verbatim *> TRANS is CHARACTER*1 -*> Specifies the operation applied to the matrix A. -*> = 'N': System is A * x = b +*> Used only if A is nonsymmetric; specifies the operation +*> applied to the matrix A. +*> = 'N': B := A * X (No transpose) *> = 'T': B := A**T * X (Transpose) *> = 'C': B := A**H * X (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/EIG/zlarhs.f b/TESTING/EIG/zlarhs.f index d9ab7cc9fd..788e033289 100644 --- a/TESTING/EIG/zlarhs.f +++ b/TESTING/EIG/zlarhs.f @@ -85,9 +85,9 @@ *> TRANS is CHARACTER*1 *> Used only if A is nonsymmetric; specifies the operation *> applied to the matrix A. -*> = 'N': B := A * X -*> = 'T': B := A**T * X -*> = 'C': B := A**H * X +*> = 'N': B := A * X (No transpose) +*> = 'T': B := A**T * X (Transpose) +*> = 'C': B := A**H * X (Conjugate transpose) *> \endverbatim *> *> \param[in] M diff --git a/TESTING/LIN/cgbt05.f b/TESTING/LIN/cgbt05.f index 859f8b3eb2..88fc405699 100644 --- a/TESTING/LIN/cgbt05.f +++ b/TESTING/LIN/cgbt05.f @@ -50,7 +50,7 @@ *> \verbatim *> TRANS is CHARACTER*1 *> Specifies the form of the system of equations. -*> = 'N': A * X = B (No transpose) +*> = 'N': A * X = B (No transpose) *> = 'T': A**T * X = B (Transpose) *> = 'C': A**H * X = B (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/cgtt02.f b/TESTING/LIN/cgtt02.f index abfdfb75a8..587c43a6e5 100644 --- a/TESTING/LIN/cgtt02.f +++ b/TESTING/LIN/cgtt02.f @@ -40,7 +40,7 @@ *> \verbatim *> TRANS is CHARACTER *> Specifies the form of the residual. -*> = 'N': B - A * X (No transpose) +*> = 'N': B - A * X (No transpose) *> = 'T': B - A**T * X (Transpose) *> = 'C': B - A**H * X (Conjugate transpose) *> \endverbatim diff --git a/TESTING/LIN/clarhs.f b/TESTING/LIN/clarhs.f index fcb196e1bf..684d4489cc 100644 --- a/TESTING/LIN/clarhs.f +++ b/TESTING/LIN/clarhs.f @@ -85,9 +85,9 @@ *> TRANS is CHARACTER*1 *> Used only if A is nonsymmetric; specifies the operation *> applied to the matrix A. -*> = 'N': B := A * X -*> = 'T': B := A**T * X -*> = 'C': B := A**H * X +*> = 'N': B := A * X (No transpose) +*> = 'T': B := A**T * X (Transpose) +*> = 'C': B := A**H * X (Conjugate transpose) *> \endverbatim *> *> \param[in] M diff --git a/TESTING/LIN/ctbt02.f b/TESTING/LIN/ctbt02.f index fc6f2dab7f..f2a00a65b4 100644 --- a/TESTING/LIN/ctbt02.f +++ b/TESTING/LIN/ctbt02.f @@ -53,9 +53,9 @@ *> \verbatim *> TRANS is CHARACTER*1 *> Specifies the operation applied to A. -*> = 'N': A *x = b (No transpose) -*> = 'T': A**T *x = b (Transpose) -*> = 'C': A**H *x = b (Conjugate transpose) +*> = 'N': A * X = B (No transpose) +*> = 'T': A**T * X = B (Transpose) +*> = 'C': A**H * X = B (Conjugate transpose) *> \endverbatim *> *> \param[in] DIAG diff --git a/TESTING/LIN/dgbt05.f b/TESTING/LIN/dgbt05.f index ba85eaf7d9..793e73951d 100644 --- a/TESTING/LIN/dgbt05.f +++ b/TESTING/LIN/dgbt05.f @@ -50,7 +50,7 @@ *> \verbatim *> TRANS is CHARACTER*1 *> Specifies the form of the system of equations. -*> = 'N': A * X = B (No transpose) +*> = 'N': A * X = B (No transpose) *> = 'T': A**T * X = B (Transpose) *> = 'C': A**H * X = B (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/dgtt02.f b/TESTING/LIN/dgtt02.f index 30de915b73..36cfef7fee 100644 --- a/TESTING/LIN/dgtt02.f +++ b/TESTING/LIN/dgtt02.f @@ -40,7 +40,7 @@ *> \verbatim *> TRANS is CHARACTER *> Specifies the form of the residual. -*> = 'N': B - A * X (No transpose) +*> = 'N': B - A * X (No transpose) *> = 'T': B - A**T * X (Transpose) *> = 'C': B - A**H * X (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/dlarhs.f b/TESTING/LIN/dlarhs.f index 8e63387e8c..48a1d54a68 100644 --- a/TESTING/LIN/dlarhs.f +++ b/TESTING/LIN/dlarhs.f @@ -79,8 +79,9 @@ *> \param[in] TRANS *> \verbatim *> TRANS is CHARACTER*1 -*> Specifies the operation applied to the matrix A. -*> = 'N': System is A * x = b +*> Used only if A is nonsymmetric; specifies the operation +*> applied to the matrix A. +*> = 'N': B := A * X (No transpose) *> = 'T': B := A**T * X (Transpose) *> = 'C': B := A**H * X (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/dtbt02.f b/TESTING/LIN/dtbt02.f index 429800c868..d570c5c137 100644 --- a/TESTING/LIN/dtbt02.f +++ b/TESTING/LIN/dtbt02.f @@ -51,7 +51,7 @@ *> \verbatim *> TRANS is CHARACTER*1 *> Specifies the operation applied to A. -*> = 'N': A *x = b (No transpose) +*> = 'N': A * X = B (No transpose) *> = 'T': A**T * X = B (Transpose) *> = 'C': A**H * X = B (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/sgbt05.f b/TESTING/LIN/sgbt05.f index cfd6dd743e..bba80b56da 100644 --- a/TESTING/LIN/sgbt05.f +++ b/TESTING/LIN/sgbt05.f @@ -50,7 +50,7 @@ *> \verbatim *> TRANS is CHARACTER*1 *> Specifies the form of the system of equations. -*> = 'N': A * X = B (No transpose) +*> = 'N': A * X = B (No transpose) *> = 'T': A**T * X = B (Transpose) *> = 'C': A**H * X = B (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/sgtt02.f b/TESTING/LIN/sgtt02.f index 0d19b70f24..4356e82fe6 100644 --- a/TESTING/LIN/sgtt02.f +++ b/TESTING/LIN/sgtt02.f @@ -40,7 +40,7 @@ *> \verbatim *> TRANS is CHARACTER *> Specifies the form of the residual. -*> = 'N': B - A * X (No transpose) +*> = 'N': B - A * X (No transpose) *> = 'T': B - A**T * X (Transpose) *> = 'C': B - A**H * X (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/slarhs.f b/TESTING/LIN/slarhs.f index 2dc3361f5e..769714ce37 100644 --- a/TESTING/LIN/slarhs.f +++ b/TESTING/LIN/slarhs.f @@ -79,8 +79,9 @@ *> \param[in] TRANS *> \verbatim *> TRANS is CHARACTER*1 -*> Specifies the operation applied to the matrix A. -*> = 'N': System is A * x = b +*> Used only if A is nonsymmetric; specifies the operation +*> applied to the matrix A. +*> = 'N': B := A * X (No transpose) *> = 'T': B := A**T * X (Transpose) *> = 'C': B := A**H * X (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/stbt02.f b/TESTING/LIN/stbt02.f index 3b602112d3..da23a9595e 100644 --- a/TESTING/LIN/stbt02.f +++ b/TESTING/LIN/stbt02.f @@ -50,7 +50,7 @@ *> \verbatim *> TRANS is CHARACTER*1 *> Specifies the operation applied to A. -*> = 'N': A *x = b (No transpose) +*> = 'N': A * X = B (No transpose) *> = 'T': A**T * X = B (Transpose) *> = 'C': A**H * X = B (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/zgbt05.f b/TESTING/LIN/zgbt05.f index f74287f421..6c39fa2a7d 100644 --- a/TESTING/LIN/zgbt05.f +++ b/TESTING/LIN/zgbt05.f @@ -50,7 +50,7 @@ *> \verbatim *> TRANS is CHARACTER*1 *> Specifies the form of the system of equations. -*> = 'N': A * X = B (No transpose) +*> = 'N': A * X = B (No transpose) *> = 'T': A**T * X = B (Transpose) *> = 'C': A**H * X = B (Conjugate transpose = Transpose) *> \endverbatim diff --git a/TESTING/LIN/zgtt02.f b/TESTING/LIN/zgtt02.f index ded2256064..ed4d340bed 100644 --- a/TESTING/LIN/zgtt02.f +++ b/TESTING/LIN/zgtt02.f @@ -40,7 +40,7 @@ *> \verbatim *> TRANS is CHARACTER *> Specifies the form of the residual. -*> = 'N': B - A * X (No transpose) +*> = 'N': B - A * X (No transpose) *> = 'T': B - A**T * X (Transpose) *> = 'C': B - A**H * X (Conjugate transpose) *> \endverbatim diff --git a/TESTING/LIN/zlarhs.f b/TESTING/LIN/zlarhs.f index 87e463c4d6..04032cacee 100644 --- a/TESTING/LIN/zlarhs.f +++ b/TESTING/LIN/zlarhs.f @@ -85,9 +85,9 @@ *> TRANS is CHARACTER*1 *> Used only if A is nonsymmetric; specifies the operation *> applied to the matrix A. -*> = 'N': B := A * X -*> = 'T': B := A**T * X -*> = 'C': B := A**H * X +*> = 'N': B := A * X (No transpose) +*> = 'T': B := A**T * X (Transpose) +*> = 'C': B := A**H * X (Conjugate transpose) *> \endverbatim *> *> \param[in] M diff --git a/TESTING/LIN/ztbt02.f b/TESTING/LIN/ztbt02.f index 2553d9eaca..8ef83cfd18 100644 --- a/TESTING/LIN/ztbt02.f +++ b/TESTING/LIN/ztbt02.f @@ -53,9 +53,9 @@ *> \verbatim *> TRANS is CHARACTER*1 *> Specifies the operation applied to A. -*> = 'N': A *x = b (No transpose) -*> = 'T': A**T *x = b (Transpose) -*> = 'C': A**H *x = b (Conjugate transpose) +*> = 'N': A * X = B (No transpose) +*> = 'T': A**T * X = B (Transpose) +*> = 'C': A**H * X = B (Conjugate transpose) *> \endverbatim *> *> \param[in] DIAG From 47c5704326813a0626a2ab6b86d8bb14c96276f6 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Wed, 2 Jun 2021 22:24:17 +1000 Subject: [PATCH 06/16] comments unified: replace a notation ||expr|| by norm(expr) in xQRT17 --- TESTING/LIN/cqrt17.f | 6 +++--- TESTING/LIN/dqrt17.f | 6 +++--- TESTING/LIN/sqrt17.f | 6 +++--- TESTING/LIN/zqrt17.f | 6 +++--- 4 files changed, 12 insertions(+), 12 deletions(-) diff --git a/TESTING/LIN/cqrt17.f b/TESTING/LIN/cqrt17.f index d38486e3fb..0f2f437fc3 100644 --- a/TESTING/LIN/cqrt17.f +++ b/TESTING/LIN/cqrt17.f @@ -28,12 +28,12 @@ *> *> CQRT17 computes the ratio *> -*> || R**H * op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) +*> norm(R**H * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), *> *> where R = op(A)*X - B, op(A) is A or A**H, and *> -*> alpha = ||B|| if IRESID = 1 (zero-residual problem) -*> alpha = ||R|| if IRESID = 2 (otherwise). +*> alpha = norm(B) if IRESID = 1 (zero-residual problem) +*> alpha = norm(R) if IRESID = 2 (otherwise). *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/dqrt17.f b/TESTING/LIN/dqrt17.f index e34ce3381b..e18c6b7788 100644 --- a/TESTING/LIN/dqrt17.f +++ b/TESTING/LIN/dqrt17.f @@ -28,12 +28,12 @@ *> *> DQRT17 computes the ratio *> -*> || R**T * op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) +*> norm(R**T * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), *> *> where R = op(A)*X - B, op(A) is A or A**T, and *> -*> alpha = ||B|| if IRESID = 1 (zero-residual problem) -*> alpha = ||R|| if IRESID = 2 (otherwise). +*> alpha = norm(B) if IRESID = 1 (zero-residual problem) +*> alpha = norm(R) if IRESID = 2 (otherwise). *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/sqrt17.f b/TESTING/LIN/sqrt17.f index fe2914aba3..9615b91e86 100644 --- a/TESTING/LIN/sqrt17.f +++ b/TESTING/LIN/sqrt17.f @@ -28,12 +28,12 @@ *> *> SQRT17 computes the ratio *> -*> || R**T * op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) +*> norm(R**T * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), *> *> where R = op(A)*X - B, op(A) is A or A**T, and *> -*> alpha = ||B|| if IRESID = 1 (zero-residual problem) -*> alpha = ||R|| if IRESID = 2 (otherwise). +*> alpha = norm(B) if IRESID = 1 (zero-residual problem) +*> alpha = norm(R) if IRESID = 2 (otherwise). *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/zqrt17.f b/TESTING/LIN/zqrt17.f index 42c658b293..136e79a40c 100644 --- a/TESTING/LIN/zqrt17.f +++ b/TESTING/LIN/zqrt17.f @@ -28,12 +28,12 @@ *> *> ZQRT17 computes the ratio *> -*> || R**H * op(A) ||/(||A||*alpha*max(M,N,NRHS)*eps) +*> norm(R**H * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), *> *> where R = op(A)*X - B, op(A) is A or A**H, and *> -*> alpha = ||B|| if IRESID = 1 (zero-residual problem) -*> alpha = ||R|| if IRESID = 2 (otherwise). +*> alpha = norm(B) if IRESID = 1 (zero-residual problem) +*> alpha = norm(R) if IRESID = 2 (otherwise). *> \endverbatim * * Arguments: From 8671e8d098f543e329f7c610b5a9c682027769bf Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Thu, 3 Jun 2021 16:38:24 +1000 Subject: [PATCH 07/16] comments unified: about op(A) in xLARHS and xQRT17 --- TESTING/EIG/clarhs.f | 3 +-- TESTING/EIG/zlarhs.f | 3 +-- TESTING/LIN/clarhs.f | 3 +-- TESTING/LIN/cqrt17.f | 2 +- TESTING/LIN/dqrt17.f | 2 +- TESTING/LIN/slarhs.f | 2 +- TESTING/LIN/sqrt17.f | 2 +- TESTING/LIN/zlarhs.f | 3 +-- TESTING/LIN/zqrt17.f | 2 +- 9 files changed, 9 insertions(+), 13 deletions(-) diff --git a/TESTING/EIG/clarhs.f b/TESTING/EIG/clarhs.f index 4505ab34bd..9842df9599 100644 --- a/TESTING/EIG/clarhs.f +++ b/TESTING/EIG/clarhs.f @@ -30,8 +30,7 @@ *> CLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system *> op(A) * X = B, -*> where op(A) = A, A**T (transpose of A) or A**H (conjugate -*> transpose of A). +*> where op(A) = A, A**T or A**H, depending on TRANS. *> \endverbatim * * Arguments: diff --git a/TESTING/EIG/zlarhs.f b/TESTING/EIG/zlarhs.f index 788e033289..1dc003f2e4 100644 --- a/TESTING/EIG/zlarhs.f +++ b/TESTING/EIG/zlarhs.f @@ -30,8 +30,7 @@ *> ZLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system *> op(A) * X = B, -*> where op(A) = A, A**T (transpose of A), or A**H (conjugate -*> transpose of A). +*> where op(A) = A, A**T, or A**H, depending on TRANS. *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/clarhs.f b/TESTING/LIN/clarhs.f index 684d4489cc..279c98f953 100644 --- a/TESTING/LIN/clarhs.f +++ b/TESTING/LIN/clarhs.f @@ -30,8 +30,7 @@ *> CLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system *> op(A) * X = B, -*> where op(A) = A, A**T (transpose of A), or A**H (conjugate -*> transpose of A). +*> where op(A) = A, A**T, or A**H, depending on TRANS. *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/cqrt17.f b/TESTING/LIN/cqrt17.f index 0f2f437fc3..bafc3123a9 100644 --- a/TESTING/LIN/cqrt17.f +++ b/TESTING/LIN/cqrt17.f @@ -30,7 +30,7 @@ *> *> norm(R**H * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), *> -*> where R = op(A)*X - B, op(A) is A or A**H, and +*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, and *> *> alpha = norm(B) if IRESID = 1 (zero-residual problem) *> alpha = norm(R) if IRESID = 2 (otherwise). diff --git a/TESTING/LIN/dqrt17.f b/TESTING/LIN/dqrt17.f index e18c6b7788..a727a10909 100644 --- a/TESTING/LIN/dqrt17.f +++ b/TESTING/LIN/dqrt17.f @@ -30,7 +30,7 @@ *> *> norm(R**T * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), *> -*> where R = op(A)*X - B, op(A) is A or A**T, and +*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, and *> *> alpha = norm(B) if IRESID = 1 (zero-residual problem) *> alpha = norm(R) if IRESID = 2 (otherwise). diff --git a/TESTING/LIN/slarhs.f b/TESTING/LIN/slarhs.f index 769714ce37..6a8a592c8c 100644 --- a/TESTING/LIN/slarhs.f +++ b/TESTING/LIN/slarhs.f @@ -30,7 +30,7 @@ *> SLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system *> op(A) * X = B, -*> where op(A) = A or A**T. +*> where op(A) = A or A**T, depending on TRANS. *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/sqrt17.f b/TESTING/LIN/sqrt17.f index 9615b91e86..4d60cafa4d 100644 --- a/TESTING/LIN/sqrt17.f +++ b/TESTING/LIN/sqrt17.f @@ -30,7 +30,7 @@ *> *> norm(R**T * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), *> -*> where R = op(A)*X - B, op(A) is A or A**T, and +*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, and *> *> alpha = norm(B) if IRESID = 1 (zero-residual problem) *> alpha = norm(R) if IRESID = 2 (otherwise). diff --git a/TESTING/LIN/zlarhs.f b/TESTING/LIN/zlarhs.f index 04032cacee..210a4a3c14 100644 --- a/TESTING/LIN/zlarhs.f +++ b/TESTING/LIN/zlarhs.f @@ -30,8 +30,7 @@ *> ZLARHS chooses a set of NRHS random solution vectors and sets *> up the right hand sides for the linear system *> op(A) * X = B, -*> where op(A) = A, A**T (transpose of A), or A**H (conjugate -*> transpose of A). +*> where op(A) = A, A**T, or A**H, depending on TRANS. *> \endverbatim * * Arguments: diff --git a/TESTING/LIN/zqrt17.f b/TESTING/LIN/zqrt17.f index 136e79a40c..5c2a90f6ee 100644 --- a/TESTING/LIN/zqrt17.f +++ b/TESTING/LIN/zqrt17.f @@ -30,7 +30,7 @@ *> *> norm(R**H * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), *> -*> where R = op(A)*X - B, op(A) is A or A**H, and +*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, and *> *> alpha = norm(B) if IRESID = 1 (zero-residual problem) *> alpha = norm(R) if IRESID = 2 (otherwise). From b1143c56a77cd83f0d53031f6ef759ac18f15943 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Thu, 3 Jun 2021 16:47:05 +1000 Subject: [PATCH 08/16] comments fixed: line length in CBDT01 and xPOT01 --- TESTING/EIG/cbdt01.f | 3 ++- TESTING/LIN/cpot01.f | 4 ++-- TESTING/LIN/dpot01.f | 4 ++-- TESTING/LIN/spot01.f | 4 ++-- TESTING/LIN/zpot01.f | 2 +- 5 files changed, 9 insertions(+), 8 deletions(-) diff --git a/TESTING/EIG/cbdt01.f b/TESTING/EIG/cbdt01.f index 2f678ef233..f5633bdc31 100644 --- a/TESTING/EIG/cbdt01.f +++ b/TESTING/EIG/cbdt01.f @@ -127,7 +127,8 @@ *> \param[out] RESID *> \verbatim *> RESID is REAL -*> The test ratio: norm(A - Q * B * P**H) / ( n * norm(A) * EPS ) +*> The test ratio: +*> norm(A - Q * B * P**H) / ( n * norm(A) * EPS ) *> \endverbatim * * Authors: diff --git a/TESTING/LIN/cpot01.f b/TESTING/LIN/cpot01.f index 30dc55d884..fc8de3f045 100644 --- a/TESTING/LIN/cpot01.f +++ b/TESTING/LIN/cpot01.f @@ -69,8 +69,8 @@ *> AFAC is COMPLEX array, dimension (LDAFAC,N) *> On entry, the factor L or U from the L * L**H or U**H * U *> factorization of A. -*> Overwritten with the reconstructed matrix, and then with the -*> difference L * L**H - A (or U**H * U - A). +*> Overwritten with the reconstructed matrix, and then with +*> the difference L * L**H - A (or U**H * U - A). *> \endverbatim *> *> \param[in] LDAFAC diff --git a/TESTING/LIN/dpot01.f b/TESTING/LIN/dpot01.f index 5c463c12c2..b48d066101 100644 --- a/TESTING/LIN/dpot01.f +++ b/TESTING/LIN/dpot01.f @@ -67,8 +67,8 @@ *> AFAC is DOUBLE PRECISION array, dimension (LDAFAC,N) *> On entry, the factor L or U from the L * L**T or U**T * U *> factorization of A. -*> Overwritten with the reconstructed matrix, and then with the -*> difference L * L**T - A (or U**T * U - A). +*> Overwritten with the reconstructed matrix, and then with +*> the difference L * L**T - A (or U**T * U - A). *> \endverbatim *> *> \param[in] LDAFAC diff --git a/TESTING/LIN/spot01.f b/TESTING/LIN/spot01.f index 05c85815b7..584bdce6d7 100644 --- a/TESTING/LIN/spot01.f +++ b/TESTING/LIN/spot01.f @@ -67,8 +67,8 @@ *> AFAC is REAL array, dimension (LDAFAC,N) *> On entry, the factor L or U from the L * L**T or U**T * U *> factorization of A. -*> Overwritten with the reconstructed matrix, and then with the -*> difference L * L**T - A (or U**T * U - A). +*> Overwritten with the reconstructed matrix, and then with +*> the difference L * L**T - A (or U**T * U - A). *> \endverbatim *> *> \param[in] LDAFAC diff --git a/TESTING/LIN/zpot01.f b/TESTING/LIN/zpot01.f index 2c48ff001f..995851a2c8 100644 --- a/TESTING/LIN/zpot01.f +++ b/TESTING/LIN/zpot01.f @@ -69,7 +69,7 @@ *> AFAC is COMPLEX*16 array, dimension (LDAFAC,N) *> On entry, the factor L or U from the L * L**H or U**H * U *> factorization of A. -*> Overwritten with the reconstructed matrix, and then with the +*> Overwritten with the reconstructed matrix, and then with *> the difference L * L**H - A (or U**H * U - A). *> \endverbatim *> From dc6719220f936957fa89af5cb6f6652fa5afc19e Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Thu, 3 Jun 2021 16:52:01 +1000 Subject: [PATCH 09/16] comments unified: compress space chains in xGBT01 and xPOT01 --- TESTING/LIN/cgbt01.f | 2 +- TESTING/LIN/cpot01.f | 2 +- TESTING/LIN/dgbt01.f | 2 +- TESTING/LIN/sgbt01.f | 2 +- TESTING/LIN/spot01.f | 2 +- TESTING/LIN/zgbt01.f | 2 +- TESTING/LIN/zpot01.f | 2 +- 7 files changed, 7 insertions(+), 7 deletions(-) diff --git a/TESTING/LIN/cgbt01.f b/TESTING/LIN/cgbt01.f index 375fbc1e3c..148076e9ac 100644 --- a/TESTING/LIN/cgbt01.f +++ b/TESTING/LIN/cgbt01.f @@ -26,7 +26,7 @@ *> *> \verbatim *> -*> CGBT01 reconstructs a band matrix A from its L*U factorization and +*> CGBT01 reconstructs a band matrix A from its L*U factorization and *> computes the residual: *> norm(L*U - A) / ( N * norm(A) * EPS ), *> where EPS is the machine epsilon. diff --git a/TESTING/LIN/cpot01.f b/TESTING/LIN/cpot01.f index fc8de3f045..3c8e9a8859 100644 --- a/TESTING/LIN/cpot01.f +++ b/TESTING/LIN/cpot01.f @@ -206,7 +206,7 @@ SUBROUTINE CPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 30 CONTINUE END IF * -* Compute the difference L * L**H - A (or U**H * U - A). +* Compute the difference L * L**H - A (or U**H * U - A). * IF( LSAME( UPLO, 'U' ) ) THEN DO 50 J = 1, N diff --git a/TESTING/LIN/dgbt01.f b/TESTING/LIN/dgbt01.f index bf0e77a4fa..6a095876ed 100644 --- a/TESTING/LIN/dgbt01.f +++ b/TESTING/LIN/dgbt01.f @@ -26,7 +26,7 @@ *> *> \verbatim *> -*> DGBT01 reconstructs a band matrix A from its L*U factorization and +*> DGBT01 reconstructs a band matrix A from its L*U factorization and *> computes the residual: *> norm(L*U - A) / ( N * norm(A) * EPS ), *> where EPS is the machine epsilon. diff --git a/TESTING/LIN/sgbt01.f b/TESTING/LIN/sgbt01.f index 9abcf56e28..c5ae6b11a5 100644 --- a/TESTING/LIN/sgbt01.f +++ b/TESTING/LIN/sgbt01.f @@ -26,7 +26,7 @@ *> *> \verbatim *> -*> SGBT01 reconstructs a band matrix A from its L*U factorization and +*> SGBT01 reconstructs a band matrix A from its L*U factorization and *> computes the residual: *> norm(L*U - A) / ( N * norm(A) * EPS ), *> where EPS is the machine epsilon. diff --git a/TESTING/LIN/spot01.f b/TESTING/LIN/spot01.f index 584bdce6d7..5f14b3f0ff 100644 --- a/TESTING/LIN/spot01.f +++ b/TESTING/LIN/spot01.f @@ -191,7 +191,7 @@ SUBROUTINE SPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 20 CONTINUE END IF * -* Compute the difference L * L**T - A (or U**T * U - A). +* Compute the difference L * L**T - A (or U**T * U - A). * IF( LSAME( UPLO, 'U' ) ) THEN DO 40 J = 1, N diff --git a/TESTING/LIN/zgbt01.f b/TESTING/LIN/zgbt01.f index 67efa58f8e..e041ab9fa0 100644 --- a/TESTING/LIN/zgbt01.f +++ b/TESTING/LIN/zgbt01.f @@ -26,7 +26,7 @@ *> *> \verbatim *> -*> ZGBT01 reconstructs a band matrix A from its L*U factorization and +*> ZGBT01 reconstructs a band matrix A from its L*U factorization and *> computes the residual: *> norm(L*U - A) / ( N * norm(A) * EPS ), *> where EPS is the machine epsilon. diff --git a/TESTING/LIN/zpot01.f b/TESTING/LIN/zpot01.f index 995851a2c8..b083856211 100644 --- a/TESTING/LIN/zpot01.f +++ b/TESTING/LIN/zpot01.f @@ -206,7 +206,7 @@ SUBROUTINE ZPOT01( UPLO, N, A, LDA, AFAC, LDAFAC, RWORK, RESID ) 30 CONTINUE END IF * -* Compute the difference L * L**H - A (or U**H * U - A). +* Compute the difference L * L**H - A (or U**H * U - A). * IF( LSAME( UPLO, 'U' ) ) THEN DO 50 J = 1, N From a98ffc640d82443837eeef547a77d72ae1fcd1b0 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Thu, 3 Jun 2021 16:54:24 +1000 Subject: [PATCH 10/16] comments unified: specify where result will be stored in xGBT01 and xPOT01 Were unified with xGET02, xGET08, xPOT06, xQRT16, xSYT02. --- TESTING/LIN/cgtt02.f | 2 +- TESTING/LIN/dgtt02.f | 2 +- TESTING/LIN/sgtt02.f | 2 +- TESTING/LIN/zgtt02.f | 2 +- 4 files changed, 4 insertions(+), 4 deletions(-) diff --git a/TESTING/LIN/cgtt02.f b/TESTING/LIN/cgtt02.f index 587c43a6e5..0d067c3378 100644 --- a/TESTING/LIN/cgtt02.f +++ b/TESTING/LIN/cgtt02.f @@ -182,7 +182,7 @@ SUBROUTINE CGTT02( TRANS, N, NRHS, DL, D, DU, X, LDX, B, LDB, RETURN END IF * -* Compute B - op(A)*X. +* Compute B - op(A)*X and store in B. * CALL CLAGTM( TRANS, N, NRHS, -ONE, DL, D, DU, X, LDX, ONE, B, $ LDB ) diff --git a/TESTING/LIN/dgtt02.f b/TESTING/LIN/dgtt02.f index 36cfef7fee..e91f002249 100644 --- a/TESTING/LIN/dgtt02.f +++ b/TESTING/LIN/dgtt02.f @@ -182,7 +182,7 @@ SUBROUTINE DGTT02( TRANS, N, NRHS, DL, D, DU, X, LDX, B, LDB, RETURN END IF * -* Compute B - op(A)*X. +* Compute B - op(A)*X and store in B. * CALL DLAGTM( TRANS, N, NRHS, -ONE, DL, D, DU, X, LDX, ONE, B, $ LDB ) diff --git a/TESTING/LIN/sgtt02.f b/TESTING/LIN/sgtt02.f index 4356e82fe6..0b03e7f592 100644 --- a/TESTING/LIN/sgtt02.f +++ b/TESTING/LIN/sgtt02.f @@ -182,7 +182,7 @@ SUBROUTINE SGTT02( TRANS, N, NRHS, DL, D, DU, X, LDX, B, LDB, RETURN END IF * -* Compute B - op(A)*X. +* Compute B - op(A)*X and store in B. * CALL SLAGTM( TRANS, N, NRHS, -ONE, DL, D, DU, X, LDX, ONE, B, $ LDB ) diff --git a/TESTING/LIN/zgtt02.f b/TESTING/LIN/zgtt02.f index ed4d340bed..72cc9ea712 100644 --- a/TESTING/LIN/zgtt02.f +++ b/TESTING/LIN/zgtt02.f @@ -182,7 +182,7 @@ SUBROUTINE ZGTT02( TRANS, N, NRHS, DL, D, DU, X, LDX, B, LDB, RETURN END IF * -* Compute B - op(A)*X. +* Compute B - op(A)*X and store in B. * CALL ZLAGTM( TRANS, N, NRHS, -ONE, DL, D, DU, X, LDX, ONE, B, $ LDB ) From 80409a52f4d5ed986f67656c026ba3da596317ba Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Thu, 3 Jun 2021 17:11:21 +1000 Subject: [PATCH 11/16] comments unified: add an EPS description in xQRT17 --- TESTING/LIN/cqrt17.f | 5 +++-- TESTING/LIN/dqrt17.f | 5 +++-- TESTING/LIN/sqrt17.f | 5 +++-- TESTING/LIN/zqrt17.f | 5 +++-- 4 files changed, 12 insertions(+), 8 deletions(-) diff --git a/TESTING/LIN/cqrt17.f b/TESTING/LIN/cqrt17.f index bafc3123a9..85443c8387 100644 --- a/TESTING/LIN/cqrt17.f +++ b/TESTING/LIN/cqrt17.f @@ -28,9 +28,10 @@ *> *> CQRT17 computes the ratio *> -*> norm(R**H * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), +*> norm(R**H * op(A)) / ( norm(A) * alpha * max(M,N,NRHS) * EPS ), *> -*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, and +*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, EPS +*> is the machine epsilon, and *> *> alpha = norm(B) if IRESID = 1 (zero-residual problem) *> alpha = norm(R) if IRESID = 2 (otherwise). diff --git a/TESTING/LIN/dqrt17.f b/TESTING/LIN/dqrt17.f index a727a10909..15479f8e52 100644 --- a/TESTING/LIN/dqrt17.f +++ b/TESTING/LIN/dqrt17.f @@ -28,9 +28,10 @@ *> *> DQRT17 computes the ratio *> -*> norm(R**T * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), +*> norm(R**T * op(A)) / ( norm(A) * alpha * max(M,N,NRHS) * EPS ), *> -*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, and +*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, EPS +*> is the machine epsilon, and *> *> alpha = norm(B) if IRESID = 1 (zero-residual problem) *> alpha = norm(R) if IRESID = 2 (otherwise). diff --git a/TESTING/LIN/sqrt17.f b/TESTING/LIN/sqrt17.f index 4d60cafa4d..83c3ca98fd 100644 --- a/TESTING/LIN/sqrt17.f +++ b/TESTING/LIN/sqrt17.f @@ -28,9 +28,10 @@ *> *> SQRT17 computes the ratio *> -*> norm(R**T * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), +*> norm(R**T * op(A)) / ( norm(A) * alpha * max(M,N,NRHS) * EPS ), *> -*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, and +*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, EPS +*> is the machine epsilon, and *> *> alpha = norm(B) if IRESID = 1 (zero-residual problem) *> alpha = norm(R) if IRESID = 2 (otherwise). diff --git a/TESTING/LIN/zqrt17.f b/TESTING/LIN/zqrt17.f index 5c2a90f6ee..060d46cb9a 100644 --- a/TESTING/LIN/zqrt17.f +++ b/TESTING/LIN/zqrt17.f @@ -28,9 +28,10 @@ *> *> ZQRT17 computes the ratio *> -*> norm(R**H * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps), +*> norm(R**H * op(A)) / ( norm(A) * alpha * max(M,N,NRHS) * EPS ), *> -*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, and +*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, EPS +*> is the machine epsilon, and *> *> alpha = norm(B) if IRESID = 1 (zero-residual problem) *> alpha = norm(R) if IRESID = 2 (otherwise). From b3166e1b7594e9ab373bd1ffbaf2822f03035f66 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Tue, 1 Jun 2021 13:37:55 +1000 Subject: [PATCH 12/16] comments fixed: NBVAL dimension --- TESTING/LIN/cchkge.f | 2 +- TESTING/LIN/cchkhe.f | 2 +- TESTING/LIN/cchkhe_aa.f | 2 +- TESTING/LIN/cchkhe_aa_2stage.f | 2 +- TESTING/LIN/cchkhe_rk.f | 2 +- TESTING/LIN/cchkhe_rook.f | 2 +- TESTING/LIN/cchklqt.f | 2 +- TESTING/LIN/cchklqtp.f | 2 +- TESTING/LIN/cchkpb.f | 2 +- TESTING/LIN/cchkpo.f | 2 +- TESTING/LIN/cchkps.f | 2 +- TESTING/LIN/cchkqrt.f | 2 +- TESTING/LIN/cchkqrtp.f | 2 +- TESTING/LIN/cchksy.f | 2 +- TESTING/LIN/cchksy_aa.f | 2 +- TESTING/LIN/cchksy_aa_2stage.f | 2 +- TESTING/LIN/cchksy_rk.f | 2 +- TESTING/LIN/cchksy_rook.f | 2 +- TESTING/LIN/cchktsqr.f | 2 +- TESTING/LIN/cchkunhr_col.f | 2 +- TESTING/LIN/dchkge.f | 2 +- TESTING/LIN/dchklqt.f | 2 +- TESTING/LIN/dchklqtp.f | 2 +- TESTING/LIN/dchkorhr_col.f | 2 +- TESTING/LIN/dchkpb.f | 2 +- TESTING/LIN/dchkpo.f | 2 +- TESTING/LIN/dchkps.f | 2 +- TESTING/LIN/dchkqrt.f | 2 +- TESTING/LIN/dchkqrtp.f | 2 +- TESTING/LIN/dchksy.f | 2 +- TESTING/LIN/dchksy_aa.f | 2 +- TESTING/LIN/dchksy_aa_2stage.f | 2 +- TESTING/LIN/dchksy_rk.f | 2 +- TESTING/LIN/dchksy_rook.f | 2 +- TESTING/LIN/dchktsqr.f | 2 +- TESTING/LIN/schkge.f | 2 +- TESTING/LIN/schklqt.f | 2 +- TESTING/LIN/schklqtp.f | 2 +- TESTING/LIN/schkorhr_col.f | 2 +- TESTING/LIN/schkpb.f | 2 +- TESTING/LIN/schkpo.f | 2 +- TESTING/LIN/schkps.f | 2 +- TESTING/LIN/schkqrt.f | 2 +- TESTING/LIN/schkqrtp.f | 2 +- TESTING/LIN/schksy.f | 2 +- TESTING/LIN/schksy_aa.f | 2 +- TESTING/LIN/schksy_aa_2stage.f | 2 +- TESTING/LIN/schksy_rk.f | 2 +- TESTING/LIN/schksy_rook.f | 2 +- TESTING/LIN/schktsqr.f | 2 +- TESTING/LIN/zchkgb.f | 2 +- TESTING/LIN/zchkge.f | 2 +- TESTING/LIN/zchkhe.f | 2 +- TESTING/LIN/zchkhe_aa.f | 2 +- TESTING/LIN/zchkhe_aa_2stage.f | 2 +- TESTING/LIN/zchkhe_rk.f | 2 +- TESTING/LIN/zchkhe_rook.f | 2 +- TESTING/LIN/zchklqt.f | 2 +- TESTING/LIN/zchklqtp.f | 2 +- TESTING/LIN/zchkpb.f | 2 +- TESTING/LIN/zchkpo.f | 2 +- TESTING/LIN/zchkps.f | 2 +- TESTING/LIN/zchkqrt.f | 2 +- TESTING/LIN/zchkqrtp.f | 2 +- TESTING/LIN/zchksy.f | 2 +- TESTING/LIN/zchksy_aa.f | 2 +- TESTING/LIN/zchksy_aa_2stage.f | 2 +- TESTING/LIN/zchksy_rk.f | 2 +- TESTING/LIN/zchksy_rook.f | 2 +- TESTING/LIN/zchktsqr.f | 2 +- TESTING/LIN/zchkunhr_col.f | 2 +- 71 files changed, 71 insertions(+), 71 deletions(-) diff --git a/TESTING/LIN/cchkge.f b/TESTING/LIN/cchkge.f index d970571ac3..0ee1c85c57 100644 --- a/TESTING/LIN/cchkge.f +++ b/TESTING/LIN/cchkge.f @@ -78,7 +78,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkhe.f b/TESTING/LIN/cchkhe.f index 875d8e3997..1ff0257579 100644 --- a/TESTING/LIN/cchkhe.f +++ b/TESTING/LIN/cchkhe.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkhe_aa.f b/TESTING/LIN/cchkhe_aa.f index 826faf7d42..487f9941e8 100644 --- a/TESTING/LIN/cchkhe_aa.f +++ b/TESTING/LIN/cchkhe_aa.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkhe_aa_2stage.f b/TESTING/LIN/cchkhe_aa_2stage.f index 6ba9c309ca..30a61261f5 100644 --- a/TESTING/LIN/cchkhe_aa_2stage.f +++ b/TESTING/LIN/cchkhe_aa_2stage.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkhe_rk.f b/TESTING/LIN/cchkhe_rk.f index 96a2fe1d69..4f9801f88d 100644 --- a/TESTING/LIN/cchkhe_rk.f +++ b/TESTING/LIN/cchkhe_rk.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkhe_rook.f b/TESTING/LIN/cchkhe_rook.f index 52ab5db702..72a4a1d4d0 100644 --- a/TESTING/LIN/cchkhe_rook.f +++ b/TESTING/LIN/cchkhe_rook.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchklqt.f b/TESTING/LIN/cchklqt.f index 57b843fbc0..7ec8b2c622 100644 --- a/TESTING/LIN/cchklqt.f +++ b/TESTING/LIN/cchklqt.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchklqtp.f b/TESTING/LIN/cchklqtp.f index 9fd21a964b..987c0c3654 100644 --- a/TESTING/LIN/cchklqtp.f +++ b/TESTING/LIN/cchklqtp.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkpb.f b/TESTING/LIN/cchkpb.f index bcfc771f77..4d1ab6f4f9 100644 --- a/TESTING/LIN/cchkpb.f +++ b/TESTING/LIN/cchkpb.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkpo.f b/TESTING/LIN/cchkpo.f index 87d885bf26..8cf85bf6e4 100644 --- a/TESTING/LIN/cchkpo.f +++ b/TESTING/LIN/cchkpo.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkps.f b/TESTING/LIN/cchkps.f index 45702fa692..d77ffc0262 100644 --- a/TESTING/LIN/cchkps.f +++ b/TESTING/LIN/cchkps.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the block size NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkqrt.f b/TESTING/LIN/cchkqrt.f index 5dc69dd034..7aeb6df4be 100644 --- a/TESTING/LIN/cchkqrt.f +++ b/TESTING/LIN/cchkqrt.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkqrtp.f b/TESTING/LIN/cchkqrtp.f index 246e761215..6820c90042 100644 --- a/TESTING/LIN/cchkqrtp.f +++ b/TESTING/LIN/cchkqrtp.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchksy.f b/TESTING/LIN/cchksy.f index 3392188237..c933f0c967 100644 --- a/TESTING/LIN/cchksy.f +++ b/TESTING/LIN/cchksy.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchksy_aa.f b/TESTING/LIN/cchksy_aa.f index 735e79bd04..015d5a68e6 100644 --- a/TESTING/LIN/cchksy_aa.f +++ b/TESTING/LIN/cchksy_aa.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchksy_aa_2stage.f b/TESTING/LIN/cchksy_aa_2stage.f index 7a0dcc5673..bede7cca5e 100644 --- a/TESTING/LIN/cchksy_aa_2stage.f +++ b/TESTING/LIN/cchksy_aa_2stage.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchksy_rk.f b/TESTING/LIN/cchksy_rk.f index 9ba653c093..f4954054f4 100644 --- a/TESTING/LIN/cchksy_rk.f +++ b/TESTING/LIN/cchksy_rk.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchksy_rook.f b/TESTING/LIN/cchksy_rook.f index bbe18ebf9c..2031089801 100644 --- a/TESTING/LIN/cchksy_rook.f +++ b/TESTING/LIN/cchksy_rook.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchktsqr.f b/TESTING/LIN/cchktsqr.f index e71656671a..0a85d7cf1e 100644 --- a/TESTING/LIN/cchktsqr.f +++ b/TESTING/LIN/cchktsqr.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/cchkunhr_col.f b/TESTING/LIN/cchkunhr_col.f index 0d6a9063db..b2f233e258 100644 --- a/TESTING/LIN/cchkunhr_col.f +++ b/TESTING/LIN/cchkunhr_col.f @@ -82,7 +82,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchkge.f b/TESTING/LIN/dchkge.f index 705a9dbbc2..467667a42b 100644 --- a/TESTING/LIN/dchkge.f +++ b/TESTING/LIN/dchkge.f @@ -77,7 +77,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchklqt.f b/TESTING/LIN/dchklqt.f index ba67011852..aa77d6a62e 100644 --- a/TESTING/LIN/dchklqt.f +++ b/TESTING/LIN/dchklqt.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchklqtp.f b/TESTING/LIN/dchklqtp.f index 8090213305..d9251b1f09 100644 --- a/TESTING/LIN/dchklqtp.f +++ b/TESTING/LIN/dchklqtp.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchkorhr_col.f b/TESTING/LIN/dchkorhr_col.f index 0e2d44d8d2..a849f39106 100644 --- a/TESTING/LIN/dchkorhr_col.f +++ b/TESTING/LIN/dchkorhr_col.f @@ -82,7 +82,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchkpb.f b/TESTING/LIN/dchkpb.f index 74b8979d65..c5d51230ed 100644 --- a/TESTING/LIN/dchkpb.f +++ b/TESTING/LIN/dchkpb.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchkpo.f b/TESTING/LIN/dchkpo.f index e5abd9c7a1..e9d5af1863 100644 --- a/TESTING/LIN/dchkpo.f +++ b/TESTING/LIN/dchkpo.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchkps.f b/TESTING/LIN/dchkps.f index 47d43d0a62..1e4a9e58de 100644 --- a/TESTING/LIN/dchkps.f +++ b/TESTING/LIN/dchkps.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the block size NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchkqrt.f b/TESTING/LIN/dchkqrt.f index 690b041df7..566f273d4e 100644 --- a/TESTING/LIN/dchkqrt.f +++ b/TESTING/LIN/dchkqrt.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchkqrtp.f b/TESTING/LIN/dchkqrtp.f index 1e53bcded7..924ba7ea37 100644 --- a/TESTING/LIN/dchkqrtp.f +++ b/TESTING/LIN/dchkqrtp.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchksy.f b/TESTING/LIN/dchksy.f index 84589d8e02..8d127cd3b6 100644 --- a/TESTING/LIN/dchksy.f +++ b/TESTING/LIN/dchksy.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchksy_aa.f b/TESTING/LIN/dchksy_aa.f index 81139440fa..4c79ebbbf1 100644 --- a/TESTING/LIN/dchksy_aa.f +++ b/TESTING/LIN/dchksy_aa.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchksy_aa_2stage.f b/TESTING/LIN/dchksy_aa_2stage.f index 405f8f12d7..bc4e77a5aa 100644 --- a/TESTING/LIN/dchksy_aa_2stage.f +++ b/TESTING/LIN/dchksy_aa_2stage.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchksy_rk.f b/TESTING/LIN/dchksy_rk.f index 3731a3c666..701cb9ce6a 100644 --- a/TESTING/LIN/dchksy_rk.f +++ b/TESTING/LIN/dchksy_rk.f @@ -63,7 +63,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchksy_rook.f b/TESTING/LIN/dchksy_rook.f index 1af46409ba..03b16b78ac 100644 --- a/TESTING/LIN/dchksy_rook.f +++ b/TESTING/LIN/dchksy_rook.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/dchktsqr.f b/TESTING/LIN/dchktsqr.f index c7f06c11c7..27076b5c26 100644 --- a/TESTING/LIN/dchktsqr.f +++ b/TESTING/LIN/dchktsqr.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schkge.f b/TESTING/LIN/schkge.f index 6a2c0b0cfd..0921f79257 100644 --- a/TESTING/LIN/schkge.f +++ b/TESTING/LIN/schkge.f @@ -77,7 +77,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schklqt.f b/TESTING/LIN/schklqt.f index 966142f8f2..2230583bae 100644 --- a/TESTING/LIN/schklqt.f +++ b/TESTING/LIN/schklqt.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schklqtp.f b/TESTING/LIN/schklqtp.f index aeb3b43710..20c0dd0199 100644 --- a/TESTING/LIN/schklqtp.f +++ b/TESTING/LIN/schklqtp.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schkorhr_col.f b/TESTING/LIN/schkorhr_col.f index f61b749028..9c91b85353 100644 --- a/TESTING/LIN/schkorhr_col.f +++ b/TESTING/LIN/schkorhr_col.f @@ -82,7 +82,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schkpb.f b/TESTING/LIN/schkpb.f index 6a9662bdeb..fdad0d6dfd 100644 --- a/TESTING/LIN/schkpb.f +++ b/TESTING/LIN/schkpb.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schkpo.f b/TESTING/LIN/schkpo.f index 266abcb43f..ccfb9b7b27 100644 --- a/TESTING/LIN/schkpo.f +++ b/TESTING/LIN/schkpo.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schkps.f b/TESTING/LIN/schkps.f index 9385a5809c..7bd70c4525 100644 --- a/TESTING/LIN/schkps.f +++ b/TESTING/LIN/schkps.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the block size NB. *> \endverbatim *> diff --git a/TESTING/LIN/schkqrt.f b/TESTING/LIN/schkqrt.f index 23515d95d4..ed25a230d4 100644 --- a/TESTING/LIN/schkqrt.f +++ b/TESTING/LIN/schkqrt.f @@ -74,7 +74,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schkqrtp.f b/TESTING/LIN/schkqrtp.f index 5be01b6cd7..eee963c472 100644 --- a/TESTING/LIN/schkqrtp.f +++ b/TESTING/LIN/schkqrtp.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schksy.f b/TESTING/LIN/schksy.f index 39353671d3..dc94b34924 100644 --- a/TESTING/LIN/schksy.f +++ b/TESTING/LIN/schksy.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schksy_aa.f b/TESTING/LIN/schksy_aa.f index 573b905aa8..75b895bd43 100644 --- a/TESTING/LIN/schksy_aa.f +++ b/TESTING/LIN/schksy_aa.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schksy_aa_2stage.f b/TESTING/LIN/schksy_aa_2stage.f index 8a5ed94dec..d3c27ae561 100644 --- a/TESTING/LIN/schksy_aa_2stage.f +++ b/TESTING/LIN/schksy_aa_2stage.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schksy_rk.f b/TESTING/LIN/schksy_rk.f index 0158ead9f5..932e6e1860 100644 --- a/TESTING/LIN/schksy_rk.f +++ b/TESTING/LIN/schksy_rk.f @@ -63,7 +63,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schksy_rook.f b/TESTING/LIN/schksy_rook.f index c0ad39ce96..eb91eba4f0 100644 --- a/TESTING/LIN/schksy_rook.f +++ b/TESTING/LIN/schksy_rook.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/schktsqr.f b/TESTING/LIN/schktsqr.f index a3994cbed3..5d7aa2d374 100644 --- a/TESTING/LIN/schktsqr.f +++ b/TESTING/LIN/schktsqr.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkgb.f b/TESTING/LIN/zchkgb.f index d4d261d93a..3b0cd707ab 100644 --- a/TESTING/LIN/zchkgb.f +++ b/TESTING/LIN/zchkgb.f @@ -78,7 +78,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkge.f b/TESTING/LIN/zchkge.f index df46a3fc8a..d87759df5c 100644 --- a/TESTING/LIN/zchkge.f +++ b/TESTING/LIN/zchkge.f @@ -78,7 +78,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkhe.f b/TESTING/LIN/zchkhe.f index f0e746890a..e443eab04f 100644 --- a/TESTING/LIN/zchkhe.f +++ b/TESTING/LIN/zchkhe.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkhe_aa.f b/TESTING/LIN/zchkhe_aa.f index acb0dd31cf..1c4ee88e0e 100644 --- a/TESTING/LIN/zchkhe_aa.f +++ b/TESTING/LIN/zchkhe_aa.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkhe_aa_2stage.f b/TESTING/LIN/zchkhe_aa_2stage.f index f4d9585f33..726e1fea6d 100644 --- a/TESTING/LIN/zchkhe_aa_2stage.f +++ b/TESTING/LIN/zchkhe_aa_2stage.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkhe_rk.f b/TESTING/LIN/zchkhe_rk.f index e670bf3279..63697d33a1 100644 --- a/TESTING/LIN/zchkhe_rk.f +++ b/TESTING/LIN/zchkhe_rk.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkhe_rook.f b/TESTING/LIN/zchkhe_rook.f index c18528b3c6..84f5a84220 100644 --- a/TESTING/LIN/zchkhe_rook.f +++ b/TESTING/LIN/zchkhe_rook.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchklqt.f b/TESTING/LIN/zchklqt.f index 1e941f0bcb..680c1135e2 100644 --- a/TESTING/LIN/zchklqt.f +++ b/TESTING/LIN/zchklqt.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchklqtp.f b/TESTING/LIN/zchklqtp.f index a3ad56ed63..75e3ddfb1d 100644 --- a/TESTING/LIN/zchklqtp.f +++ b/TESTING/LIN/zchklqtp.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkpb.f b/TESTING/LIN/zchkpb.f index a935d10c89..1ce8c3b88f 100644 --- a/TESTING/LIN/zchkpb.f +++ b/TESTING/LIN/zchkpb.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkpo.f b/TESTING/LIN/zchkpo.f index 7ceb65f7c5..93ea721633 100644 --- a/TESTING/LIN/zchkpo.f +++ b/TESTING/LIN/zchkpo.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkps.f b/TESTING/LIN/zchkps.f index 8cf05ccd26..79fc951e22 100644 --- a/TESTING/LIN/zchkps.f +++ b/TESTING/LIN/zchkps.f @@ -64,7 +64,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the block size NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkqrt.f b/TESTING/LIN/zchkqrt.f index cf096ef4eb..cb07af72ce 100644 --- a/TESTING/LIN/zchkqrt.f +++ b/TESTING/LIN/zchkqrt.f @@ -75,7 +75,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkqrtp.f b/TESTING/LIN/zchkqrtp.f index 54cb9ab725..9337a68375 100644 --- a/TESTING/LIN/zchkqrtp.f +++ b/TESTING/LIN/zchkqrtp.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchksy.f b/TESTING/LIN/zchksy.f index 2232497b6f..7877eb7e2d 100644 --- a/TESTING/LIN/zchksy.f +++ b/TESTING/LIN/zchksy.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchksy_aa.f b/TESTING/LIN/zchksy_aa.f index 3f0c52173a..9f53593d57 100644 --- a/TESTING/LIN/zchksy_aa.f +++ b/TESTING/LIN/zchksy_aa.f @@ -65,7 +65,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchksy_aa_2stage.f b/TESTING/LIN/zchksy_aa_2stage.f index 8671d09796..4c8ff94f78 100644 --- a/TESTING/LIN/zchksy_aa_2stage.f +++ b/TESTING/LIN/zchksy_aa_2stage.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchksy_rk.f b/TESTING/LIN/zchksy_rk.f index 276fa569e3..3a206b30d8 100644 --- a/TESTING/LIN/zchksy_rk.f +++ b/TESTING/LIN/zchksy_rk.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchksy_rook.f b/TESTING/LIN/zchksy_rook.f index 4daac3f951..dcb5b6edd1 100644 --- a/TESTING/LIN/zchksy_rook.f +++ b/TESTING/LIN/zchksy_rook.f @@ -66,7 +66,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchktsqr.f b/TESTING/LIN/zchktsqr.f index 151a925aef..ca7828f1ad 100644 --- a/TESTING/LIN/zchktsqr.f +++ b/TESTING/LIN/zchktsqr.f @@ -76,7 +76,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> diff --git a/TESTING/LIN/zchkunhr_col.f b/TESTING/LIN/zchkunhr_col.f index 395ea178af..a6ec2f0f81 100644 --- a/TESTING/LIN/zchkunhr_col.f +++ b/TESTING/LIN/zchkunhr_col.f @@ -82,7 +82,7 @@ *> *> \param[in] NBVAL *> \verbatim -*> NBVAL is INTEGER array, dimension (NBVAL) +*> NBVAL is INTEGER array, dimension (NNB) *> The values of the blocksize NB. *> \endverbatim *> From 1d98feca5650f55304905fb00ea74c1b71de3394 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Thu, 3 Jun 2021 17:52:46 +1000 Subject: [PATCH 13/16] comments fixed: A**H was missed in xGBT05 --- TESTING/LIN/cgbt05.f | 2 +- TESTING/LIN/zgbt05.f | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/TESTING/LIN/cgbt05.f b/TESTING/LIN/cgbt05.f index 88fc405699..34b81058fd 100644 --- a/TESTING/LIN/cgbt05.f +++ b/TESTING/LIN/cgbt05.f @@ -30,7 +30,7 @@ *> CGBT05 tests the error bounds from iterative refinement for the *> computed solution to a system of equations op(A)*X = B, where A is a *> general band matrix of order n with kl subdiagonals and ku -*> superdiagonals and op(A) = A or A**T, depending on TRANS. +*> superdiagonals and op(A) = A, A**T, or A**H, depending on TRANS. *> *> RESLTS(1) = test of the error bound *> = norm(X - XACT) / ( norm(X) * FERR ) diff --git a/TESTING/LIN/zgbt05.f b/TESTING/LIN/zgbt05.f index 6c39fa2a7d..0dc6852118 100644 --- a/TESTING/LIN/zgbt05.f +++ b/TESTING/LIN/zgbt05.f @@ -30,7 +30,7 @@ *> ZGBT05 tests the error bounds from iterative refinement for the *> computed solution to a system of equations op(A)*X = B, where A is a *> general band matrix of order n with kl subdiagonals and ku -*> superdiagonals and op(A) = A or A**T, depending on TRANS. +*> superdiagonals and op(A) = A, A**T, or A**H, depending on TRANS. *> *> RESLTS(1) = test of the error bound *> = norm(X - XACT) / ( norm(X) * FERR ) From 0277cd9362d8849d0b2d7df4359d240d41de8f3b Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Thu, 3 Jun 2021 17:53:35 +1000 Subject: [PATCH 14/16] comments fixed: expression for residual in xTBT02 --- TESTING/LIN/ctbt02.f | 4 ++-- TESTING/LIN/dtbt02.f | 4 ++-- TESTING/LIN/stbt02.f | 4 ++-- TESTING/LIN/ztbt02.f | 4 ++-- 4 files changed, 8 insertions(+), 8 deletions(-) diff --git a/TESTING/LIN/ctbt02.f b/TESTING/LIN/ctbt02.f index f2a00a65b4..359eb467f0 100644 --- a/TESTING/LIN/ctbt02.f +++ b/TESTING/LIN/ctbt02.f @@ -142,7 +142,7 @@ *> \verbatim *> RESID is REAL *> The maximum over the number of right hand sides of -*> norm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ). +*> norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ). *> \endverbatim * * Authors: @@ -221,7 +221,7 @@ SUBROUTINE CTBT02( UPLO, TRANS, DIAG, N, KD, NRHS, AB, LDAB, X, END IF * * Compute the maximum over the number of right hand sides of -* norm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ). +* norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ). * RESID = ZERO DO 10 J = 1, NRHS diff --git a/TESTING/LIN/dtbt02.f b/TESTING/LIN/dtbt02.f index d570c5c137..12f4ec2b5d 100644 --- a/TESTING/LIN/dtbt02.f +++ b/TESTING/LIN/dtbt02.f @@ -135,7 +135,7 @@ *> \verbatim *> RESID is DOUBLE PRECISION *> The maximum over the number of right hand sides of -*> norm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ). +*> norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ). *> \endverbatim * * Authors: @@ -213,7 +213,7 @@ SUBROUTINE DTBT02( UPLO, TRANS, DIAG, N, KD, NRHS, AB, LDAB, X, END IF * * Compute the maximum over the number of right hand sides of -* norm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ). +* norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ). * RESID = ZERO DO 10 J = 1, NRHS diff --git a/TESTING/LIN/stbt02.f b/TESTING/LIN/stbt02.f index da23a9595e..5c3932c711 100644 --- a/TESTING/LIN/stbt02.f +++ b/TESTING/LIN/stbt02.f @@ -134,7 +134,7 @@ *> \verbatim *> RESID is REAL *> The maximum over the number of right hand sides of -*> norm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ). +*> norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ). *> \endverbatim * * Authors: @@ -212,7 +212,7 @@ SUBROUTINE STBT02( UPLO, TRANS, DIAG, N, KD, NRHS, AB, LDAB, X, END IF * * Compute the maximum over the number of right hand sides of -* norm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ). +* norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ). * RESID = ZERO DO 10 J = 1, NRHS diff --git a/TESTING/LIN/ztbt02.f b/TESTING/LIN/ztbt02.f index 8ef83cfd18..9943d94672 100644 --- a/TESTING/LIN/ztbt02.f +++ b/TESTING/LIN/ztbt02.f @@ -142,7 +142,7 @@ *> \verbatim *> RESID is DOUBLE PRECISION *> The maximum over the number of right hand sides of -*> norm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ). +*> norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ). *> \endverbatim * * Authors: @@ -221,7 +221,7 @@ SUBROUTINE ZTBT02( UPLO, TRANS, DIAG, N, KD, NRHS, AB, LDAB, X, END IF * * Compute the maximum over the number of right hand sides of -* norm(op(A)*x - b) / ( norm(op(A)) * norm(x) * EPS ). +* norm(B - op(A)*X) / ( norm(op(A)) * norm(X) * EPS ). * RESID = ZERO DO 10 J = 1, NRHS From 39f707db0c9a45e8f6aa56e65421b64cdfd33694 Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Tue, 1 Jun 2021 18:05:11 +1000 Subject: [PATCH 15/16] comments fixed: subroutine name --- SRC/zhetf2_rk.f | 4 ++-- SRC/zhetf2_rook.f | 4 ++-- SRC/zlahef_rk.f | 8 ++++---- SRC/zlahef_rook.f | 8 ++++---- 4 files changed, 12 insertions(+), 12 deletions(-) diff --git a/SRC/zhetf2_rk.f b/SRC/zhetf2_rk.f index fe737e631b..050c7993d0 100644 --- a/SRC/zhetf2_rk.f +++ b/SRC/zhetf2_rk.f @@ -417,7 +417,7 @@ SUBROUTINE ZHETF2_RK( UPLO, N, A, LDA, E, IPIV, INFO ) * * Case(2) * Equivalent to testing for -* ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX +* ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) ) @@ -770,7 +770,7 @@ SUBROUTINE ZHETF2_RK( UPLO, N, A, LDA, E, IPIV, INFO ) * * Case(2) * Equivalent to testing for -* ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX +* ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) ) diff --git a/SRC/zhetf2_rook.f b/SRC/zhetf2_rook.f index f9fd900c34..94bb297367 100644 --- a/SRC/zhetf2_rook.f +++ b/SRC/zhetf2_rook.f @@ -357,7 +357,7 @@ SUBROUTINE ZHETF2_ROOK( UPLO, N, A, LDA, IPIV, INFO ) * * Case(2) * Equivalent to testing for -* ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX +* ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) ) @@ -669,7 +669,7 @@ SUBROUTINE ZHETF2_ROOK( UPLO, N, A, LDA, IPIV, INFO ) * * Case(2) * Equivalent to testing for -* ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX +* ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( A( IMAX, IMAX ) ) ) diff --git a/SRC/zlahef_rk.f b/SRC/zlahef_rk.f index 901fb25a4b..ed66e133da 100644 --- a/SRC/zlahef_rk.f +++ b/SRC/zlahef_rk.f @@ -460,7 +460,7 @@ SUBROUTINE ZLAHEF_RK( UPLO, N, NB, KB, A, LDA, E, IPIV, W, LDW, * * Case(2) * Equivalent to testing for -* ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX +* ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( W( IMAX,KW-1 ) ) ) @@ -599,7 +599,7 @@ SUBROUTINE ZLAHEF_RK( UPLO, N, NB, KB, A, LDA, E, IPIV, W, LDW, * A(1:k-1,k) := U(1:k-1,k) = W(1:k-1,kw)/D(k,k) * * (NOTE: No need to use for Hermitian matrix -* A( K, K ) = REAL( W( K, K) ) to separately copy diagonal +* A( K, K ) = DBLE( W( K, K) ) to separately copy diagonal * element D(k,k) from W (potentially saves only one load)) CALL ZCOPY( K, W( 1, KW ), 1, A( 1, K ), 1 ) IF( K.GT.1 ) THEN @@ -912,7 +912,7 @@ SUBROUTINE ZLAHEF_RK( UPLO, N, NB, KB, A, LDA, E, IPIV, W, LDW, * * Case(2) * Equivalent to testing for -* ABS( REAL( W( IMAX,K+1 ) ) ).GE.ALPHA*ROWMAX +* ABS( DBLE( W( IMAX,K+1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( W( IMAX,K+1 ) ) ) @@ -1042,7 +1042,7 @@ SUBROUTINE ZLAHEF_RK( UPLO, N, NB, KB, A, LDA, E, IPIV, W, LDW, * A(k+1:N,k) := L(k+1:N,k) = W(k+1:N,k)/D(k,k) * * (NOTE: No need to use for Hermitian matrix -* A( K, K ) = REAL( W( K, K) ) to separately copy diagonal +* A( K, K ) = DBLE( W( K, K) ) to separately copy diagonal * element D(k,k) from W (potentially saves only one load)) CALL ZCOPY( N-K+1, W( K, K ), 1, A( K, K ), 1 ) IF( K.LT.N ) THEN diff --git a/SRC/zlahef_rook.f b/SRC/zlahef_rook.f index 6905a95db4..65f70a5978 100644 --- a/SRC/zlahef_rook.f +++ b/SRC/zlahef_rook.f @@ -370,7 +370,7 @@ SUBROUTINE ZLAHEF_ROOK( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, * * Case(2) * Equivalent to testing for -* ABS( REAL( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX +* ABS( DBLE( W( IMAX,KW-1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( W( IMAX,KW-1 ) ) ) @@ -509,7 +509,7 @@ SUBROUTINE ZLAHEF_ROOK( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, * A(1:k-1,k) := U(1:k-1,k) = W(1:k-1,kw)/D(k,k) * * (NOTE: No need to use for Hermitian matrix -* A( K, K ) = REAL( W( K, K) ) to separately copy diagonal +* A( K, K ) = DBLE( W( K, K) ) to separately copy diagonal * element D(k,k) from W (potentially saves only one load)) CALL ZCOPY( K, W( 1, KW ), 1, A( 1, K ), 1 ) IF( K.GT.1 ) THEN @@ -833,7 +833,7 @@ SUBROUTINE ZLAHEF_ROOK( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, * * Case(2) * Equivalent to testing for -* ABS( REAL( W( IMAX,K+1 ) ) ).GE.ALPHA*ROWMAX +* ABS( DBLE( W( IMAX,K+1 ) ) ).GE.ALPHA*ROWMAX * (used to handle NaN and Inf) * IF( .NOT.( ABS( DBLE( W( IMAX,K+1 ) ) ) @@ -963,7 +963,7 @@ SUBROUTINE ZLAHEF_ROOK( UPLO, N, NB, KB, A, LDA, IPIV, W, LDW, * A(k+1:N,k) := L(k+1:N,k) = W(k+1:N,k)/D(k,k) * * (NOTE: No need to use for Hermitian matrix -* A( K, K ) = REAL( W( K, K) ) to separately copy diagonal +* A( K, K ) = DBLE( W( K, K) ) to separately copy diagonal * element D(k,k) from W (potentially saves only one load)) CALL ZCOPY( N-K+1, W( K, K ), 1, A( K, K ), 1 ) IF( K.LT.N ) THEN From 8b992fbf9942527d33446892f14f29003c02e8ce Mon Sep 17 00:00:00 2001 From: Igor Zhuravlov Date: Thu, 3 Jun 2021 21:07:14 +1000 Subject: [PATCH 16/16] code fixed: wrong intrinsic function in CABS1 definition in ZTREVC3 --- SRC/ztrevc3.f | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/SRC/ztrevc3.f b/SRC/ztrevc3.f index 93c833e62a..6300e80aec 100644 --- a/SRC/ztrevc3.f +++ b/SRC/ztrevc3.f @@ -286,13 +286,13 @@ SUBROUTINE ZTREVC3( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR, $ ZGEMM, DLABAD, ZLASET, ZLACPY * .. * .. Intrinsic Functions .. - INTRINSIC ABS, DBLE, DCMPLX, CONJG, AIMAG, MAX + INTRINSIC ABS, DBLE, DCMPLX, CONJG, DIMAG, MAX * .. * .. Statement Functions .. DOUBLE PRECISION CABS1 * .. * .. Statement Function definitions .. - CABS1( CDUM ) = ABS( DBLE( CDUM ) ) + ABS( AIMAG( CDUM ) ) + CABS1( CDUM ) = ABS( DBLE( CDUM ) ) + ABS( DIMAG( CDUM ) ) * .. * .. Executable Statements .. *