version 1.4, 2010/08/06 15:32:25
|
version 1.8, 2011/07/22 07:38:06
|
Line 26
|
Line 26
|
* of a full symmetric matrix (which was reduced to tridiagonal form) |
* of a full symmetric matrix (which was reduced to tridiagonal form) |
* are desired. |
* 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**T*u, 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 |