version 1.1, 2010/01/26 15:22:46
|
version 1.9, 2011/07/22 07:38:14
|
Line 1
|
Line 1
|
SUBROUTINE ZGEQPF( M, N, A, LDA, JPVT, TAU, WORK, RWORK, INFO ) |
SUBROUTINE ZGEQPF( M, N, A, LDA, JPVT, TAU, WORK, RWORK, INFO ) |
* |
* |
* -- LAPACK deprecated driver routine (version 3.2) -- |
* -- LAPACK deprecated computational 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 .. |
INTEGER INFO, LDA, M, N |
INTEGER INFO, LDA, M, N |
Line 69
|
Line 69
|
* |
* |
* Each H(i) has the form |
* Each H(i) has the form |
* |
* |
* H = I - tau * v * v' |
* H = I - tau * v * v**H |
* |
* |
* where tau is a complex scalar, and v is a complex vector with |
* where tau is a complex scalar, and v is a complex vector with |
* v(1:i-1) = 0 and v(i) = 1; v(i+1:m) is stored on exit in A(i+1:m,i). |
* v(1:i-1) = 0 and v(i) = 1; v(i+1:m) is stored on exit in A(i+1:m,i). |
Line 81
|
Line 81
|
* Partial column norm updating strategy modified by |
* Partial column norm updating strategy modified by |
* Z. Drmac and Z. Bujanovic, Dept. of Mathematics, |
* Z. Drmac and Z. Bujanovic, Dept. of Mathematics, |
* University of Zagreb, Croatia. |
* University of Zagreb, Croatia. |
* June 2006. |
* -- April 2011 -- |
* For more details see LAPACK Working Note 176. |
* For more details see LAPACK Working Note 176. |
* |
* |
* ===================================================================== |
* ===================================================================== |
Line 96
|
Line 96
|
COMPLEX*16 AII |
COMPLEX*16 AII |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL XERBLA, ZGEQR2, ZLARF, ZLARFP, ZSWAP, ZUNM2R |
EXTERNAL XERBLA, ZGEQR2, ZLARF, ZLARFG, ZSWAP, ZUNM2R |
* .. |
* .. |
* .. Intrinsic Functions .. |
* .. Intrinsic Functions .. |
INTRINSIC ABS, DCMPLX, DCONJG, MAX, MIN, SQRT |
INTRINSIC ABS, DCMPLX, DCONJG, MAX, MIN, SQRT |
Line 186
|
Line 186
|
* Generate elementary reflector H(i) |
* Generate elementary reflector H(i) |
* |
* |
AII = A( I, I ) |
AII = A( I, I ) |
CALL ZLARFP( M-I+1, AII, A( MIN( I+1, M ), I ), 1, |
CALL ZLARFG( M-I+1, AII, A( MIN( I+1, M ), I ), 1, |
$ TAU( I ) ) |
$ TAU( I ) ) |
A( I, I ) = AII |
A( I, I ) = AII |
* |
* |