version 1.12, 2012/08/22 09:48:37
|
version 1.18, 2016/08/27 15:35:00
|
Line 1
|
Line 1
|
*> \brief \b ZLARFT |
*> \brief \b ZLARFT forms the triangular factor T of a block reflector H = I - vtvH |
* |
* |
* =========== DOCUMENTATION =========== |
* =========== DOCUMENTATION =========== |
* |
* |
Line 130
|
Line 130
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date April 2012 |
*> \date June 2016 |
* |
* |
*> \ingroup complex16OTHERauxiliary |
*> \ingroup complex16OTHERauxiliary |
* |
* |
Line 163
|
Line 163
|
* ===================================================================== |
* ===================================================================== |
SUBROUTINE ZLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT ) |
SUBROUTINE ZLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.4.1) -- |
* -- LAPACK auxiliary routine (version 3.6.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..-- |
* April 2012 |
* June 2016 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER DIRECT, STOREV |
CHARACTER DIRECT, STOREV |
Line 187
|
Line 187
|
INTEGER I, J, PREVLASTV, LASTV |
INTEGER I, J, PREVLASTV, LASTV |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL ZGEMV, ZLACGV, ZTRMV |
EXTERNAL ZGEMV, ZLACGV, ZTRMV, ZGEMM |
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
LOGICAL LSAME |
LOGICAL LSAME |
Line 303
|
Line 303
|
* |
* |
CALL ZGEMM( 'N', 'C', K-I, 1, N-K+I-J, -TAU( I ), |
CALL ZGEMM( 'N', 'C', K-I, 1, N-K+I-J, -TAU( I ), |
$ V( I+1, J ), LDV, V( I, J ), LDV, |
$ V( I+1, J ), LDV, V( I, J ), LDV, |
$ ONE, T( I+1, I ), LDT ) |
$ ONE, T( I+1, I ), LDT ) |
END IF |
END IF |
* |
* |
* T(i+1:k,i) := T(i+1:k,i+1:k) * T(i+1:k,i) |
* T(i+1:k,i) := T(i+1:k,i+1:k) * T(i+1:k,i) |