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

Diff for /rpl/lapack/lapack/zgelq2.f between versions 1.2 and 1.9

version 1.2, 2010/04/21 13:45:28 version 1.9, 2011/07/22 07:38:14
Line 1 Line 1
       SUBROUTINE ZGELQ2( M, N, A, LDA, TAU, WORK, INFO )        SUBROUTINE ZGELQ2( M, N, A, LDA, TAU, WORK, 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, LDA, M, N        INTEGER            INFO, LDA, M, N
Line 53 Line 53
 *  *
 *  The matrix Q is represented as a product of elementary reflectors  *  The matrix Q is represented as a product of elementary reflectors
 *  *
 *     Q = H(k)' . . . H(2)' H(1)', where k = min(m,n).  *     Q = H(k)**H . . . H(2)**H H(1)**H, where k = min(m,n).
 *  *
 *  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(1:i-1) = 0 and v(i) = 1; conjg(v(i+1:n)) is stored on exit in  *  v(1:i-1) = 0 and v(i) = 1; conjg(v(i+1:n)) is stored on exit in
Line 74 Line 74
       COMPLEX*16         ALPHA        COMPLEX*16         ALPHA
 *     ..  *     ..
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFP        EXTERNAL           XERBLA, ZLACGV, ZLARF, ZLARFG
 *     ..  *     ..
 *     .. Intrinsic Functions ..  *     .. Intrinsic Functions ..
       INTRINSIC          MAX, MIN        INTRINSIC          MAX, MIN
Line 104 Line 104
 *  *
          CALL ZLACGV( N-I+1, A( I, I ), LDA )           CALL ZLACGV( N-I+1, A( I, I ), LDA )
          ALPHA = A( I, I )           ALPHA = A( I, I )
          CALL ZLARFP( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,           CALL ZLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
      $                TAU( I ) )       $                TAU( I ) )
          IF( I.LT.M ) THEN           IF( I.LT.M ) THEN
 *  *
Removed from v.1.2  
changed lines
  Added in v.1.9

CVSweb interface <joel.bertrand@systella.fr>