NAG FL Interface
f08jjf (dstebz)

1 Purpose

f08jjf computes some (or all) of the eigenvalues of a real symmetric tridiagonal matrix, by bisection.

2 Specification

Fortran Interface
Subroutine f08jjf ( range, order, n, vl, vu, il, iu, abstol, d, e, m, nsplit, w, iblock, isplit, work, iwork, info)
Integer, Intent (In) :: n, il, iu
Integer, Intent (Out) :: m, nsplit, iblock(n), isplit(n), iwork(3*n), info
Real (Kind=nag_wp), Intent (In) :: vl, vu, abstol, d(*), e(*)
Real (Kind=nag_wp), Intent (Out) :: w(n), work(4*n)
Character (1), Intent (In) :: range, order
C Header Interface
#include <nag.h>
void  f08jjf_ (const char *range, const char *order, const Integer *n, const double *vl, const double *vu, const Integer *il, const Integer *iu, const double *abstol, const double d[], const double e[], Integer *m, Integer *nsplit, double w[], Integer iblock[], Integer isplit[], double work[], Integer iwork[], Integer *info, const Charlen length_range, const Charlen length_order)
The routine may be called by the names f08jjf, nagf_lapackeig_dstebz or its LAPACK name dstebz.

3 Description

f08jjf uses bisection to compute some or all of the eigenvalues of a real symmetric tridiagonal matrix T.
It searches for zero or negligible off-diagonal elements of T to see if the matrix splits into block diagonal form:
T = T1 T2 . . . Tp .  
It performs bisection on each of the blocks Ti and returns the block index of each computed eigenvalue, so that a subsequent call to f08jkf to compute eigenvectors can also take advantage of the block structure.

4 References

Kahan W (1966) Accurate eigenvalues of a symmetric tridiagonal matrix Report CS41 Stanford University

5 Arguments

1: range Character(1) Input
On entry: indicates which eigenvalues are required.
range='A'
All the eigenvalues are required.
range='V'
All the eigenvalues in the half-open interval (vl,vu] are required.
range='I'
Eigenvalues with indices il to iu are required.
Constraint: range='A', 'V' or 'I'.
2: order Character(1) Input
On entry: indicates the order in which the eigenvalues and their block numbers are to be stored.
order='B'
The eigenvalues are to be grouped by split-off block and ordered from smallest to largest within each block.
order='E'
The eigenvalues for the entire matrix are to be ordered from smallest to largest.
Constraint: order='B' or 'E'.
3: n Integer Input
On entry: n, the order of the matrix T.
Constraint: n0.
4: vl Real (Kind=nag_wp) Input
5: vu Real (Kind=nag_wp) Input
On entry: if range='V', the lower and upper bounds, respectively, of the half-open interval
vl,vu within which the required eigenvalues lie.
If range='A' or 'I', vl is not referenced.
Constraint: if range='V', vl<vu.
6: il Integer Input
7: iu Integer Input
On entry: if range='I', the indices of the first and last eigenvalues, respectively, to be computed (assuming that the eigenvalues are in ascending order).
If range='A' or 'V', il is not referenced.
Constraint: if range='I', 1 il iu n .
8: abstol Real (Kind=nag_wp) Input
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 abstol. If abstol0.0, the tolerance is taken as machine precision×T1.
9: d* Real (Kind=nag_wp) array Input
Note: the dimension of the array d must be at least max1,n.
On entry: the diagonal elements of the tridiagonal matrix T.
10: e* Real (Kind=nag_wp) array Input
Note: the dimension of the array e must be at least max1,n-1.
On entry: the off-diagonal elements of the tridiagonal matrix T.
11: m Integer Output
On exit: m, the actual number of eigenvalues found.
12: nsplit Integer Output
On exit: the number of diagonal blocks which constitute the tridiagonal matrix T.
13: wn Real (Kind=nag_wp) array Output
On exit: the required eigenvalues of the tridiagonal matrix T stored in w1 to wm.
14: iblockn Integer array Output
On exit: at each row/column j where ej is zero or negligible, T is considered to split into a block diagonal matrix and iblocki contains the block number of the eigenvalue stored in wi, for i=1,2,,m. Note that iblocki<0 for some i whenever info=1 or 3 (see Section 6) and range='A' or 'V'.
15: isplitn Integer array Output
On exit: the leading nsplit elements contain the points at which T splits up into sub-matrices as follows. The first sub-matrix consists of rows/columns 1 to isplit1, the second sub-matrix consists of rows/columns isplit1+1 to isplit2, , and the nsplit(th) sub-matrix consists of rows/columns isplitnsplit-1+1 to isplitnsplit (=n).
16: work4×n Real (Kind=nag_wp) array Workspace
17: iwork3×n Integer array Workspace
18: info Integer Output
On exit: info=0 unless the routine detects an error (see Section 6).

6 Error Indicators and Warnings

If failures with info1 are causing persistent trouble and you have checked that the routine is being called correctly, please contact NAG.
info<0
If info=-i, argument i had an illegal value. An explanatory message is output, and execution of the program is terminated.
info=1
If range='A' or 'V', the algorithm failed to compute some (or all) of the required eigenvalues to the required accuracy. More precisely, iblockvalue<0 indicates that eigenvalue value (stored in wvalue) failed to converge.
info=2
If range='I', the algorithm failed to compute some (or all) of the required eigenvalues. Try calling the routine again with range='A'.
info=3
If range='I', see the description above for info=2.
If range='A' or 'V', see the description above for info=1.
info=4
No eigenvalues have been computed. The floating-point arithmetic on the computer is not behaving as expected.

7 Accuracy

The eigenvalues of T 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 QR method. However, the reduction to tridiagonal form (prior to calling the routine) may exclude the possibility of obtaining high relative accuracy in the small eigenvalues of the original matrix if its eigenvalues vary widely in magnitude.

8 Parallelism and Performance

f08jjf 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 routine. Please also consult the Users' Note for your implementation for any additional implementation-specific information.

9 Further Comments

There is no complex analogue of this routine.

10 Example

See Section 10 in f08fgf.