version 1.3, 2010/08/06 15:28:45
|
version 1.8, 2011/07/22 07:38:09
|
Line 1
|
Line 1
|
SUBROUTINE DORMQR( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, |
SUBROUTINE DORMQR( SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, |
$ WORK, LWORK, INFO ) |
$ WORK, LWORK, 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 .. |
CHARACTER SIDE, TRANS |
CHARACTER SIDE, TRANS |
Line 234
|
Line 234
|
$ LDA, TAU( I ), T, LDT ) |
$ LDA, TAU( I ), T, LDT ) |
IF( LEFT ) THEN |
IF( LEFT ) THEN |
* |
* |
* H or H' is applied to C(i:m,1:n) |
* H or H**T is applied to C(i:m,1:n) |
* |
* |
MI = M - I + 1 |
MI = M - I + 1 |
IC = I |
IC = I |
ELSE |
ELSE |
* |
* |
* H or H' is applied to C(1:m,i:n) |
* H or H**T is applied to C(1:m,i:n) |
* |
* |
NI = N - I + 1 |
NI = N - I + 1 |
JC = I |
JC = I |
END IF |
END IF |
* |
* |
* Apply H or H' |
* Apply H or H**T |
* |
* |
CALL DLARFB( SIDE, TRANS, 'Forward', 'Columnwise', MI, NI, |
CALL DLARFB( SIDE, TRANS, 'Forward', 'Columnwise', MI, NI, |
$ IB, A( I, I ), LDA, T, LDT, C( IC, JC ), LDC, |
$ IB, A( I, I ), LDA, T, LDT, C( IC, JC ), LDC, |