diff options
Diffstat (limited to 'TESTING/LIN/dlqt04.f')
| -rw-r--r-- | TESTING/LIN/dlqt04.f | 48 |
1 files changed, 24 insertions, 24 deletions
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| * |