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

Diff for /rpl/lapack/lapack/zungrq.f between versions 1.7 and 1.8

version 1.7, 2010/12/21 13:53:58 version 1.8, 2011/07/22 07:38:22
Line 1 Line 1
       SUBROUTINE ZUNGRQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )        SUBROUTINE ZUNGRQ( 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 19 Line 19
 *  which is defined as the last M rows of a product of K elementary  *  which is defined as the last M rows of a product of K elementary
 *  reflectors of order N  *  reflectors of order N
 *  *
 *        Q  =  H(1)' H(2)' . . . H(k)'  *        Q  =  H(1)**H H(2)**H . . . H(k)**H
 *  *
 *  as returned by ZGERQF.  *  as returned by ZGERQF.
 *  *
Line 193 Line 193
                CALL ZLARFT( 'Backward', 'Rowwise', N-K+I+IB-1, IB,                 CALL ZLARFT( '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**H to A(1:m-k+i-1,1:n-k+i+ib-1) from the right
 *  *
                CALL ZLARFB( 'Right', 'Conjugate transpose', 'Backward',                 CALL ZLARFB( 'Right', 'Conjugate transpose', 'Backward',
      $                      'Rowwise', II-1, N-K+I+IB-1, IB, A( II, 1 ),       $                      'Rowwise', II-1, N-K+I+IB-1, IB, A( II, 1 ),
Line 201 Line 201
      $                      LDWORK )       $                      LDWORK )
             END IF              END IF
 *  *
 *           Apply H' to columns 1:n-k+i+ib-1 of current block  *           Apply H**H to columns 1:n-k+i+ib-1 of current block
 *  *
             CALL ZUNGR2( IB, N-K+I+IB-1, IB, A( II, 1 ), LDA, TAU( I ),              CALL ZUNGR2( IB, N-K+I+IB-1, IB, A( II, 1 ), LDA, TAU( I ),
      $                   WORK, IINFO )       $                   WORK, IINFO )
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>