SRC/slaed4.f(3)
π NAME π
okay bestie, this page is from the Library Functions Manual and itβs all about SRC/slaed4.f π like, a whole LAPACK subroutine page, no cap.
β¨ SYNOPSIS β¨
π¦ Functions/Subroutines π¦
subroutine slaed4 (n, i, d, z, delta, rho, dlam, info)β SLAED4 used by SSTEDC. Finds a single root of the secular equation.
π Function/Subroutine Documentation π
π subroutine slaed4 π
subroutine slaed4 (integer n, integer i, real, dimension( * ) d, real,
dimension( * ) z, real, dimension( * ) delta, real rho, real dlam,
integer info)
SLAED4 used by SSTEDC. Finds a single root of the secular equation. π«Ά
β¨ Purpose β¨
okay so this subroutine computes the I-th updated eigenvalue of a symmetric rank-one modification to a diagonal matrix whose elements are given in the array d, and that:
D(i) < D(j) for i < j
and that RHO > 0. π This is arranged by the calling routine, and is no loss in generality. The rank-one modified system is thus:
diag( D ) + RHO * Z * Z_transpose.
where we assume the Euclidean norm of Z is 1. π
The method consists of approximating the rational functions in the secular equation by simpler interpolating rational functions. literally iconic π
π§ Parameters π§
NβN is INTEGER. The length of all arrays.IβI is INTEGER. The index of the eigenvalue to be computed.1 <= I <= N.DβD is REAL array, dimension (N). The original eigenvalues. It is assumed that they are in order,D(I) < D(J) for I < J.ZβZ is REAL array, dimension (N). The components of the updating vector.DELTAβDELTA is REAL array, dimension (N). IfN > 2,DELTAcontains(D(j) - lambda_I)in its j-th component. IfN = 1, thenDELTA(1) = 1. IfN = 2, seeSLAED5for detail. The vectorDELTAcontains the information necessary to construct the eigenvectors bySLAED3andSLAED9.RHOβRHO is REAL. The scalar in the symmetric updating formula.DLAMβDLAM is REAL. The computedlambda_I, theI-th updated eigenvalue.INFOβINFO is INTEGER.= 0: successful exit> 0: ifINFO = 1, the updating process failed.
π₯ Internal Parameters π₯
- Logical variable
ORGATI(origin-at-i?) is used for distinguishing whetherD(i)orD(i+1)is treated as the origin.
ORGATI = .true. origin at i
ORGATI = .false. origin at i+1
- Logical variable
SWTCH3(switch-for-3-poles?) is for noting if we are working with THREE poles! π MAXITis the maximum number of iterations allowed for each eigenvalue.
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
πΈ Contributors πΈ
Ren-Cang Li, Computer Science Division, University of California at Berkeley, USA π
β Definition β
Definition at line 144 of file slaed4.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
bye bestie π this footer moment is brought to you by LAPACK Version 3.12.0 and the source path SRC/slaed4.f(3) β slay, periodt. β¨
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms