Step 1 of xlarfp : back out the old changes (need to run testing before going on)

diff --git a/SRC/cgelq2.f b/SRC/cgelq2.f
index 529c6a3..583f5ee 100644 (file)
--- a/SRC/cgelq2.f
+++ b/SRC/cgelq2.f
@@ -74,7 +74,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLACGV, CLARF, CLARFP, XERBLA
+      EXTERNAL           CLACGV, CLARF, CLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -104,7 +104,7 @@
 *
          CALL CLACGV( N-I+1, A( I, I ), LDA )
          ALPHA = A( I, I )
-         CALL CLARFP( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
+         CALL CLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
      $                TAU( I ) )
          IF( I.LT.M ) THEN
 *

diff --git a/SRC/cgeql2.f b/SRC/cgeql2.f
index c0aea5b..bf0bf8f 100644 (file)
--- a/SRC/cgeql2.f
+++ b/SRC/cgeql2.f
@@ -76,7 +76,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLARF, CLARFP, XERBLA
+      EXTERNAL           CLARF, CLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          CONJG, MAX, MIN
@@ -106,7 +106,7 @@
 *        A(1:m-k+i-1,n-k+i)
 *
          ALPHA = A( M-K+I, N-K+I )
-         CALL CLARFP( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )
+         CALL CLARFG( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )
 *
 *        Apply H(i)' to A(1:m-k+i,1:n-k+i-1) from the left
 *

diff --git a/SRC/cgeqpf.f b/SRC/cgeqpf.f
index d63127f..8e85f31 100644 (file)
--- a/SRC/cgeqpf.f
+++ b/SRC/cgeqpf.f
@@ -96,7 +96,7 @@
       COMPLEX            AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CGEQR2, CLARF, CLARFP, CSWAP, CUNM2R, XERBLA
+      EXTERNAL           CGEQR2, CLARF, CLARFG, CSWAP, CUNM2R, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, CMPLX, CONJG, MAX, MIN, SQRT
@@ -186,7 +186,7 @@
 *           Generate elementary reflector H(i)
 *
             AII = A( I, I )
-            CALL CLARFP( M-I+1, AII, A( MIN( I+1, M ), I ), 1,
+            CALL CLARFG( M-I+1, AII, A( MIN( I+1, M ), I ), 1,
      $                   TAU( I ) )
             A( I, I ) = AII
 *

diff --git a/SRC/cgeqr2.f b/SRC/cgeqr2.f
index 135a167..7697baf 100644 (file)
--- a/SRC/cgeqr2.f
+++ b/SRC/cgeqr2.f
@@ -74,7 +74,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLARF, CLARFP, XERBLA
+      EXTERNAL           CLARF, CLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          CONJG, MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i+1:m,i)
 *
-         CALL CLARFP( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
+         CALL CLARFG( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
      $                TAU( I ) )
          IF( I.LT.N ) THEN
 *

diff --git a/SRC/cgerq2.f b/SRC/cgerq2.f
index 50b6192..14f1ecb 100644 (file)
--- a/SRC/cgerq2.f
+++ b/SRC/cgerq2.f
@@ -76,7 +76,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLACGV, CLARF, CLARFP, XERBLA
+      EXTERNAL           CLACGV, CLARF, CLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -107,7 +107,7 @@
 *
          CALL CLACGV( N-K+I, A( M-K+I, 1 ), LDA )
          ALPHA = A( M-K+I, N-K+I )
-         CALL CLARFP( N-K+I, ALPHA, A( M-K+I, 1 ), LDA,
+         CALL CLARFG( N-K+I, ALPHA, A( M-K+I, 1 ), LDA,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right

diff --git a/SRC/claqp2.f b/SRC/claqp2.f
index b65f615..9c7d99b 100644 (file)
--- a/SRC/claqp2.f
+++ b/SRC/claqp2.f
@@ -91,7 +91,7 @@
       COMPLEX            AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLARF, CLARFP, CSWAP
+      EXTERNAL           CLARF, CLARFG, CSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, CONJG, MAX, MIN, SQRT
@@ -128,10 +128,10 @@
 *        Generate elementary reflector H(i).
 *
          IF( OFFPI.LT.M ) THEN
-            CALL CLARFP( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
+            CALL CLARFG( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
      $                   TAU( I ) )
          ELSE
-            CALL CLARFP( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
+            CALL CLARFG( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
          END IF
 *
          IF( I.LT.N ) THEN

diff --git a/SRC/claqps.f b/SRC/claqps.f
index cefdf04..30daffb 100644 (file)
--- a/SRC/claqps.f
+++ b/SRC/claqps.f
@@ -108,7 +108,7 @@
       COMPLEX            AKK
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CGEMM, CGEMV, CLARFP, CSWAP
+      EXTERNAL           CGEMM, CGEMV, CLARFG, CSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, CONJG, MAX, MIN, NINT, REAL, SQRT
@@ -162,9 +162,9 @@
 *        Generate elementary reflector H(k).
 *
          IF( RK.LT.M ) THEN
-            CALL CLARFP( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
+            CALL CLARFG( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
          ELSE
-            CALL CLARFP( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
+            CALL CLARFG( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
          END IF
 *
          AKK = A( RK, K )

diff --git a/SRC/clatrz.f b/SRC/clatrz.f
index da5e42d..68d22a1 100644 (file)
--- a/SRC/clatrz.f
+++ b/SRC/clatrz.f
@@ -91,7 +91,7 @@
       COMPLEX            ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CLACGV, CLARFP, CLARZ
+      EXTERNAL           CLACGV, CLARFG, CLARZ
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          CONJG
@@ -116,7 +116,7 @@
 *
          CALL CLACGV( L, A( I, N-L+1 ), LDA )
          ALPHA = CONJG( A( I, I ) )
-         CALL CLARFP( L+1, ALPHA, A( I, N-L+1 ), LDA, TAU( I ) )
+         CALL CLARFG( L+1, ALPHA, A( I, N-L+1 ), LDA, TAU( I ) )
          TAU( I ) = CONJG( TAU( I ) )
 *
 *        Apply H(i) to A(1:i-1,i:n) from the right

diff --git a/SRC/ctzrqf.f b/SRC/ctzrqf.f
index fb2bc23..38a38a7 100644 (file)
--- a/SRC/ctzrqf.f
+++ b/SRC/ctzrqf.f
@@ -98,7 +98,7 @@
       INTRINSIC          CONJG, MAX, MIN
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           CAXPY, CCOPY, CGEMV, CGERC, CLACGV, CLARFP,
+      EXTERNAL           CAXPY, CCOPY, CGEMV, CGERC, CLACGV, CLARFG,
      $                   XERBLA
 *     ..
 *     .. Executable Statements ..
@@ -136,7 +136,7 @@
             A( K, K ) = CONJG( A( K, K ) )
             CALL CLACGV( N-M, A( K, M1 ), LDA )
             ALPHA = A( K, K )
-            CALL CLARFP( N-M+1, ALPHA, A( K, M1 ), LDA, TAU( K ) )
+            CALL CLARFG( N-M+1, ALPHA, A( K, M1 ), LDA, TAU( K ) )
             A( K, K ) = ALPHA
             TAU( K ) = CONJG( TAU( K ) )
 *

diff --git a/SRC/dgelq2.f b/SRC/dgelq2.f
index 96ef32a..0eb8ca8 100644 (file)
--- a/SRC/dgelq2.f
+++ b/SRC/dgelq2.f
@@ -74,7 +74,7 @@
       DOUBLE PRECISION   AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, XERBLA
+      EXTERNAL           DLARF, DLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i,i+1:n)
 *
-         CALL DLARFP( N-I+1, A( I, I ), A( I, MIN( I+1, N ) ), LDA,
+         CALL DLARFG( N-I+1, A( I, I ), A( I, MIN( I+1, N ) ), LDA,
      $                TAU( I ) )
          IF( I.LT.M ) THEN
 *

diff --git a/SRC/dgeql2.f b/SRC/dgeql2.f
index a684249..02756c9 100644 (file)
--- a/SRC/dgeql2.f
+++ b/SRC/dgeql2.f
@@ -76,7 +76,7 @@
       DOUBLE PRECISION   AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, XERBLA
+      EXTERNAL           DLARF, DLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -105,7 +105,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        A(1:m-k+i-1,n-k+i)
 *
-         CALL DLARFP( M-K+I, A( M-K+I, N-K+I ), A( 1, N-K+I ), 1,
+         CALL DLARFG( M-K+I, A( M-K+I, N-K+I ), A( 1, N-K+I ), 1,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i,1:n-k+i-1) from the left

diff --git a/SRC/dgeqpf.f b/SRC/dgeqpf.f
index 4f78350..2c86e5a 100644 (file)
--- a/SRC/dgeqpf.f
+++ b/SRC/dgeqpf.f
@@ -92,7 +92,7 @@
       DOUBLE PRECISION   AII, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DGEQR2, DLARF, DLARFP, DORM2R, DSWAP, XERBLA
+      EXTERNAL           DGEQR2, DLARF, DLARFG, DORM2R, DSWAP, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, SQRT
@@ -182,9 +182,9 @@
 *           Generate elementary reflector H(i)
 *
             IF( I.LT.M ) THEN
-               CALL DLARFP( M-I+1, A( I, I ), A( I+1, I ), 1, TAU( I ) )
+               CALL DLARFG( M-I+1, A( I, I ), A( I+1, I ), 1, TAU( I ) )
             ELSE
-               CALL DLARFP( 1, A( M, M ), A( M, M ), 1, TAU( M ) )
+               CALL DLARFG( 1, A( M, M ), A( M, M ), 1, TAU( M ) )
             END IF
 *
             IF( I.LT.N ) THEN

diff --git a/SRC/dgeqr2.f b/SRC/dgeqr2.f
index 78c3aac..a1b9f23 100644 (file)
--- a/SRC/dgeqr2.f
+++ b/SRC/dgeqr2.f
@@ -74,7 +74,7 @@
       DOUBLE PRECISION   AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, XERBLA
+      EXTERNAL           DLARF, DLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i+1:m,i)
 *
-         CALL DLARFP( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
+         CALL DLARFG( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
      $                TAU( I ) )
          IF( I.LT.N ) THEN
 *

diff --git a/SRC/dgerq2.f b/SRC/dgerq2.f
index e635733..4bf2340 100644 (file)
--- a/SRC/dgerq2.f
+++ b/SRC/dgerq2.f
@@ -76,7 +76,7 @@
       DOUBLE PRECISION   AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, XERBLA
+      EXTERNAL           DLARF, DLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -105,7 +105,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        A(m-k+i,1:n-k+i-1)
 *
-         CALL DLARFP( N-K+I, A( M-K+I, N-K+I ), A( M-K+I, 1 ), LDA,
+         CALL DLARFG( N-K+I, A( M-K+I, N-K+I ), A( M-K+I, 1 ), LDA,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right

diff --git a/SRC/dlaqp2.f b/SRC/dlaqp2.f
index df455d2..b579fb4 100644 (file)
--- a/SRC/dlaqp2.f
+++ b/SRC/dlaqp2.f
@@ -88,7 +88,7 @@
       DOUBLE PRECISION   AII, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARF, DLARFP, DSWAP
+      EXTERNAL           DLARF, DLARFG, DSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, SQRT
@@ -125,10 +125,10 @@
 *        Generate elementary reflector H(i).
 *
          IF( OFFPI.LT.M ) THEN
-            CALL DLARFP( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
+            CALL DLARFG( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
      $                   TAU( I ) )
          ELSE
-            CALL DLARFP( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
+            CALL DLARFG( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
          END IF
 *
          IF( I.LE.N ) THEN

diff --git a/SRC/dlaqps.f b/SRC/dlaqps.f
index 9e00472..ddb7222 100644 (file)
--- a/SRC/dlaqps.f
+++ b/SRC/dlaqps.f
@@ -104,7 +104,7 @@
       DOUBLE PRECISION   AKK, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DGEMM, DGEMV, DLARFP, DSWAP
+      EXTERNAL           DGEMM, DGEMV, DLARFG, DSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DBLE, MAX, MIN, NINT, SQRT
@@ -152,9 +152,9 @@
 *        Generate elementary reflector H(k).
 *
          IF( RK.LT.M ) THEN
-            CALL DLARFP( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
+            CALL DLARFG( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
          ELSE
-            CALL DLARFP( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
+            CALL DLARFG( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
          END IF
 *
          AKK = A( RK, K )

diff --git a/SRC/dlatrz.f b/SRC/dlatrz.f
index ca74935..4d4d489 100644 (file)
--- a/SRC/dlatrz.f
+++ b/SRC/dlatrz.f
@@ -90,7 +90,7 @@
       INTEGER            I
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DLARFP, DLARZ
+      EXTERNAL           DLARFG, DLARZ
 *     ..
 *     .. Executable Statements ..
 *
@@ -112,7 +112,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        [ A(i,i) A(i,n-l+1:n) ]
 *
-         CALL DLARFP( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )
+         CALL DLARFG( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )
 *
 *        Apply H(i) to A(1:i-1,i:n) from the right
 *

diff --git a/SRC/dtzrqf.f b/SRC/dtzrqf.f
index 9777036..c1b425b 100644 (file)
--- a/SRC/dtzrqf.f
+++ b/SRC/dtzrqf.f
@@ -96,7 +96,7 @@
       INTRINSIC          MAX, MIN
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           DAXPY, DCOPY, DGEMV, DGER, DLARFP, XERBLA
+      EXTERNAL           DAXPY, DCOPY, DGEMV, DGER, DLARFG, XERBLA
 *     ..
 *     .. Executable Statements ..
 *
@@ -130,7 +130,7 @@
 *           Use a Householder reflection to zero the kth row of A.
 *           First set up the reflection.
 *
-            CALL DLARFP( N-M+1, A( K, K ), A( K, M1 ), LDA, TAU( K ) )
+            CALL DLARFG( N-M+1, A( K, K ), A( K, M1 ), LDA, TAU( K ) )
 *
             IF( ( TAU( K ).NE.ZERO ) .AND. ( K.GT.1 ) ) THEN
 *

diff --git a/SRC/sgelq2.f b/SRC/sgelq2.f
index 16db9bb..b4519a8 100644 (file)
--- a/SRC/sgelq2.f
+++ b/SRC/sgelq2.f
@@ -74,7 +74,7 @@
       REAL               AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, XERBLA
+      EXTERNAL           SLARF, SLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i,i+1:n)
 *
-         CALL SLARFP( N-I+1, A( I, I ), A( I, MIN( I+1, N ) ), LDA,
+         CALL SLARFG( N-I+1, A( I, I ), A( I, MIN( I+1, N ) ), LDA,
      $                TAU( I ) )
          IF( I.LT.M ) THEN
 *

diff --git a/SRC/sgeql2.f b/SRC/sgeql2.f
index 6d3bb75..6e64d81 100644 (file)
--- a/SRC/sgeql2.f
+++ b/SRC/sgeql2.f
@@ -76,7 +76,7 @@
       REAL               AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, XERBLA
+      EXTERNAL           SLARF, SLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -105,7 +105,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        A(1:m-k+i-1,n-k+i)
 *
-         CALL SLARFP( M-K+I, A( M-K+I, N-K+I ), A( 1, N-K+I ), 1,
+         CALL SLARFG( M-K+I, A( M-K+I, N-K+I ), A( 1, N-K+I ), 1,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i,1:n-k+i-1) from the left

diff --git a/SRC/sgeqpf.f b/SRC/sgeqpf.f
index 1781630..b32f8c8 100644 (file)
--- a/SRC/sgeqpf.f
+++ b/SRC/sgeqpf.f
@@ -92,7 +92,7 @@
       REAL               AII, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SGEQR2, SLARF, SLARFP, SORM2R, SSWAP, XERBLA
+      EXTERNAL           SGEQR2, SLARF, SLARFG, SORM2R, SSWAP, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, SQRT
@@ -182,9 +182,9 @@
 *           Generate elementary reflector H(i)
 *
             IF( I.LT.M ) THEN
-               CALL SLARFP( M-I+1, A( I, I ), A( I+1, I ), 1, TAU( I ) )
+               CALL SLARFG( M-I+1, A( I, I ), A( I+1, I ), 1, TAU( I ) )
             ELSE
-               CALL SLARFP( 1, A( M, M ), A( M, M ), 1, TAU( M ) )
+               CALL SLARFG( 1, A( M, M ), A( M, M ), 1, TAU( M ) )
             END IF
 *
             IF( I.LT.N ) THEN

diff --git a/SRC/sgeqr2.f b/SRC/sgeqr2.f
index 4353cd6..c70fb98 100644 (file)
--- a/SRC/sgeqr2.f
+++ b/SRC/sgeqr2.f
@@ -74,7 +74,7 @@
       REAL               AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, XERBLA
+      EXTERNAL           SLARF, SLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i+1:m,i)
 *
-         CALL SLARFP( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
+         CALL SLARFG( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
      $                TAU( I ) )
          IF( I.LT.N ) THEN
 *

diff --git a/SRC/sgerq2.f b/SRC/sgerq2.f
index e527909..b8f11be 100644 (file)
--- a/SRC/sgerq2.f
+++ b/SRC/sgerq2.f
@@ -76,7 +76,7 @@
       REAL               AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, XERBLA
+      EXTERNAL           SLARF, SLARFG, XERBLA
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -105,7 +105,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        A(m-k+i,1:n-k+i-1)
 *
-         CALL SLARFP( N-K+I, A( M-K+I, N-K+I ), A( M-K+I, 1 ), LDA,
+         CALL SLARFG( N-K+I, A( M-K+I, N-K+I ), A( M-K+I, 1 ), LDA,
      $                TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right

diff --git a/SRC/slaqp2.f b/SRC/slaqp2.f
index 2f51ecd..10d5790 100644 (file)
--- a/SRC/slaqp2.f
+++ b/SRC/slaqp2.f
@@ -88,7 +88,7 @@
       REAL               AII, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARF, SLARFP, SSWAP
+      EXTERNAL           SLARF, SLARFG, SSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, SQRT
@@ -125,10 +125,10 @@
 *        Generate elementary reflector H(i).
 *
          IF( OFFPI.LT.M ) THEN
-            CALL SLARFP( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
+            CALL SLARFG( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
      $                   TAU( I ) )
          ELSE
-            CALL SLARFP( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
+            CALL SLARFG( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
          END IF
 *
          IF( I.LT.N ) THEN

diff --git a/SRC/slaqps.f b/SRC/slaqps.f
index 3ff3697..343e16d 100644 (file)
--- a/SRC/slaqps.f
+++ b/SRC/slaqps.f
@@ -104,7 +104,7 @@
       REAL               AKK, TEMP, TEMP2, TOL3Z
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SGEMM, SGEMV, SLARFP, SSWAP
+      EXTERNAL           SGEMM, SGEMV, SLARFG, SSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, MAX, MIN, NINT, REAL, SQRT
@@ -152,9 +152,9 @@
 *        Generate elementary reflector H(k).
 *
          IF( RK.LT.M ) THEN
-            CALL SLARFP( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
+            CALL SLARFG( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
          ELSE
-            CALL SLARFP( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
+            CALL SLARFG( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
          END IF
 *
          AKK = A( RK, K )

diff --git a/SRC/slatrz.f b/SRC/slatrz.f
index 3b0893a..1adefa7 100644 (file)
--- a/SRC/slatrz.f
+++ b/SRC/slatrz.f
@@ -90,7 +90,7 @@
       INTEGER            I
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SLARFP, SLARZ
+      EXTERNAL           SLARFG, SLARZ
 *     ..
 *     .. Executable Statements ..
 *
@@ -112,7 +112,7 @@
 *        Generate elementary reflector H(i) to annihilate
 *        [ A(i,i) A(i,n-l+1:n) ]
 *
-         CALL SLARFP( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )
+         CALL SLARFG( L+1, A( I, I ), A( I, N-L+1 ), LDA, TAU( I ) )
 *
 *        Apply H(i) to A(1:i-1,i:n) from the right
 *

diff --git a/SRC/stzrqf.f b/SRC/stzrqf.f
index 30b4e20..a64965b 100644 (file)
--- a/SRC/stzrqf.f
+++ b/SRC/stzrqf.f
@@ -96,7 +96,7 @@
       INTRINSIC          MAX, MIN
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           SAXPY, SCOPY, SGEMV, SGER, SLARFP, XERBLA
+      EXTERNAL           SAXPY, SCOPY, SGEMV, SGER, SLARFG, XERBLA
 *     ..
 *     .. Executable Statements ..
 *
@@ -130,7 +130,7 @@
 *           Use a Householder reflection to zero the kth row of A.
 *           First set up the reflection.
 *
-            CALL SLARFP( N-M+1, A( K, K ), A( K, M1 ), LDA, TAU( K ) )
+            CALL SLARFG( N-M+1, A( K, K ), A( K, M1 ), LDA, TAU( K ) )
 *
             IF( ( TAU( K ).NE.ZERO ) .AND. ( K.GT.1 ) ) THEN
 *

diff --git a/SRC/zgelq2.f b/SRC/zgelq2.f
index 04ad7df..8b91ff6 100644 (file)
--- a/SRC/zgelq2.f
+++ b/SRC/zgelq2.f
@@ -74,7 +74,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFP
+      EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFG
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -104,7 +104,7 @@
 *
          CALL ZLACGV( N-I+1, A( I, I ), LDA )
          ALPHA = A( I, I )
-         CALL ZLARFP( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
+         CALL ZLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
      $                TAU( I ) )
          IF( I.LT.M ) THEN
 *

diff --git a/SRC/zgeql2.f b/SRC/zgeql2.f
index b8fb863..3ae6d72 100644 (file)
--- a/SRC/zgeql2.f
+++ b/SRC/zgeql2.f
@@ -76,7 +76,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZLARF, ZLARFP
+      EXTERNAL           XERBLA, ZLARF, ZLARFG
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          DCONJG, MAX, MIN
@@ -106,7 +106,7 @@
 *        A(1:m-k+i-1,n-k+i)
 *
          ALPHA = A( M-K+I, N-K+I )
-         CALL ZLARFP( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )
+         CALL ZLARFG( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )
 *
 *        Apply H(i)' to A(1:m-k+i,1:n-k+i-1) from the left
 *

diff --git a/SRC/zgeqpf.f b/SRC/zgeqpf.f
index 2400089..51e93c3 100644 (file)
--- a/SRC/zgeqpf.f
+++ b/SRC/zgeqpf.f
@@ -96,7 +96,7 @@
       COMPLEX*16         AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZGEQR2, ZLARF, ZLARFP, ZSWAP, ZUNM2R
+      EXTERNAL           XERBLA, ZGEQR2, ZLARF, ZLARFG, ZSWAP, ZUNM2R
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DCMPLX, DCONJG, MAX, MIN, SQRT
@@ -186,7 +186,7 @@
 *           Generate elementary reflector H(i)
 *
             AII = A( I, I )
-            CALL ZLARFP( M-I+1, AII, A( MIN( I+1, M ), I ), 1,
+            CALL ZLARFG( M-I+1, AII, A( MIN( I+1, M ), I ), 1,
      $                   TAU( I ) )
             A( I, I ) = AII
 *

diff --git a/SRC/zgeqr2.f b/SRC/zgeqr2.f
index 2924f8c..76dcba4 100644 (file)
--- a/SRC/zgeqr2.f
+++ b/SRC/zgeqr2.f
@@ -74,7 +74,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZLARF, ZLARFP
+      EXTERNAL           XERBLA, ZLARF, ZLARFG
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          DCONJG, MAX, MIN
@@ -102,7 +102,7 @@
 *
 *        Generate elementary reflector H(i) to annihilate A(i+1:m,i)
 *
-         CALL ZLARFP( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
+         CALL ZLARFG( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1,
      $                TAU( I ) )
          IF( I.LT.N ) THEN
 *

diff --git a/SRC/zgerq2.f b/SRC/zgerq2.f
index 7fcccc4..989330b 100644 (file)
--- a/SRC/zgerq2.f
+++ b/SRC/zgerq2.f
@@ -76,7 +76,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFP
+      EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFG
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN
@@ -107,7 +107,7 @@
 *
          CALL ZLACGV( N-K+I, A( M-K+I, 1 ), LDA )
          ALPHA = A( M-K+I, N-K+I )
-         CALL ZLARFP( N-K+I, ALPHA, A( M-K+I, 1 ), LDA, TAU( I ) )
+         CALL ZLARFG( N-K+I, ALPHA, A( M-K+I, 1 ), LDA, TAU( I ) )
 *
 *        Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right
 *

diff --git a/SRC/zlaqp2.f b/SRC/zlaqp2.f
index 2354a00..85a6223 100644 (file)
--- a/SRC/zlaqp2.f
+++ b/SRC/zlaqp2.f
@@ -91,7 +91,7 @@
       COMPLEX*16         AII
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ZLARF, ZLARFP, ZSWAP
+      EXTERNAL           ZLARF, ZLARFG, ZSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DCONJG, MAX, MIN, SQRT
@@ -128,10 +128,10 @@
 *        Generate elementary reflector H(i).
 *
          IF( OFFPI.LT.M ) THEN
-            CALL ZLARFP( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
+            CALL ZLARFG( M-OFFPI+1, A( OFFPI, I ), A( OFFPI+1, I ), 1,
      $                   TAU( I ) )
          ELSE
-            CALL ZLARFP( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
+            CALL ZLARFG( 1, A( M, I ), A( M, I ), 1, TAU( I ) )
          END IF
 *
          IF( I.LT.N ) THEN

diff --git a/SRC/zlaqps.f b/SRC/zlaqps.f
index fdff5f3..571af2f 100644 (file)
--- a/SRC/zlaqps.f
+++ b/SRC/zlaqps.f
@@ -103,7 +103,7 @@
       COMPLEX*16         AKK
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ZGEMM, ZGEMV, ZLARFP, ZSWAP
+      EXTERNAL           ZGEMM, ZGEMV, ZLARFG, ZSWAP
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          ABS, DBLE, DCONJG, MAX, MIN, NINT, SQRT
@@ -157,9 +157,9 @@
 *        Generate elementary reflector H(k).
 *
          IF( RK.LT.M ) THEN
-            CALL ZLARFP( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
+            CALL ZLARFG( M-RK+1, A( RK, K ), A( RK+1, K ), 1, TAU( K ) )
          ELSE
-            CALL ZLARFP( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
+            CALL ZLARFG( 1, A( RK, K ), A( RK, K ), 1, TAU( K ) )
          END IF
 *
          AKK = A( RK, K )

diff --git a/SRC/zlatrz.f b/SRC/zlatrz.f
index b2b837d..06f021f 100644 (file)
--- a/SRC/zlatrz.f
+++ b/SRC/zlatrz.f
@@ -91,7 +91,7 @@
       COMPLEX*16         ALPHA
 *     ..
 *     .. External Subroutines ..
-      EXTERNAL           ZLACGV, ZLARFP, ZLARZ
+      EXTERNAL           ZLACGV, ZLARFG, ZLARZ
 *     ..
 *     .. Intrinsic Functions ..
       INTRINSIC          DCONJG
@@ -116,7 +116,7 @@
 *
          CALL ZLACGV( L, A( I, N-L+1 ), LDA )
          ALPHA = DCONJG( A( I, I ) )
-         CALL ZLARFP( L+1, ALPHA, A( I, N-L+1 ), LDA, TAU( I ) )
+         CALL ZLARFG( L+1, ALPHA, A( I, N-L+1 ), LDA, TAU( I ) )
          TAU( I ) = DCONJG( TAU( I ) )
 *
 *        Apply H(i) to A(1:i-1,i:n) from the right

diff --git a/SRC/ztzrqf.f b/SRC/ztzrqf.f
index 4571ca6..1744ac6 100644 (file)
--- a/SRC/ztzrqf.f
+++ b/SRC/ztzrqf.f
@@ -99,7 +99,7 @@
 *     ..
 *     .. External Subroutines ..
       EXTERNAL           XERBLA, ZAXPY, ZCOPY, ZGEMV, ZGERC, ZLACGV,
-     $                   ZLARFP
+     $                   ZLARFG
 *     ..
 *     .. Executable Statements ..
 *
@@ -136,7 +136,7 @@
             A( K, K ) = DCONJG( A( K, K ) )
             CALL ZLACGV( N-M, A( K, M1 ), LDA )
             ALPHA = A( K, K )
-            CALL ZLARFP( N-M+1, ALPHA, A( K, M1 ), LDA, TAU( K ) )
+            CALL ZLARFG( N-M+1, ALPHA, A( K, M1 ), LDA, TAU( K ) )
             A( K, K ) = ALPHA
             TAU( K ) = DCONJG( TAU( K ) )
 *

diff --git a/TESTING/LIN/alahd.f b/TESTING/LIN/alahd.f
index d147bc2..8a5637e 100644 (file)
--- a/TESTING/LIN/alahd.f
+++ b/TESTING/LIN/alahd.f
@@ -357,7 +357,6 @@
          WRITE( IOUNIT, FMT = 9942 )5, 'M'
          WRITE( IOUNIT, FMT = 9941 )6, 'M'
          WRITE( IOUNIT, FMT = 9960 )7
-         WRITE( IOUNIT, FMT = 6660 )8
          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
 *
       ELSE IF( LSAMEN( 2, P2, 'LQ' ) ) THEN
@@ -727,7 +726,6 @@
      $      '( N * norm(A) * norm(AINV) * EPS )' )
  9960 FORMAT( 3X, I2, ': norm( B - A * X )  / ',
      $      '( norm(A) * norm(X) * EPS )' )
- 6660 FORMAT( 3X, I2, ': diagonal is not non-negative')
  9959 FORMAT( 3X, I2, ': norm( X - XACT )   / ',
      $      '( norm(XACT) * CNDNUM * EPS )' )
  9958 FORMAT( 3X, I2, ': norm( X - XACT )   / ',

diff --git a/TESTING/LIN/cchklq.f b/TESTING/LIN/cchklq.f
index 7f8cba1..912ce78 100644 (file)
--- a/TESTING/LIN/cchklq.f
+++ b/TESTING/LIN/cchklq.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       REAL               ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       REAL               RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            CGENND
-      EXTERNAL           CGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, CERRLQ, CGELQS, CGET02,
      $                   CLACPY, CLARHS, CLATB4, CLATMS, CLQT01, CLQT02,
@@ -244,9 +240,6 @@
 *
                         CALL CLQT01( M, N, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( .NOT.CGENND( M, N, AF, LDA ) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        NT = NT + 1
                      ELSE IF( M.LE.N ) THEN
 *
 *                       Test CUNGLQ, using factorization
@@ -254,9 +247,6 @@
 *
                         CALL CLQT02( M, N, K, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                     ELSE
-                        RESULT( 1 ) = ZERO
-                        RESULT( 2 ) = ZERO
                      END IF
                      IF( M.GE.K ) THEN
 *
@@ -299,14 +289,7 @@
      $                                  LDA, X, LDA, B, LDA, RWORK,
      $                                  RESULT( 7 ) )
                            NT = NT + 1
-                        ELSE
-                           RESULT( 7 ) = ZERO
                         END IF
-                     ELSE
-                        RESULT( 3 ) = ZERO
-                        RESULT( 4 ) = ZERO
-                        RESULT( 5 ) = ZERO
-                        RESULT( 6 ) = ZERO
                      END IF
 *
 *                    Print information about the tests that did not

diff --git a/TESTING/LIN/cchkql.f b/TESTING/LIN/cchkql.f
index fccb879..1e8253e 100644 (file)
--- a/TESTING/LIN/cchkql.f
+++ b/TESTING/LIN/cchkql.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       REAL               ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       REAL               RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            CGENND
-      EXTERNAL           CGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, CERRQL, CGEQLS, CGET02,
      $                   CLACPY, CLARHS, CLATB4, CLATMS, CQLT01, CQLT02,
@@ -244,15 +240,6 @@
 *
                         CALL CQLT01( M, N, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( M.GE.N ) THEN
-*                          Check the lower-left n-by-n corner
-                           IF( .NOT.CGENND(N, N, AF(M-N+1), LDA) )
-     $                          RESULT( 8 ) = 2*THRESH
-                        ELSE
-*                          Check the (n-m)th superdiagonal
-                           IF( .NOT.CGENND(M, M, AF(1+(N-M)*LDA), LDA) )
-     $                          RESULT( 8 ) = 2*THRESH
-                        ENDIF
                      ELSE IF( M.GE.N ) THEN
 *
 *                       Test CUNGQL, using factorization
@@ -260,9 +247,6 @@
 *
                         CALL CQLT02( M, N, K, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                     ELSE
-                        RESULT( 1 ) = ZERO
-                        RESULT( 2 ) = ZERO
                      END IF
                      IF( M.GE.K ) THEN
 *
@@ -305,14 +289,7 @@
      $                                  LDA, X( M-N+1 ), LDA, B, LDA,
      $                                  RWORK, RESULT( 7 ) )
                            NT = NT + 1
-                        ELSE
-                           RESULT( 7 ) = ZERO
                         END IF
-                     ELSE
-                        RESULT( 3 ) = ZERO
-                        RESULT( 4 ) = ZERO
-                        RESULT( 5 ) = ZERO
-                        RESULT( 6 ) = ZERO
                      END IF
 *
 *                    Print information about the tests that did not

diff --git a/TESTING/LIN/cchkqr.f b/TESTING/LIN/cchkqr.f
index f1bdeba..0a3f4ac 100644 (file)
--- a/TESTING/LIN/cchkqr.f
+++ b/TESTING/LIN/cchkqr.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       REAL               ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       REAL               RESULT( NTESTS )
 *     ..
-*     .. External Fuinctions ..
-      LOGICAL            CGENND
-      EXTERNAL           CGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, CERRQR, CGEQRS, CGET02,
      $                   CLACPY, CLARHS, CLATB4, CLATMS, CQRT01, CQRT02,
@@ -244,9 +240,6 @@
 *
                         CALL CQRT01( M, N, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( .NOT.CGENND( M, N, AF, LDA ) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        NT = NT + 1
                      ELSE IF( M.GE.N ) THEN
 *
 *                       Test CUNGQR, using factorization

diff --git a/TESTING/LIN/cchkrq.f b/TESTING/LIN/cchkrq.f
index b672d35..00349eb 100644 (file)
--- a/TESTING/LIN/cchkrq.f
+++ b/TESTING/LIN/cchkrq.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       REAL               ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       REAL               RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            CGENND
-      EXTERNAL           CGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, CERRRQ, CGERQS, CGET02,
      $                   CLACPY, CLARHS, CLATB4, CLATMS, CRQT01, CRQT02,
@@ -244,16 +240,6 @@
 *
                         CALL CRQT01( M, N, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( M.LE.N ) THEN
-*                          Check the upper-right m-by-m corner
-                           IF( .NOT.CGENND(M, M, AF(1+LDA*(N-M)), LDA) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        ELSE
-*                          Check the (m-n)th subdiagonal
-                           I = M - N
-                           IF( .NOT.CGENND(N, N, AF(I+1), LDA) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        END IF
                      ELSE IF( M.LE.N ) THEN
 *
 *                       Test CUNGRQ, using factorization
@@ -261,9 +247,6 @@
 *
                         CALL CRQT02( M, N, K, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                     ELSE
-                        RESULT( 1 ) = ZERO
-                        RESULT( 2 ) = ZERO
                      END IF
                      IF( M.GE.K ) THEN
 *
@@ -306,14 +289,7 @@
      $                                  LDA, X, LDA, B, LDA, RWORK,
      $                                  RESULT( 7 ) )
                            NT = NT + 1
-                        ELSE
-                           RESULT( 7 ) = ZERO
                         END IF
-                     ELSE
-                        RESULT( 3 ) = ZERO
-                        RESULT( 4 ) = ZERO
-                        RESULT( 5 ) = ZERO
-                        RESULT( 6 ) = ZERO
                      END IF
 *
 *                    Print information about the tests that did not

diff --git a/TESTING/LIN/dchklq.f b/TESTING/LIN/dchklq.f
index 350c7ce..5569725 100644 (file)
--- a/TESTING/LIN/dchklq.f
+++ b/TESTING/LIN/dchklq.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       DOUBLE PRECISION   ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       DOUBLE PRECISION   RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            DGENND
-      EXTERNAL           DGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, DERRLQ, DGELQS, DGET02,
      $                   DLACPY, DLARHS, DLATB4, DLATMS, DLQT01, DLQT02,
@@ -244,9 +240,6 @@
 *
                         CALL DLQT01( M, N, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( .NOT.DGENND( M, N, AF, LDA ) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        NT = NT + 1
                      ELSE IF( M.LE.N ) THEN
 *
 *                       Test DORGLQ, using factorization

diff --git a/TESTING/LIN/dchkql.f b/TESTING/LIN/dchkql.f
index 857bc2e..0e04723 100644 (file)
--- a/TESTING/LIN/dchkql.f
+++ b/TESTING/LIN/dchkql.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       DOUBLE PRECISION   ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       DOUBLE PRECISION   RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            DGENND
-      EXTERNAL           DGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, DERRQL, DGEQLS, DGET02,
      $                   DLACPY, DLARHS, DLATB4, DLATMS, DQLT01, DQLT02,
@@ -244,15 +240,6 @@
 *
                         CALL DQLT01( M, N, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( M.GE.N ) THEN
-*                          Check the lower-left n-by-n corner
-                           IF( .NOT.DGENND(N, N, AF(M-N+1), LDA) )
-     $                          RESULT( 8 ) = 2*THRESH
-                        ELSE
-*                          Check the (n-m)th superdiagonal
-                           IF( .NOT.DGENND(M, M, AF(1+(N-M)*LDA), LDA) )
-     $                          RESULT( 8 ) = 2*THRESH
-                        ENDIF
                      ELSE IF( M.GE.N ) THEN
 *
 *                       Test DORGQL, using factorization
@@ -260,9 +247,6 @@
 *
                         CALL DQLT02( M, N, K, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                     ELSE
-                        RESULT( 1 ) = ZERO
-                        RESULT( 2 ) = ZERO
                      END IF
                      IF( M.GE.K ) THEN
 *
@@ -305,14 +289,7 @@
      $                                  LDA, X( M-N+1 ), LDA, B, LDA,
      $                                  RWORK, RESULT( 7 ) )
                            NT = NT + 1
-                        ELSE
-                           RESULT( 7 ) = ZERO
                         END IF
-                     ELSE
-                        RESULT( 3 ) = ZERO
-                        RESULT( 4 ) = ZERO
-                        RESULT( 5 ) = ZERO
-                        RESULT( 6 ) = ZERO
                      END IF
 *
 *                    Print information about the tests that did not

diff --git a/TESTING/LIN/dchkqr.f b/TESTING/LIN/dchkqr.f
index d9a9b40..bab6359 100644 (file)
--- a/TESTING/LIN/dchkqr.f
+++ b/TESTING/LIN/dchkqr.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       DOUBLE PRECISION   ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       DOUBLE PRECISION   RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            DGENND
-      EXTERNAL           DGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, DERRQR, DGEQRS, DGET02,
      $                   DLACPY, DLARHS, DLATB4, DLATMS, DQRT01, DQRT02,
@@ -244,9 +240,6 @@
 *
                         CALL DQRT01( M, N, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( .NOT.DGENND( M, N, AF, LDA ) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        NT = NT + 1
                      ELSE IF( M.GE.N ) THEN
 *
 *                       Test DORGQR, using factorization

diff --git a/TESTING/LIN/dchkrq.f b/TESTING/LIN/dchkrq.f
index 6f23b41..3193ccc 100644 (file)
--- a/TESTING/LIN/dchkrq.f
+++ b/TESTING/LIN/dchkrq.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       DOUBLE PRECISION   ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       DOUBLE PRECISION   RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            DGENND
-      EXTERNAL           DGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, DERRRQ, DGERQS, DGET02,
      $                   DLACPY, DLARHS, DLATB4, DLATMS, DRQT01, DRQT02,
@@ -244,16 +240,6 @@
 *
                         CALL DRQT01( M, N, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( M.LE.N ) THEN
-*                          Check the upper-right m-by-m corner
-                           IF( .NOT.DGENND(M, M, AF(1+LDA*(N-M)), LDA) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        ELSE
-*                          Check the (m-n)th subdiagonal
-                           I = M - N
-                           IF( .NOT.DGENND(N, N, AF(I+1), LDA) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        END IF
                      ELSE IF( M.LE.N ) THEN
 *
 *                       Test DORGRQ, using factorization
@@ -261,9 +247,7 @@
 *
                         CALL DRQT02( M, N, K, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                     ELSE
-                        RESULT( 1 ) = ZERO
-                        RESULT( 2 ) = ZERO
+
                      END IF
                      IF( M.GE.K ) THEN
 *
@@ -306,14 +290,7 @@
      $                                  LDA, X, LDA, B, LDA, RWORK,
      $                                  RESULT( 7 ) )
                            NT = NT + 1
-                        ELSE
-                           RESULT( 7 ) = ZERO
                         END IF
-                     ELSE
-                        RESULT( 3 ) = ZERO
-                        RESULT( 4 ) = ZERO
-                        RESULT( 5 ) = ZERO
-                        RESULT( 6 ) = ZERO
                      END IF
 *
 *                    Print information about the tests that did not

diff --git a/TESTING/LIN/schklq.f b/TESTING/LIN/schklq.f
index 183a156..4867765 100644 (file)
--- a/TESTING/LIN/schklq.f
+++ b/TESTING/LIN/schklq.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       REAL               ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       REAL               RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            SGENND
-      EXTERNAL           SGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, SERRLQ, SGELQS, SGET02,
      $                   SLACPY, SLARHS, SLATB4, SLATMS, SLQT01, SLQT02,
@@ -244,9 +240,6 @@
 *
                         CALL SLQT01( M, N, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( .NOT.SGENND( M, N, AF, LDA ) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        NT = NT + 1
                      ELSE IF( M.LE.N ) THEN
 *
 *                       Test SORGLQ, using factorization

diff --git a/TESTING/LIN/schkql.f b/TESTING/LIN/schkql.f
index 9683eef..634ad87 100644 (file)
--- a/TESTING/LIN/schkql.f
+++ b/TESTING/LIN/schkql.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       REAL               ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       REAL               RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            SGENND
-      EXTERNAL           SGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, SERRQL, SGEQLS, SGET02,
      $                   SLACPY, SLARHS, SLATB4, SLATMS, SQLT01, SQLT02,
@@ -244,16 +240,6 @@
 *
                         CALL SQLT01( M, N, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( M.GE.N ) THEN
-*                          Check the lower-left n-by-n corner
-                           IF( .NOT.SGENND(N, N, AF(M-N+1), LDA) )
-     $                          RESULT( 8 ) = 2*THRESH
-                        ELSE
-*                          Check the (n-m)th superdiagonal
-                           IF( .NOT.SGENND(M, M, AF(1+(N-M)*LDA), LDA) )
-     $                          RESULT( 8 ) = 2*THRESH
-                        ENDIF
-                        NT = NT + 1
                      ELSE IF( M.GE.N ) THEN
 *
 *                       Test SORGQL, using factorization

diff --git a/TESTING/LIN/schkqr.f b/TESTING/LIN/schkqr.f
index cb40a13..c97ea39 100644 (file)
--- a/TESTING/LIN/schkqr.f
+++ b/TESTING/LIN/schkqr.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       REAL               ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       REAL               RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            SGENND
-      EXTERNAL           SGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, SERRQR, SGEQRS, SGET02,
      $                   SLACPY, SLARHS, SLATB4, SLATMS, SQRT01, SQRT02,
@@ -244,9 +240,6 @@
 *
                         CALL SQRT01( M, N, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( .NOT.SGENND( M, N, AF, LDA ) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        NT = NT + 1
                      ELSE IF( M.GE.N ) THEN
 *
 *                       Test SORGQR, using factorization

diff --git a/TESTING/LIN/schkrq.f b/TESTING/LIN/schkrq.f
index c996340..f9b7c50 100644 (file)
--- a/TESTING/LIN/schkrq.f
+++ b/TESTING/LIN/schkrq.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       REAL               ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       REAL               RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            SGENND
-      EXTERNAL           SGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, SERRRQ, SGERQS, SGET02,
      $                   SLACPY, SLARHS, SLATB4, SLATMS, SRQT01, SRQT02,
@@ -244,17 +240,6 @@
 *
                         CALL SRQT01( M, N, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( M.LE.N ) THEN
-*                          Check the upper-right m-by-m corner
-                           IF( .NOT.SGENND(M, M, AF(1+LDA*(N-M)), LDA) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        ELSE
-*                          Check the (m-n)th subdiagonal
-                           I = M - N
-                           IF( .NOT.SGENND(N, N, AF(I+1), LDA) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        END IF
-                        NT = NT + 1
                      ELSE IF( M.LE.N ) THEN
 *
 *                       Test SORGRQ, using factorization

diff --git a/TESTING/LIN/zchklq.f b/TESTING/LIN/zchklq.f
index 3fcf9c7..f72b4c0 100644 (file)
--- a/TESTING/LIN/zchklq.f
+++ b/TESTING/LIN/zchklq.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       DOUBLE PRECISION   ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       DOUBLE PRECISION   RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            ZGENND
-      EXTERNAL           ZGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, XLAENV, ZERRLQ, ZGELQS,
      $                   ZGET02, ZLACPY, ZLARHS, ZLATB4, ZLATMS, ZLQT01,
@@ -244,9 +240,6 @@
 *
                         CALL ZLQT01( M, N, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( .NOT.ZGENND( M, N, AF, LDA ) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        NT = NT + 1
                      ELSE IF( M.LE.N ) THEN
 *
 *                       Test ZUNGLQ, using factorization
@@ -254,9 +247,6 @@
 *
                         CALL ZLQT02( M, N, K, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                     ELSE
-                        RESULT( 1 ) = ZERO
-                        RESULT( 2 ) = ZERO
                      END IF
                      IF( M.GE.K ) THEN
 *
@@ -299,14 +289,7 @@
      $                                  LDA, X, LDA, B, LDA, RWORK,
      $                                  RESULT( 7 ) )
                            NT = NT + 1
-                        ELSE
-                           RESULT( 7 ) = ZERO
                         END IF
-                     ELSE
-                        RESULT( 3 ) = ZERO
-                        RESULT( 4 ) = ZERO
-                        RESULT( 5 ) = ZERO
-                        RESULT( 6 ) = ZERO
                      END IF
 *
 *                    Print information about the tests that did not

diff --git a/TESTING/LIN/zchkql.f b/TESTING/LIN/zchkql.f
index 8e08b1f..42861d4 100644 (file)
--- a/TESTING/LIN/zchkql.f
+++ b/TESTING/LIN/zchkql.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       DOUBLE PRECISION   ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       DOUBLE PRECISION   RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            ZGENND
-      EXTERNAL           ZGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, XLAENV, ZERRQL, ZGEQLS,
      $                   ZGET02, ZLACPY, ZLARHS, ZLATB4, ZLATMS, ZQLT01,
@@ -244,15 +240,6 @@
 *
                         CALL ZQLT01( M, N, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( M.GE.N ) THEN
-*                          Check the lower-left n-by-n corner
-                           IF( .NOT.ZGENND(N, N, AF(M-N+1), LDA) )
-     $                          RESULT( 8 ) = 2*THRESH
-                        ELSE
-*                          Check the (n-m)th superdiagonal
-                           IF( .NOT.ZGENND(M, M, AF(1+(N-M)*LDA), LDA) )
-     $                          RESULT( 8 ) = 2*THRESH
-                        ENDIF
                      ELSE IF( M.GE.N ) THEN
 *
 *                       Test ZUNGQL, using factorization
@@ -260,9 +247,6 @@
 *
                         CALL ZQLT02( M, N, K, A, AF, AQ, AL, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                     ELSE
-                        RESULT( 1 ) = ZERO
-                        RESULT( 2 ) = ZERO
                      END IF
                      IF( M.GE.K ) THEN
 *
@@ -305,14 +289,7 @@
      $                                  LDA, X( M-N+1 ), LDA, B, LDA,
      $                                  RWORK, RESULT( 7 ) )
                            NT = NT + 1
-                        ELSE
-                           RESULT( 7 ) = ZERO
                         END IF
-                     ELSE
-                        RESULT( 3 ) = ZERO
-                        RESULT( 4 ) = ZERO
-                        RESULT( 5 ) = ZERO
-                        RESULT( 6 ) = ZERO
                      END IF
 *
 *                    Print information about the tests that did not

diff --git a/TESTING/LIN/zchkqr.f b/TESTING/LIN/zchkqr.f
index 7dee29b..126b605 100644 (file)
--- a/TESTING/LIN/zchkqr.f
+++ b/TESTING/LIN/zchkqr.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       DOUBLE PRECISION   ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       DOUBLE PRECISION   RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            ZGENND
-      EXTERNAL           ZGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, XLAENV, ZERRQR, ZGEQRS,
      $                   ZGET02, ZLACPY, ZLARHS, ZLATB4, ZLATMS, ZQRT01,
@@ -244,9 +240,6 @@
 *
                         CALL ZQRT01( M, N, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( .NOT.ZGENND( M, N, AF, LDA ) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        NT = NT + 1
                      ELSE IF( M.GE.N ) THEN
 *
 *                       Test ZUNGQR, using factorization

diff --git a/TESTING/LIN/zchkrq.f b/TESTING/LIN/zchkrq.f
index 2916fa8..f9bd1f1 100644 (file)
--- a/TESTING/LIN/zchkrq.f
+++ b/TESTING/LIN/zchkrq.f
@@ -103,7 +103,7 @@
 *
 *     .. Parameters ..
       INTEGER            NTESTS
-      PARAMETER          ( NTESTS = 8 )
+      PARAMETER          ( NTESTS = 7 )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 8 )
       DOUBLE PRECISION   ZERO
@@ -121,10 +121,6 @@
       INTEGER            ISEED( 4 ), ISEEDY( 4 ), KVAL( 4 )
       DOUBLE PRECISION   RESULT( NTESTS )
 *     ..
-*     .. External Functions ..
-      LOGICAL            ZGENND
-      EXTERNAL           ZGENND
-*     ..
 *     .. External Subroutines ..
       EXTERNAL           ALAERH, ALAHD, ALASUM, XLAENV, ZERRRQ, ZGERQS,
      $                   ZGET02, ZLACPY, ZLARHS, ZLATB4, ZLATMS, ZRQT01,
@@ -244,16 +240,6 @@
 *
                         CALL ZRQT01( M, N, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                        IF( M.LE.N ) THEN
-*                          Check the upper-right m-by-m corner
-                           IF( .NOT.ZGENND(M, M, AF(1+LDA*(N-M)), LDA) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        ELSE
-*                          Check the (m-n)th subdiagonal
-                           I = M - N
-                           IF( .NOT.ZGENND(N, N, AF(I+1), LDA) )
-     $                       RESULT( 8 ) = 2*THRESH
-                        END IF
                      ELSE IF( M.LE.N ) THEN
 *
 *                       Test ZUNGRQ, using factorization
@@ -261,9 +247,6 @@
 *
                         CALL ZRQT02( M, N, K, A, AF, AQ, AR, LDA, TAU,
      $                               WORK, LWORK, RWORK, RESULT( 1 ) )
-                     ELSE
-                        RESULT( 1 ) = ZERO
-                        RESULT( 2 ) = ZERO
                      END IF
                      IF( M.GE.K ) THEN
 *
@@ -306,15 +289,8 @@
      $                                  LDA, X, LDA, B, LDA, RWORK,
      $                                  RESULT( 7 ) )
                            NT = NT + 1
-                        ELSE
-                           RESULT( 7 ) = ZERO
                         END IF
-                     ELSE
-                        RESULT( 3 ) = ZERO
-                        RESULT( 4 ) = ZERO
-                        RESULT( 5 ) = ZERO
-                        RESULT( 6 ) = ZERO
-                     END IF
+                    END IF
 *
 *                    Print information about the tests that did not
 *                    pass the threshold.