Visible to Intel only — GUID: GUID-0D4734FE-AF54-4D69-80D2-583704918540
Visible to Intel only — GUID: GUID-0D4734FE-AF54-4D69-80D2-583704918540
?stegr2a
Computes selected eigenvalues and initial representations needed for eigenvector computations.
Syntax
call sstegr2a( jobz, range, n, d, e, vl, vu, il, iu, m, w, z, ldz, nzc, work, lwork, iwork, liwork, dol, dou, needil, neediu, inderr, nsplit, pivmin, scale, wl, wu, info )
call dstegr2a( jobz, range, n, d, e, vl, vu, il, iu, m, w, z, ldz, nzc, work, lwork, iwork, liwork, dol, dou, needil, neediu, inderr, nsplit, pivmin, scale, wl, wu, info )
Description
?stegr2a computes selected eigenvalues and initial representations needed for eigenvector computations in ?stegr2b. It is invoked in the ScaLAPACK MRRR driver p?syevr and the corresponding Hermitian version when both eigenvalues and eigenvectors are computed in parallel on multiple processors. For this case, ?stegr2a implements the first part of the MRRR algorithm, parallel eigenvalue computation and finding the root RRR. At the end of ?stegr2a, other processors might have a part of the spectrum that is needed to continue the computation locally. Once this eigenvalue information has been received by the processor, the computation can then proceed by calling the second part of the parallel MRRR algorithm, ?stegr2b.
Please note:
The calling sequence has two additional integer parameters, (compared to LAPACK's stegr), these are dol and dou and should satisfy m≥dou≥dol≥1. These parameters are only relevant for the case jobz = 'V'.
Globally invoked over all processors, ?stegr2a computes all the eigenvalues specified by range.
?stegr2a locally only computes the eigenvalues corresponding to eigenvalues dol through dou in w. (That is, instead of computing the eigenvectors belonging to w(1) through w(m), only the eigenvectors belonging to eigenvalues w(dol) through w(dou) are computed. In this case, only the eigenvalues dol:dou are guaranteed to be fully accurate.
m is not the number of eigenvalues specified by range, but it is m = dou - dol + 1. Instead, m refers to the number of eigenvalues computed on this processor.
While no eigenvectors are computed in ?stegr2a itself (this is done later in ?stegr2b), the interface
If jobz = 'V' then, depending on range and dol, dou, ?stegr2a might need more workspace in z then the original ?stegr. In particular, the arrays w and z might not contain all the wanted eigenpairs locally, instead this information is distributed over other processors.
Product and Performance Information |
---|
Performance varies by use, configuration and other factors. Learn more at www.Intel.com/PerformanceIndex. Notice revision #20201201 |
Input Parameters
- jobz
-
CHARACTER*1
= 'N': Compute eigenvalues only;
= 'V': Compute eigenvalues and eigenvectors.
- range
-
CHARACTER*1
= 'A': all eigenvalues will be found.
= 'V': all eigenvalues in the half-open interval (vl,vu] will be found.
= 'I': the il-th through iu-th eigenvalues (of the entire matrix) will be found.
- n
-
INTEGER
The order of the matrix. n≥ 0.
- d
-
REAL for sstegr2a
DOUBLE PRECISION for dstegr2a
Array of size n
The n diagonal elements of the tridiagonal matrix T. Overwritten on exit.
- e
-
REAL for sstegr2a
DOUBLE PRECISION for dstegr2a
Array of size n
On entry, the (n-1) subdiagonal elements of the tridiagonal matrix T in elements 1 to n-1 of e. e(n) need not be set on input, but is used internally as workspace. Overwritten on exit.
- vl, vu
-
REAL for sstegr2a
DOUBLE PRECISION for dstegr2a
If range='V', the lower and upper bounds of the interval to be searched for eigenvalues. vl < vu.
Not referenced if range = 'A' or 'I'.
- il, iu
-
INTEGER
If range='I', the indices (in ascending order) of the smallest and largest eigenvalues to be returned. 1 ≤il≤iu≤n, if n > 0.
Not referenced if range = 'A' or 'V'.
- ldz
-
INTEGER
The leading dimension of the array z. ldz≥ 1, and if jobz = 'V', then ldz≥ max(1,n).
- nzc
-
INTEGER
The number of eigenvectors to be held in the array z.
If range = 'A', then nzc≥ max(1,n).
If range = 'V', then nzc≥ the number of eigenvalues in (vl,vu].
If range = 'I', then nzc≥iu-il+1.
If nzc = -1, then a workspace query is assumed; the routine calculates the number of columns of the array z that are needed to hold the eigenvectors. This value is returned as the first entry of the z array, and no error message related to nzc is issued.
- lwork
-
INTEGER
The size of the array work. lwork≥ max(1,18*n) if jobz = 'V', and lwork≥ max(1,12*n) if jobz = 'N'.
If lwork = -1, then a workspace query is assumed; the routine only calculates the optimal size of the work array, returns this value as the first entry of the work array, and no error message related to lwork is issued.
- liwork
-
INTEGER
The size of the array iwork. liwork≥ max(1,10*n) if the eigenvectors are desired, and liwork≥ max(1,8*n) if only the eigenvalues are to be computed.
If liwork = -1, then a workspace query is assumed; the routine only calculates the optimal size of the iwork array, returns this value as the first entry of the iwork array, and no error message related to liwork is issued.
- dol, dou
-
INTEGER
From all the eigenvalues w(1:m), only eigenvalues w(dol:dou) are computed.
OUTPUT Parameters
- m
-
INTEGER
Globally summed over all processors, m equals the total number of eigenvalues found. 0 ≤m≤n.
If range = 'A', m = n, and if range = 'I', m = iu-il+1.
The local output equals m = dou - dol + 1.
- w
-
REAL for sstegr2a
DOUBLE PRECISION for dstegr2a
Array of size n
The first m elements contain approximations to the selected eigenvalues in ascending order. Note that immediately after exiting this routine, only the eigenvalues from position dol:dou are reliable on this processor because the eigenvalue computation is done in parallel. The other entries outside dol:douare very crude preliminary approximations. Other processors hold reliable information on these other parts of the w array.
This information is communicated in the ScaLAPACK driver.
- z
-
REAL for sstegr2a
DOUBLE PRECISION for dstegr2a
Array of size (ldz, max(1,m) ).
?stegr2a does not compute eigenvectors, this is done in ?stegr2b. The argument z as well as all related other arguments only appear to keep the interface consistent and to signal to the user that this subroutine is meant to be used when eigenvectors are computed.
- work
-
On exit, if info = 0, work(1) returns the optimal (and minimal) lwork.
- iwork
-
On exit, if info = 0, iwork(1) returns the optimal liwork.
- needil, neediu
-
INTEGER
The indices of the leftmost and rightmost eigenvalues needed to accurately compute the relevant part of the representation tree. This information can be used to find out which processors have the relevant eigenvalue information needed so that it can be communicated.
- inderr
-
INTEGER
inderr points to the place in the work space where the eigenvalue uncertainties (errors) are stored.
- nsplit
-
INTEGER
The number of blocks into which T splits. 1 ≤nsplit≤n.
- pivmin
-
REAL for sstegr2a
DOUBLE PRECISION for dstegr2a
The minimum pivot in the sturm sequence for T.
- scale
-
REAL for sstegr2a
DOUBLE PRECISION for dstegr2a
The scaling factor for the tridiagonal T.
- wl, wu
-
REAL for sstegr2a
DOUBLE PRECISION for dstegr2a
The interval (wl, wu] contains all the wanted eigenvalues.
It is either given by the user or computed in ?larre2a.
- info
-
INTEGER
On exit, info = 0: successful exit
other: if info = -i, the i-th argument had an illegal value
if info = 10x, internal error in ?larre2a,
Here, the digit x = abs( iinfo ) < 10, where iinfo is the nonzero error code returned by ?larre2a.