version 1.7, 2010/12/21 13:53:24
|
version 1.8, 2011/07/22 07:38:04
|
Line 1
|
Line 1
|
SUBROUTINE DGEBRD( M, N, A, LDA, D, E, TAUQ, TAUP, WORK, LWORK, |
SUBROUTINE DGEBRD( M, N, A, LDA, D, E, TAUQ, TAUP, WORK, LWORK, |
$ INFO ) |
$ 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, LWORK, M, N |
INTEGER INFO, LDA, LWORK, M, N |
Line 99
|
Line 99
|
* |
* |
* Each H(i) and G(i) has the form: |
* Each H(i) and G(i) has the form: |
* |
* |
* H(i) = I - tauq * v * v' and G(i) = I - taup * u * u' |
* H(i) = I - tauq * v * v**T and G(i) = I - taup * u * u**T |
* |
* |
* where tauq and taup are real scalars, and v and u are real vectors; |
* where tauq and taup are real scalars, and v and u are real vectors; |
* v(1:i-1) = 0, v(i) = 1, and v(i+1:m) is stored on exit in A(i+1:m,i); |
* v(1:i-1) = 0, v(i) = 1, and v(i+1:m) is stored on exit in A(i+1:m,i); |
Line 112
|
Line 112
|
* |
* |
* Each H(i) and G(i) has the form: |
* Each H(i) and G(i) has the form: |
* |
* |
* H(i) = I - tauq * v * v' and G(i) = I - taup * u * u' |
* H(i) = I - tauq * v * v**T and G(i) = I - taup * u * u**T |
* |
* |
* where tauq and taup are real scalars, and v and u are real vectors; |
* where tauq and taup are real scalars, and v and u are real vectors; |
* v(1:i) = 0, v(i+1) = 1, and v(i+2:m) is stored on exit in A(i+2:m,i); |
* v(1:i) = 0, v(i+1) = 1, and v(i+2:m) is stored on exit in A(i+2:m,i); |
Line 232
|
Line 232
|
$ WORK( LDWRKX*NB+1 ), LDWRKY ) |
$ WORK( LDWRKX*NB+1 ), LDWRKY ) |
* |
* |
* Update the trailing submatrix A(i+nb:m,i+nb:n), using an update |
* Update the trailing submatrix A(i+nb:m,i+nb:n), using an update |
* of the form A := A - V*Y' - X*U' |
* of the form A := A - V*Y**T - X*U**T |
* |
* |
CALL DGEMM( 'No transpose', 'Transpose', M-I-NB+1, N-I-NB+1, |
CALL DGEMM( 'No transpose', 'Transpose', M-I-NB+1, N-I-NB+1, |
$ NB, -ONE, A( I+NB, I ), LDA, |
$ NB, -ONE, A( I+NB, I ), LDA, |