--- rpl/lapack/lapack/zlabrd.f 2012/12/14 14:22:49 1.13
+++ rpl/lapack/lapack/zlabrd.f 2023/08/07 08:39:28 1.20
@@ -2,25 +2,25 @@
*
* =========== DOCUMENTATION ===========
*
-* Online html documentation available at
-* http://www.netlib.org/lapack/explore-html/
+* Online html documentation available at
+* http://www.netlib.org/lapack/explore-html/
*
*> \htmlonly
-*> Download ZLABRD + dependencies
-*>
-*> [TGZ]
-*>
-*> [ZIP]
-*>
+*> Download ZLABRD + dependencies
+*>
+*> [TGZ]
+*>
+*> [ZIP]
+*>
*> [TXT]
-*> \endhtmlonly
+*> \endhtmlonly
*
* Definition:
* ===========
*
* SUBROUTINE ZLABRD( M, N, NB, A, LDA, D, E, TAUQ, TAUP, X, LDX, Y,
* LDY )
-*
+*
* .. Scalar Arguments ..
* INTEGER LDA, LDX, LDY, M, N, NB
* ..
@@ -29,7 +29,7 @@
* COMPLEX*16 A( LDA, * ), TAUP( * ), TAUQ( * ), X( LDX, * ),
* $ Y( LDY, * )
* ..
-*
+*
*
*> \par Purpose:
* =============
@@ -111,7 +111,7 @@
*>
*> \param[out] TAUQ
*> \verbatim
-*> TAUQ is COMPLEX*16 array dimension (NB)
+*> TAUQ is COMPLEX*16 array, dimension (NB)
*> The scalar factors of the elementary reflectors which
*> represent the unitary matrix Q. See Further Details.
*> \endverbatim
@@ -152,12 +152,10 @@
* Authors:
* ========
*
-*> \author Univ. of Tennessee
-*> \author Univ. of California Berkeley
-*> \author Univ. of Colorado Denver
-*> \author NAG Ltd.
-*
-*> \date September 2012
+*> \author Univ. of Tennessee
+*> \author Univ. of California Berkeley
+*> \author Univ. of Colorado Denver
+*> \author NAG Ltd.
*
*> \ingroup complex16OTHERauxiliary
*
@@ -212,10 +210,9 @@
SUBROUTINE ZLABRD( M, N, NB, A, LDA, D, E, TAUQ, TAUP, X, LDX, Y,
$ LDY )
*
-* -- LAPACK auxiliary routine (version 3.4.2) --
+* -- LAPACK auxiliary routine --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* September 2012
*
* .. Scalar Arguments ..
INTEGER LDA, LDX, LDY, M, N, NB
@@ -270,7 +267,7 @@
ALPHA = A( I, I )
CALL ZLARFG( M-I+1, ALPHA, A( MIN( I+1, M ), I ), 1,
$ TAUQ( I ) )
- D( I ) = ALPHA
+ D( I ) = DBLE( ALPHA )
IF( I.LT.N ) THEN
A( I, I ) = ONE
*
@@ -310,7 +307,7 @@
ALPHA = A( I, I+1 )
CALL ZLARFG( N-I, ALPHA, A( I, MIN( I+2, N ) ), LDA,
$ TAUP( I ) )
- E( I ) = ALPHA
+ E( I ) = DBLE( ALPHA )
A( I, I+1 ) = ONE
*
* Compute X(i+1:m,i)
@@ -354,7 +351,7 @@
ALPHA = A( I, I )
CALL ZLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
$ TAUP( I ) )
- D( I ) = ALPHA
+ D( I ) = DBLE( ALPHA )
IF( I.LT.M ) THEN
A( I, I ) = ONE
*
@@ -388,7 +385,7 @@
ALPHA = A( I+1, I )
CALL ZLARFG( M-I, ALPHA, A( MIN( I+2, M ), I ), 1,
$ TAUQ( I ) )
- E( I ) = ALPHA
+ E( I ) = DBLE( ALPHA )
A( I+1, I ) = ONE
*
* Compute Y(i+1:n,i)