version 1.3, 2010/08/06 15:28:42
|
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, |