--- rpl/lapack/lapack/ztzrqf.f 2010/12/21 13:53:58 1.8 +++ rpl/lapack/lapack/ztzrqf.f 2011/07/22 07:38:21 1.9 @@ -1,9 +1,9 @@ SUBROUTINE ZTZRQF( M, N, A, LDA, TAU, INFO ) * -* -- LAPACK routine (version 3.2.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..-- -* June 2010 +* -- April 2011 -- * * .. Scalar Arguments .. INTEGER INFO, LDA, M, N @@ -66,9 +66,9 @@ * * where * -* T( k ) = I - tau*u( k )*u( k )', u( k ) = ( 1 ), -* ( 0 ) -* ( z( k ) ) +* T( k ) = I - tau*u( k )*u( k )**H, u( k ) = ( 1 ), +* ( 0 ) +* ( z( k ) ) * * tau is a scalar and z( k ) is an ( n - m ) element vector. * tau and z( k ) are chosen to annihilate the elements of the kth row @@ -142,7 +142,7 @@ * IF( TAU( K ).NE.CZERO .AND. K.GT.1 ) THEN * -* We now perform the operation A := A*P( k )'. +* We now perform the operation A := A*P( k )**H. * * Use the first ( k - 1 ) elements of TAU to store a( k ), * where a( k ) consists of the first ( k - 1 ) elements of @@ -157,7 +157,7 @@ $ LDA, A( K, M1 ), LDA, CONE, TAU, 1 ) * * Now form a( k ) := a( k ) - conjg(tau)*w -* and B := B - conjg(tau)*w*z( k )'. +* and B := B - conjg(tau)*w*z( k )**H. * CALL ZAXPY( K-1, -DCONJG( TAU( K ) ), TAU, 1, A( 1, K ), $ 1 )