Diff for /rpl/lapack/lapack/dgeqrfp.f between versions 1.4 and 1.5

version 1.4, 2010/12/21 13:53:26 version 1.5, 2011/07/22 07:38:05
Line 1 Line 1
       SUBROUTINE DGEQRFP( M, N, A, LDA, TAU, WORK, LWORK, INFO )        SUBROUTINE DGEQRFP( M, N, A, LDA, TAU, WORK, LWORK, INFO )
 *  *
 *  -- LAPACK routine (version 3.2.2) --  *  -- LAPACK routine (version 3.3.1) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --  *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--  *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
 *     June 2010  *  -- April 2011                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            INFO, LDA, LWORK, M, N        INTEGER            INFO, LDA, LWORK, M, N
Line 69 Line 69
 *  *
 *  Each H(i) has the form  *  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  *  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),  *  v(1:i-1) = 0 and v(i) = 1; v(i+1:m) is stored on exit in A(i+1:m,i),
Line 171 Line 171
                CALL DLARFT( 'Forward', 'Columnwise', M-I+1, IB,                 CALL DLARFT( 'Forward', 'Columnwise', M-I+1, IB,
      $                      A( I, I ), LDA, TAU( I ), WORK, LDWORK )       $                      A( I, I ), LDA, TAU( I ), WORK, LDWORK )
 *  *
 *              Apply H' to A(i:m,i+ib:n) from the left  *              Apply H**T to A(i:m,i+ib:n) from the left
 *  *
                CALL DLARFB( 'Left', 'Transpose', 'Forward',                 CALL DLARFB( 'Left', 'Transpose', 'Forward',
      $                      'Columnwise', M-I+1, N-I-IB+1, IB,       $                      'Columnwise', M-I+1, N-I-IB+1, IB,

Removed from v.1.4  
changed lines
  Added in v.1.5


CVSweb interface <joel.bertrand@systella.fr>