1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathMatrixXd X = MatrixXd::Random(5,5); 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathMatrixXd A = X * X.transpose(); 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathX = MatrixXd::Random(5,5); 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathMatrixXd B = X * X.transpose(); 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathGeneralizedSelfAdjointEigenSolver<MatrixXd> es(A,B,EigenvaluesOnly); 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathcout << "The eigenvalues of the pencil (A,B) are:" << endl << es.eigenvalues() << endl; 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathes.compute(B,A,false); 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathcout << "The eigenvalues of the pencil (B,A) are:" << endl << es.eigenvalues() << endl; 10