diff options
author | Julien Langou <julien.langou@ucdenver.edu> | 2016-11-03 08:48:54 +0100 |
---|---|---|
committer | Julien Langou <julien.langou@ucdenver.edu> | 2016-11-03 08:48:54 +0100 |
commit | bbff7393714b29a6ff70e8c1565784cb16a0e746 (patch) | |
tree | 146d4bc89db9148abc4b19e2e453623f1a6c7bbe /TESTING | |
parent | bd47060bcb3a470520622de69ac1426ca4186f5e (diff) | |
download | lapack-bbff7393714b29a6ff70e8c1565784cb16a0e746.tar.gz lapack-bbff7393714b29a6ff70e8c1565784cb16a0e746.tar.bz2 lapack-bbff7393714b29a6ff70e8c1565784cb16a0e746.zip |
Lots of trailing whitespaces in the files of Syd. Cleaning this. No big deal.
Diffstat (limited to 'TESTING')
57 files changed, 795 insertions, 795 deletions
diff --git a/TESTING/LIN/aladhd.f b/TESTING/LIN/aladhd.f index 3a53e0bd..a45a56f3 100644 --- a/TESTING/LIN/aladhd.f +++ b/TESTING/LIN/aladhd.f @@ -279,7 +279,7 @@ * ELSE IF( LSAMEN( 2, P2, 'HA' ) ) THEN * -* HA: Hermitian +* HA: Hermitian * Aasen algorithm WRITE( IOUNIT, FMT = 9971 )PATH, 'Hermitian' * @@ -294,9 +294,9 @@ WRITE( IOUNIT, FMT = 9978 )5 WRITE( IOUNIT, FMT = 9976 )6 WRITE( IOUNIT, FMT = '( '' Messages:'' )' ) - - - ELSE IF( LSAMEN( 2, P2, 'HE' ) .OR. + + + ELSE IF( LSAMEN( 2, P2, 'HE' ) .OR. $ LSAMEN( 2, P2, 'HP' ) ) THEN * * HE: Hermitian indefinite full diff --git a/TESTING/LIN/cchkaa.f b/TESTING/LIN/cchkaa.f index f448e824..b3c6c0d8 100644 --- a/TESTING/LIN/cchkaa.f +++ b/TESTING/LIN/cchkaa.f @@ -652,7 +652,7 @@ * IF( TSTCHK ) THEN CALL CCHKHE_AASEN( DOTYPE, NN, NVAL, NNB2, NBVAL2, NNS, - $ NSVAL, THRESH, TSTERR, LDA, + $ NSVAL, THRESH, TSTERR, LDA, $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), $ WORK, RWORK, IWORK, NOUT ) @@ -661,9 +661,9 @@ END IF * IF( TSTDRV ) THEN - CALL CDRVHE_AASEN( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, - $ LDA, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), - $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), + CALL CDRVHE_AASEN( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, + $ LDA, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), + $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), $ WORK, RWORK, IWORK, NOUT ) ELSE WRITE( NOUT, FMT = 9988 )PATH @@ -987,7 +987,7 @@ * QX: QRT routines for triangular-pentagonal matrices * IF( TSTCHK ) THEN - CALL CCHKQRTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL CCHKQRTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -998,7 +998,7 @@ * TQ: LQT routines for general matrices * IF( TSTCHK ) THEN - CALL CCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL CCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -1009,7 +1009,7 @@ * XQ: LQT routines for triangular-pentagonal matrices * IF( TSTCHK ) THEN - CALL CCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL CCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -1020,7 +1020,7 @@ * TS: QR routines for tall-skinny matrices * IF( TSTCHK ) THEN - CALL CCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL CCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH diff --git a/TESTING/LIN/cchkhe_aasen.f b/TESTING/LIN/cchkhe_aasen.f index 355dd25b..33691128 100644 --- a/TESTING/LIN/cchkhe_aasen.f +++ b/TESTING/LIN/cchkhe_aasen.f @@ -225,7 +225,7 @@ * .. External Subroutines .. EXTERNAL ALAERH, ALAHD, ALASUM, XLAENV, CERRHE, CGET04, $ ZHECON, CHERFS, CHET01, CHETRF_AASEN, ZHETRI2, - $ CHETRS_AASEN, CLACPY, CLAIPD, CLARHS, CLATB4, + $ CHETRS_AASEN, CLACPY, CLAIPD, CLARHS, CLATB4, $ CLATMS, CPOT02, ZPOT03, ZPOT05 * .. * .. Intrinsic Functions .. @@ -433,7 +433,7 @@ * LWORK = ( NB+1 )*LDA SRNAMT = 'CHETRF_AASEN' - CALL CHETRF_AASEN( UPLO, N, AFAC, LDA, IWORK, AINV, + CALL CHETRF_AASEN( UPLO, N, AFAC, LDA, IWORK, AINV, $ LWORK, INFO ) * * Adjust the expected value of INFO to account for @@ -459,8 +459,8 @@ * Check error code from CHETRF and handle error. * IF( INFO.NE.K ) THEN - CALL ALAERH( PATH, 'CHETRF_AASEN', INFO, K, UPLO, - $ N, N, -1, -1, NB, IMAT, NFAIL, NERRS, + CALL ALAERH( PATH, 'CHETRF_AASEN', INFO, K, UPLO, + $ N, N, -1, -1, NB, IMAT, NFAIL, NERRS, $ NOUT ) END IF * diff --git a/TESTING/LIN/cchklqt.f b/TESTING/LIN/cchklqt.f index d6c4f7e1..04f3cbfc 100644 --- a/TESTING/LIN/cchklqt.f +++ b/TESTING/LIN/cchklqt.f @@ -2,15 +2,15 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE CCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE CCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) -* +* * .. Scalar Arguments .. * LOGICAL TSTERR * INTEGER NM, NN, NNB, NOUT @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE CCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE CCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -175,7 +175,7 @@ NB = NBVAL( K ) * * Test CGELQT and CUNMLQT -* +* IF( (NB.LE.MINMN).AND.(NB.GT.0) ) THEN CALL CLQT04( M, N, NB, RESULT ) * diff --git a/TESTING/LIN/cchklqtp.f b/TESTING/LIN/cchklqtp.f index 5e573e4c..46aa97b1 100644 --- a/TESTING/LIN/cchklqtp.f +++ b/TESTING/LIN/cchklqtp.f @@ -2,13 +2,13 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE CCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE CCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) * * .. Scalar Arguments .. @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE CCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE CCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -172,14 +172,14 @@ * MINMN = MIN( M, N ) DO L = 0, MINMN, MAX( MINMN, 1 ) -* +* * Do for each possible value of NB * DO K = 1, NNB NB = NBVAL( K ) * * Test DTPLQT and DTPMLQT -* +* IF( (NB.LE.M).AND.(NB.GT.0) ) THEN CALL CLQT05( M, N, L, NB, RESULT ) * @@ -212,4 +212,4 @@ * * End of CCHKLQTP * - END
\ No newline at end of file + END diff --git a/TESTING/LIN/cchktsqr.f b/TESTING/LIN/cchktsqr.f index 8c55f399..9a761348 100644 --- a/TESTING/LIN/cchktsqr.f +++ b/TESTING/LIN/cchktsqr.f @@ -2,15 +2,15 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE CCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE CCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) -* +* * .. Scalar Arguments .. * LOGICAL TSTERR * INTEGER NM, NN, NNB, NOUT @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE CCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE CCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -132,11 +132,11 @@ REAL RESULT( NTESTS ) * .. * .. External Subroutines .. - EXTERNAL ALAERH, ALAHD, ALASUM, CERRTSQR, + EXTERNAL ALAERH, ALAHD, ALASUM, CERRTSQR, $ CTSQR01, XLAENV * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. * .. Scalars in Common .. LOGICAL LERR, OK @@ -172,7 +172,7 @@ DO J = 1, NN N = NVAL( J ) IF (MIN(M,N).NE.0) THEN - DO INB = 1, NNB + DO INB = 1, NNB MB = NBVAL( INB ) CALL XLAENV( 1, MB ) DO IMB = 1, NNB @@ -180,7 +180,7 @@ CALL XLAENV( 2, NB ) * * Test DGEQR and DGEMQR -* +* CALL CTSQR01( 'TS', M, N, MB, NB, RESULT ) * * Print information about the tests that did not @@ -196,9 +196,9 @@ END IF END DO NRUN = NRUN + NTESTS - END DO - END DO - END IF + END DO + END DO + END IF END DO END DO * @@ -212,7 +212,7 @@ DO J = 1, NN N = NVAL( J ) IF (MIN(M,N).NE.0) THEN - DO INB = 1, NNB + DO INB = 1, NNB MB = NBVAL( INB ) CALL XLAENV( 1, MB ) DO IMB = 1, NNB @@ -220,7 +220,7 @@ CALL XLAENV( 2, NB ) * * Test DGEQR and DGEMQR -* +* CALL CTSQR01( 'SW', M, N, MB, NB, RESULT ) * * Print information about the tests that did not @@ -236,9 +236,9 @@ END IF END DO NRUN = NRUN + NTESTS - END DO - END DO - END IF + END DO + END DO + END IF END DO END DO * diff --git a/TESTING/LIN/cdrvhe_aasen.f b/TESTING/LIN/cdrvhe_aasen.f index 617f307a..ed645076 100644 --- a/TESTING/LIN/cdrvhe_aasen.f +++ b/TESTING/LIN/cdrvhe_aasen.f @@ -465,10 +465,10 @@ c END IF END IF * * Check error code from CHESV_AASEN . -* +* IF( INFO.NE.K ) THEN - CALL ALAERH( PATH, 'CHESV_AASEN', INFO, K, - $ UPLO, N, N, -1, -1, NRHS, + CALL ALAERH( PATH, 'CHESV_AASEN', INFO, K, + $ UPLO, N, N, -1, -1, NRHS, $ IMAT, NFAIL, NERRS, NOUT ) GO TO 120 ELSE IF( INFO.NE.0 ) THEN @@ -479,7 +479,7 @@ c END IF * residual. * CALL CHET01_AASEN( UPLO, N, A, LDA, AFAC, LDA, - $ IWORK, AINV, LDA, RWORK, + $ IWORK, AINV, LDA, RWORK, $ RESULT( 1 ) ) * * Compute residual of the computed solution. diff --git a/TESTING/LIN/cdrvls.f b/TESTING/LIN/cdrvls.f index af55dbaa..de9890f4 100644 --- a/TESTING/LIN/cdrvls.f +++ b/TESTING/LIN/cdrvls.f @@ -262,8 +262,8 @@ * .. * .. External Subroutines .. EXTERNAL ALAERH, ALAHD, ALASVM, CERRLS, CGELS, CGELSD, - $ CGELSS, CGELSY, CGEMM, CGETSLS, CLACPY, - $ CLARNV, CQRT13, CQRT15, CQRT16, CSSCAL, + $ CGELSS, CGELSY, CGEMM, CGETSLS, CLACPY, + $ CLARNV, CQRT13, CQRT15, CQRT16, CSSCAL, $ SAXPY, XLAENV * .. * .. Intrinsic Functions .. @@ -490,7 +490,7 @@ $ COPYB, LDB, B, LDB ) END IF SRNAMT = 'DGETSLS ' - CALL CGETSLS( TRANS, M, N, NRHS, A, + CALL CGETSLS( TRANS, M, N, NRHS, A, $ LDA, B, LDB, WORK, LWORK, INFO ) IF( INFO.NE.0 ) $ CALL ALAERH( PATH, 'CGETSLS ', INFO, 0, @@ -773,8 +773,8 @@ $ ', NB=', I4, ', type', I2, ', test(', I2, ')=', G12.5 ) 9998 FORMAT( ' M=', I5, ', N=', I5, ', NRHS=', I4, ', NB=', I4, $ ', type', I2, ', test(', I2, ')=', G12.5 ) - 9997 FORMAT( ' TRANS=''', A1,' M=', I5, ', N=', I5, ', NRHS=', I4, - $ ', MB=', I4,', NB=', I4,', type', I2, + 9997 FORMAT( ' TRANS=''', A1,' M=', I5, ', N=', I5, ', NRHS=', I4, + $ ', MB=', I4,', NB=', I4,', type', I2, $ ', test(', I2, ')=', G12.5 ) RETURN * diff --git a/TESTING/LIN/cerrlqt.f b/TESTING/LIN/cerrlqt.f index 008cb0a9..8308d109 100644 --- a/TESTING/LIN/cerrlqt.f +++ b/TESTING/LIN/cerrlqt.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE CERRLQT( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * @@ -81,7 +81,7 @@ * .. * .. External Subroutines .. EXTERNAL ALAESM, CHKXER, CGELQT3, CGELQT, - $ CGEMLQT + $ CGEMLQT * .. * .. Scalars in Common .. LOGICAL LERR, OK diff --git a/TESTING/LIN/cerrlqtp.f b/TESTING/LIN/cerrlqtp.f index 45797ddb..04dffb84 100644 --- a/TESTING/LIN/cerrlqtp.f +++ b/TESTING/LIN/cerrlqtp.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE CERRLQTP( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * @@ -81,7 +81,7 @@ * .. * .. External Subroutines .. EXTERNAL ALAESM, CHKXER, CTPLQT2, CTPLQT, - $ CTPMLQT + $ CTPMLQT * .. * .. Scalars in Common .. LOGICAL LERR, OK @@ -171,46 +171,46 @@ * SRNAMT = 'CTPMLQT' INFOT = 1 - CALL CTPMLQT( '/', 'N', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL CTPMLQT( '/', 'N', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 2 - CALL CTPMLQT( 'L', '/', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL CTPMLQT( 'L', '/', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 3 - CALL CTPMLQT( 'L', 'N', -1, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL CTPMLQT( 'L', 'N', -1, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 4 - CALL CTPMLQT( 'L', 'N', 0, -1, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL CTPMLQT( 'L', 'N', 0, -1, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 5 - CALL CTPMLQT( 'L', 'N', 0, 0, -1, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL CTPMLQT( 'L', 'N', 0, 0, -1, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) INFOT = 6 - CALL CTPMLQT( 'L', 'N', 0, 0, 0, -1, 1, A, 1, T, 1, B, 1, C, 1, + CALL CTPMLQT( 'L', 'N', 0, 0, 0, -1, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 7 - CALL CTPMLQT( 'L', 'N', 0, 0, 0, 0, 0, A, 1, T, 1, B, 1, C, 1, + CALL CTPMLQT( 'L', 'N', 0, 0, 0, 0, 0, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 9 - CALL CTPMLQT( 'R', 'N', 2, 2, 2, 1, 1, A, 1, T, 1, B, 1, C, 1, + CALL CTPMLQT( 'R', 'N', 2, 2, 2, 1, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 11 - CALL CTPMLQT( 'R', 'N', 1, 1, 1, 1, 1, A, 1, T, 0, B, 1, C, 1, + CALL CTPMLQT( 'R', 'N', 1, 1, 1, 1, 1, A, 1, T, 0, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 13 - CALL CTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 0, C, 1, + CALL CTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 0, C, 1, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 15 - CALL CTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 1, C, 0, + CALL CTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 1, C, 0, $ W, INFO ) CALL CHKXER( 'CTPMLQT', INFOT, NOUT, LERR, OK ) * diff --git a/TESTING/LIN/cerrtsqr.f b/TESTING/LIN/cerrtsqr.f index 3ca8b379..b8b42dcc 100644 --- a/TESTING/LIN/cerrtsqr.f +++ b/TESTING/LIN/cerrtsqr.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE CERRTSQR( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Zenver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Zenver +*> \author NAG Ltd. * *> \date November 2011 * diff --git a/TESTING/LIN/chet01_aasen.f b/TESTING/LIN/chet01_aasen.f index d87a6102..71abfc80 100644 --- a/TESTING/LIN/chet01_aasen.f +++ b/TESTING/LIN/chet01_aasen.f @@ -8,7 +8,7 @@ * Definition: * =========== * -* SUBROUTINE CHET01_AASEN( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, +* SUBROUTINE CHET01_AASEN( UPLO, N, A, LDA, AFAC, LDAFAC, IPIV, * C, LDC, RWORK, RESID ) * * .. Scalar Arguments .. @@ -145,7 +145,7 @@ * * .. Parameters .. COMPLEX CZERO, CONE - PARAMETER ( CZERO = ( 0.0E+0, 0.0E+0 ), + PARAMETER ( CZERO = ( 0.0E+0, 0.0E+0 ), $ CONE = ( 1.0E+0, 0.0E+0 ) ) REAL ZERO, ONE PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 ) diff --git a/TESTING/LIN/clqt04.f b/TESTING/LIN/clqt04.f index cdab2dfd..f1b722b0 100644 --- a/TESTING/LIN/clqt04.f +++ b/TESTING/LIN/clqt04.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE CLQT04(M,N,NB,RESULT) -* +* * .. Scalar Arguments .. * INTEGER M, N, NB * .. Return values .. * REAL RESULT(6) -* +* * *> \par Purpose: * ============= @@ -54,17 +54,17 @@ *> RESULT(2) = | I - Q Q^H | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -87,9 +87,9 @@ * ===================================================================== * * .. -* .. Local allocatable arrays +* .. Local allocatable arrays COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:), - $ L(:,:), RWORK(:), WORK( : ), T(:,:), + $ L(:,:), RWORK(:), WORK( : ), T(:,:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:) * * .. Parameters .. @@ -112,11 +112,11 @@ EXTERNAL SLAMCH, CLANGE, CLANSY, LSAME * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. * .. Data statements .. - DATA ISEED / 1988, 1989, 1990, 1991 / -* + DATA ISEED / 1988, 1989, 1990, 1991 / +* EPS = SLAMCH( 'Epsilon' ) K = MIN(M,N) LL = MAX(M,N) @@ -124,8 +124,8 @@ * * Dynamically allocate local arrays * - ALLOCATE ( A(M,N), AF(M,N), Q(N,N), L(LL,N), RWORK(LL), - $ WORK(LWORK), T(NB,N), C(M,N), CF(M,N), + ALLOCATE ( A(M,N), AF(M,N), Q(N,N), L(LL,N), RWORK(LL), + $ WORK(LWORK), T(NB,N), C(M,N), CF(M,N), $ D(N,M), DF(N,M) ) * * Put random numbers into A and copy to AF @@ -143,7 +143,7 @@ * Generate the n-by-n matrix Q * CALL CLASET( 'Full', N, N, CZERO, ONE, Q, N ) - CALL CGEMLQT( 'R', 'N', N, N, K, NB, AF, M, T, LDT, Q, N, + CALL CGEMLQT( 'R', 'N', N, N, K, NB, AF, M, T, LDT, Q, N, $ WORK, INFO ) * * Copy L @@ -179,7 +179,7 @@ * * Apply Q to C as Q*C * - CALL CGEMLQT( 'L', 'N', N, M, K, NB, AF, M, T, NB, DF, N, + CALL CGEMLQT( 'L', 'N', N, M, K, NB, AF, M, T, NB, DF, N, $ WORK, INFO) * * Compute |Q*D - Q*D| / |D| @@ -198,7 +198,7 @@ * * Apply Q to D as QT*D * - CALL CGEMLQT( 'L', 'C', N, M, K, NB, AF, M, T, NB, DF, N, + CALL CGEMLQT( 'L', 'C', N, M, K, NB, AF, M, T, NB, DF, N, $ WORK, INFO) * * Compute |QT*D - QT*D| / |D| @@ -209,7 +209,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,M)*DNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -221,8 +221,8 @@ * * Apply Q to C as C*Q * - CALL CGEMLQT( 'R', 'N', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, INFO) + CALL CGEMLQT( 'R', 'N', M, N, K, NB, AF, M, T, NB, CF, M, + $ WORK, INFO) * * Compute |C*Q - C*Q| / |C| * @@ -240,8 +240,8 @@ * * Apply Q to D as D*QT * - CALL CGEMLQT( 'R', 'C', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, INFO) + CALL CGEMLQT( 'R', 'C', M, N, K, NB, AF, M, T, NB, CF, M, + $ WORK, INFO) * * Compute |C*QT - C*QT| / |C| * diff --git a/TESTING/LIN/clqt05.f b/TESTING/LIN/clqt05.f index 22ffcc05..3eed9f3a 100644 --- a/TESTING/LIN/clqt05.f +++ b/TESTING/LIN/clqt05.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE CLQT05(M,N,L,NB,RESULT) -* +* * .. Scalar Arguments .. * INTEGER LWORK, M, N, L, NB, LDT * .. Return values .. * DOUBLE PRECISION RESULT(6) -* +* * *> \par Purpose: * ============= @@ -45,7 +45,7 @@ *> The number of rows of the upper trapezoidal part the *> lower test matrix. 0 <= L <= M. *> \endverbatim -*> +*> *> \param[in] NB *> \verbatim *> NB is INTEGER @@ -61,17 +61,17 @@ *> RESULT(2) = | I - Q^H Q | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -92,11 +92,11 @@ REAL RESULT(6) * * ===================================================================== -* +* * .. -* .. Local allocatable arrays +* .. Local allocatable arrays COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:), - $ R(:,:), RWORK(:), WORK( : ), T(:,:), + $ R(:,:), RWORK(:), WORK( : ), T(:,:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:) * * .. Parameters .. @@ -119,7 +119,7 @@ * .. * .. Data statements .. DATA ISEED / 1988, 1989, 1990, 1991 / -* +* EPS = SLAMCH( 'Epsilon' ) K = M N2 = M+N @@ -133,7 +133,7 @@ * Dynamically allocate all arrays * ALLOCATE(A(M,N2),AF(M,N2),Q(N2,N2),R(N2,N2),RWORK(N2), - $ WORK(LWORK),T(NB,M),C(N2,M),CF(N2,M), + $ WORK(LWORK),T(NB,M),C(N2,M),CF(N2,M), $ D(M,N2),DF(M,N2) ) * * Put random stuff into A @@ -151,7 +151,7 @@ END IF IF( L.GT.0 ) THEN DO J=1,L - CALL CLARNV( 2, ISEED, M-J+1, A( J, MIN(N+M,N+M-L+1) + CALL CLARNV( 2, ISEED, M-J+1, A( J, MIN(N+M,N+M-L+1) $ + J - 1 ) ) END DO END IF @@ -204,7 +204,7 @@ CALL CLACPY( 'Full', N2, M, C, N2, CF, N2 ) * * Apply Q to C as Q*C -* +* CALL CTPMLQT( 'L','N', N,M,K,L,NB,AF(1, NP1),M,T,LDT,CF,N2, $ CF(NP1,1),N2,WORK,INFO) * @@ -226,18 +226,18 @@ * Apply Q to C as QT*C * CALL CTPMLQT( 'L','C',N,M,K,L,NB,AF(1,NP1),M,T,LDT,CF,N2, - $ CF(NP1,1),N2,WORK,INFO) + $ CF(NP1,1),N2,WORK,INFO) * * Compute |QT*C - QT*C| / |C| * CALL CGEMM('C','N',N2,M,N2,-ONE,Q,N2,C,N2,ONE,CF,N2) RESID = CLANGE( '1', N2, M, CF, N2, RWORK ) - + IF( CNORM.GT.ZERO ) THEN RESULT( 4 ) = RESID / (EPS*MAX(1,N2)*CNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random m-by-n matrix D and a copy DF * @@ -269,8 +269,8 @@ * Apply Q to D as D*QT * CALL CTPMLQT('R','C',M,N,K,L,NB,AF(1,NP1),M,T,LDT,DF,M, - $ DF(1,NP1),M,WORK,INFO) - + $ DF(1,NP1),M,WORK,INFO) + * * Compute |D*QT - D*QT| / |D| * @@ -286,4 +286,4 @@ * DEALLOCATE ( A, AF, Q, R, RWORK, WORK, T, C, D, CF, DF) RETURN - END
\ No newline at end of file + END diff --git a/TESTING/LIN/ctsqr01.f b/TESTING/LIN/ctsqr01.f index a94f89f2..a437386b 100644 --- a/TESTING/LIN/ctsqr01.f +++ b/TESTING/LIN/ctsqr01.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE CTSQR01(TSSW, M,N, MB, NB, RESULT) -* +* * .. Scalar Arguments .. * INTEGER M, N, MB * .. Return values .. * REAL RESULT(6) -* +* * *> \par Purpose: * ============= @@ -65,17 +65,17 @@ *> RESULT(2) = | I - Q^H Q | or | I - Q Q^H | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -97,9 +97,9 @@ * ===================================================================== * * .. -* .. Local allocatable arrays +* .. Local allocatable arrays COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:), - $ R(:,:), RWORK(:), WORK( : ), T(:), + $ R(:,:), RWORK(:), WORK( : ), T(:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:), LQ(:,:) * * .. Parameters .. @@ -122,24 +122,24 @@ EXTERNAL SLAMCH, CLANGE, CLANSY, LSAME, ILAENV * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. Scalars in Common .. CHARACTER*32 srnamt * .. * .. Common blocks .. - COMMON / srnamc / srnamt + COMMON / srnamc / srnamt * .. * .. Data statements .. - DATA ISEED / 1988, 1989, 1990, 1991 / + DATA ISEED / 1988, 1989, 1990, 1991 / * * TEST TALL SKINNY OR SHORT WIDE * - TS = LSAME(TSSW, 'TS') -* + TS = LSAME(TSSW, 'TS') +* * TEST MATRICES WITH HALF OF MATRIX BEING ZEROS * TESTZEROS = .FALSE. -* +* EPS = SLAMCH( 'Epsilon' ) K = MIN(M,N) L = MAX(M,N,1) @@ -148,14 +148,14 @@ IF((K.GE.MNB).OR.(MNB.GE.L))THEN LT=MAX(1,L)*MNB+5 ELSE - LT=MAX(1,(L-K)/(MNB-K)+1)*L*MNB+5 + LT=MAX(1,(L-K)/(MNB-K)+1)*L*MNB+5 END IF * * Dynamically allocate local arrays * - ALLOCATE ( A(M,N), AF(M,N), Q(L,L), R(M,L), RWORK(L), - $ WORK(LWORK), T(LT), C(M,N), CF(M,N), + ALLOCATE ( A(M,N), AF(M,N), Q(L,L), R(M,L), RWORK(L), + $ WORK(LWORK), T(LT), C(M,N), CF(M,N), $ D(N,M), DF(N,M), LQ(L,N) ) * * Put random numbers into A and copy to AF @@ -183,7 +183,7 @@ * CALL CLASET( 'Full', M, M, CZERO, ONE, Q, M ) srnamt = 'CGEMQR' - CALL CGEMQR( 'L', 'N', M, M, K, AF, M, T, LT, Q, M, + CALL CGEMQR( 'L', 'N', M, M, K, AF, M, T, LT, Q, M, $ WORK, LWORK, INFO ) * * Copy R @@ -220,7 +220,7 @@ * Apply Q to C as Q*C * srnamt = 'CGEMQR' - CALL CGEMQR( 'L', 'N', M, N, K, AF, M, T, LT, CF, M, + CALL CGEMQR( 'L', 'N', M, N, K, AF, M, T, LT, CF, M, $ WORK, LWORK, INFO) * * Compute |Q*C - Q*C| / |C| @@ -240,7 +240,7 @@ * Apply Q to C as QT*C * srnamt = 'CGEMQR' - CALL CGEMQR( 'L', 'C', M, N, K, AF, M, T, LT, CF, M, + CALL CGEMQR( 'L', 'C', M, N, K, AF, M, T, LT, CF, M, $ WORK, LWORK, INFO) * * Compute |QT*C - QT*C| / |C| @@ -251,7 +251,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,M)*CNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -264,8 +264,8 @@ * Apply Q to D as D*Q * srnamt = 'CGEMQR' - CALL CGEMQR( 'R', 'N', N, M, K, AF, M, T, LT, DF, N, - $ WORK, LWORK, INFO) + CALL CGEMQR( 'R', 'N', N, M, K, AF, M, T, LT, DF, N, + $ WORK, LWORK, INFO) * * Compute |D*Q - D*Q| / |D| * @@ -283,8 +283,8 @@ * * Apply Q to D as D*QT * - CALL CGEMQR( 'R', 'C', N, M, K, AF, M, T, LT, DF, N, - $ WORK, LWORK, INFO) + CALL CGEMQR( 'R', 'C', N, M, K, AF, M, T, LT, DF, N, + $ WORK, LWORK, INFO) * * Compute |D*QT - D*QT| / |D| * @@ -307,7 +307,7 @@ * CALL CLASET( 'Full', N, N, CZERO, ONE, Q, N ) srnamt = 'CGEMLQ' - CALL CGEMLQ( 'R', 'N', N, N, K, AF, M, T, LT, Q, N, + CALL CGEMLQ( 'R', 'N', N, N, K, AF, M, T, LT, Q, N, $ WORK, LWORK, INFO ) * * Copy R @@ -343,7 +343,7 @@ * * Apply Q to C as Q*C * - CALL CGEMLQ( 'L', 'N', N, M, K, AF, M, T, LT, DF, N, + CALL CGEMLQ( 'L', 'N', N, M, K, AF, M, T, LT, DF, N, $ WORK, LWORK, INFO) * * Compute |Q*D - Q*D| / |D| @@ -362,7 +362,7 @@ * * Apply Q to D as QT*D * - CALL CGEMLQ( 'L', 'C', N, M, K, AF, M, T, LT, DF, N, + CALL CGEMLQ( 'L', 'C', N, M, K, AF, M, T, LT, DF, N, $ WORK, LWORK, INFO) * * Compute |QT*D - QT*D| / |D| @@ -373,7 +373,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,N)*DNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -385,8 +385,8 @@ * * Apply Q to C as C*Q * - CALL CGEMLQ( 'R', 'N', M, N, K, AF, M, T, LT, CF, M, - $ WORK, LWORK, INFO) + CALL CGEMLQ( 'R', 'N', M, N, K, AF, M, T, LT, CF, M, + $ WORK, LWORK, INFO) * * Compute |C*Q - C*Q| / |C| * @@ -404,8 +404,8 @@ * * Apply Q to D as D*QT * - CALL CGEMLQ( 'R', 'C', M, N, K, AF, M, T, LT, CF, M, - $ WORK, LWORK, INFO) + CALL CGEMLQ( 'R', 'C', M, N, K, AF, M, T, LT, CF, M, + $ WORK, LWORK, INFO) * * Compute |C*QT - C*QT| / |C| * @@ -424,4 +424,4 @@ DEALLOCATE ( A, AF, Q, R, RWORK, WORK, T, C, D, CF, DF) * RETURN - END
\ No newline at end of file + END diff --git a/TESTING/LIN/dchkaa.f b/TESTING/LIN/dchkaa.f index 2e38bb51..68a748c9 100644 --- a/TESTING/LIN/dchkaa.f +++ b/TESTING/LIN/dchkaa.f @@ -676,9 +676,9 @@ CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) * IF( TSTCHK ) THEN - CALL DCHKSY_AASEN( DOTYPE, NN, NVAL, NNB2, NBVAL2, NNS, + CALL DCHKSY_AASEN( DOTYPE, NN, NVAL, NNB2, NBVAL2, NNS, $ NSVAL, THRESH, TSTERR, LDA, - $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), + $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), $ WORK, RWORK, IWORK, NOUT ) ELSE @@ -910,7 +910,7 @@ * QX: QRT routines for triangular-pentagonal matrices * IF( TSTCHK ) THEN - CALL DCHKQRTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL DCHKQRTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -921,7 +921,7 @@ * TQ: LQT routines for general matrices * IF( TSTCHK ) THEN - CALL DCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL DCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -932,7 +932,7 @@ * XQ: LQT routines for triangular-pentagonal matrices * IF( TSTCHK ) THEN - CALL DCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL DCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -943,7 +943,7 @@ * TS: QR routines for tall-skinny matrices * IF( TSTCHK ) THEN - CALL DCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL DCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH diff --git a/TESTING/LIN/dchklqt.f b/TESTING/LIN/dchklqt.f index 1726090e..9a9ba65a 100644 --- a/TESTING/LIN/dchklqt.f +++ b/TESTING/LIN/dchklqt.f @@ -2,15 +2,15 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE DCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE DCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) -* +* * .. Scalar Arguments .. * LOGICAL TSTERR * INTEGER NM, NN, NNB, NOUT @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE DCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE DCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -175,7 +175,7 @@ NB = NBVAL( K ) * * Test DGELQT and DGEMLQT -* +* IF( (NB.LE.MINMN).AND.(NB.GT.0) ) THEN CALL DLQT04( M, N, NB, RESULT ) * diff --git a/TESTING/LIN/dchklqtp.f b/TESTING/LIN/dchklqtp.f index 1cc82ec5..d4b486ca 100644 --- a/TESTING/LIN/dchklqtp.f +++ b/TESTING/LIN/dchklqtp.f @@ -2,13 +2,13 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE DCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE DCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) * * .. Scalar Arguments .. @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE DCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE DCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -172,14 +172,14 @@ * MINMN = MIN( M, N ) DO L = 0, MINMN, MAX( MINMN, 1 ) -* +* * Do for each possible value of NB * DO K = 1, NNB NB = NBVAL( K ) * * Test DTPLQT and DTPMLQT -* +* IF( (NB.LE.M).AND.(NB.GT.0) ) THEN CALL DLQT05( M, N, L, NB, RESULT ) * diff --git a/TESTING/LIN/dchksy_aasen.f b/TESTING/LIN/dchksy_aasen.f index 041ef754..bd8e4fba 100644 --- a/TESTING/LIN/dchksy_aasen.f +++ b/TESTING/LIN/dchksy_aasen.f @@ -430,7 +430,7 @@ * SRNAMT = 'DSYTRF_AASEN' LWORK = N*NB + N - CALL DSYTRF_AASEN( UPLO, N, AFAC, LDA, IWORK, AINV, + CALL DSYTRF_AASEN( UPLO, N, AFAC, LDA, IWORK, AINV, $ LWORK, INFO ) * * Adjust the expected value of INFO to account for @@ -456,8 +456,8 @@ * Check error code from DSYTRF and handle error. * IF( INFO.NE.K ) THEN - CALL ALAERH( PATH, 'DSYTRF_AASEN', INFO, K, UPLO, - $ N, N, -1, -1, NB, IMAT, NFAIL, NERRS, + CALL ALAERH( PATH, 'DSYTRF_AASEN', INFO, K, UPLO, + $ N, N, -1, -1, NB, IMAT, NFAIL, NERRS, $ NOUT ) END IF * @@ -517,7 +517,7 @@ * SRNAMT = 'DSYTRS_AASEN' LWORK = 3*N-2 - CALL DSYTRS_AASEN( UPLO, N, NRHS, AFAC, LDA, + CALL DSYTRS_AASEN( UPLO, N, NRHS, AFAC, LDA, $ IWORK, X, LDA, WORK, LWORK, $ INFO ) * diff --git a/TESTING/LIN/dchktsqr.f b/TESTING/LIN/dchktsqr.f index 0c3de46e..0e2d0ef2 100644 --- a/TESTING/LIN/dchktsqr.f +++ b/TESTING/LIN/dchktsqr.f @@ -2,15 +2,15 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE DCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE DCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) -* +* * .. Scalar Arguments .. * LOGICAL TSTERR * INTEGER NM, NN, NNB, NOUT @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE DCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE DCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -132,11 +132,11 @@ DOUBLE PRECISION RESULT( NTESTS ) * .. * .. External Subroutines .. - EXTERNAL ALAERH, ALAHD, ALASUM, DERRTSQR, + EXTERNAL ALAERH, ALAHD, ALASUM, DERRTSQR, $ DTSQR01, XLAENV * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. * .. Scalars in Common .. LOGICAL LERR, OK @@ -172,7 +172,7 @@ DO J = 1, NN N = NVAL( J ) IF (MIN(M,N).NE.0) THEN - DO INB = 1, NNB + DO INB = 1, NNB MB = NBVAL( INB ) CALL XLAENV( 1, MB ) DO IMB = 1, NNB @@ -180,7 +180,7 @@ CALL XLAENV( 2, NB ) * * Test DGEQR and DGEMQR -* +* CALL DTSQR01( 'TS', M, N, MB, NB, RESULT ) * * Print information about the tests that did not @@ -196,9 +196,9 @@ END IF END DO NRUN = NRUN + NTESTS - END DO - END DO - END IF + END DO + END DO + END IF END DO END DO * @@ -212,7 +212,7 @@ DO J = 1, NN N = NVAL( J ) IF (MIN(M,N).NE.0) THEN - DO INB = 1, NNB + DO INB = 1, NNB MB = NBVAL( INB ) CALL XLAENV( 1, MB ) DO IMB = 1, NNB @@ -220,7 +220,7 @@ CALL XLAENV( 2, NB ) * * Test DGEQR and DGEMQR -* +* CALL DTSQR01( 'SW', M, N, MB, NB, RESULT ) * * Print information about the tests that did not @@ -236,9 +236,9 @@ END IF END DO NRUN = NRUN + NTESTS - END DO - END DO - END IF + END DO + END DO + END IF END DO END DO * @@ -254,4 +254,4 @@ * * End of DCHKQRT * - END
\ No newline at end of file + END diff --git a/TESTING/LIN/ddrvls.f b/TESTING/LIN/ddrvls.f index 2a194551..d11f910e 100644 --- a/TESTING/LIN/ddrvls.f +++ b/TESTING/LIN/ddrvls.f @@ -234,9 +234,9 @@ * .. Local Scalars .. CHARACTER TRANS CHARACTER*3 PATH - INTEGER CRANK, I, IM, IMB, IN, INB, INFO, INS, IRANK, - $ ISCALE, ITRAN, ITYPE, J, K, LDA, LDB, LDWORK, - $ LWLSY, LWORK, M, MNMIN, N, NB, NCOLS, NERRS, + INTEGER CRANK, I, IM, IMB, IN, INB, INFO, INS, IRANK, + $ ISCALE, ITRAN, ITYPE, J, K, LDA, LDB, LDWORK, + $ LWLSY, LWORK, M, MNMIN, N, NB, NCOLS, NERRS, $ NFAIL, NLVL, NRHS, NROWS, NRUN, RANK, MB, LWTS DOUBLE PRECISION EPS, NORMA, NORMB, RCOND * .. @@ -324,7 +324,7 @@ LWORK = MAX( 1, ( M+NRHS )*( N+2 ), ( N+NRHS )*( M+2 ), $ M*N+4*MNMIN+MAX( M, N ), 12*MNMIN+2*MNMIN*SMLSIZ+ $ 8*MNMIN*NLVL+MNMIN*NRHS+(SMLSIZ+1)**2,LWTS) - $ + $ * DO 120 IRANK = 1, 2 DO 110 ISCALE = 1, 3 @@ -484,7 +484,7 @@ $ COPYB, LDB, B, LDB ) END IF SRNAMT = 'DGETSLS ' - CALL DGETSLS( TRANS, M, N, NRHS, A, + CALL DGETSLS( TRANS, M, N, NRHS, A, $ LDA, B, LDB, WORK, LWORK, INFO ) IF( INFO.NE.0 ) $ CALL ALAERH( PATH, 'DGETSLS ', INFO, 0, @@ -765,8 +765,8 @@ $ ', NB=', I4, ', type', I2, ', test(', I2, ')=', G12.5 ) 9998 FORMAT( ' M=', I5, ', N=', I5, ', NRHS=', I4, ', NB=', I4, $ ', type', I2, ', test(', I2, ')=', G12.5 ) - 9997 FORMAT( ' TRANS=''', A1,' M=', I5, ', N=', I5, ', NRHS=', I4, - $ ', MB=', I4,', NB=', I4,', type', I2, + 9997 FORMAT( ' TRANS=''', A1,' M=', I5, ', N=', I5, ', NRHS=', I4, + $ ', MB=', I4,', NB=', I4,', type', I2, $ ', test(', I2, ')=', G12.5 ) RETURN * diff --git a/TESTING/LIN/ddrvsy_aasen.f b/TESTING/LIN/ddrvsy_aasen.f index a3520cb1..ab8b136c 100644 --- a/TESTING/LIN/ddrvsy_aasen.f +++ b/TESTING/LIN/ddrvsy_aasen.f @@ -467,7 +467,7 @@ c END IF * residual. * CALL DSYT01_AASEN( UPLO, N, A, LDA, AFAC, LDA, - $ IWORK, AINV, LDA, RWORK, + $ IWORK, AINV, LDA, RWORK, $ RESULT( 1 ) ) * * Compute residual of the computed solution. diff --git a/TESTING/LIN/derrlqt.f b/TESTING/LIN/derrlqt.f index 5a768f01..926f4197 100644 --- a/TESTING/LIN/derrlqt.f +++ b/TESTING/LIN/derrlqt.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE DERRLQT( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * @@ -81,7 +81,7 @@ * .. * .. External Subroutines .. EXTERNAL ALAESM, CHKXER, DGELQT3, DGELQT, - $ DGEMLQT + $ DGEMLQT * .. * .. Scalars in Common .. LOGICAL LERR, OK diff --git a/TESTING/LIN/derrlqtp.f b/TESTING/LIN/derrlqtp.f index ae118af9..b642c349 100644 --- a/TESTING/LIN/derrlqtp.f +++ b/TESTING/LIN/derrlqtp.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE DERRLQTP( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * @@ -81,7 +81,7 @@ * .. * .. External Subroutines .. EXTERNAL ALAESM, CHKXER, DTPLQT2, DTPLQT, - $ DTPMLQT + $ DTPMLQT * .. * .. Scalars in Common .. LOGICAL LERR, OK @@ -171,46 +171,46 @@ * SRNAMT = 'DTPMLQT' INFOT = 1 - CALL DTPMLQT( '/', 'N', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL DTPMLQT( '/', 'N', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 2 - CALL DTPMLQT( 'L', '/', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL DTPMLQT( 'L', '/', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 3 - CALL DTPMLQT( 'L', 'N', -1, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL DTPMLQT( 'L', 'N', -1, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 4 - CALL DTPMLQT( 'L', 'N', 0, -1, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL DTPMLQT( 'L', 'N', 0, -1, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 5 - CALL DTPMLQT( 'L', 'N', 0, 0, -1, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL DTPMLQT( 'L', 'N', 0, 0, -1, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) INFOT = 6 - CALL DTPMLQT( 'L', 'N', 0, 0, 0, -1, 1, A, 1, T, 1, B, 1, C, 1, + CALL DTPMLQT( 'L', 'N', 0, 0, 0, -1, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 7 - CALL DTPMLQT( 'L', 'N', 0, 0, 0, 0, 0, A, 1, T, 1, B, 1, C, 1, + CALL DTPMLQT( 'L', 'N', 0, 0, 0, 0, 0, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 9 - CALL DTPMLQT( 'R', 'N', 2, 2, 2, 1, 1, A, 1, T, 1, B, 1, C, 1, + CALL DTPMLQT( 'R', 'N', 2, 2, 2, 1, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 11 - CALL DTPMLQT( 'R', 'N', 1, 1, 1, 1, 1, A, 1, T, 0, B, 1, C, 1, + CALL DTPMLQT( 'R', 'N', 1, 1, 1, 1, 1, A, 1, T, 0, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 13 - CALL DTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 0, C, 1, + CALL DTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 0, C, 1, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 15 - CALL DTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 1, C, 0, + CALL DTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 1, C, 0, $ W, INFO ) CALL CHKXER( 'DTPMLQT', INFOT, NOUT, LERR, OK ) * diff --git a/TESTING/LIN/derrtsqr.f b/TESTING/LIN/derrtsqr.f index aa9f3674..4a5ad5e6 100644 --- a/TESTING/LIN/derrtsqr.f +++ b/TESTING/LIN/derrtsqr.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE DERRTSQR( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * diff --git a/TESTING/LIN/dlqt04.f b/TESTING/LIN/dlqt04.f index 216ef3ea..9e6e11cf 100644 --- a/TESTING/LIN/dlqt04.f +++ b/TESTING/LIN/dlqt04.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE DLQT04(M,N,NB,RESULT) -* +* * .. Scalar Arguments .. * INTEGER M, N, NB, LDT * .. Return values .. * DOUBLE PRECISION RESULT(6) -* +* * *> \par Purpose: * ============= @@ -54,17 +54,17 @@ *> RESULT(2) = | I - Q Q^H | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -87,9 +87,9 @@ * ===================================================================== * * .. -* .. Local allocatable arrays +* .. Local allocatable arrays DOUBLE PRECISION, ALLOCATABLE :: AF(:,:), Q(:,:), - $ L(:,:), RWORK(:), WORK( : ), T(:,:), + $ L(:,:), RWORK(:), WORK( : ), T(:,:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:) * * .. Parameters .. @@ -109,11 +109,11 @@ EXTERNAL DLAMCH, DLANGE, DLANSY, LSAME * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. * .. Data statements .. - DATA ISEED / 1988, 1989, 1990, 1991 / -* + DATA ISEED / 1988, 1989, 1990, 1991 / +* EPS = DLAMCH( 'Epsilon' ) K = MIN(M,N) LL = MAX(M,N) @@ -121,8 +121,8 @@ * * Dynamically allocate local arrays * - ALLOCATE ( A(M,N), AF(M,N), Q(N,N), L(LL,N), RWORK(LL), - $ WORK(LWORK), T(NB,N), C(M,N), CF(M,N), + ALLOCATE ( A(M,N), AF(M,N), Q(N,N), L(LL,N), RWORK(LL), + $ WORK(LWORK), T(NB,N), C(M,N), CF(M,N), $ D(N,M), DF(N,M) ) * * Put random numbers into A and copy to AF @@ -140,7 +140,7 @@ * Generate the n-by-n matrix Q * CALL DLASET( 'Full', N, N, ZERO, ONE, Q, N ) - CALL DGEMLQT( 'R', 'N', N, N, K, NB, AF, M, T, LDT, Q, N, + CALL DGEMLQT( 'R', 'N', N, N, K, NB, AF, M, T, LDT, Q, N, $ WORK, INFO ) * * Copy R @@ -176,7 +176,7 @@ * * Apply Q to C as Q*C * - CALL DGEMLQT( 'L', 'N', N, M, K, NB, AF, M, T, NB, DF, N, + CALL DGEMLQT( 'L', 'N', N, M, K, NB, AF, M, T, NB, DF, N, $ WORK, INFO) * * Compute |Q*D - Q*D| / |D| @@ -195,7 +195,7 @@ * * Apply Q to D as QT*D * - CALL DGEMLQT( 'L', 'T', N, M, K, NB, AF, M, T, NB, DF, N, + CALL DGEMLQT( 'L', 'T', N, M, K, NB, AF, M, T, NB, DF, N, $ WORK, INFO) * * Compute |QT*D - QT*D| / |D| @@ -206,7 +206,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,M)*DNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -218,8 +218,8 @@ * * Apply Q to C as C*Q * - CALL DGEMLQT( 'R', 'N', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, INFO) + CALL DGEMLQT( 'R', 'N', M, N, K, NB, AF, M, T, NB, CF, M, + $ WORK, INFO) * * Compute |C*Q - C*Q| / |C| * @@ -237,8 +237,8 @@ * * Apply Q to D as D*QT * - CALL DGEMLQT( 'R', 'T', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, INFO) + CALL DGEMLQT( 'R', 'T', M, N, K, NB, AF, M, T, NB, CF, M, + $ WORK, INFO) * * Compute |C*QT - C*QT| / |C| * diff --git a/TESTING/LIN/dlqt05.f b/TESTING/LIN/dlqt05.f index b357dcb5..88681eba 100644 --- a/TESTING/LIN/dlqt05.f +++ b/TESTING/LIN/dlqt05.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE DLQT05(M,N,L,NB,RESULT) -* +* * .. Scalar Arguments .. * INTEGER LWORK, M, N, L, NB, LDT * .. Return values .. * DOUBLE PRECISION RESULT(6) -* +* * *> \par Purpose: * ============= @@ -45,7 +45,7 @@ *> The number of rows of the upper trapezoidal part the *> lower test matrix. 0 <= L <= M. *> \endverbatim -*> +*> *> \param[in] NB *> \verbatim *> NB is INTEGER @@ -61,17 +61,17 @@ *> RESULT(2) = | I - Q^H Q | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -92,11 +92,11 @@ DOUBLE PRECISION RESULT(6) * * ===================================================================== -* +* * .. -* .. Local allocatable arrays +* .. Local allocatable arrays DOUBLE PRECISION, ALLOCATABLE :: AF(:,:), Q(:,:), - $ R(:,:), RWORK(:), WORK( : ), T(:,:), + $ R(:,:), RWORK(:), WORK( : ), T(:,:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:) * * .. Parameters .. @@ -117,7 +117,7 @@ * .. * .. Data statements .. DATA ISEED / 1988, 1989, 1990, 1991 / -* +* EPS = DLAMCH( 'Epsilon' ) K = M N2 = M+N @@ -131,7 +131,7 @@ * Dynamically allocate all arrays * ALLOCATE(A(M,N2),AF(M,N2),Q(N2,N2),R(N2,N2),RWORK(N2), - $ WORK(LWORK),T(NB,M),C(N2,M),CF(N2,M), + $ WORK(LWORK),T(NB,M),C(N2,M),CF(N2,M), $ D(M,N2),DF(M,N2) ) * * Put random stuff into A @@ -149,7 +149,7 @@ END IF IF( L.GT.0 ) THEN DO J=1,L - CALL DLARNV( 2, ISEED, M-J+1, A( J, MIN(N+M,N+M-L+1) + CALL DLARNV( 2, ISEED, M-J+1, A( J, MIN(N+M,N+M-L+1) $ + J - 1 ) ) END DO END IF @@ -201,7 +201,7 @@ CALL DLACPY( 'Full', N2, M, C, N2, CF, N2 ) * * Apply Q to C as Q*C -* +* CALL DTPMLQT( 'L','N', N,M,K,L,NB,AF(1, NP1),M,T,LDT,CF,N2, $ CF(NP1,1),N2,WORK,INFO) * @@ -223,18 +223,18 @@ * Apply Q to C as QT*C * CALL DTPMLQT( 'L','T',N,M,K,L,NB,AF(1,NP1),M,T,LDT,CF,N2, - $ CF(NP1,1),N2,WORK,INFO) + $ CF(NP1,1),N2,WORK,INFO) * * Compute |QT*C - QT*C| / |C| * CALL DGEMM('T','N',N2,M,N2,-ONE,Q,N2,C,N2,ONE,CF,N2) RESID = DLANGE( '1', N2, M, CF, N2, RWORK ) - + IF( CNORM.GT.ZERO ) THEN RESULT( 4 ) = RESID / (EPS*MAX(1,N2)*CNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random m-by-n matrix D and a copy DF * @@ -266,8 +266,8 @@ * Apply Q to D as D*QT * CALL DTPMLQT('R','T',M,N,K,L,NB,AF(1,NP1),M,T,LDT,DF,M, - $ DF(1,NP1),M,WORK,INFO) - + $ DF(1,NP1),M,WORK,INFO) + * * Compute |D*QT - D*QT| / |D| * @@ -283,4 +283,4 @@ * DEALLOCATE ( A, AF, Q, R, RWORK, WORK, T, C, D, CF, DF) RETURN - END
\ No newline at end of file + END diff --git a/TESTING/LIN/dtplqt.f b/TESTING/LIN/dtplqt.f index 27965442..a233eb18 100644 --- a/TESTING/LIN/dtplqt.f +++ b/TESTING/LIN/dtplqt.f @@ -2,41 +2,41 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * *> \htmlonly -*> Download DTPQRT + dependencies -*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dtplqt.f"> -*> [TGZ]</a> -*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dtplqt.f"> -*> [ZIP]</a> -*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dtplqt.f"> +*> Download DTPQRT + dependencies +*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dtplqt.f"> +*> [TGZ]</a> +*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dtplqt.f"> +*> [ZIP]</a> +*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dtplqt.f"> *> [TXT]</a> -*> \endhtmlonly +*> \endhtmlonly * * Definition: * =========== * * SUBROUTINE DTPLQT( M, N, L, MB, A, LDA, B, LDB, T, LDT, WORK, * INFO ) -* +* * .. Scalar Arguments .. * INTEGER INFO, LDA, LDB, LDT, N, M, L, MB * .. * .. Array Arguments .. * DOUBLE PRECISION A( LDA, * ), B( LDB, * ), T( LDT, * ), WORK( * ) * .. -* +* * *> \par Purpose: * ============= *> *> \verbatim *> -*> DTPLQT computes a blocked LQ factorization of a real -*> "triangular-pentagonal" matrix C, which is composed of a -*> triangular block A and pentagonal block B, using the compact +*> DTPLQT computes a blocked LQ factorization of a real +*> "triangular-pentagonal" matrix C, which is composed of a +*> triangular block A and pentagonal block B, using the compact *> WY representation for Q. *> \endverbatim * @@ -47,7 +47,7 @@ *> \verbatim *> M is INTEGER *> The number of rows of the matrix B, and the order of the -*> triangular matrix A. +*> triangular matrix A. *> M >= 0. *> \endverbatim *> @@ -88,7 +88,7 @@ *> \param[in,out] B *> \verbatim *> B is DOUBLE PRECISION array, dimension (LDB,N) -*> On entry, the pentagonal M-by-N matrix B. The first N-L columns +*> On entry, the pentagonal M-by-N matrix B. The first N-L columns *> are rectangular, and the last L columns are lower trapezoidal. *> On exit, B contains the pentagonal matrix V. See Further Details. *> \endverbatim @@ -105,7 +105,7 @@ *> The lower triangular block reflectors stored in compact form *> as a sequence of upper triangular blocks. See Further Details. *> \endverbatim -*> +*> *> \param[in] LDT *> \verbatim *> LDT is INTEGER @@ -127,10 +127,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2013 * @@ -141,45 +141,45 @@ *> *> \verbatim *> -*> The input matrix C is a M-by-(M+N) matrix +*> The input matrix C is a M-by-(M+N) matrix *> *> C = [ A ] [ B ] -*> +*> *> *> where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal *> matrix consisting of a M-by-(N-L) rectangular matrix B1 on left of a M-by-L *> upper trapezoidal matrix B2: -*> [ B ] = [ B1 ] [ B2 ] +*> [ B ] = [ B1 ] [ B2 ] *> [ B1 ] <- M-by-(N-L) rectangular *> [ B2 ] <- M-by-L upper trapezoidal. *> *> The lower trapezoidal matrix B2 consists of the first L columns of a -*> N-by-N lower triangular matrix, where 0 <= L <= MIN(M,N). If L=0, -*> B is rectangular M-by-N; if M=L=N, B is lower triangular. +*> N-by-N lower triangular matrix, where 0 <= L <= MIN(M,N). If L=0, +*> B is rectangular M-by-N; if M=L=N, B is lower triangular. *> *> The matrix W stores the elementary reflectors H(i) in the i-th row *> above the diagonal (of A) in the M-by-(M+N) input matrix C -*> [ C ] = [ A ] [ B ] +*> [ C ] = [ A ] [ B ] *> [ A ] <- lower triangular N-by-N *> [ B ] <- M-by-N pentagonal *> *> so that W can be represented as -*> [ W ] = [ I ] [ V ] +*> [ W ] = [ I ] [ V ] *> [ I ] <- identity, N-by-N *> [ V ] <- M-by-N, same form as B. *> *> Thus, all of information needed for W is contained on exit in B, which -*> we call V above. Note that V has the same form as B; that is, -*> [ V ] = [ V1 ] [ V2 ] +*> we call V above. Note that V has the same form as B; that is, +*> [ V ] = [ V1 ] [ V2 ] *> [ V1 ] <- M-by-(N-L) rectangular *> [ V2 ] <- M-by-L lower trapezoidal. *> -*> The rows of V represent the vectors which define the H(i)'s. +*> The rows of V represent the vectors which define the H(i)'s. *> *> The number of blocks is B = ceiling(M/MB), where each -*> block is of order MB except for the last block, which is of order +*> block is of order MB except for the last block, which is of order *> IB = M - (M-1)*MB. For each of the B blocks, a upper triangular block -*> reflector factor is computed: T1, T2, ..., TB. The MB-by-MB (and IB-by-IB +*> reflector factor is computed: T1, T2, ..., TB. The MB-by-MB (and IB-by-IB *> for the last block) T's are stored in the MB-by-N matrix T as *> *> T = [T1 T2 ... TB]. @@ -240,7 +240,7 @@ IF( M.EQ.0 .OR. N.EQ.0 ) RETURN * DO I = 1, M, MB -* +* * Compute the QR factorization of the current block * IB = MIN( M-I+1, MB ) @@ -251,20 +251,20 @@ LB = NB-N+L-I+1 END IF * - CALL DTPLQT2( IB, NB, LB, A(I,I), LDA, B( I, 1 ), LDB, + CALL DTPLQT2( IB, NB, LB, A(I,I), LDA, B( I, 1 ), LDB, $ T(1, I ), LDT, IINFO ) * * Update by applying H**T to B(I+IB:M,:) from the right * IF( I+IB.LE.M ) THEN CALL DTPRFB( 'R', 'N', 'F', 'R', M-I-IB+1, NB, IB, LB, - $ B( I, 1 ), LDB, T( 1, I ), LDT, - $ A( I+IB, I ), LDA, B( I+IB, 1 ), LDB, + $ B( I, 1 ), LDB, T( 1, I ), LDT, + $ A( I+IB, I ), LDA, B( I+IB, 1 ), LDB, $ WORK, M-I-IB+1) END IF END DO RETURN -* +* * End of DTPLQT * END diff --git a/TESTING/LIN/dtsqr01.f b/TESTING/LIN/dtsqr01.f index 29d4b63e..a9ac1635 100644 --- a/TESTING/LIN/dtsqr01.f +++ b/TESTING/LIN/dtsqr01.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE DTSQR01(TSSW, M,N, MB, NB, RESULT) -* +* * .. Scalar Arguments .. * INTEGER M, N, MB * .. Return values .. * DOUBLE PRECISION RESULT(6) -* +* * *> \par Purpose: * ============= @@ -65,17 +65,17 @@ *> RESULT(2) = | I - Q^H Q | or | I - Q Q^H | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -99,9 +99,9 @@ * ===================================================================== * * .. -* .. Local allocatable arrays +* .. Local allocatable arrays DOUBLE PRECISION, ALLOCATABLE :: AF(:,:), Q(:,:), - $ R(:,:), RWORK(:), WORK( : ), T(:), + $ R(:,:), RWORK(:), WORK( : ), T(:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:), LQ(:,:) * * .. Parameters .. @@ -123,24 +123,24 @@ EXTERNAL DLAMCH, DLANGE, DLANSY, LSAME, ILAENV * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. Scalars in Common .. CHARACTER*32 srnamt * .. * .. Common blocks .. - COMMON / srnamc / srnamt + COMMON / srnamc / srnamt * .. * .. Data statements .. - DATA ISEED / 1988, 1989, 1990, 1991 / + DATA ISEED / 1988, 1989, 1990, 1991 / * * TEST TALL SKINNY OR SHORT WIDE * TS = LSAME(TSSW, 'TS') -* +* * TEST MATRICES WITH HALF OF MATRIX BEING ZEROS * TESTZEROS = .FALSE. -* +* EPS = DLAMCH( 'Epsilon' ) K = MIN(M,N) L = MAX(M,N,1) @@ -149,14 +149,14 @@ IF((K.GE.MNB).OR.(MNB.GE.L))THEN LT=MAX(1,L)*MNB+5 ELSE - LT=MAX(1,(L-K)/(MNB-K)+1)*L*MNB+5 + LT=MAX(1,(L-K)/(MNB-K)+1)*L*MNB+5 END IF * * Dynamically allocate local arrays * - ALLOCATE ( A(M,N), AF(M,N), Q(L,L), R(M,L), RWORK(L), - $ WORK(LWORK), T(LT), C(M,N), CF(M,N), + ALLOCATE ( A(M,N), AF(M,N), Q(L,L), R(M,L), RWORK(L), + $ WORK(LWORK), T(LT), C(M,N), CF(M,N), $ D(N,M), DF(N,M), LQ(L,N) ) * * Put random numbers into A and copy to AF @@ -184,7 +184,7 @@ * CALL DLASET( 'Full', M, M, ZERO, ONE, Q, M ) srnamt = 'DGEMQR' - CALL DGEMQR( 'L', 'N', M, M, K, AF, M, T, LT, Q, M, + CALL DGEMQR( 'L', 'N', M, M, K, AF, M, T, LT, Q, M, $ WORK, LWORK, INFO ) * * Copy R @@ -221,7 +221,7 @@ * Apply Q to C as Q*C * srnamt = 'DGEMQR' - CALL DGEMQR( 'L', 'N', M, N, K, AF, M, T, LT, CF, M, + CALL DGEMQR( 'L', 'N', M, N, K, AF, M, T, LT, CF, M, $ WORK, LWORK, INFO) * * Compute |Q*C - Q*C| / |C| @@ -241,7 +241,7 @@ * Apply Q to C as QT*C * srnamt = 'DGEMQR' - CALL DGEMQR( 'L', 'T', M, N, K, AF, M, T, LT, CF, M, + CALL DGEMQR( 'L', 'T', M, N, K, AF, M, T, LT, CF, M, $ WORK, LWORK, INFO) * * Compute |QT*C - QT*C| / |C| @@ -252,7 +252,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,M)*CNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -265,8 +265,8 @@ * Apply Q to D as D*Q * srnamt = 'DGEMQR' - CALL DGEMQR( 'R', 'N', N, M, K, AF, M, T, LT, DF, N, - $ WORK, LWORK, INFO) + CALL DGEMQR( 'R', 'N', N, M, K, AF, M, T, LT, DF, N, + $ WORK, LWORK, INFO) * * Compute |D*Q - D*Q| / |D| * @@ -284,8 +284,8 @@ * * Apply Q to D as D*QT * - CALL DGEMQR( 'R', 'T', N, M, K, AF, M, T, LT, DF, N, - $ WORK, LWORK, INFO) + CALL DGEMQR( 'R', 'T', N, M, K, AF, M, T, LT, DF, N, + $ WORK, LWORK, INFO) * * Compute |D*QT - D*QT| / |D| * @@ -308,7 +308,7 @@ * CALL DLASET( 'Full', N, N, ZERO, ONE, Q, N ) srnamt = 'DGEMLQ' - CALL DGEMLQ( 'R', 'N', N, N, K, AF, M, T, LT, Q, N, + CALL DGEMLQ( 'R', 'N', N, N, K, AF, M, T, LT, Q, N, $ WORK, LWORK, INFO ) * * Copy R @@ -344,7 +344,7 @@ * * Apply Q to C as Q*C * - CALL DGEMLQ( 'L', 'N', N, M, K, AF, M, T, LT, DF, N, + CALL DGEMLQ( 'L', 'N', N, M, K, AF, M, T, LT, DF, N, $ WORK, LWORK, INFO) * * Compute |Q*D - Q*D| / |D| @@ -363,7 +363,7 @@ * * Apply Q to D as QT*D * - CALL DGEMLQ( 'L', 'T', N, M, K, AF, M, T, LT, DF, N, + CALL DGEMLQ( 'L', 'T', N, M, K, AF, M, T, LT, DF, N, $ WORK, LWORK, INFO) * * Compute |QT*D - QT*D| / |D| @@ -374,7 +374,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,N)*DNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -386,8 +386,8 @@ * * Apply Q to C as C*Q * - CALL DGEMLQ( 'R', 'N', M, N, K, AF, M, T, LT, CF, M, - $ WORK, LWORK, INFO) + CALL DGEMLQ( 'R', 'N', M, N, K, AF, M, T, LT, CF, M, + $ WORK, LWORK, INFO) * * Compute |C*Q - C*Q| / |C| * @@ -405,8 +405,8 @@ * * Apply Q to D as D*QT * - CALL DGEMLQ( 'R', 'T', M, N, K, AF, M, T, LT, CF, M, - $ WORK, LWORK, INFO) + CALL DGEMLQ( 'R', 'T', M, N, K, AF, M, T, LT, CF, M, + $ WORK, LWORK, INFO) * * Compute |C*QT - C*QT| / |C| * @@ -425,4 +425,4 @@ DEALLOCATE ( A, AF, Q, R, RWORK, WORK, T, C, D, CF, DF) * RETURN - END
\ No newline at end of file + END diff --git a/TESTING/LIN/schkaa.f b/TESTING/LIN/schkaa.f index b58a3a66..de4efb60 100644 --- a/TESTING/LIN/schkaa.f +++ b/TESTING/LIN/schkaa.f @@ -906,7 +906,7 @@ * QX: QRT routines for triangular-pentagonal matrices * IF( TSTCHK ) THEN - CALL SCHKQRTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL SCHKQRTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -917,7 +917,7 @@ * TQ: LQT routines for general matrices * IF( TSTCHK ) THEN - CALL SCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL SCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -928,7 +928,7 @@ * XQ: LQT routines for triangular-pentagonal matrices * IF( TSTCHK ) THEN - CALL SCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL SCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -939,7 +939,7 @@ * TS: QR routines for tall-skinny matrices * IF( TSTCHK ) THEN - CALL SCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL SCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH diff --git a/TESTING/LIN/schklqt.f b/TESTING/LIN/schklqt.f index fd449b1a..5e96546d 100644 --- a/TESTING/LIN/schklqt.f +++ b/TESTING/LIN/schklqt.f @@ -2,15 +2,15 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE SCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE SCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) -* +* * .. Scalar Arguments .. * LOGICAL TSTERR * INTEGER NM, NN, NNB, NOUT @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE SCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE SCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -175,7 +175,7 @@ NB = NBVAL( K ) * * Test DGELQT and DGEMLQT -* +* IF( (NB.LE.MINMN).AND.(NB.GT.0) ) THEN CALL SLQT04( M, N, NB, RESULT ) * diff --git a/TESTING/LIN/schklqtp.f b/TESTING/LIN/schklqtp.f index d85ef8d1..be6b84ba 100644 --- a/TESTING/LIN/schklqtp.f +++ b/TESTING/LIN/schklqtp.f @@ -2,13 +2,13 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE SCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE SCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) * * .. Scalar Arguments .. @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE SCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE SCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -172,14 +172,14 @@ * MINMN = MIN( M, N ) DO L = 0, MINMN, MAX( MINMN, 1 ) -* +* * Do for each possible value of NB * DO K = 1, NNB NB = NBVAL( K ) * * Test DTPLQT and DTPMLQT -* +* IF( (NB.LE.M).AND.(NB.GT.0) ) THEN CALL SLQT05( M, N, L, NB, RESULT ) * diff --git a/TESTING/LIN/schksy_aasen.f b/TESTING/LIN/schksy_aasen.f index 659e3fd3..949f8118 100644 --- a/TESTING/LIN/schksy_aasen.f +++ b/TESTING/LIN/schksy_aasen.f @@ -431,7 +431,7 @@ * SRNAMT = 'SSYTRF_AASEN' LWORK = N*NB + N - CALL SSYTRF_AASEN( UPLO, N, AFAC, LDA, IWORK, AINV, + CALL SSYTRF_AASEN( UPLO, N, AFAC, LDA, IWORK, AINV, $ LWORK, INFO ) * * Adjust the expected value of INFO to account for @@ -457,8 +457,8 @@ * Check error code from SSYTRF and handle error. * IF( INFO.NE.K ) THEN - CALL ALAERH( PATH, 'SSYTRF_AASEN', INFO, K, UPLO, - $ N, N, -1, -1, NB, IMAT, NFAIL, NERRS, + CALL ALAERH( PATH, 'SSYTRF_AASEN', INFO, K, UPLO, + $ N, N, -1, -1, NB, IMAT, NFAIL, NERRS, $ NOUT ) END IF * @@ -518,7 +518,7 @@ * SRNAMT = 'SSYTRS_AASEN' LWORK = 3*N-2 - CALL SSYTRS_AASEN( UPLO, N, NRHS, AFAC, LDA, + CALL SSYTRS_AASEN( UPLO, N, NRHS, AFAC, LDA, $ IWORK, X, LDA, WORK, LWORK, $ INFO ) * @@ -526,7 +526,7 @@ * IF( INFO.NE.0 ) THEN CALL ALAERH( PATH, 'SSYTRS_AASEN', INFO, 0, - $ UPLO, N, N, -1, -1, NRHS, IMAT, + $ UPLO, N, N, -1, -1, NRHS, IMAT, $ NFAIL, NERRS, NOUT ) END IF * diff --git a/TESTING/LIN/schktsqr.f b/TESTING/LIN/schktsqr.f index a4303143..3bb238ff 100644 --- a/TESTING/LIN/schktsqr.f +++ b/TESTING/LIN/schktsqr.f @@ -2,15 +2,15 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE SCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE SCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) -* +* * .. Scalar Arguments .. * LOGICAL TSTERR * INTEGER NM, NN, NNB, NOUT @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE SCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE SCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -132,11 +132,11 @@ REAL RESULT( NTESTS ) * .. * .. External Subroutines .. - EXTERNAL ALAERH, ALAHD, ALASUM, SERRTSQR, + EXTERNAL ALAERH, ALAHD, ALASUM, SERRTSQR, $ STSQR01, XLAENV * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. * .. Scalars in Common .. LOGICAL LERR, OK @@ -172,7 +172,7 @@ DO J = 1, NN N = NVAL( J ) IF (MIN(M,N).NE.0) THEN - DO INB = 1, NNB + DO INB = 1, NNB MB = NBVAL( INB ) CALL XLAENV( 1, MB ) DO IMB = 1, NNB @@ -180,7 +180,7 @@ CALL XLAENV( 2, NB ) * * Test SGEQR and SGEMQR -* +* CALL STSQR01('TS', M, N, MB, NB, RESULT ) * * Print information about the tests that did not @@ -197,7 +197,7 @@ END DO NRUN = NRUN + NTESTS END DO - END DO + END DO END IF END DO END DO @@ -212,7 +212,7 @@ DO J = 1, NN N = NVAL( J ) IF (MIN(M,N).NE.0) THEN - DO INB = 1, NNB + DO INB = 1, NNB MB = NBVAL( INB ) CALL XLAENV( 1, MB ) DO IMB = 1, NNB @@ -220,7 +220,7 @@ CALL XLAENV( 2, NB ) * * Test SGEQR and SGEMQR -* +* CALL STSQR01('SW', M, N, MB, NB, RESULT ) * * Print information about the tests that did not @@ -236,9 +236,9 @@ END IF END DO NRUN = NRUN + NTESTS - END DO - END DO - END IF + END DO + END DO + END IF END DO END DO * diff --git a/TESTING/LIN/sdrvls.f b/TESTING/LIN/sdrvls.f index 372ec9fc..03598937 100644 --- a/TESTING/LIN/sdrvls.f +++ b/TESTING/LIN/sdrvls.f @@ -234,9 +234,9 @@ * .. Local Scalars .. CHARACTER TRANS CHARACTER*3 PATH - INTEGER CRANK, I, IM, IMB, IN, INB, INFO, INS, IRANK, - $ ISCALE, ITRAN, ITYPE, J, K, LDA, LDB, LDWORK, - $ LWLSY, LWORK, M, MNMIN, N, NB, NCOLS, NERRS, + INTEGER CRANK, I, IM, IMB, IN, INB, INFO, INS, IRANK, + $ ISCALE, ITRAN, ITYPE, J, K, LDA, LDB, LDWORK, + $ LWLSY, LWORK, M, MNMIN, N, NB, NCOLS, NERRS, $ NFAIL, NLVL, NRHS, NROWS, NRUN, RANK, MB, LWTS REAL EPS, NORMA, NORMB, RCOND * .. @@ -324,7 +324,7 @@ LWORK = MAX( 1, ( M+NRHS )*( N+2 ), ( N+NRHS )*( M+2 ), $ M*N+4*MNMIN+MAX( M, N ), 12*MNMIN+2*MNMIN*SMLSIZ+ $ 8*MNMIN*NLVL+MNMIN*NRHS+(SMLSIZ+1)**2,LWTS) - $ + $ * DO 120 IRANK = 1, 2 DO 110 ISCALE = 1, 3 @@ -446,7 +446,7 @@ CALL XLAENV( 1, MB ) DO 62 IMB = 1, NNB NB = NBVAL( IMB ) - CALL XLAENV( 2, NB ) + CALL XLAENV( 2, NB ) * DO 60 ITRAN = 1, 2 IF( ITRAN.EQ.1 ) THEN @@ -484,7 +484,7 @@ $ COPYB, LDB, B, LDB ) END IF SRNAMT = 'SGETSLS ' - CALL SGETSLS( TRANS, M, N, NRHS, A, + CALL SGETSLS( TRANS, M, N, NRHS, A, $ LDA, B, LDB, WORK, LWORK, INFO ) IF( INFO.NE.0 ) $ CALL ALAERH( PATH, 'SGETSLS ', INFO, 0, @@ -765,8 +765,8 @@ $ ', NB=', I4, ', type', I2, ', test(', I2, ')=', G12.5 ) 9998 FORMAT( ' M=', I5, ', N=', I5, ', NRHS=', I4, ', NB=', I4, $ ', type', I2, ', test(', I2, ')=', G12.5 ) - 9997 FORMAT( ' TRANS=''', A1,' M=', I5, ', N=', I5, ', NRHS=', I4, - $ ', MB=', I4,', NB=', I4,', type', I2, + 9997 FORMAT( ' TRANS=''', A1,' M=', I5, ', N=', I5, ', NRHS=', I4, + $ ', MB=', I4,', NB=', I4,', type', I2, $ ', test(', I2, ')=', G12.5 ) RETURN * diff --git a/TESTING/LIN/sdrvsy_aasen.f b/TESTING/LIN/sdrvsy_aasen.f index afa5dc49..6de06fda 100644 --- a/TESTING/LIN/sdrvsy_aasen.f +++ b/TESTING/LIN/sdrvsy_aasen.f @@ -455,8 +455,8 @@ c END IF * Check error code from SSYSV_AASEN . * IF( INFO.NE.K ) THEN - CALL ALAERH( PATH, 'SSYSV_AASEN ', INFO, K, - $ UPLO, N, N, -1, -1, NRHS, + CALL ALAERH( PATH, 'SSYSV_AASEN ', INFO, K, + $ UPLO, N, N, -1, -1, NRHS, $ IMAT, NFAIL, NERRS, NOUT ) GO TO 120 ELSE IF( INFO.NE.0 ) THEN @@ -467,7 +467,7 @@ c END IF * residual. * CALL SSYT01_AASEN( UPLO, N, A, LDA, AFAC, LDA, - $ IWORK, AINV, LDA, RWORK, + $ IWORK, AINV, LDA, RWORK, $ RESULT( 1 ) ) * * Compute residual of the computed solution. diff --git a/TESTING/LIN/serrlqt.f b/TESTING/LIN/serrlqt.f index 2c2c575b..d73f990d 100644 --- a/TESTING/LIN/serrlqt.f +++ b/TESTING/LIN/serrlqt.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE SERRLQT( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * @@ -81,7 +81,7 @@ * .. * .. External Subroutines .. EXTERNAL ALAESM, CHKXER, SGELQT3, SGELQT, - $ SGEMLQT + $ SGEMLQT * .. * .. Scalars in Common .. LOGICAL LERR, OK diff --git a/TESTING/LIN/serrlqtp.f b/TESTING/LIN/serrlqtp.f index 319ee91c..782de113 100644 --- a/TESTING/LIN/serrlqtp.f +++ b/TESTING/LIN/serrlqtp.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE SERRLQTP( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * @@ -81,7 +81,7 @@ * .. * .. External Subroutines .. EXTERNAL ALAESM, CHKXER, STPLQT2, STPLQT, - $ STPMLQT + $ STPMLQT * .. * .. Scalars in Common .. LOGICAL LERR, OK @@ -171,46 +171,46 @@ * SRNAMT = 'STPMLQT' INFOT = 1 - CALL STPMLQT( '/', 'N', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL STPMLQT( '/', 'N', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 2 - CALL STPMLQT( 'L', '/', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL STPMLQT( 'L', '/', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 3 - CALL STPMLQT( 'L', 'N', -1, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL STPMLQT( 'L', 'N', -1, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 4 - CALL STPMLQT( 'L', 'N', 0, -1, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL STPMLQT( 'L', 'N', 0, -1, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 5 - CALL STPMLQT( 'L', 'N', 0, 0, -1, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL STPMLQT( 'L', 'N', 0, 0, -1, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) INFOT = 6 - CALL STPMLQT( 'L', 'N', 0, 0, 0, -1, 1, A, 1, T, 1, B, 1, C, 1, + CALL STPMLQT( 'L', 'N', 0, 0, 0, -1, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 7 - CALL STPMLQT( 'L', 'N', 0, 0, 0, 0, 0, A, 1, T, 1, B, 1, C, 1, + CALL STPMLQT( 'L', 'N', 0, 0, 0, 0, 0, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 9 - CALL STPMLQT( 'R', 'N', 2, 2, 2, 1, 1, A, 1, T, 1, B, 1, C, 1, + CALL STPMLQT( 'R', 'N', 2, 2, 2, 1, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 11 - CALL STPMLQT( 'R', 'N', 1, 1, 1, 1, 1, A, 1, T, 0, B, 1, C, 1, + CALL STPMLQT( 'R', 'N', 1, 1, 1, 1, 1, A, 1, T, 0, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 13 - CALL STPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 0, C, 1, + CALL STPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 0, C, 1, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 15 - CALL STPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 1, C, 0, + CALL STPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 1, C, 0, $ W, INFO ) CALL CHKXER( 'STPMLQT', INFOT, NOUT, LERR, OK ) * diff --git a/TESTING/LIN/serrtsqr.f b/TESTING/LIN/serrtsqr.f index 0ba37978..eddadbee 100644 --- a/TESTING/LIN/serrtsqr.f +++ b/TESTING/LIN/serrtsqr.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE SERRTSQR( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * diff --git a/TESTING/LIN/slqt04.f b/TESTING/LIN/slqt04.f index debae5ca..2f4637b9 100644 --- a/TESTING/LIN/slqt04.f +++ b/TESTING/LIN/slqt04.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE SLQT04(M,N,NB,RESULT) -* +* * .. Scalar Arguments .. * INTEGER M, N, NB, LDT * .. Return values .. * REAL RESULT(6) -* +* * *> \par Purpose: * ============= @@ -54,17 +54,17 @@ *> RESULT(2) = | I - Q Q^H | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -87,9 +87,9 @@ * ===================================================================== * * .. -* .. Local allocatable arrays +* .. Local allocatable arrays REAL, ALLOCATABLE :: AF(:,:), Q(:,:), - $ L(:,:), RWORK(:), WORK( : ), T(:,:), + $ L(:,:), RWORK(:), WORK( : ), T(:,:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:) * * .. Parameters .. @@ -109,11 +109,11 @@ EXTERNAL SLAMCH, SLANGE, SLANSY, LSAME * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. * .. Data statements .. - DATA ISEED / 1988, 1989, 1990, 1991 / -* + DATA ISEED / 1988, 1989, 1990, 1991 / +* EPS = SLAMCH( 'Epsilon' ) K = MIN(M,N) LL = MAX(M,N) @@ -121,8 +121,8 @@ * * Dynamically allocate local arrays * - ALLOCATE ( A(M,N), AF(M,N), Q(N,N), L(LL,N), RWORK(LL), - $ WORK(LWORK), T(NB,N), C(M,N), CF(M,N), + ALLOCATE ( A(M,N), AF(M,N), Q(N,N), L(LL,N), RWORK(LL), + $ WORK(LWORK), T(NB,N), C(M,N), CF(M,N), $ D(N,M), DF(N,M) ) * * Put random numbers into A and copy to AF @@ -140,7 +140,7 @@ * Generate the n-by-n matrix Q * CALL SLASET( 'Full', N, N, ZERO, ONE, Q, N ) - CALL SGEMLQT( 'R', 'N', N, N, K, NB, AF, M, T, LDT, Q, N, + CALL SGEMLQT( 'R', 'N', N, N, K, NB, AF, M, T, LDT, Q, N, $ WORK, INFO ) * * Copy R @@ -176,7 +176,7 @@ * * Apply Q to C as Q*C * - CALL SGEMLQT( 'L', 'N', N, M, K, NB, AF, M, T, NB, DF, N, + CALL SGEMLQT( 'L', 'N', N, M, K, NB, AF, M, T, NB, DF, N, $ WORK, INFO) * * Compute |Q*D - Q*D| / |D| @@ -195,7 +195,7 @@ * * Apply Q to D as QT*D * - CALL SGEMLQT( 'L', 'T', N, M, K, NB, AF, M, T, NB, DF, N, + CALL SGEMLQT( 'L', 'T', N, M, K, NB, AF, M, T, NB, DF, N, $ WORK, INFO) * * Compute |QT*D - QT*D| / |D| @@ -206,7 +206,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,M)*DNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -218,8 +218,8 @@ * * Apply Q to C as C*Q * - CALL SGEMLQT( 'R', 'N', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, INFO) + CALL SGEMLQT( 'R', 'N', M, N, K, NB, AF, M, T, NB, CF, M, + $ WORK, INFO) * * Compute |C*Q - C*Q| / |C| * @@ -237,8 +237,8 @@ * * Apply Q to D as D*QT * - CALL SGEMLQT( 'R', 'T', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, INFO) + CALL SGEMLQT( 'R', 'T', M, N, K, NB, AF, M, T, NB, CF, M, + $ WORK, INFO) * * Compute |C*QT - C*QT| / |C| * diff --git a/TESTING/LIN/slqt05.f b/TESTING/LIN/slqt05.f index 5ad3a4b2..7ce36c75 100644 --- a/TESTING/LIN/slqt05.f +++ b/TESTING/LIN/slqt05.f @@ -2,12 +2,12 @@ * =========== * * SUBROUTINE SLQT05(M,N,L,NB,RESULT) -* +* * .. Scalar Arguments .. * INTEGER LWORK, M, N, L, NB, LDT * .. Return values .. * REAL RESULT(6) -* +* * *> \par Purpose: * ============= @@ -38,7 +38,7 @@ *> The number of rows of the upper trapezoidal part the *> lower test matrix. 0 <= L <= M. *> \endverbatim -*> +*> *> \param[in] NB *> \verbatim *> NB is INTEGER @@ -54,17 +54,17 @@ *> RESULT(2) = | I - Q^H Q | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -85,11 +85,11 @@ REAL RESULT(6) * * ===================================================================== -* +* * .. -* .. Local allocatable arrays +* .. Local allocatable arrays REAL, ALLOCATABLE :: AF(:,:), Q(:,:), - $ R(:,:), RWORK(:), WORK( : ), T(:,:), + $ R(:,:), RWORK(:), WORK( : ), T(:,:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:) * * .. Parameters .. @@ -110,7 +110,7 @@ * .. * .. Data statements .. DATA ISEED / 1988, 1989, 1990, 1991 / -* +* EPS = SLAMCH( 'Epsilon' ) K = M N2 = M+N @@ -124,7 +124,7 @@ * Dynamically allocate all arrays * ALLOCATE(A(M,N2),AF(M,N2),Q(N2,N2),R(N2,N2),RWORK(N2), - $ WORK(LWORK),T(NB,M),C(N2,M),CF(N2,M), + $ WORK(LWORK),T(NB,M),C(N2,M),CF(N2,M), $ D(M,N2),DF(M,N2) ) * * Put random stuff into A @@ -142,7 +142,7 @@ END IF IF( L.GT.0 ) THEN DO J=1,L - CALL SLARNV( 2, ISEED, M-J+1, A( J, MIN(N+M,N+M-L+1) + CALL SLARNV( 2, ISEED, M-J+1, A( J, MIN(N+M,N+M-L+1) $ + J - 1 ) ) END DO END IF @@ -194,7 +194,7 @@ CALL SLACPY( 'Full', N2, M, C, N2, CF, N2 ) * * Apply Q to C as Q*C -* +* CALL STPMLQT( 'L','N', N,M,K,L,NB,AF(1, NP1),M,T,LDT,CF,N2, $ CF(NP1,1),N2,WORK,INFO) * @@ -216,18 +216,18 @@ * Apply Q to C as QT*C * CALL STPMLQT( 'L','T',N,M,K,L,NB,AF(1,NP1),M,T,LDT,CF,N2, - $ CF(NP1,1),N2,WORK,INFO) + $ CF(NP1,1),N2,WORK,INFO) * * Compute |QT*C - QT*C| / |C| * CALL SGEMM('T','N',N2,M,N2,-ONE,Q,N2,C,N2,ONE,CF,N2) RESID = SLANGE( '1', N2, M, CF, N2, RWORK ) - + IF( CNORM.GT.ZERO ) THEN RESULT( 4 ) = RESID / (EPS*MAX(1,N2)*CNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random m-by-n matrix D and a copy DF * @@ -259,8 +259,8 @@ * Apply Q to D as D*QT * CALL STPMLQT('R','T',M,N,K,L,NB,AF(1,NP1),M,T,LDT,DF,M, - $ DF(1,NP1),M,WORK,INFO) - + $ DF(1,NP1),M,WORK,INFO) + * * Compute |D*QT - D*QT| / |D| * @@ -276,4 +276,4 @@ * DEALLOCATE ( A, AF, Q, R, RWORK, WORK, T, C, D, CF, DF) RETURN - END
\ No newline at end of file + END diff --git a/TESTING/LIN/stplqt.f b/TESTING/LIN/stplqt.f index adbbfe8b..4e03ae65 100644 --- a/TESTING/LIN/stplqt.f +++ b/TESTING/LIN/stplqt.f @@ -3,23 +3,23 @@ * * SUBROUTINE STPLQT( M, N, L, MB, A, LDA, B, LDB, T, LDT, WORK, * INFO ) -* +* * .. Scalar Arguments .. * INTEGER INFO, LDA, LDB, LDT, N, M, L, MB * .. * .. Array Arguments .. * REAL A( LDA, * ), B( LDB, * ), T( LDT, * ), WORK( * ) * .. -* +* * *> \par Purpose: * ============= *> *> \verbatim *> -*> STPLQT computes a blocked LQ factorization of a real -*> "triangular-pentagonal" matrix C, which is composed of a -*> triangular block A and pentagonal block B, using the compact +*> STPLQT computes a blocked LQ factorization of a real +*> "triangular-pentagonal" matrix C, which is composed of a +*> triangular block A and pentagonal block B, using the compact *> WY representation for Q. *> \endverbatim * @@ -30,7 +30,7 @@ *> \verbatim *> M is INTEGER *> The number of rows of the matrix B, and the order of the -*> triangular matrix A. +*> triangular matrix A. *> M >= 0. *> \endverbatim *> @@ -71,7 +71,7 @@ *> \param[in,out] B *> \verbatim *> B is REAL array, dimension (LDB,N) -*> On entry, the pentagonal M-by-N matrix B. The first N-L columns +*> On entry, the pentagonal M-by-N matrix B. The first N-L columns *> are rectangular, and the last L columns are lower trapezoidal. *> On exit, B contains the pentagonal matrix V. See Further Details. *> \endverbatim @@ -88,7 +88,7 @@ *> The lower triangular block reflectors stored in compact form *> as a sequence of upper triangular blocks. See Further Details. *> \endverbatim -*> +*> *> \param[in] LDT *> \verbatim *> LDT is INTEGER @@ -110,10 +110,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2013 * @@ -124,45 +124,45 @@ *> *> \verbatim *> -*> The input matrix C is a M-by-(M+N) matrix +*> The input matrix C is a M-by-(M+N) matrix *> *> C = [ A ] [ B ] -*> +*> *> *> where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal *> matrix consisting of a M-by-(N-L) rectangular matrix B1 on left of a M-by-L *> upper trapezoidal matrix B2: -*> [ B ] = [ B1 ] [ B2 ] +*> [ B ] = [ B1 ] [ B2 ] *> [ B1 ] <- M-by-(N-L) rectangular *> [ B2 ] <- M-by-L upper trapezoidal. *> *> The lower trapezoidal matrix B2 consists of the first L columns of a -*> N-by-N lower triangular matrix, where 0 <= L <= MIN(M,N). If L=0, -*> B is rectangular M-by-N; if M=L=N, B is lower triangular. +*> N-by-N lower triangular matrix, where 0 <= L <= MIN(M,N). If L=0, +*> B is rectangular M-by-N; if M=L=N, B is lower triangular. *> *> The matrix W stores the elementary reflectors H(i) in the i-th row *> above the diagonal (of A) in the M-by-(M+N) input matrix C -*> [ C ] = [ A ] [ B ] +*> [ C ] = [ A ] [ B ] *> [ A ] <- lower triangular N-by-N *> [ B ] <- M-by-N pentagonal *> *> so that W can be represented as -*> [ W ] = [ I ] [ V ] +*> [ W ] = [ I ] [ V ] *> [ I ] <- identity, N-by-N *> [ V ] <- M-by-N, same form as B. *> *> Thus, all of information needed for W is contained on exit in B, which -*> we call V above. Note that V has the same form as B; that is, -*> [ V ] = [ V1 ] [ V2 ] +*> we call V above. Note that V has the same form as B; that is, +*> [ V ] = [ V1 ] [ V2 ] *> [ V1 ] <- M-by-(N-L) rectangular *> [ V2 ] <- M-by-L lower trapezoidal. *> -*> The rows of V represent the vectors which define the H(i)'s. +*> The rows of V represent the vectors which define the H(i)'s. *> *> The number of blocks is B = ceiling(M/MB), where each -*> block is of order MB except for the last block, which is of order +*> block is of order MB except for the last block, which is of order *> IB = M - (M-1)*MB. For each of the B blocks, a upper triangular block -*> reflector factor is computed: T1, T2, ..., TB. The MB-by-MB (and IB-by-IB +*> reflector factor is computed: T1, T2, ..., TB. The MB-by-MB (and IB-by-IB *> for the last block) T's are stored in the MB-by-N matrix T as *> *> T = [T1 T2 ... TB]. @@ -223,7 +223,7 @@ IF( M.EQ.0 .OR. N.EQ.0 ) RETURN * DO I = 1, M, MB -* +* * Compute the QR factorization of the current block * IB = MIN( M-I+1, MB ) @@ -234,20 +234,20 @@ LB = NB-N+L-I+1 END IF * - CALL STPLQT2( IB, NB, LB, A(I,I), LDA, B( I, 1 ), LDB, + CALL STPLQT2( IB, NB, LB, A(I,I), LDA, B( I, 1 ), LDB, $ T(1, I ), LDT, IINFO ) * * Update by applying H**T to B(I+IB:M,:) from the right * IF( I+IB.LE.M ) THEN CALL STPRFB( 'R', 'N', 'F', 'R', M-I-IB+1, NB, IB, LB, - $ B( I, 1 ), LDB, T( 1, I ), LDT, - $ A( I+IB, I ), LDA, B( I+IB, 1 ), LDB, + $ B( I, 1 ), LDB, T( 1, I ), LDT, + $ A( I+IB, I ), LDA, B( I+IB, 1 ), LDB, $ WORK, M-I-IB+1) END IF END DO RETURN -* +* * End of STPLQT * END diff --git a/TESTING/LIN/stsqr01.f b/TESTING/LIN/stsqr01.f index dbaf3aac..4cebfc88 100644 --- a/TESTING/LIN/stsqr01.f +++ b/TESTING/LIN/stsqr01.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE STSQR01(TSSW, M,N, MB, NB, RESULT) -* +* * .. Scalar Arguments .. * INTEGER M, N, MB * .. Return values .. * REAL RESULT(6) -* +* * *> \par Purpose: * ============= @@ -65,17 +65,17 @@ *> RESULT(2) = | I - Q^H Q | or | I - Q Q^H | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -99,9 +99,9 @@ * ===================================================================== * * .. -* .. Local allocatable arrays +* .. Local allocatable arrays REAL, ALLOCATABLE :: AF(:,:), Q(:,:), - $ R(:,:), RWORK(:), WORK( : ), T(:), + $ R(:,:), RWORK(:), WORK( : ), T(:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:), LQ(:,:) * * .. Parameters .. @@ -123,24 +123,24 @@ EXTERNAL SLAMCH, SLARNV, SLANGE, SLANSY, LSAME, ILAENV * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. Scalars in Common .. CHARACTER*32 srnamt * .. * .. Common blocks .. - COMMON / srnamc / srnamt + COMMON / srnamc / srnamt * .. * .. Data statements .. - DATA ISEED / 1988, 1989, 1990, 1991 / + DATA ISEED / 1988, 1989, 1990, 1991 / * * TEST TALL SKINNY OR SHORT WIDE * - TS = LSAME(TSSW, 'TS') -* + TS = LSAME(TSSW, 'TS') +* * TEST MATRICES WITH HALF OF MATRIX BEING ZEROS * TESTZEROS = .FALSE. -* +* EPS = SLAMCH( 'Epsilon' ) K = MIN(M,N) L = MAX(M,N,1) @@ -149,14 +149,14 @@ IF((K.GE.MNB).OR.(MNB.GE.L))THEN LT=MAX(1,L)*MNB+5 ELSE - LT=MAX(1,(L-K)/(MNB-K)+1)*L*MNB+5 + LT=MAX(1,(L-K)/(MNB-K)+1)*L*MNB+5 END IF * * Dynamically allocate local arrays * - ALLOCATE ( A(M,N), AF(M,N), Q(L,L), R(M,L), RWORK(L), - $ WORK(LWORK), T(LT), C(M,N), CF(M,N), + ALLOCATE ( A(M,N), AF(M,N), Q(L,L), R(M,L), RWORK(L), + $ WORK(LWORK), T(LT), C(M,N), CF(M,N), $ D(N,M), DF(N,M), LQ(L,N) ) * * Put random numbers into A and copy to AF @@ -184,7 +184,7 @@ * CALL SLASET( 'Full', M, M, ZERO, ONE, Q, M ) srnamt = 'SGEMQR' - CALL SGEMQR( 'L', 'N', M, M, K, AF, M, T, LT, Q, M, + CALL SGEMQR( 'L', 'N', M, M, K, AF, M, T, LT, Q, M, $ WORK, LWORK, INFO ) * * Copy R @@ -221,7 +221,7 @@ * Apply Q to C as Q*C * srnamt = 'DGEQR' - CALL SGEMQR( 'L', 'N', M, N, K, AF, M, T, LT, CF, M, + CALL SGEMQR( 'L', 'N', M, N, K, AF, M, T, LT, CF, M, $ WORK, LWORK, INFO) * * Compute |Q*C - Q*C| / |C| @@ -241,7 +241,7 @@ * Apply Q to C as QT*C * srnamt = 'DGEQR' - CALL SGEMQR( 'L', 'T', M, N, K, AF, M, T, LT, CF, M, + CALL SGEMQR( 'L', 'T', M, N, K, AF, M, T, LT, CF, M, $ WORK, LWORK, INFO) * * Compute |QT*C - QT*C| / |C| @@ -252,7 +252,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,M)*CNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -265,8 +265,8 @@ * Apply Q to D as D*Q * srnamt = 'DGEQR' - CALL SGEMQR( 'R', 'N', N, M, K, AF, M, T, LT, DF, N, - $ WORK, LWORK, INFO) + CALL SGEMQR( 'R', 'N', N, M, K, AF, M, T, LT, DF, N, + $ WORK, LWORK, INFO) * * Compute |D*Q - D*Q| / |D| * @@ -284,8 +284,8 @@ * * Apply Q to D as D*QT * - CALL SGEMQR( 'R', 'T', N, M, K, AF, M, T, LT, DF, N, - $ WORK, LWORK, INFO) + CALL SGEMQR( 'R', 'T', N, M, K, AF, M, T, LT, DF, N, + $ WORK, LWORK, INFO) * * Compute |D*QT - D*QT| / |D| * @@ -308,7 +308,7 @@ * CALL SLASET( 'Full', N, N, ZERO, ONE, Q, N ) srnamt = 'SGEMQR' - CALL SGEMLQ( 'R', 'N', N, N, K, AF, M, T, LT, Q, N, + CALL SGEMLQ( 'R', 'N', N, N, K, AF, M, T, LT, Q, N, $ WORK, LWORK, INFO ) * * Copy R @@ -344,7 +344,7 @@ * * Apply Q to C as Q*C * - CALL SGEMLQ( 'L', 'N', N, M, K, AF, M, T, LT, DF, N, + CALL SGEMLQ( 'L', 'N', N, M, K, AF, M, T, LT, DF, N, $ WORK, LWORK, INFO) * * Compute |Q*D - Q*D| / |D| @@ -363,7 +363,7 @@ * * Apply Q to D as QT*D * - CALL SGEMLQ( 'L', 'T', N, M, K, AF, M, T, LT, DF, N, + CALL SGEMLQ( 'L', 'T', N, M, K, AF, M, T, LT, DF, N, $ WORK, LWORK, INFO) * * Compute |QT*D - QT*D| / |D| @@ -374,7 +374,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,N)*DNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -386,8 +386,8 @@ * * Apply Q to C as C*Q * - CALL SGEMLQ( 'R', 'N', M, N, K, AF, M, T, LT, CF, M, - $ WORK, LWORK, INFO) + CALL SGEMLQ( 'R', 'N', M, N, K, AF, M, T, LT, CF, M, + $ WORK, LWORK, INFO) * * Compute |C*Q - C*Q| / |C| * @@ -405,8 +405,8 @@ * * Apply Q to D as D*QT * - CALL SGEMLQ( 'R', 'T', M, N, K, AF, M, T, LT, CF, M, - $ WORK, LWORK, INFO) + CALL SGEMLQ( 'R', 'T', M, N, K, AF, M, T, LT, CF, M, + $ WORK, LWORK, INFO) * * Compute |C*QT - C*QT| / |C| * @@ -425,4 +425,4 @@ DEALLOCATE ( A, AF, Q, R, RWORK, WORK, T, C, D, CF, DF) * RETURN - END
\ No newline at end of file + END diff --git a/TESTING/LIN/zchkaa.f b/TESTING/LIN/zchkaa.f index b5830619..97872490 100644 --- a/TESTING/LIN/zchkaa.f +++ b/TESTING/LIN/zchkaa.f @@ -651,7 +651,7 @@ * IF( TSTCHK ) THEN CALL ZCHKHE_AASEN( DOTYPE, NN, NVAL, NNB2, NBVAL2, NNS, - $ NSVAL, THRESH, TSTERR, LDA, + $ NSVAL, THRESH, TSTERR, LDA, $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), $ WORK, RWORK, IWORK, NOUT ) @@ -661,8 +661,8 @@ * IF( TSTDRV ) THEN CALL ZDRVHE_AASEN( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, - $ LDA, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), - $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), + $ LDA, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), + $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), $ WORK, RWORK, IWORK, NOUT ) ELSE WRITE( NOUT, FMT = 9988 )PATH @@ -1042,7 +1042,7 @@ * XQ: LQT routines for triangular-pentagonal matrices * IF( TSTCHK ) THEN - CALL ZCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL ZCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH @@ -1053,7 +1053,7 @@ * TS: QR routines for tall-skinny matrices * IF( TSTCHK ) THEN - CALL ZCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + CALL ZCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) ELSE WRITE( NOUT, FMT = 9989 )PATH diff --git a/TESTING/LIN/zchkhe_aasen.f b/TESTING/LIN/zchkhe_aasen.f index 1ebe5ef4..64423d89 100644 --- a/TESTING/LIN/zchkhe_aasen.f +++ b/TESTING/LIN/zchkhe_aasen.f @@ -225,7 +225,7 @@ * .. External Subroutines .. EXTERNAL ALAERH, ALAHD, ALASUM, XLAENV, ZERRHE, ZGET04, $ ZHECON, ZHERFS, ZHET01, ZHETRF_AASEN, ZHETRI2, - $ ZHETRS_AASEN, ZLACPY, ZLAIPD, ZLARHS, ZLATB4, + $ ZHETRS_AASEN, ZLACPY, ZLAIPD, ZLARHS, ZLATB4, $ ZLATMS, ZPOT02, ZPOT03, ZPOT05 * .. * .. Intrinsic Functions .. @@ -432,7 +432,7 @@ * LWORK = ( NB+1 )*LDA SRNAMT = 'ZHETRF_AASEN' - CALL ZHETRF_AASEN( UPLO, N, AFAC, LDA, IWORK, AINV, + CALL ZHETRF_AASEN( UPLO, N, AFAC, LDA, IWORK, AINV, $ LWORK, INFO ) * * Adjust the expected value of INFO to account for @@ -458,8 +458,8 @@ * Check error code from ZHETRF and handle error. * IF( INFO.NE.K ) THEN - CALL ALAERH( PATH, 'ZHETRF_AASEN', INFO, K, UPLO, - $ N, N, -1, -1, NB, IMAT, NFAIL, NERRS, + CALL ALAERH( PATH, 'ZHETRF_AASEN', INFO, K, UPLO, + $ N, N, -1, -1, NB, IMAT, NFAIL, NERRS, $ NOUT ) END IF * @@ -513,7 +513,7 @@ * SRNAMT = 'ZLARHS' CALL ZLARHS( MATPATH, XTYPE, UPLO, ' ', N, N, - $ KL, KU, NRHS, A, LDA, XACT, LDA, + $ KL, KU, NRHS, A, LDA, XACT, LDA, $ B, LDA, ISEED, INFO ) CALL ZLACPY( 'Full', N, NRHS, B, LDA, X, LDA ) * diff --git a/TESTING/LIN/zchklqt.f b/TESTING/LIN/zchklqt.f index e15793be..31c885bc 100644 --- a/TESTING/LIN/zchklqt.f +++ b/TESTING/LIN/zchklqt.f @@ -2,15 +2,15 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE ZCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE ZCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) -* +* * .. Scalar Arguments .. * LOGICAL TSTERR * INTEGER NM, NN, NNB, NOUT @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE ZCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE ZCHKLQT( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -175,7 +175,7 @@ NB = NBVAL( K ) * * Test ZGELQT and ZUNMLQT -* +* IF( (NB.LE.MINMN).AND.(NB.GT.0) ) THEN CALL ZLQT04( M, N, NB, RESULT ) * diff --git a/TESTING/LIN/zchklqtp.f b/TESTING/LIN/zchklqtp.f index 10f7363a..efde54bd 100644 --- a/TESTING/LIN/zchklqtp.f +++ b/TESTING/LIN/zchklqtp.f @@ -2,13 +2,13 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE ZCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE ZCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) * * .. Scalar Arguments .. @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE ZCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE ZCHKLQTP( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -172,14 +172,14 @@ * MINMN = MIN( M, N ) DO L = 0, MINMN, MAX( MINMN, 1 ) -* +* * Do for each possible value of NB * DO K = 1, NNB NB = NBVAL( K ) * * Test DTPLQT and DTPMLQT -* +* IF( (NB.LE.M).AND.(NB.GT.0) ) THEN CALL ZLQT05( M, N, L, NB, RESULT ) * @@ -212,4 +212,4 @@ * * End of ZCHKLQTP * - END
\ No newline at end of file + END diff --git a/TESTING/LIN/zchktsqr.f b/TESTING/LIN/zchktsqr.f index c79a92b7..6d8fc99b 100644 --- a/TESTING/LIN/zchktsqr.f +++ b/TESTING/LIN/zchktsqr.f @@ -2,15 +2,15 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * -* SUBROUTINE ZCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, +* SUBROUTINE ZCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, * NBVAL, NOUT ) -* +* * .. Scalar Arguments .. * LOGICAL TSTERR * INTEGER NM, NN, NNB, NOUT @@ -89,17 +89,17 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * *> \ingroup double_lin * * ===================================================================== - SUBROUTINE ZCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, + SUBROUTINE ZCHKTSQR( THRESH, TSTERR, NM, MVAL, NN, NVAL, NNB, $ NBVAL, NOUT ) IMPLICIT NONE * @@ -132,11 +132,11 @@ DOUBLE PRECISION RESULT( NTESTS ) * .. * .. External Subroutines .. - EXTERNAL ALAERH, ALAHD, ALASUM, DERRTSQR, + EXTERNAL ALAERH, ALAHD, ALASUM, DERRTSQR, $ DTSQR01, XLAENV * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. * .. Scalars in Common .. LOGICAL LERR, OK @@ -172,7 +172,7 @@ DO J = 1, NN N = NVAL( J ) IF (MIN(M,N).NE.0) THEN - DO INB = 1, NNB + DO INB = 1, NNB MB = NBVAL( INB ) CALL XLAENV( 1, MB ) DO IMB = 1, NNB @@ -180,7 +180,7 @@ CALL XLAENV( 2, NB ) * * Test ZGEQR and ZGEMQR -* +* CALL ZTSQR01( 'TS', M, N, MB, NB, RESULT ) * * Print information about the tests that did not @@ -196,9 +196,9 @@ END IF END DO NRUN = NRUN + NTESTS - END DO - END DO - END IF + END DO + END DO + END IF END DO END DO * @@ -212,7 +212,7 @@ DO J = 1, NN N = NVAL( J ) IF (MIN(M,N).NE.0) THEN - DO INB = 1, NNB + DO INB = 1, NNB MB = NBVAL( INB ) CALL XLAENV( 1, MB ) DO IMB = 1, NNB @@ -220,7 +220,7 @@ CALL XLAENV( 2, NB ) * * Test ZGELQ and ZGEMLQ -* +* CALL ZTSQR01( 'SW', M, N, MB, NB, RESULT ) * * Print information about the tests that did not @@ -236,9 +236,9 @@ END IF END DO NRUN = NRUN + NTESTS - END DO - END DO - END IF + END DO + END DO + END IF END DO END DO * diff --git a/TESTING/LIN/zdrvhe_aasen.f b/TESTING/LIN/zdrvhe_aasen.f index 0170873e..4e283564 100644 --- a/TESTING/LIN/zdrvhe_aasen.f +++ b/TESTING/LIN/zdrvhe_aasen.f @@ -150,7 +150,7 @@ * * ===================================================================== SUBROUTINE ZDRVHE_AASEN( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, - $ NMAX, A, AFAC, AINV, B, X, XACT, WORK, + $ NMAX, A, AFAC, AINV, B, X, XACT, WORK, $ RWORK, IWORK, NOUT ) * * -- LAPACK test routine (version 3.7.0) -- @@ -201,7 +201,7 @@ * .. * .. External Subroutines .. EXTERNAL ALADHD, ALAERH, ALASVM, XLAENV, ZERRVX, ZGET04, - $ ZHESV_AASEN, ZHET01_AASEN, ZHETRF_AASEN, + $ ZHESV_AASEN, ZHET01_AASEN, ZHETRF_AASEN, $ ZHETRI2, ZLACPY, ZLAIPD, ZLARHS, ZLATB4, ZLATMS, $ ZPOT02 * .. @@ -475,7 +475,7 @@ c END IF * residual. * CALL ZHET01_AASEN( UPLO, N, A, LDA, AFAC, LDA, - $ IWORK, AINV, LDA, RWORK, + $ IWORK, AINV, LDA, RWORK, $ RESULT( 1 ) ) * * Compute residual of the computed solution. diff --git a/TESTING/LIN/zdrvls.f b/TESTING/LIN/zdrvls.f index 1ffa585f..63fcc69e 100644 --- a/TESTING/LIN/zdrvls.f +++ b/TESTING/LIN/zdrvls.f @@ -488,7 +488,7 @@ $ COPYB, LDB, B, LDB ) END IF SRNAMT = 'DGETSLS ' - CALL ZGETSLS( TRANS, M, N, NRHS, A, + CALL ZGETSLS( TRANS, M, N, NRHS, A, $ LDA, B, LDB, WORK, LWORK, INFO ) IF( INFO.NE.0 ) $ CALL ALAERH( PATH, 'ZGETSLS ', INFO, 0, @@ -771,8 +771,8 @@ $ ', NB=', I4, ', type', I2, ', test(', I2, ')=', G12.5 ) 9998 FORMAT( ' M=', I5, ', N=', I5, ', NRHS=', I4, ', NB=', I4, $ ', type', I2, ', test(', I2, ')=', G12.5 ) - 9997 FORMAT( ' TRANS=''', A1,' M=', I5, ', N=', I5, ', NRHS=', I4, - $ ', MB=', I4,', NB=', I4,', type', I2, + 9997 FORMAT( ' TRANS=''', A1,' M=', I5, ', N=', I5, ', NRHS=', I4, + $ ', MB=', I4,', NB=', I4,', type', I2, $ ', test(', I2, ')=', G12.5 ) RETURN * diff --git a/TESTING/LIN/zerrlqt.f b/TESTING/LIN/zerrlqt.f index fd6b4527..f65f8951 100644 --- a/TESTING/LIN/zerrlqt.f +++ b/TESTING/LIN/zerrlqt.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE ZERRLQT( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * @@ -81,7 +81,7 @@ * .. * .. External Subroutines .. EXTERNAL ALAESM, CHKXER, ZGELQT3, ZGELQT, - $ ZGEMLQT + $ ZGEMLQT * .. * .. Scalars in Common .. LOGICAL LERR, OK diff --git a/TESTING/LIN/zerrlqtp.f b/TESTING/LIN/zerrlqtp.f index 25a079ec..6e00e4f9 100644 --- a/TESTING/LIN/zerrlqtp.f +++ b/TESTING/LIN/zerrlqtp.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE ZERRLQTP( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2011 * @@ -81,7 +81,7 @@ * .. * .. External Subroutines .. EXTERNAL ALAESM, CHKXER, ZTPLQT2, ZTPLQT, - $ ZTPMLQT + $ ZTPMLQT * .. * .. Scalars in Common .. LOGICAL LERR, OK @@ -171,46 +171,46 @@ * SRNAMT = 'ZTPMLQT' INFOT = 1 - CALL ZTPMLQT( '/', 'N', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL ZTPMLQT( '/', 'N', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 2 - CALL ZTPMLQT( 'L', '/', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL ZTPMLQT( 'L', '/', 0, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 3 - CALL ZTPMLQT( 'L', 'N', -1, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL ZTPMLQT( 'L', 'N', -1, 0, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 4 - CALL ZTPMLQT( 'L', 'N', 0, -1, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL ZTPMLQT( 'L', 'N', 0, -1, 0, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 5 - CALL ZTPMLQT( 'L', 'N', 0, 0, -1, 0, 1, A, 1, T, 1, B, 1, C, 1, + CALL ZTPMLQT( 'L', 'N', 0, 0, -1, 0, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) INFOT = 6 - CALL ZTPMLQT( 'L', 'N', 0, 0, 0, -1, 1, A, 1, T, 1, B, 1, C, 1, + CALL ZTPMLQT( 'L', 'N', 0, 0, 0, -1, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 7 - CALL ZTPMLQT( 'L', 'N', 0, 0, 0, 0, 0, A, 1, T, 1, B, 1, C, 1, + CALL ZTPMLQT( 'L', 'N', 0, 0, 0, 0, 0, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 9 - CALL ZTPMLQT( 'R', 'N', 2, 2, 2, 1, 1, A, 1, T, 1, B, 1, C, 1, + CALL ZTPMLQT( 'R', 'N', 2, 2, 2, 1, 1, A, 1, T, 1, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 11 - CALL ZTPMLQT( 'R', 'N', 1, 1, 1, 1, 1, A, 1, T, 0, B, 1, C, 1, + CALL ZTPMLQT( 'R', 'N', 1, 1, 1, 1, 1, A, 1, T, 0, B, 1, C, 1, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 13 - CALL ZTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 0, C, 1, + CALL ZTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 0, C, 1, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) INFOT = 15 - CALL ZTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 1, C, 0, + CALL ZTPMLQT( 'L', 'N', 1, 1, 1, 1, 1, A, 1, T, 1, B, 1, C, 0, $ W, INFO ) CALL CHKXER( 'ZTPMLQT', INFOT, NOUT, LERR, OK ) * diff --git a/TESTING/LIN/zerrtsqr.f b/TESTING/LIN/zerrtsqr.f index 19c99805..3aa3e4a5 100644 --- a/TESTING/LIN/zerrtsqr.f +++ b/TESTING/LIN/zerrtsqr.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE ZERRTSQR( PATH, NUNIT ) -* +* * .. Scalar Arguments .. * CHARACTER*3 PATH * INTEGER NUNIT * .. -* +* * *> \par Purpose: * ============= @@ -43,10 +43,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Zenver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Zenver +*> \author NAG Ltd. * *> \date November 2011 * diff --git a/TESTING/LIN/zhet01_aasen.f b/TESTING/LIN/zhet01_aasen.f index 89b87990..f658029b 100644 --- a/TESTING/LIN/zhet01_aasen.f +++ b/TESTING/LIN/zhet01_aasen.f @@ -145,7 +145,7 @@ * * .. Parameters .. COMPLEX*16 CZERO, CONE - PARAMETER ( CZERO = ( 0.0D+0, 0.0D+0 ), + PARAMETER ( CZERO = ( 0.0D+0, 0.0D+0 ), $ CONE = ( 1.0D+0, 0.0D+0 ) ) DOUBLE PRECISION ZERO, ONE PARAMETER ( ZERO = 0.0D+0, ONE = 1.0D+0 ) diff --git a/TESTING/LIN/zlqt04.f b/TESTING/LIN/zlqt04.f index a1aff90e..4571215e 100644 --- a/TESTING/LIN/zlqt04.f +++ b/TESTING/LIN/zlqt04.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE ZLQT04(M,N,NB,RESULT) -* +* * .. Scalar Arguments .. * INTEGER M, N, NB * .. Return values .. * DOUBLE PRECISION RESULT(6) -* +* * *> \par Purpose: * ============= @@ -54,17 +54,17 @@ *> RESULT(2) = | I - Q Q^H | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -87,9 +87,9 @@ * ===================================================================== * * .. -* .. Local allocatable arrays +* .. Local allocatable arrays COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:), - $ L(:,:), RWORK(:), WORK( : ), T(:,:), + $ L(:,:), RWORK(:), WORK( : ), T(:,:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:) * * .. Parameters .. @@ -112,11 +112,11 @@ EXTERNAL DLAMCH, ZLANGE, ZLANSY, LSAME * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. * .. Data statements .. - DATA ISEED / 1988, 1989, 1990, 1991 / -* + DATA ISEED / 1988, 1989, 1990, 1991 / +* EPS = DLAMCH( 'Epsilon' ) K = MIN(M,N) LL = MAX(M,N) @@ -124,8 +124,8 @@ * * Dynamically allocate local arrays * - ALLOCATE ( A(M,N), AF(M,N), Q(N,N), L(LL,N), RWORK(LL), - $ WORK(LWORK), T(NB,N), C(M,N), CF(M,N), + ALLOCATE ( A(M,N), AF(M,N), Q(N,N), L(LL,N), RWORK(LL), + $ WORK(LWORK), T(NB,N), C(M,N), CF(M,N), $ D(N,M), DF(N,M) ) * * Put random numbers into A and copy to AF @@ -143,7 +143,7 @@ * Generate the n-by-n matrix Q * CALL ZLASET( 'Full', N, N, CZERO, ONE, Q, N ) - CALL ZGEMLQT( 'R', 'N', N, N, K, NB, AF, M, T, LDT, Q, N, + CALL ZGEMLQT( 'R', 'N', N, N, K, NB, AF, M, T, LDT, Q, N, $ WORK, INFO ) * * Copy L @@ -179,7 +179,7 @@ * * Apply Q to C as Q*C * - CALL ZGEMLQT( 'L', 'N', N, M, K, NB, AF, M, T, NB, DF, N, + CALL ZGEMLQT( 'L', 'N', N, M, K, NB, AF, M, T, NB, DF, N, $ WORK, INFO) * * Compute |Q*D - Q*D| / |D| @@ -198,7 +198,7 @@ * * Apply Q to D as QT*D * - CALL ZGEMLQT( 'L', 'C', N, M, K, NB, AF, M, T, NB, DF, N, + CALL ZGEMLQT( 'L', 'C', N, M, K, NB, AF, M, T, NB, DF, N, $ WORK, INFO) * * Compute |QT*D - QT*D| / |D| @@ -209,7 +209,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,M)*DNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -221,8 +221,8 @@ * * Apply Q to C as C*Q * - CALL ZGEMLQT( 'R', 'N', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, INFO) + CALL ZGEMLQT( 'R', 'N', M, N, K, NB, AF, M, T, NB, CF, M, + $ WORK, INFO) * * Compute |C*Q - C*Q| / |C| * @@ -240,8 +240,8 @@ * * Apply Q to D as D*QT * - CALL ZGEMLQT( 'R', 'C', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, INFO) + CALL ZGEMLQT( 'R', 'C', M, N, K, NB, AF, M, T, NB, CF, M, + $ WORK, INFO) * * Compute |C*QT - C*QT| / |C| * diff --git a/TESTING/LIN/zlqt05.f b/TESTING/LIN/zlqt05.f index 676c95b8..196750fd 100644 --- a/TESTING/LIN/zlqt05.f +++ b/TESTING/LIN/zlqt05.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE ZLQT05(M,N,L,NB,RESULT) -* +* * .. Scalar Arguments .. * INTEGER LWORK, M, N, L, NB, LDT * .. Return values .. * DOUBLE PRECISION RESULT(6) -* +* * *> \par Purpose: * ============= @@ -45,7 +45,7 @@ *> The number of rows of the upper trapezoidal part the *> lower test matrix. 0 <= L <= M. *> \endverbatim -*> +*> *> \param[in] NB *> \verbatim *> NB is INTEGER @@ -61,17 +61,17 @@ *> RESULT(2) = | I - Q^H Q | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -92,11 +92,11 @@ DOUBLE PRECISION RESULT(6) * * ===================================================================== -* +* * .. -* .. Local allocatable arrays +* .. Local allocatable arrays COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:), - $ R(:,:), RWORK(:), WORK( : ), T(:,:), + $ R(:,:), RWORK(:), WORK( : ), T(:,:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:) * * .. Parameters .. @@ -119,7 +119,7 @@ * .. * .. Data statements .. DATA ISEED / 1988, 1989, 1990, 1991 / -* +* EPS = DLAMCH( 'Epsilon' ) K = M N2 = M+N @@ -133,7 +133,7 @@ * Dynamically allocate all arrays * ALLOCATE(A(M,N2),AF(M,N2),Q(N2,N2),R(N2,N2),RWORK(N2), - $ WORK(LWORK),T(NB,M),C(N2,M),CF(N2,M), + $ WORK(LWORK),T(NB,M),C(N2,M),CF(N2,M), $ D(M,N2),DF(M,N2) ) * * Put random stuff into A @@ -151,7 +151,7 @@ END IF IF( L.GT.0 ) THEN DO J=1,L - CALL ZLARNV( 2, ISEED, M-J+1, A( J, MIN(N+M,N+M-L+1) + CALL ZLARNV( 2, ISEED, M-J+1, A( J, MIN(N+M,N+M-L+1) $ + J - 1 ) ) END DO END IF @@ -204,7 +204,7 @@ CALL ZLACPY( 'Full', N2, M, C, N2, CF, N2 ) * * Apply Q to C as Q*C -* +* CALL ZTPMLQT( 'L','N', N,M,K,L,NB,AF(1, NP1),M,T,LDT,CF,N2, $ CF(NP1,1),N2,WORK,INFO) * @@ -226,18 +226,18 @@ * Apply Q to C as QT*C * CALL ZTPMLQT( 'L','C',N,M,K,L,NB,AF(1,NP1),M,T,LDT,CF,N2, - $ CF(NP1,1),N2,WORK,INFO) + $ CF(NP1,1),N2,WORK,INFO) * * Compute |QT*C - QT*C| / |C| * CALL ZGEMM('C','N',N2,M,N2,-ONE,Q,N2,C,N2,ONE,CF,N2) RESID = ZLANGE( '1', N2, M, CF, N2, RWORK ) - + IF( CNORM.GT.ZERO ) THEN RESULT( 4 ) = RESID / (EPS*MAX(1,N2)*CNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random m-by-n matrix D and a copy DF * @@ -269,8 +269,8 @@ * Apply Q to D as D*QT * CALL ZTPMLQT('R','C',M,N,K,L,NB,AF(1,NP1),M,T,LDT,DF,M, - $ DF(1,NP1),M,WORK,INFO) - + $ DF(1,NP1),M,WORK,INFO) + * * Compute |D*QT - D*QT| / |D| * @@ -286,4 +286,4 @@ * DEALLOCATE ( A, AF, Q, R, RWORK, WORK, T, C, D, CF, DF) RETURN - END
\ No newline at end of file + END diff --git a/TESTING/LIN/ztsqr01.f b/TESTING/LIN/ztsqr01.f index 5f39ae7e..38ace9c8 100644 --- a/TESTING/LIN/ztsqr01.f +++ b/TESTING/LIN/ztsqr01.f @@ -2,19 +2,19 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE ZTSQR01(TSSW, M,N, MB, NB, RESULT) -* +* * .. Scalar Arguments .. * INTEGER M, N, MB * .. Return values .. * DOUBLE PRECISION RESULT(6) -* +* * *> \par Purpose: * ============= @@ -65,17 +65,17 @@ *> RESULT(2) = | I - Q^H Q | or | I - Q Q^H | *> RESULT(3) = | Q C - Q C | *> RESULT(4) = | Q^H C - Q^H C | -*> RESULT(5) = | C Q - C Q | +*> RESULT(5) = | C Q - C Q | *> RESULT(6) = | C Q^H - C Q^H | *> \endverbatim * * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date April 2012 * @@ -97,9 +97,9 @@ * ===================================================================== * * .. -* .. Local allocatable arrays +* .. Local allocatable arrays COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:), - $ R(:,:), RWORK(:), WORK( : ), T(:), + $ R(:,:), RWORK(:), WORK( : ), T(:), $ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:), LQ(:,:) * * .. Parameters .. @@ -122,24 +122,24 @@ EXTERNAL DLAMCH, ZLANGE, ZLANSY, LSAME, ILAENV * .. * .. Intrinsic Functions .. - INTRINSIC MAX, MIN + INTRINSIC MAX, MIN * .. Scalars in Common .. CHARACTER*32 srnamt * .. * .. Common blocks .. - COMMON / srnamc / srnamt + COMMON / srnamc / srnamt * .. * .. Data statements .. - DATA ISEED / 1988, 1989, 1990, 1991 / + DATA ISEED / 1988, 1989, 1990, 1991 / * * TEST TALL SKINNY OR SHORT WIDE * - TS = LSAME(TSSW, 'TS') -* + TS = LSAME(TSSW, 'TS') +* * TEST MATRICES WITH HALF OF MATRIX BEING ZEROS * TESTZEROS = .FALSE. -* +* EPS = DLAMCH( 'Epsilon' ) K = MIN(M,N) L = MAX(M,N,1) @@ -148,14 +148,14 @@ IF((K.GE.MNB).OR.(MNB.GE.L))THEN LT=MAX(1,L)*MNB+5 ELSE - LT=MAX(1,(L-K)/(MNB-K)+1)*L*MNB+5 + LT=MAX(1,(L-K)/(MNB-K)+1)*L*MNB+5 END IF * * Dynamically allocate local arrays * - ALLOCATE ( A(M,N), AF(M,N), Q(L,L), R(M,L), RWORK(L), - $ WORK(LWORK), T(LT), C(M,N), CF(M,N), + ALLOCATE ( A(M,N), AF(M,N), Q(L,L), R(M,L), RWORK(L), + $ WORK(LWORK), T(LT), C(M,N), CF(M,N), $ D(N,M), DF(N,M), LQ(L,N) ) * * Put random numbers into A and copy to AF @@ -183,7 +183,7 @@ * CALL ZLASET( 'Full', M, M, CZERO, ONE, Q, M ) srnamt = 'ZGEMQR' - CALL ZGEMQR( 'L', 'N', M, M, K, AF, M, T, LT, Q, M, + CALL ZGEMQR( 'L', 'N', M, M, K, AF, M, T, LT, Q, M, $ WORK, LWORK, INFO ) * * Copy R @@ -220,7 +220,7 @@ * Apply Q to C as Q*C * srnamt = 'ZGEMQR' - CALL ZGEMQR( 'L', 'N', M, N, K, AF, M, T, LT, CF, M, + CALL ZGEMQR( 'L', 'N', M, N, K, AF, M, T, LT, CF, M, $ WORK, LWORK, INFO) * * Compute |Q*C - Q*C| / |C| @@ -240,7 +240,7 @@ * Apply Q to C as QT*C * srnamt = 'ZGEMQR' - CALL ZGEMQR( 'L', 'C', M, N, K, AF, M, T, LT, CF, M, + CALL ZGEMQR( 'L', 'C', M, N, K, AF, M, T, LT, CF, M, $ WORK, LWORK, INFO) * * Compute |QT*C - QT*C| / |C| @@ -251,7 +251,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,M)*CNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -264,8 +264,8 @@ * Apply Q to D as D*Q * srnamt = 'ZGEMQR' - CALL ZGEMQR( 'R', 'N', N, M, K, AF, M, T, LT, DF, N, - $ WORK, LWORK, INFO) + CALL ZGEMQR( 'R', 'N', N, M, K, AF, M, T, LT, DF, N, + $ WORK, LWORK, INFO) * * Compute |D*Q - D*Q| / |D| * @@ -283,8 +283,8 @@ * * Apply Q to D as D*QT * - CALL ZGEMQR( 'R', 'C', N, M, K, AF, M, T, LT, DF, N, - $ WORK, LWORK, INFO) + CALL ZGEMQR( 'R', 'C', N, M, K, AF, M, T, LT, DF, N, + $ WORK, LWORK, INFO) * * Compute |D*QT - D*QT| / |D| * @@ -307,7 +307,7 @@ * CALL ZLASET( 'Full', N, N, CZERO, ONE, Q, N ) srnamt = 'ZGEMLQ' - CALL ZGEMLQ( 'R', 'N', N, N, K, AF, M, T, LT, Q, N, + CALL ZGEMLQ( 'R', 'N', N, N, K, AF, M, T, LT, Q, N, $ WORK, LWORK, INFO ) * * Copy R @@ -343,7 +343,7 @@ * * Apply Q to C as Q*C * - CALL ZGEMLQ( 'L', 'N', N, M, K, AF, M, T, LT, DF, N, + CALL ZGEMLQ( 'L', 'N', N, M, K, AF, M, T, LT, DF, N, $ WORK, LWORK, INFO) * * Compute |Q*D - Q*D| / |D| @@ -362,7 +362,7 @@ * * Apply Q to D as QT*D * - CALL ZGEMLQ( 'L', 'C', N, M, K, AF, M, T, LT, DF, N, + CALL ZGEMLQ( 'L', 'C', N, M, K, AF, M, T, LT, DF, N, $ WORK, LWORK, INFO) * * Compute |QT*D - QT*D| / |D| @@ -373,7 +373,7 @@ RESULT( 4 ) = RESID / (EPS*MAX(1,N)*DNORM) ELSE RESULT( 4 ) = ZERO - END IF + END IF * * Generate random n-by-m matrix D and a copy DF * @@ -385,8 +385,8 @@ * * Apply Q to C as C*Q * - CALL ZGEMLQ( 'R', 'N', M, N, K, AF, M, T, LT, CF, M, - $ WORK, LWORK, INFO) + CALL ZGEMLQ( 'R', 'N', M, N, K, AF, M, T, LT, CF, M, + $ WORK, LWORK, INFO) * * Compute |C*Q - C*Q| / |C| * @@ -404,8 +404,8 @@ * * Apply Q to D as D*QT * - CALL ZGEMLQ( 'R', 'C', M, N, K, AF, M, T, LT, CF, M, - $ WORK, LWORK, INFO) + CALL ZGEMLQ( 'R', 'C', M, N, K, AF, M, T, LT, CF, M, + $ WORK, LWORK, INFO) * * Compute |C*QT - C*QT| / |C| * @@ -424,4 +424,4 @@ DEALLOCATE ( A, AF, Q, R, RWORK, WORK, T, C, D, CF, DF) * RETURN - END
\ No newline at end of file + END |