--- rpl/lapack/lapack/zlatzm.f 2010/12/21 13:53:52 1.7 +++ rpl/lapack/lapack/zlatzm.f 2011/07/22 07:38:18 1.8 @@ -1,9 +1,9 @@ SUBROUTINE ZLATZM( SIDE, M, N, V, INCV, TAU, C1, C2, LDC, WORK ) * -* -- 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 .. CHARACTER SIDE @@ -21,8 +21,8 @@ * * ZLATZM applies a Householder matrix generated by ZTZRQF to a matrix. * -* Let P = I - tau*u*u', u = ( 1 ), -* ( v ) +* Let P = I - tau*u*u**H, u = ( 1 ), +* ( v ) * where v is an (m-1) vector if SIDE = 'L', or a (n-1) vector if * SIDE = 'R'. * @@ -112,14 +112,14 @@ * IF( LSAME( SIDE, 'L' ) ) THEN * -* w := conjg( C1 + v' * C2 ) +* w := ( C1 + v**H * C2 )**H * CALL ZCOPY( N, C1, LDC, WORK, 1 ) CALL ZLACGV( N, WORK, 1 ) CALL ZGEMV( 'Conjugate transpose', M-1, N, ONE, C2, LDC, V, $ INCV, ONE, WORK, 1 ) * -* [ C1 ] := [ C1 ] - tau* [ 1 ] * w' +* [ C1 ] := [ C1 ] - tau* [ 1 ] * w**H * [ C2 ] [ C2 ] [ v ] * CALL ZLACGV( N, WORK, 1 ) @@ -134,7 +134,7 @@ CALL ZGEMV( 'No transpose', M, N-1, ONE, C2, LDC, V, INCV, ONE, $ WORK, 1 ) * -* [ C1, C2 ] := [ C1, C2 ] - tau* w * [ 1 , v'] +* [ C1, C2 ] := [ C1, C2 ] - tau* w * [ 1 , v**H] * CALL ZAXPY( M, -TAU, WORK, 1, C1, 1 ) CALL ZGERC( M, N-1, -TAU, WORK, 1, V, INCV, C2, LDC )