version 1.8, 2010/12/21 13:53:58
|
version 1.9, 2011/07/22 07:38:21
|
Line 1
|
Line 1
|
SUBROUTINE ZTZRQF( M, N, A, LDA, TAU, INFO ) |
SUBROUTINE ZTZRQF( M, N, A, LDA, TAU, 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 66
|
Line 66
|
* |
* |
* where |
* where |
* |
* |
* T( k ) = I - tau*u( k )*u( k )', u( k ) = ( 1 ), |
* T( k ) = I - tau*u( k )*u( k )**H, u( k ) = ( 1 ), |
* ( 0 ) |
* ( 0 ) |
* ( z( k ) ) |
* ( z( k ) ) |
* |
* |
* tau is a scalar and z( k ) is an ( n - m ) element vector. |
* tau is a scalar and z( k ) is an ( n - m ) element vector. |
* tau and z( k ) are chosen to annihilate the elements of the kth row |
* tau and z( k ) are chosen to annihilate the elements of the kth row |
Line 142
|
Line 142
|
* |
* |
IF( TAU( K ).NE.CZERO .AND. K.GT.1 ) THEN |
IF( TAU( K ).NE.CZERO .AND. K.GT.1 ) THEN |
* |
* |
* We now perform the operation A := A*P( k )'. |
* We now perform the operation A := A*P( k )**H. |
* |
* |
* Use the first ( k - 1 ) elements of TAU to store a( k ), |
* Use the first ( k - 1 ) elements of TAU to store a( k ), |
* where a( k ) consists of the first ( k - 1 ) elements of |
* where a( k ) consists of the first ( k - 1 ) elements of |
Line 157
|
Line 157
|
$ LDA, A( K, M1 ), LDA, CONE, TAU, 1 ) |
$ LDA, A( K, M1 ), LDA, CONE, TAU, 1 ) |
* |
* |
* Now form a( k ) := a( k ) - conjg(tau)*w |
* Now form a( k ) := a( k ) - conjg(tau)*w |
* and B := B - conjg(tau)*w*z( k )'. |
* and B := B - conjg(tau)*w*z( k )**H. |
* |
* |
CALL ZAXPY( K-1, -DCONJG( TAU( K ) ), TAU, 1, A( 1, K ), |
CALL ZAXPY( K-1, -DCONJG( TAU( K ) ), TAU, 1, A( 1, K ), |
$ 1 ) |
$ 1 ) |