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

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

version 1.7, 2010/12/21 13:53:32 version 1.8, 2011/07/22 07:38:07
Line 1 Line 1
       SUBROUTINE DLARZT( DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT )        SUBROUTINE DLARZT( DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT )
 *  *
 *  -- 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 ..
       CHARACTER          DIRECT, STOREV        CHARACTER          DIRECT, STOREV
Line 27 Line 27
 *  If STOREV = 'C', the vector which defines the elementary reflector  *  If STOREV = 'C', the vector which defines the elementary reflector
 *  H(i) is stored in the i-th column of the array V, and  *  H(i) is stored in the i-th column of the array V, and
 *  *
 *     H  =  I - V * T * V'  *     H  =  I - V * T * V**T
 *  *
 *  If STOREV = 'R', the vector which defines the elementary reflector  *  If STOREV = 'R', the vector which defines the elementary reflector
 *  H(i) is stored in the i-th row of the array V, and  *  H(i) is stored in the i-th row of the array V, and
 *  *
 *     H  =  I - V' * T * V  *     H  =  I - V**T * T * V
 *  *
 *  Currently, only STOREV = 'R' and DIRECT = 'B' are supported.  *  Currently, only STOREV = 'R' and DIRECT = 'B' are supported.
 *  *
Line 164 Line 164
 *  *
             IF( I.LT.K ) THEN              IF( I.LT.K ) THEN
 *  *
 *              T(i+1:k,i) = - tau(i) * V(i+1:k,1:n) * V(i,1:n)'  *              T(i+1:k,i) = - tau(i) * V(i+1:k,1:n) * V(i,1:n)**T
 *  *
                CALL DGEMV( 'No transpose', K-I, N, -TAU( I ),                 CALL DGEMV( 'No transpose', K-I, N, -TAU( I ),
      $                     V( I+1, 1 ), LDV, V( I, 1 ), LDV, ZERO,       $                     V( I+1, 1 ), LDV, V( I, 1 ), LDV, ZERO,
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>