1 EigenSolver
<MatrixXf
> es
;
2 MatrixXf A
= MatrixXf::Random(4,4);
3 es
.compute(A
, /* computeEigenvectors = */ false);
4 cout
<< "The eigenvalues of A are: " << es
.eigenvalues().transpose() << endl
;
5 es
.compute(A
+ MatrixXf::Identity(4,4), false); // re-use es to compute eigenvalues of A+I
6 cout
<< "The eigenvalues of A+I are: " << es
.eigenvalues().transpose() << endl
;
