From 1b5afc4229b27adcfb5a1f761e308710ff4c7d56 Mon Sep 17 00:00:00 2001 From: julie Date: Tue, 10 Apr 2012 20:59:35 +0000 Subject: Fix call to GEMQRT (no lwork) and lwork declarations --- TESTING/LIN/cqrt04.f | 14 +++++++------- TESTING/LIN/dqrt04.f | 14 +++++++------- TESTING/LIN/sqrt04.f | 14 +++++++------- TESTING/LIN/zqrt04.f | 16 ++++++++-------- 4 files changed, 29 insertions(+), 29 deletions(-) diff --git a/TESTING/LIN/cqrt04.f b/TESTING/LIN/cqrt04.f index 68c9b879..cd8a3181 100644 --- a/TESTING/LIN/cqrt04.f +++ b/TESTING/LIN/cqrt04.f @@ -11,7 +11,7 @@ * SUBROUTINE CQRT04(M,N,NB,RESULT) * * .. Scalar Arguments .. -* INTEGER LWORK, M, N, NB, LDT +* INTEGER M, N, NB, LDT * .. Return values .. * REAL RESULT(6) * @@ -80,7 +80,7 @@ * November 2011 * * .. Scalar Arguments .. - INTEGER LWORK, M, N, NB, LDT + INTEGER M, N, NB, LDT * .. Return values .. REAL RESULT(6) * @@ -98,7 +98,7 @@ PARAMETER( ZERO = 0.0, ONE = (1.0,0.0), CZERO=(0.0,0.0) ) * .. * .. Local Scalars .. - INTEGER INFO, J, K, L + INTEGER INFO, J, K, L, LWORK REAL ANORM, EPS, RESID, CNORM, DNORM * .. * .. Local Arrays .. @@ -179,7 +179,7 @@ * Apply Q to C as Q*C * CALL CGEMQRT( 'L', 'N', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |Q*C - Q*C| / |C| * @@ -198,7 +198,7 @@ * Apply Q to C as QT*C * CALL CGEMQRT( 'L', 'C', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |QT*C - QT*C| / |C| * @@ -221,7 +221,7 @@ * Apply Q to D as D*Q * CALL CGEMQRT( 'R', 'N', N, M, K, NB, AF, M, T, NB, DF, N, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |D*Q - D*Q| / |D| * @@ -240,7 +240,7 @@ * Apply Q to D as D*QT * CALL CGEMQRT( 'R', 'C', N, M, K, NB, AF, M, T, NB, DF, N, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |D*QT - D*QT| / |D| * diff --git a/TESTING/LIN/dqrt04.f b/TESTING/LIN/dqrt04.f index 48d96d59..795e9b7e 100644 --- a/TESTING/LIN/dqrt04.f +++ b/TESTING/LIN/dqrt04.f @@ -11,7 +11,7 @@ * SUBROUTINE DQRT04(M,N,NB,RESULT) * * .. Scalar Arguments .. -* INTEGER LWORK, M, N, NB, LDT +* INTEGER M, N, NB, LDT * .. Return values .. * DOUBLE PRECISION RESULT(6) * @@ -80,7 +80,7 @@ * November 2011 * * .. Scalar Arguments .. - INTEGER LWORK, M, N, NB, LDT + INTEGER M, N, NB, LDT * .. Return values .. DOUBLE PRECISION RESULT(6) * @@ -97,7 +97,7 @@ PARAMETER( ZERO = 0.0, ONE = 1.0 ) * .. * .. Local Scalars .. - INTEGER INFO, J, K, L + INTEGER INFO, J, K, L, LWORK DOUBLE PRECISION ANORM, EPS, RESID, CNORM, DNORM * .. * .. Local Arrays .. @@ -177,7 +177,7 @@ * Apply Q to C as Q*C * CALL DGEMQRT( 'L', 'N', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |Q*C - Q*C| / |C| * @@ -196,7 +196,7 @@ * Apply Q to C as QT*C * CALL DGEMQRT( 'L', 'T', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |QT*C - QT*C| / |C| * @@ -219,7 +219,7 @@ * Apply Q to D as D*Q * CALL DGEMQRT( 'R', 'N', N, M, K, NB, AF, M, T, NB, DF, N, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |D*Q - D*Q| / |D| * @@ -238,7 +238,7 @@ * Apply Q to D as D*QT * CALL DGEMQRT( 'R', 'T', N, M, K, NB, AF, M, T, NB, DF, N, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |D*QT - D*QT| / |D| * diff --git a/TESTING/LIN/sqrt04.f b/TESTING/LIN/sqrt04.f index ac1841e6..9f7ec408 100644 --- a/TESTING/LIN/sqrt04.f +++ b/TESTING/LIN/sqrt04.f @@ -11,7 +11,7 @@ * SUBROUTINE SQRT04(M,N,NB,RESULT) * * .. Scalar Arguments .. -* INTEGER LWORK, M, N, NB, LDT +* INTEGER M, N, NB, LDT * .. Return values .. * REAL RESULT(6) * @@ -80,7 +80,7 @@ * November 2011 * * .. Scalar Arguments .. - INTEGER LWORK, M, N, NB, LDT + INTEGER M, N, NB, LDT * .. Return values .. REAL RESULT(6) * @@ -97,7 +97,7 @@ PARAMETER( ZERO = 0.0, ONE = 1.0 ) * .. * .. Local Scalars .. - INTEGER INFO, J, K, L + INTEGER INFO, J, K, L, LWORK REAL ANORM, EPS, RESID, CNORM, DNORM * .. * .. Local Arrays .. @@ -178,7 +178,7 @@ * Apply Q to C as Q*C * CALL SGEMQRT( 'L', 'N', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |Q*C - Q*C| / |C| * @@ -197,7 +197,7 @@ * Apply Q to C as QT*C * CALL SGEMQRT( 'L', 'T', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |QT*C - QT*C| / |C| * @@ -220,7 +220,7 @@ * Apply Q to D as D*Q * CALL SGEMQRT( 'R', 'N', N, M, K, NB, AF, M, T, NB, DF, N, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |D*Q - D*Q| / |D| * @@ -239,7 +239,7 @@ * Apply Q to D as D*QT * CALL SGEMQRT( 'R', 'T', N, M, K, NB, AF, M, T, NB, DF, N, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |D*QT - D*QT| / |D| * diff --git a/TESTING/LIN/zqrt04.f b/TESTING/LIN/zqrt04.f index e453c594..cec97122 100644 --- a/TESTING/LIN/zqrt04.f +++ b/TESTING/LIN/zqrt04.f @@ -11,7 +11,7 @@ * SUBROUTINE ZQRT04(M,N,NB,RESULT) * * .. Scalar Arguments .. -* INTEGER LWORK, M, N, NB, LDT +* INTEGER M, N, NB, LDT * .. Return values .. * DOUBLE PRECISION RESULT(6) * @@ -80,7 +80,7 @@ * November 2011 * * .. Scalar Arguments .. - INTEGER LWORK, M, N, NB, LDT + INTEGER M, N, NB, LDT * .. Return values .. DOUBLE PRECISION RESULT(6) * @@ -98,7 +98,7 @@ PARAMETER( ZERO = 0.0, ONE = (1.0,0.0), CZERO=(0.0,0.0) ) * .. * .. Local Scalars .. - INTEGER INFO, J, K, L + INTEGER INFO, J, K, L, LWORK DOUBLE PRECISION ANORM, EPS, RESID, CNORM, DNORM * .. * .. Local Arrays .. @@ -142,7 +142,7 @@ * Generate the m-by-m matrix Q * CALL ZLASET( 'Full', M, M, CZERO, ONE, Q, M ) - CALL ZGEMQRT( 'R', 'N', M, M, K, NB, AF, M, T, LDT, Q, M, + CALL ZGEMQRT( 'R', 'N', M, M, K, NB, AF, M, T, LDT, Q, M, $ WORK, INFO ) * * Copy R @@ -179,7 +179,7 @@ * Apply Q to C as Q*C * CALL ZGEMQRT( 'L', 'N', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |Q*C - Q*C| / |C| * @@ -198,7 +198,7 @@ * Apply Q to C as QT*C * CALL ZGEMQRT( 'L', 'C', M, N, K, NB, AF, M, T, NB, CF, M, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |QT*C - QT*C| / |C| * @@ -221,7 +221,7 @@ * Apply Q to D as D*Q * CALL ZGEMQRT( 'R', 'N', N, M, K, NB, AF, M, T, NB, DF, N, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |D*Q - D*Q| / |D| * @@ -240,7 +240,7 @@ * Apply Q to D as D*QT * CALL ZGEMQRT( 'R', 'C', N, M, K, NB, AF, M, T, NB, DF, N, - $ WORK, LWORK, INFO) + $ WORK, INFO) * * Compute |D*QT - D*QT| / |D| * -- cgit v1.2.3