--- rpl/lapack/lapack/dgels.f	2010/12/21 13:53:25	1.7
+++ rpl/lapack/lapack/dgels.f	2011/07/22 07:38:04	1.8
@@ -1,10 +1,10 @@
       SUBROUTINE DGELS( TRANS, M, N, NRHS, A, LDA, B, LDB, 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 ..
       CHARACTER          TRANS
@@ -277,7 +277,7 @@
 *
 *           Least-Squares Problem min || A * X - B ||
 *
-*           B(1:M,1:NRHS) := Q' * B(1:M,1:NRHS)
+*           B(1:M,1:NRHS) := Q**T * B(1:M,1:NRHS)
 *
             CALL DORMQR( 'Left', 'Transpose', M, NRHS, N, A, LDA,
      $                   WORK( 1 ), B, LDB, WORK( MN+1 ), LWORK-MN,
@@ -298,9 +298,9 @@
 *
          ELSE
 *
-*           Overdetermined system of equations A' * X = B
+*           Overdetermined system of equations A**T * X = B
 *
-*           B(1:N,1:NRHS) := inv(R') * B(1:N,1:NRHS)
+*           B(1:N,1:NRHS) := inv(R**T) * B(1:N,1:NRHS)
 *
             CALL DTRTRS( 'Upper', 'Transpose', 'Non-unit', N, NRHS,
      $                   A, LDA, B, LDB, INFO )
@@ -359,7 +359,7 @@
    30          CONTINUE
    40       CONTINUE
 *
-*           B(1:N,1:NRHS) := Q(1:N,:)' * B(1:M,1:NRHS)
+*           B(1:N,1:NRHS) := Q(1:N,:)**T * B(1:M,1:NRHS)
 *
             CALL DORMLQ( 'Left', 'Transpose', N, NRHS, M, A, LDA,
      $                   WORK( 1 ), B, LDB, WORK( MN+1 ), LWORK-MN,
@@ -371,7 +371,7 @@
 *
          ELSE
 *
-*           overdetermined system min || A' * X - B ||
+*           overdetermined system min || A**T * X - B ||
 *
 *           B(1:N,1:NRHS) := Q * B(1:N,1:NRHS)
 *
@@ -381,7 +381,7 @@
 *
 *           workspace at least NRHS, optimally NRHS*NB
 *
-*           B(1:M,1:NRHS) := inv(L') * B(1:M,1:NRHS)
+*           B(1:M,1:NRHS) := inv(L**T) * B(1:M,1:NRHS)
 *
             CALL DTRTRS( 'Lower', 'Transpose', 'Non-unit', M, NRHS,
      $                   A, LDA, B, LDB, INFO )