Diff for /rpl/lapack/lapack/zgeql2.f between versions 1.8 and 1.9

version 1.8, 2010/12/21 13:53:43 version 1.9, 2011/07/22 07:38:14
Line 1 Line 1
       SUBROUTINE ZGEQL2( M, N, A, LDA, TAU, WORK, INFO )        SUBROUTINE ZGEQL2( M, N, A, LDA, TAU, WORK, 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, M, N        INTEGER            INFO, LDA, M, N
Line 59 Line 59
 *  *
 *  Each H(i) has the form  *  Each H(i) has the form
 *  *
 *     H(i) = I - tau * v * v'  *     H(i) = I - tau * v * v**H
 *  *
 *  where tau is a complex scalar, and v is a complex vector with  *  where tau is a complex scalar, and v is a complex vector with
 *  v(m-k+i+1:m) = 0 and v(m-k+i) = 1; v(1:m-k+i-1) is stored on exit in  *  v(m-k+i+1:m) = 0 and v(m-k+i) = 1; v(1:m-k+i-1) is stored on exit in
Line 108 Line 108
          ALPHA = A( M-K+I, N-K+I )           ALPHA = A( M-K+I, N-K+I )
          CALL ZLARFG( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )           CALL ZLARFG( M-K+I, ALPHA, A( 1, N-K+I ), 1, TAU( I ) )
 *  *
 *        Apply H(i)' to A(1:m-k+i,1:n-k+i-1) from the left  *        Apply H(i)**H to A(1:m-k+i,1:n-k+i-1) from the left
 *  *
          A( M-K+I, N-K+I ) = ONE           A( M-K+I, N-K+I ) = ONE
          CALL ZLARF( 'Left', M-K+I, N-K+I-1, A( 1, N-K+I ), 1,           CALL ZLARF( 'Left', M-K+I, N-K+I-1, A( 1, N-K+I ), 1,

Removed from v.1.8  
changed lines
  Added in v.1.9


CVSweb interface <joel.bertrand@systella.fr>