SRC/slarrf.f(3)

๐ŸŒธ SRC/slarrf.f ๐ŸŒธ

From the Library Functions Manual โ€” LAPACK edition, bestie. This is the page for the slarrf source file, and it's giving subroutine with main character energy fr fr. ๐Ÿ’…โœจ


โœจ SYNOPSIS โœจ

๐Ÿ’ซ Functions/Subroutines ๐Ÿ’ซ

subroutine slarrf (n, d, l, ld, clstrt, clend, w, wgap, werr, spdiam,
    clgapl, clgapr, pivmin, sigma, dplus, lplus, work, info)

SLARRF finds a new relatively robust representation such that at least one of the eigenvalues is relatively isolated.

๐ŸŽ€ Function/Subroutine Documentation ๐ŸŽ€

subroutine slarrf (integer n, real, dimension( * ) d, real, dimension( * )
    l, real, dimension( * ) ld, integer clstrt, integer clend, real,
    dimension( * ) w, real, dimension( * ) wgap, real, dimension( * ) werr,
    real spdiam, real clgapl, real clgapr, real pivmin, real sigma, real,
    dimension( * ) dplus, real, dimension( * ) lplus, real, dimension( * )
    work, integer info)

SLARRF finds a new relatively robust representation such that at least one of the eigenvalues is relatively isolated.

๐ŸŽฏ Purpose:

Okay so here's the vibe: you hand SLARRF the initial representation L D L^T and its cluster of close eigenvalues (in a relative measure), namely W( CLSTRT ), W( CLSTRT+1 ), ... W( CLEND ). Then SLARRF finds a new relatively robust representation L D L^T - SIGMA I = L(+) D(+) L(+)^T such that at least one of the eigenvalues of L(+) D(+) L(+)^T is relatively isolated. That's the whole mission, and it ate. ๐Ÿ’…๐Ÿ”ฅ

๐Ÿงฎ Parameters

  • N โ€” N is INTEGER. The order of the matrix (subblock, if the matrix split). ๐Ÿซถ
  • D โ€” D is REAL array, dimension (N). The N diagonal elements of the diagonal matrix D.
  • L โ€” L is REAL array, dimension (N-1). The (N-1) subdiagonal elements of the unit bidiagonal matrix L.
  • LD โ€” LD is REAL array, dimension (N-1). The (N-1) elements L(i)*D(i).
  • CLSTRT โ€” CLSTRT is INTEGER. The index of the first eigenvalue in the cluster.
  • CLEND โ€” CLEND is INTEGER. The index of the last eigenvalue in the cluster.
  • W โ€” W is REAL array, dimension; dimension is >= (CLEND-CLSTRT+1). The eigenvalue APPROXIMATIONS of L D L^T in ascending order. W( CLSTRT ) through W( CLEND ) form the cluster of relatively close eigenalues. (Their spelling, not mine, but we keep it ๐Ÿ’€)
  • WGAP โ€” WGAP is REAL array, dimension; dimension is >= (CLEND-CLSTRT+1). The separation from the right neighbor eigenvalue in W.
  • WERR โ€” WERR is REAL array, dimension; dimension is >= (CLEND-CLSTRT+1). WERR contain the semiwidth of the uncertainty interval of the corresponding eigenvalue APPROXIMATION in W.
  • SPDIAM โ€” SPDIAM is REAL: estimate of the spectral diameter obtained from the Gerschgorin intervals.
  • CLGAPL โ€” CLGAPL is REAL. Short and sweet, bestie. ๐Ÿ’…
  • CLGAPR โ€” CLGAPR is REAL: absolute gap on each end of the cluster. Set by the calling routine to protect against shifts too close to eigenvalues outside the cluster.
  • PIVMIN โ€” PIVMIN is REAL. The minimum pivot allowed in the Sturm sequence.
  • SIGMA โ€” SIGMA is REAL. The shift used to form L(+) D(+) L(+)^T.
  • DPLUS โ€” DPLUS is REAL array, dimension (N). The N diagonal elements of the diagonal matrix D(+).
  • LPLUS โ€” LPLUS is REAL array, dimension (N-1). The first (N-1) elements of LPLUS contain the subdiagonal elements of the unit bidiagonal matrix L(+).
  • WORK โ€” WORK is REAL array, dimension (2*N). Workspace.
  • INFO โ€” INFO is INTEGER. Signals processing OK (=0) or failure (=1).

๐Ÿง‘โ€๐Ÿ”ฌ Author

  • Univ. of Tennessee
  • Univ. of California Berkeley
  • Univ. of Colorado Denver
  • NAG Ltd.

๐ŸŒŸ Contributors:

  • Beresford Parlett, University of California, Berkeley, USA
  • Jim Demmel, University of California, Berkeley, USA
  • Inderjit Dhillon, University of Texas, Austin, USA
  • Osni Marques, LBNL/NERSC, USA
  • Christof Voemel, University of California, Berkeley, USA
Definition at line 189 of file slarrf.f.

๐Ÿ“œ Author

Generated automatically by Doxygen for LAPACK from the source code.


Sign-off: LAPACK, Version 3.12.0, page SRC/slarrf.f(3) โ€” that's the footer, bestie (no extra date line on this one)! ๐Ÿ’–โœจ๐Ÿฆ‹

๐Ÿ“ฆ where this page comes from