--- rpl/lapack/lapack/zhetrd_he2hb.f 2017/06/17 11:06:49 1.2
+++ rpl/lapack/lapack/zhetrd_he2hb.f 2023/08/07 08:39:24 1.6
@@ -9,11 +9,11 @@
*
*> \htmlonly
*> Download ZHETRD_HE2HB + dependencies
-*>
+*>
*> [TGZ]
-*>
+*>
*> [ZIP]
-*>
+*>
*> [TXT]
*> \endhtmlonly
*
@@ -123,7 +123,7 @@
*>
*> \param[out] WORK
*> \verbatim
-*> WORK is COMPLEX*16 array, dimension LWORK.
+*> WORK is COMPLEX*16 array, dimension (LWORK)
*> On exit, if INFO = 0, or if LWORK=-1,
*> WORK(1) returns the size of LWORK.
*> \endverbatim
@@ -132,7 +132,7 @@
*> \verbatim
*> LWORK is INTEGER
*> The dimension of the array WORK which should be calculated
-* by a workspace query. LWORK = MAX(1, LWORK_QUERY)
+*> by a workspace query. LWORK = MAX(1, LWORK_QUERY)
*> If LWORK = -1, then a workspace query is assumed; the routine
*> only calculates the optimal size of the WORK array, returns
*> this value as the first entry of the WORK array, and no error
@@ -158,8 +158,6 @@
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date December 2016
-*
*> \ingroup complex16HEcomputational
*
*> \par Further Details:
@@ -222,7 +220,7 @@
*>
*> where tau is a complex scalar, and v is a complex vector with
*> v(kd+1:i) = 0 and v(i+kd+1) = 1; v(i+kd+2:n) is stored on exit in
-* A(i+kd+2:n,i), and tau in TAU(i).
+*> A(i+kd+2:n,i), and tau in TAU(i).
*>
*> The contents of A on exit are illustrated by the following examples
*> with n = 5:
@@ -245,10 +243,9 @@
*
IMPLICIT NONE
*
-* -- 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 ..
CHARACTER UPLO
@@ -277,7 +274,7 @@
$ TPOS, WPOS, S2POS, S1POS
* ..
* .. External Subroutines ..
- EXTERNAL XERBLA, ZHER2K, ZHEMM, ZGEMM,
+ EXTERNAL XERBLA, ZHER2K, ZHEMM, ZGEMM, ZCOPY,
$ ZLARFT, ZGELQF, ZGEQRF, ZLASET
* ..
* .. Intrinsic Functions ..
@@ -285,8 +282,8 @@
* ..
* .. External Functions ..
LOGICAL LSAME
- INTEGER ILAENV
- EXTERNAL LSAME, ILAENV
+ INTEGER ILAENV2STAGE
+ EXTERNAL LSAME, ILAENV2STAGE
* ..
* .. Executable Statements ..
*
@@ -296,7 +293,7 @@
INFO = 0
UPPER = LSAME( UPLO, 'U' )
LQUERY = ( LWORK.EQ.-1 )
- LWMIN = ILAENV( 20, 'ZHETRD_HE2HB', '', N, KD, -1, -1 )
+ LWMIN = ILAENV2STAGE( 4, 'ZHETRD_HE2HB', ' ', N, KD, -1, -1 )
IF( .NOT.UPPER .AND. .NOT.LSAME( UPLO, 'L' ) ) THEN
INFO = -1
@@ -363,7 +360,7 @@
*
*
* Set the workspace of the triangular matrix T to zero once such a
-* way everytime T is generated the upper/lower portion will be always zero
+* way every time T is generated the upper/lower portion will be always zero
*
CALL ZLASET( "A", LDT, KD, ZERO, ZERO, WORK( TPOS ), LDT )
*