*tplqt.f, *tplqt2.f: incorrect dimensions

8 files changed, 52 insertions, 52 deletions

diff --git a/SRC/ctplqt.f b/SRC/ctplqt.f
index 322b6dc5..4aa96499 100644
--- a/SRC/ctplqt.f
+++ b/SRC/ctplqt.f
@@ -56,8 +56,8 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is COMPLEX array, dimension (LDA,N)
-*>          On entry, the lower triangular N-by-N matrix A.
+*>          A is COMPLEX array, dimension (LDA,M)
+*>          On entry, the lower triangular M-by-M matrix A.
 *>          On exit, the elements on and below the diagonal of the array
 *>          contain the lower triangular matrix L.
 *> \endverbatim
@@ -65,7 +65,7 @@
 *> \param[in] LDA
 *> \verbatim
 *>          LDA is INTEGER
-*>          The leading dimension of the array A.  LDA >= max(1,N).
+*>          The leading dimension of the array A.  LDA >= max(1,M).
 *> \endverbatim
 *>
 *> \param[in,out] B
@@ -129,26 +129,26 @@
 *>               C = [ A ] [ B ]
 *>
 *>
-*>  where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal
+*>  where A is an lower triangular M-by-M matrix, and B is M-by-N pentagonal
 *>  matrix consisting of a M-by-(N-L) rectangular matrix B1 on left of a M-by-L
 *>  upper trapezoidal matrix B2:
 *>          [ B ] = [ B1 ] [ B2 ]
 *>                   [ B1 ]  <- M-by-(N-L) rectangular
-*>                   [ B2 ]  <-     M-by-L upper trapezoidal.
+*>                   [ B2 ]  <-     M-by-L lower trapezoidal.
 *>
 *>  The lower trapezoidal matrix B2 consists of the first L columns of a
-*>  N-by-N lower triangular matrix, where 0 <= L <= MIN(M,N).  If L=0,
+*>  M-by-M lower triangular matrix, where 0 <= L <= MIN(M,N).  If L=0,
 *>  B is rectangular M-by-N; if M=L=N, B is lower triangular.
 *>
 *>  The matrix W stores the elementary reflectors H(i) in the i-th row
 *>  above the diagonal (of A) in the M-by-(M+N) input matrix C
 *>            [ C ] = [ A ] [ B ]
-*>                   [ A ]  <- lower triangular N-by-N
+*>                   [ A ]  <- lower triangular M-by-M
 *>                   [ B ]  <- M-by-N pentagonal
 *>
 *>  so that W can be represented as
 *>            [ W ] = [ I ] [ V ]
-*>                   [ I ]  <- identity, N-by-N
+*>                   [ I ]  <- identity, M-by-M
 *>                   [ V ]  <- M-by-N, same form as B.
 *>
 *>  Thus, all of information needed for W is contained on exit in B, which
diff --git a/SRC/ctplqt2.f b/SRC/ctplqt2.f
index 1c9b128e..d1e28849 100644
--- a/SRC/ctplqt2.f
+++ b/SRC/ctplqt2.f
@@ -48,7 +48,7 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is COMPLEX array, dimension (LDA,N)
+*>          A is COMPLEX array, dimension (LDA,M)
 *>          On entry, the lower triangular M-by-M matrix A.
 *>          On exit, the elements on and below the diagonal of the array
 *>          contain the lower triangular matrix L.
@@ -57,7 +57,7 @@
 *> \param[in] LDA
 *> \verbatim
 *>          LDA is INTEGER
-*>          The leading dimension of the array A.  LDA >= max(1,N).
+*>          The leading dimension of the array A.  LDA >= max(1,M).
 *> \endverbatim
 *>
 *> \param[in,out] B
@@ -116,7 +116,7 @@
 *>               C = [ A ][ B ]
 *>
 *>
-*>  where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal
+*>  where A is an lower triangular M-by-M matrix, and B is M-by-N pentagonal
 *>  matrix consisting of a M-by-(N-L) rectangular matrix B1 left of a M-by-L
 *>  upper trapezoidal matrix B2:
 *>
@@ -132,13 +132,13 @@
 *>  above the diagonal (of A) in the M-by-(M+N) input matrix C
 *>
 *>               C = [ A ][ B ]
-*>                   [ A ]  <- lower triangular N-by-N
+*>                   [ A ]  <- lower triangular M-by-M
 *>                   [ B ]  <- M-by-N pentagonal
 *>
 *>  so that W can be represented as
 *>
 *>               W = [ I ][ V ]
-*>                   [ I ]  <- identity, N-by-N
+*>                   [ I ]  <- identity, M-by-M
 *>                   [ V ]  <- M-by-N, same form as B.
 *>
 *>  Thus, all of information needed for W is contained on exit in B, which
diff --git a/SRC/dtplqt.f b/SRC/dtplqt.f
index b312c501..a5a92744 100644
--- a/SRC/dtplqt.f
+++ b/SRC/dtplqt.f
@@ -73,8 +73,8 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is DOUBLE PRECISION array, dimension (LDA,N)
-*>          On entry, the lower triangular N-by-N matrix A.
+*>          A is DOUBLE PRECISION array, dimension (LDA,M)
+*>          On entry, the lower triangular M-by-M matrix A.
 *>          On exit, the elements on and below the diagonal of the array
 *>          contain the lower triangular matrix L.
 *> \endverbatim
@@ -82,7 +82,7 @@
 *> \param[in] LDA
 *> \verbatim
 *>          LDA is INTEGER
-*>          The leading dimension of the array A.  LDA >= max(1,N).
+*>          The leading dimension of the array A.  LDA >= max(1,M).
 *> \endverbatim
 *>
 *> \param[in,out] B
@@ -146,26 +146,26 @@
 *>               C = [ A ] [ B ]
 *>
 *>
-*>  where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal
+*>  where A is an lower triangular M-by-M matrix, and B is M-by-N pentagonal
 *>  matrix consisting of a M-by-(N-L) rectangular matrix B1 on left of a M-by-L
 *>  upper trapezoidal matrix B2:
 *>          [ B ] = [ B1 ] [ B2 ]
 *>                   [ B1 ]  <- M-by-(N-L) rectangular
-*>                   [ B2 ]  <-     M-by-L upper trapezoidal.
+*>                   [ B2 ]  <-     M-by-L lower trapezoidal.
 *>
 *>  The lower trapezoidal matrix B2 consists of the first L columns of a
-*>  N-by-N lower triangular matrix, where 0 <= L <= MIN(M,N).  If L=0,
+*>  M-by-M lower triangular matrix, where 0 <= L <= MIN(M,N).  If L=0,
 *>  B is rectangular M-by-N; if M=L=N, B is lower triangular.
 *>
 *>  The matrix W stores the elementary reflectors H(i) in the i-th row
 *>  above the diagonal (of A) in the M-by-(M+N) input matrix C
 *>            [ C ] = [ A ] [ B ]
-*>                   [ A ]  <- lower triangular N-by-N
+*>                   [ A ]  <- lower triangular M-by-M
 *>                   [ B ]  <- M-by-N pentagonal
 *>
 *>  so that W can be represented as
 *>            [ W ] = [ I ] [ V ]
-*>                   [ I ]  <- identity, N-by-N
+*>                   [ I ]  <- identity, M-by-M
 *>                   [ V ]  <- M-by-N, same form as B.
 *>
 *>  Thus, all of information needed for W is contained on exit in B, which
diff --git a/SRC/dtplqt2.f b/SRC/dtplqt2.f
index 7e87e6c5..beb5ad6b 100644
--- a/SRC/dtplqt2.f
+++ b/SRC/dtplqt2.f
@@ -65,7 +65,7 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is DOUBLE PRECISION array, dimension (LDA,N)
+*>          A is DOUBLE PRECISION array, dimension (LDA,M)
 *>          On entry, the lower triangular M-by-M matrix A.
 *>          On exit, the elements on and below the diagonal of the array
 *>          contain the lower triangular matrix L.
@@ -74,7 +74,7 @@
 *> \param[in] LDA
 *> \verbatim
 *>          LDA is INTEGER
-*>          The leading dimension of the array A.  LDA >= max(1,N).
+*>          The leading dimension of the array A.  LDA >= max(1,M).
 *> \endverbatim
 *>
 *> \param[in,out] B
@@ -133,7 +133,7 @@
 *>               C = [ A ][ B ]
 *>
 *>
-*>  where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal
+*>  where A is an lower triangular M-by-M matrix, and B is M-by-N pentagonal
 *>  matrix consisting of a M-by-(N-L) rectangular matrix B1 left of a M-by-L
 *>  upper trapezoidal matrix B2:
 *>
@@ -149,13 +149,13 @@
 *>  above the diagonal (of A) in the M-by-(M+N) input matrix C
 *>
 *>               C = [ A ][ B ]
-*>                   [ A ]  <- lower triangular N-by-N
+*>                   [ A ]  <- lower triangular M-by-M
 *>                   [ B ]  <- M-by-N pentagonal
 *>
 *>  so that W can be represented as
 *>
 *>               W = [ I ][ V ]
-*>                   [ I ]  <- identity, N-by-N
+*>                   [ I ]  <- identity, M-by-M
 *>                   [ V ]  <- M-by-N, same form as B.
 *>
 *>  Thus, all of information needed for W is contained on exit in B, which
diff --git a/SRC/stplqt.f b/SRC/stplqt.f
index e3c37abf..482ce568 100644
--- a/SRC/stplqt.f
+++ b/SRC/stplqt.f
@@ -73,8 +73,8 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is REAL array, dimension (LDA,N)
-*>          On entry, the lower triangular N-by-N matrix A.
+*>          A is REAL array, dimension (LDA,M)
+*>          On entry, the lower triangular M-by-M matrix A.
 *>          On exit, the elements on and below the diagonal of the array
 *>          contain the lower triangular matrix L.
 *> \endverbatim
@@ -82,7 +82,7 @@
 *> \param[in] LDA
 *> \verbatim
 *>          LDA is INTEGER
-*>          The leading dimension of the array A.  LDA >= max(1,N).
+*>          The leading dimension of the array A.  LDA >= max(1,M).
 *> \endverbatim
 *>
 *> \param[in,out] B
@@ -146,26 +146,26 @@
 *>               C = [ A ] [ B ]
 *>
 *>
-*>  where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal
+*>  where A is an lower triangular M-by-M matrix, and B is M-by-N pentagonal
 *>  matrix consisting of a M-by-(N-L) rectangular matrix B1 on left of a M-by-L
 *>  upper trapezoidal matrix B2:
 *>          [ B ] = [ B1 ] [ B2 ]
 *>                   [ B1 ]  <- M-by-(N-L) rectangular
-*>                   [ B2 ]  <-     M-by-L upper trapezoidal.
+*>                   [ B2 ]  <-     M-by-L lower trapezoidal.
 *>
 *>  The lower trapezoidal matrix B2 consists of the first L columns of a
-*>  N-by-N lower triangular matrix, where 0 <= L <= MIN(M,N).  If L=0,
+*>  M-by-M lower triangular matrix, where 0 <= L <= MIN(M,N).  If L=0,
 *>  B is rectangular M-by-N; if M=L=N, B is lower triangular.
 *>
 *>  The matrix W stores the elementary reflectors H(i) in the i-th row
 *>  above the diagonal (of A) in the M-by-(M+N) input matrix C
 *>            [ C ] = [ A ] [ B ]
-*>                   [ A ]  <- lower triangular N-by-N
+*>                   [ A ]  <- lower triangular M-by-M
 *>                   [ B ]  <- M-by-N pentagonal
 *>
 *>  so that W can be represented as
 *>            [ W ] = [ I ] [ V ]
-*>                   [ I ]  <- identity, N-by-N
+*>                   [ I ]  <- identity, M-by-M
 *>                   [ V ]  <- M-by-N, same form as B.
 *>
 *>  Thus, all of information needed for W is contained on exit in B, which
diff --git a/SRC/stplqt2.f b/SRC/stplqt2.f
index f1b8e030..b16738fb 100644
--- a/SRC/stplqt2.f
+++ b/SRC/stplqt2.f
@@ -65,7 +65,7 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is REAL array, dimension (LDA,N)
+*>          A is REAL array, dimension (LDA,M)
 *>          On entry, the lower triangular M-by-M matrix A.
 *>          On exit, the elements on and below the diagonal of the array
 *>          contain the lower triangular matrix L.
@@ -74,7 +74,7 @@
 *> \param[in] LDA
 *> \verbatim
 *>          LDA is INTEGER
-*>          The leading dimension of the array A.  LDA >= max(1,N).
+*>          The leading dimension of the array A.  LDA >= max(1,M).
 *> \endverbatim
 *>
 *> \param[in,out] B
@@ -133,7 +133,7 @@
 *>               C = [ A ][ B ]
 *>
 *>
-*>  where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal
+*>  where A is an lower triangular M-by-M matrix, and B is M-by-N pentagonal
 *>  matrix consisting of a M-by-(N-L) rectangular matrix B1 left of a M-by-L
 *>  upper trapezoidal matrix B2:
 *>
@@ -149,13 +149,13 @@
 *>  above the diagonal (of A) in the M-by-(M+N) input matrix C
 *>
 *>               C = [ A ][ B ]
-*>                   [ A ]  <- lower triangular N-by-N
+*>                   [ A ]  <- lower triangular M-by-M
 *>                   [ B ]  <- M-by-N pentagonal
 *>
 *>  so that W can be represented as
 *>
 *>               W = [ I ][ V ]
-*>                   [ I ]  <- identity, N-by-N
+*>                   [ I ]  <- identity, M-by-M
 *>                   [ V ]  <- M-by-N, same form as B.
 *>
 *>  Thus, all of information needed for W is contained on exit in B, which
diff --git a/SRC/ztplqt.f b/SRC/ztplqt.f
index 28740208..a9a035ef 100644
--- a/SRC/ztplqt.f
+++ b/SRC/ztplqt.f
@@ -73,8 +73,8 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is COMPLEX*16 array, dimension (LDA,N)
-*>          On entry, the lower triangular N-by-N matrix A.
+*>          A is COMPLEX*16 array, dimension (LDA,M)
+*>          On entry, the lower triangular M-by-M matrix A.
 *>          On exit, the elements on and below the diagonal of the array
 *>          contain the lower triangular matrix L.
 *> \endverbatim
@@ -82,7 +82,7 @@
 *> \param[in] LDA
 *> \verbatim
 *>          LDA is INTEGER
-*>          The leading dimension of the array A.  LDA >= max(1,N).
+*>          The leading dimension of the array A.  LDA >= max(1,M).
 *> \endverbatim
 *>
 *> \param[in,out] B
@@ -146,26 +146,26 @@
 *>               C = [ A ] [ B ]
 *>
 *>
-*>  where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal
+*>  where A is an lower triangular M-by-M matrix, and B is M-by-N pentagonal
 *>  matrix consisting of a M-by-(N-L) rectangular matrix B1 on left of a M-by-L
 *>  upper trapezoidal matrix B2:
 *>          [ B ] = [ B1 ] [ B2 ]
 *>                   [ B1 ]  <- M-by-(N-L) rectangular
-*>                   [ B2 ]  <-     M-by-L upper trapezoidal.
+*>                   [ B2 ]  <-     M-by-L lower trapezoidal.
 *>
 *>  The lower trapezoidal matrix B2 consists of the first L columns of a
-*>  N-by-N lower triangular matrix, where 0 <= L <= MIN(M,N).  If L=0,
+*>  M-by-M lower triangular matrix, where 0 <= L <= MIN(M,N).  If L=0,
 *>  B is rectangular M-by-N; if M=L=N, B is lower triangular.
 *>
 *>  The matrix W stores the elementary reflectors H(i) in the i-th row
 *>  above the diagonal (of A) in the M-by-(M+N) input matrix C
 *>            [ C ] = [ A ] [ B ]
-*>                   [ A ]  <- lower triangular N-by-N
+*>                   [ A ]  <- lower triangular M-by-M
 *>                   [ B ]  <- M-by-N pentagonal
 *>
 *>  so that W can be represented as
 *>            [ W ] = [ I ] [ V ]
-*>                   [ I ]  <- identity, N-by-N
+*>                   [ I ]  <- identity, M-by-M
 *>                   [ V ]  <- M-by-N, same form as B.
 *>
 *>  Thus, all of information needed for W is contained on exit in B, which
diff --git a/SRC/ztplqt2.f b/SRC/ztplqt2.f
index 733f9dcc..a1d96dba 100644
--- a/SRC/ztplqt2.f
+++ b/SRC/ztplqt2.f
@@ -65,7 +65,7 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is COMPLEX*16 array, dimension (LDA,N)
+*>          A is COMPLEX*16 array, dimension (LDA,M)
 *>          On entry, the lower triangular M-by-M matrix A.
 *>          On exit, the elements on and below the diagonal of the array
 *>          contain the lower triangular matrix L.
@@ -74,7 +74,7 @@
 *> \param[in] LDA
 *> \verbatim
 *>          LDA is INTEGER
-*>          The leading dimension of the array A.  LDA >= max(1,N).
+*>          The leading dimension of the array A.  LDA >= max(1,M).
 *> \endverbatim
 *>
 *> \param[in,out] B
@@ -133,7 +133,7 @@
 *>               C = [ A ][ B ]
 *>
 *>
-*>  where A is an lower triangular N-by-N matrix, and B is M-by-N pentagonal
+*>  where A is an lower triangular M-by-M matrix, and B is M-by-N pentagonal
 *>  matrix consisting of a M-by-(N-L) rectangular matrix B1 left of a M-by-L
 *>  upper trapezoidal matrix B2:
 *>
@@ -149,13 +149,13 @@
 *>  above the diagonal (of A) in the M-by-(M+N) input matrix C
 *>
 *>               C = [ A ][ B ]
-*>                   [ A ]  <- lower triangular N-by-N
+*>                   [ A ]  <- lower triangular M-by-M
 *>                   [ B ]  <- M-by-N pentagonal
 *>
 *>  so that W can be represented as
 *>
 *>               W = [ I ][ V ]
-*>                   [ I ]  <- identity, N-by-N
+*>                   [ I ]  <- identity, M-by-M
 *>                   [ V ]  <- M-by-N, same form as B.
 *>
 *>  Thus, all of information needed for W is contained on exit in B, which