--- rpl/lapack/lapack/zgehd2.f 2010/12/21 13:53:43 1.7 +++ rpl/lapack/lapack/zgehd2.f 2011/07/22 07:38:14 1.8 @@ -1,9 +1,9 @@ SUBROUTINE ZGEHD2( N, ILO, IHI, A, LDA, TAU, WORK, INFO ) * -* -- LAPACK routine (version 3.2) -- +* -- LAPACK routine (version 3.3.1) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* November 2006 +* -- April 2011 -- * * .. Scalar Arguments .. INTEGER IHI, ILO, INFO, LDA, N @@ -16,7 +16,7 @@ * ======= * * ZGEHD2 reduces a complex general matrix A to upper Hessenberg form H -* by a unitary similarity transformation: Q' * A * Q = H . +* by a unitary similarity transformation: Q**H * A * Q = H . * * Arguments * ========= @@ -63,7 +63,7 @@ * * Each H(i) has the form * -* H(i) = I - tau * v * v' +* H(i) = I - tau * v * v**H * * where tau is a complex scalar, and v is a complex vector with * v(1:i) = 0, v(i+1) = 1 and v(ihi+1:n) = 0; v(i+2:ihi) is stored on @@ -134,7 +134,7 @@ CALL ZLARF( 'Right', IHI, IHI-I, A( I+1, I ), 1, TAU( I ), $ A( 1, I+1 ), LDA, WORK ) * -* Apply H(i)' to A(i+1:ihi,i+1:n) from the left +* Apply H(i)**H to A(i+1:ihi,i+1:n) from the left * CALL ZLARF( 'Left', IHI-I, N-I, A( I+1, I ), 1, $ DCONJG( TAU( I ) ), A( I+1, I+1 ), LDA, WORK )