SRC/dsytri2.f(3)
okay bestie, welcome to SRC/dsytri2.f(3) from the Library Functions Manual β no cap, we're about to make this LAPACK routine make sense. π
πΈ NAME πΈ
SRC/dsytri2.f
β¨ SYNOPSIS β¨
π«Ά Functions/Subroutines
Here's the group-chat receipt for the callable stuff:
subroutine dsytri2 (uplo, n, a, lda, ipiv, work, lwork, info)
DSYTRI2
π» Function/Subroutine Documentation π»
π subroutine dsytri2 π
subroutine dsytri2 (character uplo, integer n, double precision, dimension(
lda, * ) a, integer lda, integer, dimension( * ) ipiv, double precision,
dimension( * ) work, integer lwork, integer info)
DSYTRI2
π Purpose π
okay so the whole point: DSYTRI2 computes the inverse of a DOUBLE PRECISION symmetric indefinite matrix A using the factorization A = U*D*U**T or A = L*D*L**T computed by DSYTRF. DSYTRI2 sets the LEADING DIMENSION of the workspace before calling DSYTRI2X that actually computes the inverse. That's the serve. π
π Parameters π
One bullet per parameter, bestie β every value stays exactly as it was:
UPLOβ isCHARACTER*1. It tells you whether the details of the factorization are stored as an upper or lower triangular matrix.= 'U': Upper triangular, form is A = U*D*U**T;= 'L': Lower triangular, form is A = L*D*L**T.Nβ isINTEGER. The order of the matrixA.N >= 0.Aβ isDOUBLE PRECISION array, dimension (LDA,N). On entry, it's the block diagonal matrixDand the multipliers used to obtain the factorUorLas computed byDSYTRF. On exit, ifINFO = 0, it's the (symmetric) inverse of the original matrix. IfUPLO = 'U', the upper triangular part of the inverse is formed and the part ofAbelow the diagonal is not referenced; ifUPLO = 'L', the lower triangular part of the inverse is formed and the part ofAabove the diagonal is not referenced. This is the main character transformation arc, literally. β¨LDAβ isINTEGER. The leading dimension of the arrayA.LDA >= max(1,N).IPIVβ isINTEGER array, dimension (N). Details of the interchanges and the block structure ofDas determined byDSYTRF.WORKβ isDOUBLE PRECISION array, dimension (MAX(1,LWORK)).LWORKβ isINTEGER. The dimension of the arrayWORK. IfN = 0,LWORK >= 1, elseLWORK >= (N+NB+1)*(NB+3). IfLWORK = -1, then a workspace query is assumed; the routine calculates:- the optimal size of the
WORKarray, returns this value as the first entry of theWORKarray, - and no error message related to
LWORKis issued byXERBLA.
- the optimal size of the
INFOβ isINTEGER.= 0: successful exit.< 0: ifINFO = -i, the i-th argument had an illegal value.> 0: ifINFO = i,D(i,i) = 0; the matrix is singular and its inverse could not be computed.
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Definition at line π
Definition at line 126 of file dsytri2.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
mwah π β this page is from LAPACK, Version 3.12.0, for SRC/dsytri2.f(3). That's the footer, bestie β go compute that inverse and slay. π«Άβ¨
π¦ 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