|
version 1.1.1.1, 2010/01/26 15:22:46
|
version 1.9, 2011/07/22 07:38:14
|
|
Line 1
|
Line 1
|
| SUBROUTINE ZGERQ2( M, N, A, LDA, TAU, WORK, INFO ) |
SUBROUTINE ZGERQ2( M, N, A, LDA, TAU, WORK, INFO ) |
| * |
* |
| * -- LAPACK routine (version 3.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..-- |
| * November 2006 |
* -- April 2011 -- |
| * |
* |
| * .. Scalar Arguments .. |
* .. Scalar Arguments .. |
| INTEGER INFO, LDA, M, N |
INTEGER INFO, LDA, M, N |
|
Line 55
|
Line 55
|
| * |
* |
| * The matrix Q is represented as a product of elementary reflectors |
* The matrix Q is represented as a product of elementary reflectors |
| * |
* |
| * Q = H(1)' H(2)' . . . H(k)', where k = min(m,n). |
* Q = H(1)**H H(2)**H . . . H(k)**H, where k = min(m,n). |
| * |
* |
| * Each H(i) has the form |
* Each H(i) has the form |
| * |
* |
| * H(i) = I - tau * v * v' |
* H(i) = 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(n-k+i+1:n) = 0 and v(n-k+i) = 1; conjg(v(1:n-k+i-1)) is stored on |
* v(n-k+i+1:n) = 0 and v(n-k+i) = 1; conjg(v(1:n-k+i-1)) is stored on |
|
Line 76
|
Line 76
|
| COMPLEX*16 ALPHA |
COMPLEX*16 ALPHA |
| * .. |
* .. |
| * .. External Subroutines .. |
* .. External Subroutines .. |
| EXTERNAL XERBLA, ZLACGV, ZLARF, ZLARFP |
EXTERNAL XERBLA, ZLACGV, ZLARF, ZLARFG |
| * .. |
* .. |
| * .. Intrinsic Functions .. |
* .. Intrinsic Functions .. |
| INTRINSIC MAX, MIN |
INTRINSIC MAX, MIN |
|
Line 107
|
Line 107
|
| * |
* |
| CALL ZLACGV( N-K+I, A( M-K+I, 1 ), LDA ) |
CALL ZLACGV( N-K+I, A( M-K+I, 1 ), LDA ) |
| ALPHA = A( M-K+I, N-K+I ) |
ALPHA = A( M-K+I, N-K+I ) |
| CALL ZLARFP( N-K+I, ALPHA, A( M-K+I, 1 ), LDA, TAU( I ) ) |
CALL ZLARFG( N-K+I, ALPHA, A( M-K+I, 1 ), LDA, TAU( I ) ) |
| * |
* |
| * Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right |
* Apply H(i) to A(1:m-k+i-1,1:n-k+i) from the right |
| * |
* |
![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to zgerq2.f CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [local] / rpl / lapack / lapack