| 281 | } |
| 282 | |
| 283 | void EigenValues(const SymmetricMatrix& A, DiagonalMatrix& D, Matrix& Z) |
| 284 | { REPORT DiagonalMatrix E; tred2(A, D, E, Z); tql2(D, E, Z); SortSV(D,Z,true); } |
| 285 | |
| 286 | void EigenValues(const SymmetricMatrix& X, DiagonalMatrix& D) |
| 287 | { REPORT DiagonalMatrix E; SymmetricMatrix A; tred3(X,D,E,A); tql1(D,E); } |