--- rpl/lapack/lapack/dgelsx.f 2010/12/21 13:53:25 1.7 +++ rpl/lapack/lapack/dgelsx.f 2011/07/22 07:38:04 1.8 @@ -1,10 +1,10 @@ SUBROUTINE DGELSX( M, N, NRHS, A, LDA, B, LDB, JPVT, RCOND, RANK, $ WORK, 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, M, N, NRHS, RANK @@ -44,8 +44,8 @@ * A * P = Q * [ T11 0 ] * Z * [ 0 0 ] * The minimum-norm solution is then -* X = P * Z' [ inv(T11)*Q1'*B ] -* [ 0 ] +* X = P * Z**T [ inv(T11)*Q1**T*B ] +* [ 0 ] * where Q1 consists of the first RANK columns of Q. * * Arguments @@ -267,7 +267,7 @@ * * Details of Householder rotations stored in WORK(MN+1:2*MN) * -* B(1:M,1:NRHS) := Q' * B(1:M,1:NRHS) +* B(1:M,1:NRHS) := Q**T * B(1:M,1:NRHS) * CALL DORM2R( 'Left', 'Transpose', M, NRHS, MN, A, LDA, WORK( 1 ), $ B, LDB, WORK( 2*MN+1 ), INFO ) @@ -285,7 +285,7 @@ 30 CONTINUE 40 CONTINUE * -* B(1:N,1:NRHS) := Y' * B(1:N,1:NRHS) +* B(1:N,1:NRHS) := Y**T * B(1:N,1:NRHS) * IF( RANK.LT.N ) THEN DO 50 I = 1, RANK