Developer Reference for Intel® oneAPI Math Kernel Library for Fortran
?larrj
Performs refinement of the initial estimates of the eigenvalues of the matrix T .
Syntax
call slarrj ( n , d , e2 , ifirst , ilast , rtol , offset , w , werr , work , iwork , pivmin, spdiam , info )
call dlarrj ( n , d , e2 , ifirst , ilast , rtol , offset , w , werr , work , iwork , pivmin, spdiam , info )
Include Files
mkl.fi
Description
slarrj dlarrj
Given the initial eigenvalue approximations of T , this routine does bisection to refine the eigenvalues of T , w(ifirst-offset) through w(ilast-offset) , to more accuracy. Initial guesses for these eigenvalues are input in w , the corresponding estimate of the error in these guesses in werr . During bisection, intervals [a,b] are maintained by storing their mid-points and semi-widths in the arrays w and werr respectively.
Input Parameters
- n
-
INTEGER . The order of the matrix T .
- d
-
REAL for slarrj DOUBLE PRECISION for dlarrj Array, DIMENSION ( n ). Contains n diagonal elements of the matrix T .
- e2
-
REAL for slarrj DOUBLE PRECISION for dlarrj Array, DIMENSION ( n-1 ). Contains (n-1) squared sub-diagonal elements of the T .
- ifirst
-
INTEGER .
- ilast
-
INTEGER .
- rtol
-
REAL for slarrj DOUBLE PRECISION for dlarrj Tolerance for the convergence of the bisection intervals. An interval [a,b] is considered to be converged if (b-a) ≤ rtol*max(|a|,|b|) .
- offset
-
INTEGER .
- w
-
REAL for slarrj DOUBLE PRECISION for dlarrj Array, DIMENSION ( n ). On input, w(ifirst-offset) through w(ilast-offset) are estimates of the eigenvalues of L*D*L^{T} indexed ifirst through ilast .
- werr
-
REAL for slarrj DOUBLE PRECISION for dlarrj Array, DIMENSION ( n ). On input, werr(ifirst-offset) through werr(ilast-offset) are the errors in the estimates of the corresponding elements in w .
- work
-
REAL for slarrj DOUBLE PRECISION for dlarrj Workspace array, DIMENSION ( 2*n ).
- iwork
-
INTEGER . Workspace array, DIMENSION ( 2*n ).
- pivmin
-
REAL for slarrj DOUBLE PRECISION for dlarrj The minimum pivot in the Sturm sequence for the matrix T .
- spdiam
-
REAL for slarrj DOUBLE PRECISION for dlarrj The spectral diameter of the matrix T .
Output Parameters
- w
-
On exit, contains the refined estimates of the eigenvalues.
- werr
-
On exit, contains the refined errors in the estimates of the corresponding elements in w .
info
INTEGER .
Now it is not used and always is set to 0.