6.4 Computing Selected Eigenvalues
Often you require only some of the eigenvalues: either all the eigenvalues in a range of real numbers (a,b), or the ith through jth eigenvalues. You can solve these problems using the server objects RWSymRangeEigServer<double> and RWeSymSomeEigServer<double>. In addition to showing how to solve this specific problem, this section demonstrates how to use a server object to adjust an algorithm's parameters.
This program computes the eigenvalues of the matrix A in the range (a,b). Eigenvectors are not computed.
#include <iostream>
#include <rw/lapack/symmat.h>
#include <rw/lapack/seigsrv.h>
int main()
{
RWSymMat<double> A; // 1
double a,b;
std::cin >> A >> a >> b;
RWSymRangeEigServer<double> server; // 2
server.computeEigenVectors(true); // 3
server.setRange(a,b); // 4
RWSymEigDecomp<double> decomp = server(A); // 5
std::cout << decomp.eigenValues() << std::endl;
return 0;
}
| //1 | First, the matrix and the range are read in from standard input. |
| //2 | When the eigenvalue server is constructed, the algorithm parameters are all set to default values. |
| //3 | Here we change the setting of a parameter so that only eigenvalues, and not eigenvectors, are computed. |
| //4 | Here we set the range in which we are interested and... |
| //5 | We compute the decomposition. |
If instead of needing all the eigenvalues in a certain range, you need the ith through jth eigenvalues, you use a RWSymSomeEigServer<double> object like this:
#include <iostream>
#include <rw/lapack/symmat.h>
#include <rw/lapack/seigsrv.h>
int main()
{
RWSymMat<double> A;
int i,j;
std::cin >> A >> i >> j;
RWSymSomeEigServer<double> server;
server.computeEigenVectors(true);
server.setRange(RWRange(i,j)); // 1
RWSymEigDecomp<double> decomp = server(A);
std::cout << decomp.eigenValues() << std::endl;
return 0;
}
The program is the same as the previous program except for //1, where we set the range of interest. The range is set using one of the subscripting objects RWSlice, RWRange, or RWToEnd. See the Essential Math Module User's Guide and the class descriptions in the Essential Math Module Reference Guide for more information on the subscripting objects.
© Copyright Rogue Wave Software, Inc. All Rights Reserved.
Rogue Wave and SourcePro are registered trademarks of Rogue Wave Software, Inc. in the United States and other countries. All other trademarks are the property of their respective owners.
Contact Rogue Wave about documentation or support issues.