--- rpl/lapack/lapack/zlarz.f	2010/12/21 13:53:52	1.7
+++ rpl/lapack/lapack/zlarz.f	2011/07/22 07:38:18	1.8
@@ -1,9 +1,9 @@
       SUBROUTINE ZLARZ( SIDE, M, N, L, V, INCV, TAU, C, 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,13 +21,13 @@
 *  M-by-N matrix C, from either the left or the right. H is represented
 *  in the form
 *
-*        H = I - tau * v * v'
+*        H = I - tau * v * v**H
 *
 *  where tau is a complex scalar and v is a complex vector.
 *
 *  If tau = 0, then H is taken to be the unit matrix.
 *
-*  To apply H' (the conjugate transpose of H), supply conjg(tau) instead
+*  To apply H**H (the conjugate transpose of H), supply conjg(tau) instead
 *  tau.
 *
 *  H is a product of k elementary reflectors as returned by ZTZRZF.
@@ -105,7 +105,7 @@
             CALL ZCOPY( N, C, LDC, WORK, 1 )
             CALL ZLACGV( N, WORK, 1 )
 *
-*           w( 1:n ) = conjg( w( 1:n ) + C( m-l+1:m, 1:n )' * v( 1:l ) )
+*           w( 1:n ) = conjg( w( 1:n ) + C( m-l+1:m, 1:n )**H * v( 1:l ) )
 *
             CALL ZGEMV( 'Conjugate transpose', L, N, ONE, C( M-L+1, 1 ),
      $                  LDC, V, INCV, ONE, WORK, 1 )
@@ -116,7 +116,7 @@
             CALL ZAXPY( N, -TAU, WORK, 1, C, LDC )
 *
 *           C( m-l+1:m, 1:n ) = C( m-l+1:m, 1:n ) - ...
-*                               tau * v( 1:l ) * conjg( w( 1:n )' )
+*                               tau * v( 1:l ) * w( 1:n )**H
 *
             CALL ZGERU( L, N, -TAU, V, INCV, WORK, 1, C( M-L+1, 1 ),
      $                  LDC )
@@ -142,7 +142,7 @@
             CALL ZAXPY( M, -TAU, WORK, 1, C, 1 )
 *
 *           C( 1:m, n-l+1:n ) = C( 1:m, n-l+1:n ) - ...
-*                               tau * w( 1:m ) * v( 1:l )'
+*                               tau * w( 1:m ) * v( 1:l )**H
 *
             CALL ZGERC( M, L, -TAU, WORK, 1, V, INCV, C( 1, N-L+1 ),
      $                  LDC )