version 1.8, 2010/12/21 13:53:41
|
version 1.9, 2011/07/22 07:38:13
|
Line 1
|
Line 1
|
SUBROUTINE DTZRQF( M, N, A, LDA, TAU, INFO ) |
SUBROUTINE DTZRQF( 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 )**T, 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 149
|
Line 149
|
$ LDA, A( K, M1 ), LDA, ONE, TAU, 1 ) |
$ LDA, A( K, M1 ), LDA, ONE, TAU, 1 ) |
* |
* |
* Now form a( k ) := a( k ) - tau*w |
* Now form a( k ) := a( k ) - tau*w |
* and B := B - tau*w*z( k )'. |
* and B := B - tau*w*z( k )**T. |
* |
* |
CALL DAXPY( K-1, -TAU( K ), TAU, 1, A( 1, K ), 1 ) |
CALL DAXPY( K-1, -TAU( K ), TAU, 1, A( 1, K ), 1 ) |
CALL DGER( K-1, N-M, -TAU( K ), TAU, 1, A( K, M1 ), LDA, |
CALL DGER( K-1, N-M, -TAU( K ), TAU, 1, A( K, M1 ), LDA, |