SRC/dlasd3.f(3)
π
SRC/dlasd3.f π
Okay bestie, you've wandered into the Library Functions Manual and this page is literally all about SRC/dlasd3.f. It's a LAPACK subroutine page, and no cap, it's giving "math girlie updates singular vectors" energy. Let's go πβ¨
π NAME π
SRC/dlasd3.f
π¦ SYNOPSIS π¦
Functions/Subroutines
subroutine dlasd3 (nl, nr, sqre, k, d, q, ldq, dsigma, u, ldu, u2, ldu2,
vt, ldvt, vt2, ldvt2, idxc, ctot, z, info)
DLASD3 finds all square roots of the roots of the secular equation, as defined by the values in D and Z, and then updates the singular vectors by matrix multiplication. Used by sbdsdc.
πΈ Function/Subroutine Documentation πΈ
subroutine dlasd3 (integer nl, integer nr, integer sqre, integer k, double
precision, dimension( * ) d, double precision, dimension( ldq, * ) q,
integer ldq, double precision, dimension( * ) dsigma, double precision,
dimension( ldu, * ) u, integer ldu, double precision, dimension( ldu2, *
) u2, integer ldu2, double precision, dimension( ldvt, * ) vt, integer
ldvt, double precision, dimension( ldvt2, * ) vt2, integer ldvt2,
integer, dimension( * ) idxc, integer, dimension( * ) ctot, double
precision, dimension( * ) z, integer info)
DLASD3 finds all square roots of the roots of the secular equation, as defined by the values in D and Z, and then updates the singular vectors by matrix multiplication. Used by sbdsdc.
β¨ Purpose β¨
DLASD3finds all the square roots of the roots of the secular equation, as defined by the values inDandZ. It makes the appropriate calls toDLASD4and then updates the singular vectors by matrix multiplication.DLASD3is called fromDLASD1.
π«Ά Parameters π«Ά
NLβNL is INTEGER; the row dimension of the upper block.NL >= 1. πNRβNR is INTEGER; the row dimension of the lower block.NR >= 1. β¨SQREβSQRE is INTEGER= 0: the lower block is anNR-by-NRsquare matrix.= 1: the lower block is anNR-by-(NR+1)rectangular matrix.
The bidiagonal matrix has N = NL + NR + 1 rows and M = N + SQRE >= N columns. π¦
KβK is INTEGER; the size of the secular equation,1 =< K = < N. πDβD is DOUBLE PRECISION array, dimension(K). On exit the square roots of the roots of the secular equation, in ascending order. πΈQβQ is DOUBLE PRECISION array, dimension (LDQ,K). πLDQβLDQ is INTEGER; the leading dimension of the arrayQ.LDQ >= K. β¨DSIGMAβDSIGMA is DOUBLE PRECISION array, dimension(K). The firstKelements of this array contain the old roots of the deflated updating problem. These are the poles of the secular equation. πUβU is DOUBLE PRECISION array, dimension (LDU, N). The lastN - Kcolumns of this matrix contain the deflated left singular vectors. π«ΆLDUβLDU is INTEGER; the leading dimension of the arrayU.LDU >= N. πU2βU2 is DOUBLE PRECISION array, dimension (LDU2, N). The firstKcolumns of this matrix contain the non-deflated left singular vectors for the split problem. πLDU2βLDU2 is INTEGER; the leading dimension of the arrayU2.LDU2 >= N. βVTβVT is DOUBLE PRECISION array, dimension (LDVT, M). The lastM - Kcolumns ofVT**Tcontain the deflated right singular vectors. π¦LDVTβLDVT is INTEGER; the leading dimension of the arrayVT.LDVT >= N. πVT2βVT2 is DOUBLE PRECISION array, dimension (LDVT2, N). The firstKcolumns ofVT2**Tcontain the non-deflated right singular vectors for the split problem. πLDVT2βLDVT2 is INTEGER; the leading dimension of the arrayVT2.LDVT2 >= N. β¨IDXCβIDXC is INTEGER array, dimension ( N ). The permutation used to arrange the columns ofU(and rows ofVT) into three groups: the first group contains non-zero entries only at and above (or before)NL +1; the second contains non-zero entries only at and below (or after)NL+2; and the third is dense. The first column ofUand the row ofVTare treated separately, however.
The rows of the singular vectors found by DLASD4 must be likewise permuted before the matrix multiplies can take place. π
CTOTβCTOT is INTEGER array, dimension ( 4 ). A count of the total number of the various types of columns inU(or rows inVT), as described inIDXC. The fourth column type is any column which has been deflated. π₯ZβZ is DOUBLE PRECISION array, dimension (K). The firstKelements of this array contain the components of the deflation-adjusted updating row vector. πINFOβINFO is INTEGER= 0: successful exit.< 0: ifINFO = -i, thei-th argument had an illegal value.> 0: ifINFO = 1, a singular value did not converge. π
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
πΈ Contributors: πΈ
Ming Gu and Huan Ren, Computer Science Division, University of California at Berkeley, USA.
π Definition at line 214 of file dlasd3.f. π
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
That's the whole page, bestie β source/version sign-off: SRC/dlasd3.f(3) and LAPACK Version 3.12.0 π«Άπβ¨ 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