--- rpl/lapack/lapack/zggglm.f	2010/12/21 13:53:45	1.7
+++ rpl/lapack/lapack/zggglm.f	2011/07/22 07:38:14	1.8
@@ -1,10 +1,10 @@
       SUBROUTINE ZGGGLM( N, M, P, A, LDA, B, LDB, D, X, Y, WORK, LWORK,
      $                   INFO )
 *
-*  -- LAPACK driver routine (version 3.2) --
+*  -- LAPACK driver 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, LDB, LWORK, M, N, P
@@ -187,9 +187,9 @@
 *
 *     Compute the GQR factorization of matrices A and B:
 *
-*            Q'*A = ( R11 ) M,    Q'*B*Z' = ( T11   T12 ) M
-*                   (  0  ) N-M             (  0    T22 ) N-M
-*                      M                     M+P-N  N-M
+*          Q**H*A = ( R11 ) M,    Q**H*B*Z**H = ( T11   T12 ) M
+*                   (  0  ) N-M                 (  0    T22 ) N-M
+*                      M                         M+P-N  N-M
 *
 *     where R11 and T22 are upper triangular, and Q and Z are
 *     unitary.
@@ -198,8 +198,8 @@
      $             WORK( M+NP+1 ), LWORK-M-NP, INFO )
       LOPT = WORK( M+NP+1 )
 *
-*     Update left-hand-side vector d = Q'*d = ( d1 ) M
-*                                             ( d2 ) N-M
+*     Update left-hand-side vector d = Q**H*d = ( d1 ) M
+*                                               ( d2 ) N-M
 *
       CALL ZUNMQR( 'Left', 'Conjugate transpose', N, 1, M, A, LDA, WORK,
      $             D, MAX( 1, N ), WORK( M+NP+1 ), LWORK-M-NP, INFO )
@@ -246,7 +246,7 @@
          CALL ZCOPY( M, D, 1, X, 1 )
       END IF
 *
-*     Backward transformation y = Z'*y
+*     Backward transformation y = Z**H *y
 *
       CALL ZUNMRQ( 'Left', 'Conjugate transpose', P, 1, NP,
      $             B( MAX( 1, N-P+1 ), 1 ), LDB, WORK( M+1 ), Y,