Fix call to GEMQRT (no lwork) and lwork declarations

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|
 *