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