mirror of
https://gitlab.com/libeigen/eigen.git
synced 2024-12-21 07:19:46 +08:00
a594ac3966
Add an assert to guard against using eigenvalues that have not converged. Add call to info() in tutorial example to cover non-convergence.
19 lines
534 B
C++
19 lines
534 B
C++
#include <iostream>
|
|
#include <Eigen/Dense>
|
|
|
|
using namespace std;
|
|
using namespace Eigen;
|
|
|
|
int main()
|
|
{
|
|
Matrix2f A;
|
|
A << 1, 2, 2, 3;
|
|
cout << "Here is the matrix A:\n" << A << endl;
|
|
SelfAdjointEigenSolver<Matrix2f> eigensolver(A);
|
|
if (eigensolver.info() != Success) abort();
|
|
cout << "The eigenvalues of A are:\n" << eigensolver.eigenvalues() << endl;
|
|
cout << "Here's a matrix whose columns are eigenvectors of A \n"
|
|
<< "corresponding to these eigenvalues:\n"
|
|
<< eigensolver.eigenvectors() << endl;
|
|
}
|