The following subroutines are used to solve non-symmetric generalized eigenvalue problems in real arithmetic. These routines are appropriate when is a general non-symmetric matrix and is symmetric and positive semi-definite. The reverse communication interface routine for the non-symmetric double precision eigenvalue problem is dnaupd. The routine is called as shown in Figure 3.6. The specification of which nev eigenvalues is controlled by the character*2 argument which. Table 3.4 lists the choices available.
Eigenvalues of largest magnitude.
Eigenvalues of smallest magnitude.
Eigenvalues of largest real part.
Eigenvalues of smallest real part.
Eigenvalues of largest imaginary part.
Eigenvalues of smallest imaginary part.
There are three different shift-invert modes for non-symmetric eigenvalue problems. These modes are specified by setting the parameter entry iparam(7) = mode where mode = 1,2,3, or 4.
In the following list, the specification of OP and are given for the various modes. Also, the iparam(7) and bmat settings are listed along with the name of the sample driver for the given mode. Sample drivers for the following modes may be found in the EXAMPLES/NONSYM subdirectory.
Note that there are two shift-invert modes with complex shifts (See dndrv5 and dndrv6). Since is complex, these both require the factorization of the matrix in complex arithmetic even though both and are real. The only advantage of using this option instead of using the standard shift-invert mode in complex arithmetic with the routine znaupd is that all of the internal operations in the IRAM are executed in real arithmetic. This results in a factor of two savings in storage and a factor of four savings in arithmetic. There is additional post-processing that is somewhat more complicated than the other modes in order to get the eigenvalues and eigenvectors of the original problem. These modes are only recommended if storage is extremely critical.