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

Diff for /rpl/lapack/lapack/zunglq.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 ZUNGLQ( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )        SUBROUTINE ZUNGLQ( 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 first M rows of a product of K elementary  *  which is defined as the first M rows of a product of K elementary
 *  reflectors of order N  *  reflectors of order N
 *  *
 *        Q  =  H(k)' . . . H(2)' H(1)'  *        Q  =  H(k)**H . . . H(2)**H H(1)**H
 *  *
 *  as returned by ZGELQF.  *  as returned by ZGELQF.
 *  *
Line 185 Line 185
                CALL ZLARFT( 'Forward', 'Rowwise', N-I+1, IB, A( I, I ),                 CALL ZLARFT( 'Forward', 'Rowwise', N-I+1, IB, A( I, I ),
      $                      LDA, TAU( I ), WORK, LDWORK )       $                      LDA, TAU( I ), WORK, LDWORK )
 *  *
 *              Apply H' to A(i+ib:m,i:n) from the right  *              Apply H**H to A(i+ib:m,i:n) from the right
 *  *
                CALL ZLARFB( 'Right', 'Conjugate transpose', 'Forward',                 CALL ZLARFB( 'Right', 'Conjugate transpose', 'Forward',
      $                      'Rowwise', M-I-IB+1, N-I+1, IB, A( I, I ),       $                      'Rowwise', M-I-IB+1, N-I+1, IB, A( I, I ),
Line 193 Line 193
      $                      WORK( IB+1 ), LDWORK )       $                      WORK( IB+1 ), LDWORK )
             END IF              END IF
 *  *
 *           Apply H' to columns i:n of current block  *           Apply H**H to columns i:n of current block
 *  *
             CALL ZUNGL2( IB, N-I+1, IB, A( I, I ), LDA, TAU( I ), WORK,              CALL ZUNGL2( IB, N-I+1, IB, A( I, I ), LDA, TAU( I ), WORK,
      $                   IINFO )       $                   IINFO )
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>