--- rpl/lapack/lapack/zgelqt3.f 2017/06/17 11:02:52 1.1 +++ rpl/lapack/lapack/zgelqt3.f 2023/08/07 08:39:17 1.5 @@ -6,7 +6,7 @@ * http://www.netlib.org/lapack/explore-html/ * *> \htmlonly -*> Download DGEQRT3 + dependencies +*> Download ZGELQT3 + dependencies *> *> [TGZ] *> @@ -33,7 +33,7 @@ *> *> \verbatim *> -*> DGELQT3 recursively computes a LQ factorization of a complex M-by-N +*> ZGELQT3 recursively computes a LQ factorization of a complex M-by-N *> matrix A, using the compact WY representation of Q. *> *> Based on the algorithm of Elmroth and Gustavson, @@ -58,7 +58,7 @@ *> \param[in,out] A *> \verbatim *> A is COMPLEX*16 array, dimension (LDA,N) -*> On entry, the real M-by-N matrix A. On exit, the elements on and +*> On entry, the complex M-by-N matrix A. On exit, the elements on and *> below the diagonal contain the N-by-N lower triangular matrix L; the *> elements above the diagonal are the rows of V. See below for *> further details. @@ -100,8 +100,6 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date December 2016 -* *> \ingroup doubleGEcomputational * *> \par Further Details: @@ -109,8 +107,8 @@ *> *> \verbatim *> -*> The matrix V stores the elementary reflectors H(i) in the i-th column -*> below the diagonal. For example, if M=5 and N=3, the matrix V is +*> The matrix V stores the elementary reflectors H(i) in the i-th row +*> above the diagonal. For example, if M=5 and N=3, the matrix V is *> *> V = ( 1 v1 v1 v1 v1 ) *> ( 1 v2 v2 v2 ) @@ -131,10 +129,9 @@ * ===================================================================== RECURSIVE SUBROUTINE ZGELQT3( M, N, A, LDA, T, LDT, INFO ) * -* -- LAPACK computational routine (version 3.7.0) -- +* -- LAPACK computational routine -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* December 2016 * * .. Scalar Arguments .. INTEGER INFO, LDA, M, N, LDT @@ -151,7 +148,7 @@ PARAMETER ( ZERO = (0.0D+00,0.0D+00)) * .. * .. Local Scalars .. - INTEGER I, I1, J, J1, M1, M2, N1, N2, IINFO + INTEGER I, I1, J, J1, M1, M2, IINFO * .. * .. External Subroutines .. EXTERNAL ZLARFG, ZTRMM, ZGEMM, XERBLA @@ -175,7 +172,7 @@ * IF( M.EQ.1 ) THEN * -* Compute Householder transform when N=1 +* Compute Householder transform when M=1 * CALL ZLARFG( N, A, A( 1, MIN( 2, N ) ), LDA, T ) T(1,1)=CONJG(T(1,1))