--- rpl/lapack/lapack/zungrq.f	2010/12/21 13:53:58	1.7
+++ rpl/lapack/lapack/zungrq.f	2011/07/22 07:38:22	1.8
@@ -1,9 +1,9 @@
       SUBROUTINE ZUNGRQ( M, N, K, A, LDA, TAU, 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, K, LDA, LWORK, M, N
@@ -19,7 +19,7 @@
 *  which is defined as the last M rows of a product of K elementary
 *  reflectors of order N
 *
-*        Q  =  H(1)' H(2)' . . . H(k)'
+*        Q  =  H(1)**H H(2)**H . . . H(k)**H
 *
 *  as returned by ZGERQF.
 *
@@ -193,7 +193,7 @@
                CALL ZLARFT( 'Backward', 'Rowwise', N-K+I+IB-1, IB,
      $                      A( II, 1 ), LDA, TAU( I ), WORK, LDWORK )
 *
-*              Apply H' to A(1:m-k+i-1,1:n-k+i+ib-1) from the right
+*              Apply H**H to A(1:m-k+i-1,1:n-k+i+ib-1) from the right
 *
                CALL ZLARFB( 'Right', 'Conjugate transpose', 'Backward',
      $                      'Rowwise', II-1, N-K+I+IB-1, IB, A( II, 1 ),
@@ -201,7 +201,7 @@
      $                      LDWORK )
             END IF
 *
-*           Apply H' to columns 1:n-k+i+ib-1 of current block
+*           Apply H**H to columns 1:n-k+i+ib-1 of current block
 *
             CALL ZUNGR2( IB, N-K+I+IB-1, IB, A( II, 1 ), LDA, TAU( I ),
      $                   WORK, IINFO )