The basic object in the Linear Algebra Module nonsymmetric eigenvalue classes is the RWEigDecomp<double> object, which represents the eigenvalues and eigenvectors of a matrix. Here is a short program that reads a matrix from standard input, computes its eigenvalue decomposition, and prints out its eigenvalues and eigenvectors:
#include <iostream>
#include <rw/lapack/eig.h>
int main()
{
RWGenMat<double> A; // 1
std::cin >> A;
RWEigDecomp<double> decomp(A); // 2
std::cout << decomp.eigenValues() << std::endl; // 3
std::cout << decomp.leftEigenVectors() << std::endl; // 4
std::cout << decomp.rightEigenVectors() << std::endl;
return 0;
}
|
| //1 | A is defined to be a general dense matrix, and is read in from standard input. |
| //2 | Here we construct an eigenvalue decomposition object decomp from the matrix. The decomposition object contains all the eigenvalues and eigenvectors, as well as information about whether these quantities are fully accurate and were successfully computed. See the SourcePro C++ API Reference Guide for a complete description of the available member functions. |
| //3 | We can access a vector containing all the eigenvalues as shown in this line. It is also possible to access one eigenvalue at a time, using the member function eigenValue(int). Note that the eigenvalues are of type DComplex, since in general they are complex. |
| //4 | The next two lines print out matrices whose columns are the eigenvectors. You can also access a vector containing just one eigenvector using the member functions leftEigenVector(int) and rightEigenVector(int). |
Copyright © Rogue Wave Software, Inc. All Rights Reserved.
The Rogue Wave name and logo, and SourcePro, are registered trademarks of Rogue Wave Software. All other trademarks are the property of their respective owners.
Provide feedback to Rogue Wave about its documentation.