/home/abui...t01_rook.f(3)
π NAME π
okay so bestie, this is the Library Functions Manual page for /home/abui...t01_rook.f(3), and the NAME line is literally the full path below (yes, the whole source path, no cap):
/home/abuild/rpmbuild/BUILD/lapack-
static-3.12.0-build/lapack-3.12.0/TESTING/LIN/dsyt01_rook.f
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine DSYT01_ROOK (uplo, n, a, lda, afac, ldafac, ipiv, c, ldc,
rwork, resid)
and then it literally says DSYT01_ROOK again underneath, because iconic behavior π
π¦ Function/Subroutine Documentation π¦
DSYT01_ROOK is the main character here, and here's the full signature, copied exactly because we do not play with code:
subroutine DSYT01_ROOK (character uplo, integer n, double precision,
dimension( lda, * ) a, integer lda, double precision, dimension( ldafac,
* ) afac, integer ldafac, integer, dimension( * ) ipiv, double
precision, dimension( ldc, * ) c, integer ldc, double precision,
dimension( * ) rwork, double precision resid)
then, yes, DSYT01_ROOK appears again under the signature. slay β¨
π Purpose π
okay so here's the tea: DSYT01_ROOK reconstructs a symmetric indefinite matrix A from its block L*D*L' or U*D*U' factorization and computes the residual
norm( C - A ) / ( N * norm(A) * EPS )
where C is the reconstructed matrix and EPS is the machine epsilon. it's giving factorization vibe check π₯
π Parameters π
UPLOβ isCHARACTER*1. Specifies whether the upper or lower triangular part of the symmetric matrixAis stored:= 'U': Upper triangular;= 'L': Lower triangular. pick your triangle, bestie πNβ isINTEGER. The number of rows and columns of the matrixA.N >= 0.Aβ isDOUBLE PRECISION array, dimension (LDA,N). The original symmetric matrixA.LDAβ isINTEGER. The leading dimension of the arrayA.LDA >= max(1,N)AFACβ isDOUBLE PRECISION array, dimension (LDAFAC,N). The factored form of the matrixA.AFACcontains the block diagonal matrixDand the multipliers used to obtain the factorLorUfrom the blockL*D*L'orU*D*U'factorization as computed byDSYTRF_ROOK.LDAFACβ isINTEGER. The leading dimension of the arrayAFAC.LDAFAC >= max(1,N).IPIVβ isINTEGER array, dimension (N). The pivot indices fromDSYTRF_ROOK.Cβ isDOUBLE PRECISION array, dimension (LDC,N).LDCβ isINTEGER. The leading dimension of the arrayC.LDC >= max(1,N).RWORKβ isDOUBLE PRECISION array, dimension (N).RESIDβ isDOUBLE PRECISION. IfUPLO = 'L',norm(L*D*L' - A) / ( N * norm(A) * EPS ); IfUPLO = 'U',norm(U*D*U' - A) / ( N * norm(A) * EPS ). residual said what it said π
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 122 of file dsyt01_rook.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
cute sign-off bestie: this footer is source/version LAPACK Version 3.12.0 and title /home/abui...t01_rook.f(3) β thereβs no separate date printed on that footer line, so Iβm not inventing one π
π«ΆπΈ
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms