--- rpl/lapack/lapack/dgeqr2.f 2010/01/26 15:22:46 1.1 +++ rpl/lapack/lapack/dgeqr2.f 2011/07/22 07:38:05 1.9 @@ -1,9 +1,9 @@ SUBROUTINE DGEQR2( M, N, A, LDA, TAU, WORK, 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, LDA, M, N @@ -57,7 +57,7 @@ * * Each H(i) has the form * -* H(i) = I - tau * v * v' +* H(i) = I - tau * v * v**T * * where tau is a real scalar, and v is a real vector with * v(1:i-1) = 0 and v(i) = 1; v(i+1:m) is stored on exit in A(i+1:m,i), @@ -74,7 +74,7 @@ DOUBLE PRECISION AII * .. * .. External Subroutines .. - EXTERNAL DLARF, DLARFP, XERBLA + EXTERNAL DLARF, DLARFG, XERBLA * .. * .. Intrinsic Functions .. INTRINSIC MAX, MIN @@ -102,7 +102,7 @@ * * Generate elementary reflector H(i) to annihilate A(i+1:m,i) * - CALL DLARFP( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1, + CALL DLARFG( M-I+1, A( I, I ), A( MIN( I+1, M ), I ), 1, $ TAU( I ) ) IF( I.LT.N ) THEN *