projects / platform / upstream / lapack.git / commitdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | patch | inline | side by side (parent: 85f43a2)
Fix call to GEMQRT (no lwork) and lwork declarations
authorjulie <julielangou@users.noreply.github.com>
Tue, 10 Apr 2012 20:59:35 +0000 (20:59 +0000)
committerjulie <julielangou@users.noreply.github.com>
Tue, 10 Apr 2012 20:59:35 +0000 (20:59 +0000)
TESTING/LIN/cqrt04.f patch | blob | history
TESTING/LIN/dqrt04.f patch | blob | history
TESTING/LIN/sqrt04.f patch | blob | history
TESTING/LIN/zqrt04.f patch | blob | history

diff --git a/TESTING/LIN/cqrt04.f b/TESTING/LIN/cqrt04.f
index 68c9b879ea0c34ce882c5468a6635fb6a33c74f1..cd8a31812bb2a46b09822bbf88fa8bf8727f6fd0 100644 (file)
--- 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 48d96d5962e6c0850b3db8718c1a5c8ec0137cd2..795e9b7e6fbc744dc44d26f202cba5952677f372 100644 (file)
--- 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 ac1841e6bbcc3c7564d8104855d7dfd7af266c70..9f7ec4081538df909229cec22a3acfda549cb8ad 100644 (file)
--- 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 e453c594beb836086469671d26e02c28a8d217f8..cec97122297258127da8cc3234fd2a9909466a17 100644 (file)
--- 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|
 *
Domain: System / Base;