The function may be called by the names: f08jjc, nag_lapackeig_dstebz or nag_dstebz.
f08jjc uses bisection to compute some or all of the eigenvalues of a real symmetric tridiagonal matrix .
It searches for zero or negligible off-diagonal elements of to see if the matrix splits into block diagonal form:
It performs bisection on each of the blocks and returns the block index of each computed eigenvalue, so that a subsequent call to f08jkc to compute eigenvectors can also take advantage of the block structure.
Kahan W (1966) Accurate eigenvalues of a symmetric tridiagonal matrix Report CS41 Stanford University
1: – Nag_RangeTypeInput
On entry: indicates which eigenvalues are required.
All the eigenvalues are required.
All the eigenvalues in the half-open interval (vl,vu] are required.
On entry: the absolute tolerance to which each eigenvalue is required. An eigenvalue (or cluster) is considered to have converged if it lies in an interval of width . If , the tolerance is taken as .
9: – const doubleInput
Note: the dimension, dim, of the array d
must be at least
On entry: the diagonal elements of the tridiagonal matrix .
10: – const doubleInput
Note: the dimension, dim, of the array e
must be at least
On entry: the off-diagonal elements of the tridiagonal matrix .
11: – Integer *Output
On exit: , the actual number of eigenvalues found.
12: – Integer *Output
On exit: the number of diagonal blocks which constitute the tridiagonal matrix .
13: – doubleOutput
On exit: the required eigenvalues of the tridiagonal matrix stored in to .
14: – IntegerOutput
On exit: at each row/column where is zero or negligible, is considered to split into a block diagonal matrix and
contains the block number of the eigenvalue stored in , for . Note that for some whenever NE_CONVERGENCE (see Section 6) and or .
15: – IntegerOutput
On exit: the leading nsplit elements contain the points at which splits up into sub-matrices as follows. The first sub-matrix consists of rows/columns to , the second sub-matrix consists of rows/columns to , , and the nsplit(th) sub-matrix consists of rows/columns to ().
16: – NagError *Input/Output
The NAG error argument (see Section 7 in the Introduction to the NAG Library CL Interface).
6Error Indicators and Warnings
Dynamic memory allocation failed.
See Section 3.1.2 in the Introduction to the NAG Library CL Interface for further information.
On entry, argument had an illegal value.
If or , the algorithm failed to compute some (or all) of the required eigenvalues to the required accuracy. More precisely, indicates that eigenvalue (stored in ) failed to converge.
If , the algorithm failed to compute some (or all) of the required eigenvalues. Try calling the function again with .
No eigenvalues have been computed. The floating-point arithmetic on the computer is not behaving as expected.
On entry, , , and .
Constraint: if , .
On entry, , and .
Constraint: if , .
On entry, .
An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact NAG for assistance.
See Section 7.5 in the Introduction to the NAG Library CL Interface for further information.
Your licence key may have expired or may not have been installed correctly.
See Section 8 in the Introduction to the NAG Library CL Interface for further information.
The eigenvalues of are computed to high relative accuracy which means that if they vary widely in magnitude, then any small eigenvalues will be computed more accurately than, for example, with the standard method. However, the reduction to tridiagonal form (prior to calling the function) may exclude the possibility of obtaining high relative accuracy in the small eigenvalues of the original matrix if its eigenvalues vary widely in magnitude.
8Parallelism and Performance
Background information to multithreading can be found in the Multithreading documentation.
f08jjc is threaded by NAG for parallel execution in multithreaded implementations of the NAG Library.
Please consult the X06 Chapter Introduction for information on how to control and interrogate the OpenMP environment used within this function. Please also consult the Users' Note for your implementation for any additional implementation-specific information.