SRC/dlar1v.f(3)

okay so bestie, this is SRC/dlar1v.f(3) from the Library Functions Manual — manual name folded in right here like a little name tag. 🦋 It’s a LAPACK Fortran page, and it’s obsessed with exactly one subroutine: dlar1v.

NAME

SRC/dlar1v.f

SYNOPSIS

Functions/Subroutines

subroutine dlar1v (n, b1, bn, lambda, d, l, ld, lld, pivmin, gaptol, z,
    wantnc, negcnt, ztz, mingma, r, isuppz, nrminv, resid, rqcorr, work)
    DLAR1V computes the (scaled) r-th column of the inverse of the
    submatrix in rows b1 through bn of the tridiagonal matrix LDLT - λI.

Function/Subroutine Documentation

Here’s the full typed-out signature, no cap:

subroutine dlar1v (integer n, integer b1, integer bn, double precision
    lambda, double precision, dimension( * ) d, double precision, dimension(
    * ) l, double precision, dimension( * ) ld, double precision, dimension(
    * ) lld, double precision pivmin, double precision gaptol, double
    precision, dimension( * ) z, logical wantnc, integer negcnt, double
    precision ztz, double precision mingma, integer r, integer, dimension( *
    ) isuppz, double precision nrminv, double precision resid, double
    precision rqcorr, double precision, dimension( * ) work)
    DLAR1V computes the (scaled) r-th column of the inverse of the submatrix
    in rows b1 through bn of the tridiagonal matrix LDLT - λI.

Purpose 🎯

DLAR1V computes the (scaled) r-th column of the inverse of the submatrix in rows B1 through BN of the tridiagonal matrix L D L**T - sigma I. When sigma is close to an eigenvalue, the computed vector is an accurate eigenvector. Usually, r corresponds to the index where the eigenvector is largest in magnitude.

The computation does these steps, in order:

• (a) Stationary qd transform: L D L**T - sigma I = L(+) D(+) L(+)**T.
• (b) Progressive qd transform: L D L**T - sigma I = U(-) D(-) U(-)**T.
• (c) Computation of the diagonal elements of the inverse of L D L**T - sigma I by combining the above transforms, and choosing r as the index where the diagonal of the inverse is (one of the) largest in magnitude.
• (d) Computation of the (scaled) r-th column of the inverse using the twisted factorization obtained by combining the top part of the stationary and the bottom part of the progressive transform.

Parameters 🎀

• N — N is INTEGER. The order of the matrix L D L**T.
• B1 — B1 is INTEGER. First index of the submatrix of L D L**T.
• BN — BN is INTEGER. Last index of the submatrix of L D L**T.
• LAMBDA — LAMBDA is DOUBLE PRECISION. The shift. In order to compute an accurate eigenvector, LAMBDA should be a good approximation to an eigenvalue of L D L**T.
• L — L is DOUBLE PRECISION array, dimension (N-1). The (n-1) subdiagonal elements of the unit bidiagonal matrix L, in elements 1 to N-1.
• D — D is DOUBLE PRECISION array, dimension (N). The n diagonal elements of the diagonal matrix D.
• LD — LD is DOUBLE PRECISION array, dimension (N-1). The n-1 elements L(i)*D(i).
• LLD — LLD is DOUBLE PRECISION array, dimension (N-1). The n-1 elements L(i)*L(i)*D(i).
• PIVMIN — PIVMIN is DOUBLE PRECISION. The minimum pivot in the Sturm sequence.
• GAPTOL — GAPTOL is DOUBLE PRECISION. Tolerance that indicates when eigenvector entries are negligible w.r.t. their contribution to the residual.
• Z — Z is DOUBLE PRECISION array, dimension (N). On input, all entries of Z must be set to 0. On output, Z contains the (scaled) r-th column of the inverse. The scaling is such that Z(R) equals 1.
• WANTNC — WANTNC is LOGICAL. Specifies whether NEGCNT has to be computed.
• NEGCNT — NEGCNT is INTEGER. If WANTNC is .TRUE. then NEGCNT = the number of pivots < pivmin in the matrix factorization L D L**T, and NEGCNT = -1 otherwise.
• ZTZ — ZTZ is DOUBLE PRECISION. The square of the 2-norm of Z.
• MINGMA — MINGMA is DOUBLE PRECISION. The reciprocal of the largest (in magnitude) diagonal element of the inverse of L D L**T - sigma I.
• R — R is INTEGER. The twist index for the twisted factorization used to compute Z. On input, 0 <= R <= N. If R is input as 0, R is set to the index where (L D L**T - sigma I)^{-1} is largest in magnitude. If 1 <= R <= N, R is unchanged. On output, R contains the twist index used to compute Z. Ideally, R designates the position of the maximum entry in the eigenvector.
• ISUPPZ — ISUPPZ is INTEGER array, dimension (2). The support of the vector in Z, i.e., the vector Z is nonzero only in elements ISUPPZ(1) through ISUPPZ( 2 ).
• NRMINV — NRMINV is DOUBLE PRECISION. NRMINV = 1/SQRT( ZTZ ).
• RESID — RESID is DOUBLE PRECISION. The residual of the FP vector. RESID = ABS( MINGMA )/SQRT( ZTZ ).
• RQCORR — RQCORR is DOUBLE PRECISION. The Rayleigh Quotient correction to LAMBDA. RQCORR = MINGMA*TMP.
• WORK — WORK is DOUBLE PRECISION array, dimension (4*N).

🚨 real talk, bestie: on input, every entry of Z must be set to 0. That’s not a vibe suggestion, that’s the spec.

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 227 of file dlar1v.f.

Author 🌟 (the automatic one)

Generated automatically by Doxygen for LAPACK from the source code.

---

💖 That’s the whole vibe: LAPACK Version 3.12.0, SRC/dlar1v.f(3) — date not shown on the page, so no made-up date, periodt. 💖