[BACK]Return to dorgrq.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dorgrq.f between versions 1.4 and 1.8

version 1.4, 2010/08/06 15:32:31 version 1.8, 2011/07/22 07:38:08
Line 1 Line 1
       SUBROUTINE DORGRQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )        SUBROUTINE DORGRQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
 *  *
 *  -- LAPACK routine (version 3.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..--
 *     November 2006  *  -- April 2011                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            INFO, K, LDA, LWORK, M, N        INTEGER            INFO, K, LDA, LWORK, M, N
Line 193 Line 193
                CALL DLARFT( 'Backward', 'Rowwise', N-K+I+IB-1, IB,                 CALL DLARFT( 'Backward', 'Rowwise', N-K+I+IB-1, IB,
      $                      A( II, 1 ), LDA, TAU( I ), WORK, LDWORK )       $                      A( II, 1 ), LDA, TAU( I ), WORK, LDWORK )
 *  *
 *              Apply H' to A(1:m-k+i-1,1:n-k+i+ib-1) from the right  *              Apply H**T to A(1:m-k+i-1,1:n-k+i+ib-1) from the right
 *  *
                CALL DLARFB( 'Right', 'Transpose', 'Backward', 'Rowwise',                 CALL DLARFB( 'Right', 'Transpose', 'Backward', 'Rowwise',
      $                      II-1, N-K+I+IB-1, IB, A( II, 1 ), LDA, WORK,       $                      II-1, N-K+I+IB-1, IB, A( II, 1 ), LDA, WORK,
      $                      LDWORK, A, LDA, WORK( IB+1 ), LDWORK )       $                      LDWORK, A, LDA, WORK( IB+1 ), LDWORK )
             END IF              END IF
 *  *
 *           Apply H' to columns 1:n-k+i+ib-1 of current block  *           Apply H**T to columns 1:n-k+i+ib-1 of current block
 *  *
             CALL DORGR2( IB, N-K+I+IB-1, IB, A( II, 1 ), LDA, TAU( I ),              CALL DORGR2( IB, N-K+I+IB-1, IB, A( II, 1 ), LDA, TAU( I ),
      $                   WORK, IINFO )       $                   WORK, IINFO )
Removed from v.1.4  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>