--- rpl/lapack/lapack/zgebrd.f 2010/12/21 13:53:43 1.7 +++ rpl/lapack/lapack/zgebrd.f 2011/07/22 07:38:13 1.8 @@ -1,10 +1,10 @@ SUBROUTINE ZGEBRD( M, N, A, LDA, D, E, TAUQ, TAUP, WORK, LWORK, $ 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 INFO, LDA, LWORK, M, N @@ -99,7 +99,7 @@ * * Each H(i) and G(i) has the form: * -* H(i) = I - tauq * v * v' and G(i) = I - taup * u * u' +* H(i) = I - tauq * v * v**H and G(i) = I - taup * u * u**H * * where tauq and taup are complex scalars, and v and u are complex * vectors; v(1:i-1) = 0, v(i) = 1, and v(i+1:m) is stored on exit in @@ -112,7 +112,7 @@ * * Each H(i) and G(i) has the form: * -* H(i) = I - tauq * v * v' and G(i) = I - taup * u * u' +* H(i) = I - tauq * v * v**H and G(i) = I - taup * u * u**H * * where tauq and taup are complex scalars, and v and u are complex * vectors; v(1:i) = 0, v(i+1) = 1, and v(i+2:m) is stored on exit in @@ -232,7 +232,7 @@ $ WORK( LDWRKX*NB+1 ), LDWRKY ) * * Update the trailing submatrix A(i+ib:m,i+ib:n), using -* an update of the form A := A - V*Y' - X*U' +* an update of the form A := A - V*Y**H - X*U**H * CALL ZGEMM( 'No transpose', 'Conjugate transpose', M-I-NB+1, $ N-I-NB+1, NB, -ONE, A( I+NB, I ), LDA,