version 1.7, 2010/12/21 13:53:29
|
version 1.8, 2011/07/22 07:38:06
|
Line 3
|
Line 3
|
$ PERM, GIVPTR, GIVCOL, GIVNUM, WORK, IWORK, |
$ PERM, GIVPTR, GIVCOL, GIVNUM, WORK, IWORK, |
$ 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 CURLVL, CURPBM, CUTPNT, ICOMPQ, INFO, LDQ, N, |
INTEGER CURLVL, CURPBM, CUTPNT, ICOMPQ, INFO, LDQ, N, |
Line 31
|
Line 31
|
* the case in which all eigenvalues and eigenvectors of a symmetric |
* the case in which all eigenvalues and eigenvectors of a symmetric |
* tridiagonal matrix are desired. |
* tridiagonal matrix are desired. |
* |
* |
* T = Q(in) ( D(in) + RHO * Z*Z' ) Q'(in) = Q(out) * D(out) * Q'(out) |
* T = Q(in) ( D(in) + RHO * Z*Z**T ) Q**T(in) = Q(out) * D(out) * Q**T(out) |
* |
* |
* where Z = Q'u, u is a vector of length N with ones in the |
* where Z = Q**Tu, u is a vector of length N with ones in the |
* CUTPNT and CUTPNT + 1 th elements and zeros elsewhere. |
* CUTPNT and CUTPNT + 1 th elements and zeros elsewhere. |
* |
* |
* The eigenvectors of the original matrix are stored in Q, and the |
* The eigenvectors of the original matrix are stored in Q, and the |