--- rpl/lapack/lapack/zgebd2.f 2018/05/29 07:18:13 1.19 +++ rpl/lapack/lapack/zgebd2.f 2023/08/07 08:39:16 1.20 @@ -132,8 +132,6 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date June 2017 -* *> \ingroup complex16GEcomputational * *> \par Further Details: @@ -189,10 +187,9 @@ * ===================================================================== SUBROUTINE ZGEBD2( M, N, A, LDA, D, E, TAUQ, TAUP, WORK, INFO ) * -* -- LAPACK computational routine (version 3.7.1) -- +* -- LAPACK computational routine -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* June 2017 * * .. Scalar Arguments .. INTEGER INFO, LDA, M, N @@ -247,7 +244,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 ) A( I, I ) = ONE * * Apply H(i)**H to A(i:m,i+1:n) from the left @@ -266,7 +263,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 * * Apply G(i) to A(i+1:m,i+1:n) from the right @@ -291,7 +288,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 ) A( I, I ) = ONE * * Apply G(i) to A(i+1:m,i:n) from the right @@ -310,7 +307,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 * * Apply H(i)**H to A(i+1:m,i+1:n) from the left