--- 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))