From 9aeb8a17f0752bdec5559217709e2966e50fbaf5 Mon Sep 17 00:00:00 2001 From: "eugene.chereshnev" Date: Thu, 16 Mar 2017 11:25:21 -0700 Subject: *tplqt.f, *tplqt2.f: incorrect dimensions --- SRC/ctplqt2.f | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) (limited to 'SRC/ctplqt2.f') 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 -- cgit v1.2.3