version 1.1, 2010/01/26 15:22:46
|
version 1.8, 2010/12/21 13:53:43
|
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.2.2) -- |
* -- 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 |
* June 2010 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER INFO, LDA, M, N |
INTEGER INFO, LDA, M, N |
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. |
* June 2010 |
* 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 |
* |
* |