doc/snippets/SelfAdjointEigenSolver_compute_MatrixType2.cpp


1
2
3
4
5
6
7
8
9

MatrixXd X = MatrixXd::Random(5,5);
MatrixXd A = X * X.transpose();
X = MatrixXd::Random(5,5);
MatrixXd B = X * X.transpose();

SelfAdjointEigenSolver<MatrixXd> es(A,B,false);
cout << "The eigenvalues of the pencil (A,B) are:" << endl << es.eigenvalues() << endl;
es.compute(B,A,false);
cout << "The eigenvalues of the pencil (B,A) are:" << endl << es.eigenvalues() << endl;